[R] Error in nlme ranef plot()
dyang at nrcan.gc.ca
Mon Apr 8 21:51:11 CEST 2002
Dear R list members;
I have a 10 x 423 data frame which consisting of response, time, subject,
site, plot and covariates (continueous and categorical) measured at the plot
level. When the data frame was converted into a groupedData object, a
> A <- groupedData(ht ~ time | Subject, data = tt, outer = ~ site * plot,
+ labels=list(y = "Height", x = "Time", units = list( y = "(m)", x =
argument lengths differ in: split(x, f).
The warning did not affect the estimation of a nlme model with parameters
B1, B2, and B3. However, when plotting the estimated ranef versus
covariates, an error resulted
> A1Mod.nlmeRE <- ranef(A1Mod.nlme, aug=T)
> plot(A1Mod.nlmeRE, form = B1 ~ Ps + Dr + El + As + Sl) results in
Error in max(length(x0), length(x1), length(y0), length(y1)) :
Argument "x0" is missing, with no default
The same data frame and commands works fine with Splus 6 on W2K. Is this a
bug in nlme 3.3.x? I use R 1.4.1.
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch
More information about the R-help