ages <- seq(0,20,1) # sigma is ignored here
oldpar <- par(no.readonly=TRUE)
pars <- c(a=26.0,b=2.8,c=-0.65,sigma=1.0) # a, b, c, sigma;
plot1(ages,Gz(pars,ages),xlab="Age",ylab="Length",lwd=2)
par(oldpar)
Run the code above in your browser using DataLab