[R] Modelling compound logistic growth curves
Dan Bebber
danbebber at forestecology.co.uk
Wed Jul 28 17:24:31 CEST 2004
Motivated by the discovery of 'loglet analysis'
(http://phe.rockefeller.edu/LogletLab/) that allows one to decompose growth
curves into a series of logistic equations, I attempted to do the same thing
in R.
SIMULATED DATA
Time <- 1:200
pop.size <- SSlogis(Time,10,20,5) + SSlogis(Time,20,100,20) +
rnorm(length(Time))
MY ANALYSIS
results <- nls(size ~ SSlogis(Time, Asym1, xmid1, scal1) + SSlogis(Time,
Asym2, xmid2, scal2),
start = list(Asym1=5, xmid1=15, scal1=30, Asym2=25, xmid2=67, scal2=25))
THE RESULT
I get the error message:
Error in nls(size ~ SSlogis(Time, Asym1, xmid1, scal1) + SSlogis(Time, :
step factor 0.000488281 reduced below `minFactor' of 0.000976563
Assistance in doing this analysis would be much appreciated.
Dan Bebber
____________________________
Dr. Daniel P. Bebber
Department of Plant Sciences
University of Oxford
South Parks Road
Oxford
OX1 3RB
Tel. 01865 275060
Web. http://www.forestecology.co.uk/
"Data, data, data!" he cried impatiently. "I can't make bricks without
clay"
- Sherlock Holmes, The Adventure of the Copper Beeches, 1892
More information about the R-help
mailing list