# NOT RUN {
data(blackisland)
oldpar <- par(no.readonly=TRUE)
plot(blackisland$l1,blackisland$dl,type="p",pch=16,
xlab="Initial Length mm",ylab="Growth Increment mm",panel.first=grid())
abline(h=0)
param <- c(25, 130, 35, 3) # MaxDL, L50, delta, sigma
predDL <- invl(param,blackisland,initL="l1",delT="dt")
lines(blackisland$l1,predDL,col=2,lwd=2)
par(oldpar)
# }
Run the code above in your browser using DataLab