[R] NLME augPred Error
Yang, Richard
dyang at NRCan.gc.ca
Wed Apr 7 17:20:38 CEST 2004
Dear all;
This error has been posted to the list a couple of times, but no
answer / walk around solution is available. An example on P. 377 MEMSS
>plot(augPred(fm5CO2.nlme, levels=0:1), layout =c(6,2))
generates an error: "Error in predict.nlme(object, value[1:(nrow(value)/nL),
, drop = FALSE], : Levels Quebec, Mississippi not allowed for Type.
for nlme ver. 3.1-45 in R.1.8.1, but it works fine in Splus6.2 or earlier
versions.
The source code for augPred() is very cryptica and identical in both S+ and
R:
function (object, primary = NULL, minimum = min(primary), maximum =
max(primary),
length.out = 51, ...)
UseMethod("augPred").
The augPred() creates an augPred class object in S+ but generates the above
error in R. Is there any difference in nlme3 in S+ and nlme versions in R?
Any ideas?
Richard
More information about the R-help
mailing list