[R-sig-ME] Error when using models with lmer
Douglas Bates
bates at stat.wisc.edu
Wed Feb 8 17:51:30 CET 2012
On Tue, Feb 7, 2012 at 2:03 PM, Olga Vitek <ovitek at stat.purdue.edu> wrote:
> Hello,
>
> I am executing the following code
>
>> library("lme4")
>> library(faraway)
>> data(penicillin)
>> fit4 <- lmer(yield ~ treat + (1|blend), penicillin)
>
> and the lines below generate error messages
>
>> library(gmodels)
>> ci(fit4)
Yes. Why would you want to try to create a vector from an S4 fitted
model object?
It almost never makes sense to use
c(foo)
Most of the time when people use that idiom what they really mean is
foo
(I say "almost never" because in the old days many of us would use
c(myMatrix) to create a vector from a matrix. The preferred form is
now as.vector(myMatrix).)
> Error in as.vector(data) :
> no method for coercing this S4 class to a vector
>
>> contrast.function <- c(0, -1, 0, 0)
>> test.result <- estimable(fit4, contrast.function)
> Error in as.vector(data) :
I assume that the "estimable" function is in the gmodels package, so
you should contact the author of that package about this.
> no method for coercing this S4 class to a vector
>
>
> This worked with with lme4_0.999375-40 and R version 2.13.1, but not with lme4_0.999375-42 and R version 2.14.1.
>
> Would you have suggestions on how to make this work?
> Thank you in advance
> Olga Vitek
>
>> sessionInfo()
> R version 2.14.1 (2011-12-22)
> Platform: i386-apple-darwin9.8.0/i386 (32-bit)
>
> locale:
> [1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
>
> attached base packages:
> [1] stats graphics grDevices utils datasets methods base
>
> other attached packages:
> [1] gmodels_2.15.1 faraway_1.0.5 lme4_0.999375-42 Matrix_1.0-3 lattice_0.20-0
>
> loaded via a namespace (and not attached):
> [1] gdata_2.8.2 grid_2.14.1 gtools_2.6.2 MASS_7.3-16 nlme_3.1-102 stats4_2.14.1 tools_2.14.1
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
More information about the R-sig-mixed-models
mailing list