3 agomp200 bgomp308 cgomp03

Info iconThis preview shows page 1. Sign up to view the full content.

View Full Document Right Arrow Icon
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: am.tree*c.gomp))+rnorm(200,0,2) plot(diam.tree,height.gomp,type="n",xlab="Tree Diameter",ylab="Tree Height", ylim=c(130,250)) points(diam.tree,height.gomp,pch=15) lines(sort(diam.tree),sort(pred.gomp),col=15,lwd=3) range(height.gomp) title("Gompertz gamma1 original") # Gompertz model # parameters: a=200,b=3,c=0.3 # a.gomp = 200 b.gomp = 3 * 0.8 c.gomp = 0.3 # pred.gomp = a.gomp*exp(­b.gomp*exp(­diam.tree*c.gomp)) height.gomp = a.gomp*exp(­b.gomp*exp(­diam.tree*c.gomp))+rnorm(200,0,2) plot(diam.tree,height.gomp,type="n",xlab="Tree Diameter",ylab="Tree Height", ylim=c(130,250)) points(diam.tree,height.gomp,pch=15) lines(sort(diam.tree),sort(pred.gomp),col=15,lwd=3) range(height.gomp) title("Gompertz gamma1 20% smaller") # Gompertz model # parameters: a=200,b=3,c=0.3 # a.gomp = 200 b.gomp = 3 * 1.2 c.gomp = 0.3 # pred.gomp = a.gomp*exp(­b.gomp*exp(­diam.tree*c.gomp)) height.gomp = a.gomp*exp(­b.gomp*exp(­diam.tree*c.gomp))+rnorm(200,0,2) plot(diam.tree,h...
View Full Document

Ask a homework question - tutors are online