[R] Bootstrapping nlme models

lynnland lynn.landriault at ontario.ca
Fri Jan 6 20:31:40 CET 2012


Hi, Let me start my saying that I am new to R hence my grasp of the
appropriate used of R coding is undoubtedly way behind many on this forum.

I am trying to use boostrapping to derive errors around my parameter
estimate for the fixed effects in the following model. It is simply
estimating the number of times an animal might cross a road based on the
road's distance from a stream. I have tried several different methods but
have been unsuccessful. I have tried both using the boot package and simply
developing a bit of code that resamples my data and drops the parameter
estimates into a new vector. I would be interested in running both a
parametric and non-parametric version.  

model2<-nlme(nwfcross~a*exp(-b*nwfdist), 
fixed=a+b~1, 
random=a+b~1|nwfid,
start=c(a=300,b=0.016))

Any help would be greatly appreciated. If it is helpful I can add in some of
the code I have tried and the resulting error messages.

Lynn




--
View this message in context: http://r.789695.n4.nabble.com/Bootstrapping-nlme-models-tp4270557p4270557.html
Sent from the R help mailing list archive at Nabble.com.



More information about the R-help mailing list