Displaying 1 result from an estimated 1 matches for "tmmean".
Did you mean:
tmean
2010 Oct 20
0
bootstrapping nonlinear mixed effects models
...ot <- nls(weight ~ SSlogis(Time, a, b, c), data=Soybean[indices, ])
return(coef(mboot))
}
require(boot)
boot1 <- boot(Soybean, bstat, 2000)
boot1
For nonlinear mixed effects models it doesn't work.
> stat.nlme <- function(A, indices) {
+ boot.nlme <- nlme(Gdiff/0.5 ~ SSlogis(tmmean,a,b,c),
+ fixed=list(a+b+c~1),
+ random = list(a+b+c~1), data = Soybean[indices,])
+ return(fixef(boot.nlme))
+ }
>
> require(boot)
> bb <- boot(Soybean, stat.nlme, 2000)
Error in parse(text = paste("~", paste(nVal, collapse = "/"))) :
unexpected end of input...