Question

我正在运行这样的模型：

model<-lme(Juice~1+Fruit ,random=list(Place=pdIdent(~1)), method="REML",
                           control=lmeControl(niterEM=150 ,msMaxIter=200), na.action=na.omit,
                           data=customData, keep.data=FALSE)

我想为每种水果提取果汁的调整平均值。我在lsmean()软件包上找到了一些涉及lmerTest的解决方案，但它引发了错误说：

该模型不是线性混合效应模型

所以，这里有点无能为力。

对于可重现的示例（这是我的数据看起来如何）：

Fruit=c(sort(rep(1:40,5)),sort(rep(1:40,5)))
Juice=c(rnorm(200, mean=15:25),rnorm(200, mean=13:23))
Place=c(rep("A",200), rep("B",200))

customData=data.frame(Fruit,Juice,Place)

提前致谢。

Answer 1

由于@ eipi10已经声明lsmeans（）不适用于lme4包生成的对象。 lsmeans库可能会给你你想要的东西：

library(lsmeans)
lsmeans(model, "Fruit")

提取lme函数的调整方法

1 个答案: