## Data from the Lerna population.
nonadult <- subset(lerna, lerna$age <= 10)
adult <- subset(lerna, lerna$age > 17)
female <- subset(adult, adult$sex == "f")
## Simple optimization without probability.
(opt.lerna <- warnOptim(
age = nonadult$age,
d15N = nonadult$d15N,
female.mean = mean(female$d15N),
control = list(maxit = 10000, ndeps = 1e-2, reltol = 1e-7)))
Run the code above in your browser using DataLab