Displaying 3 results from an estimated 3 matches for "nlme1".
Did you mean:
nlme
2008 Mar 24
0
Error using augPred with nlme grouped data
...the vonBertelanffy
growth curve using:
vanfemfunc<-function(age.years, A, lk, t0)
A*((1-exp(-exp(lk)*(age.years-t0)))^3)
vf.list<-nlsList(mass~vanfemfunc(age.years, A, lk, t0), data=vanfem,
start=c("A"=201.8, "lk"=-1.17, "t0"=-2.71))
and then my nlme with:
vf.nlme1<-nlme(vf.list)
Following the methods of Pinheiro and Bates (2000), I try and plot my
predicted values against my actual values by "name" using:
plot(augPred(vf.nlme1, level=0:1))
But it returns this error:
Error in `[[<-.data.frame`(`*tmp*`, nm, value = c(3L, 3L, 3L, 3L, 3L, :...
2006 Apr 10
3
SE estimates for treatment groups from nlme
I am wondering how to obtain SE estimates for fixed effects from a nonlinear mixed effects model?
I have fixed effects corresponding to three factors A, B and C with 2, 3 and 3 levels respectively. I have fit a model of the following general form:
nlme1<-nlme(y~ SasympOrig(x, Asym, lrc), data=df, fixed=list(Asym~A*B*C, lrc~A*B*C),
start=c(fixef(ETR.nlme)[1], rep(0,17), fixef(ETR.nlme)[2], rep(0,17)))
I am using the default ("contr.treatment" and "contr.poly"). The summary table (tTable) gives me the baseline coefficients...
2005 Apr 25
4
panel ordering in nlme and augPred plots
Dear all
I am trying nlme together with Pinheiro/Bates book. I constructed
grouped data object with suitable plotting layout (according to
some common factor, panels from bottom to top are in increasing
order).
When I do nlme(... some stuff...) I get fitted object which I can plot
with
plot(augPred(fit.nlme6, level=0:1))
but it results in completely different ordering. Is there any way