Anova lmer从方差分析中获得p值

Anova lmer从方差分析中获得p值,anova,mixed-models,p-value,Anova,Mixed Models,P Value,lmerTest被设计为一个包装器,允许使用分母自由度(ddf)的Satterthwaite估计值,从lmer混合模型分析中估计p值。但是lmerTest现在似乎被破坏了。目前,它返回一条消息,指出存在内部计算错误,并仅返回lmer结果(无p值)。我已经能够通过summary()函数计算p值,使用Dan Mirman的优秀代码计算Kenward Rogers的估计值。但在lmer模型的anova调用中,我找不到计算p值的等效代码。我怀疑一个人只需要输入一个ddf,但我不知道怎么做 提前感谢任何能

lmerTest
被设计为一个包装器,允许使用分母自由度(ddf)的Satterthwaite估计值,从lmer混合模型分析中估计
p值。但是
lmerTest
现在似乎被破坏了。目前,它返回一条消息,指出存在内部计算错误,并仅返回lmer结果(无
p值
)。我已经能够通过summary()函数计算
p值
,使用Dan Mirman的优秀代码计算Kenward Rogers的估计值。但在lmer模型的
anova
调用中,我找不到计算
p值的等效代码。我怀疑一个人只需要输入一个ddf,但我不知道怎么做

提前感谢任何能为这个问题提出解决方案的人


Larry Hunsicker

lmerTest在获得Satterthwaite近似值时(例如,在计算渐近方差协方差矩阵时)出现计算错误时,返回lme4包的方差分析输出。lmerTest并没有被打破,只是有一些例子表明Satterthwaite的近似值无法计算。根据我的经验,这种情况并不经常发生。

lmerTest在获得Satterthwaite近似值时(例如,在计算渐近方差协方差矩阵时)出现计算错误时,会返回lme4包的方差分析输出。lmerTest并没有被打破,只是有一些例子表明Satterthwaite的近似值无法计算。根据我的经验,这种情况并不经常发生。

对不起,亚历山德拉。我现在明白了,问题在于Satterthwaite近似的计算,而不是lmerTest。我现在也找到了上述问题的答案。lmerTest包中的anova函数提供了请求ddf的Kenwood-Roger近似值(ddf='Ken')的选项,而anova调用正好给出了我想要的结果。我还没有找到在lmer函数中请求K-R ddf近似值的方法。在lmer中也有使用K-R的方法吗?如果不是,这可能是一个很好的补充到您的包。感谢lmerTest。要使用Kenward-Roger的方差分析近似值,您可以编写方差分析(model,ddf=“kenw”),并在摘要函数中使用它,您可以类似地编写:summary(model,ddf=“kenw”)。lmerTest包使用pbkrtest包实现Kenward Roger的近似值My抱歉,Alexandra。我现在明白了,问题在于Satterthwaite近似的计算,而不是lmerTest。我现在也找到了上述问题的答案。lmerTest包中的anova函数提供了请求ddf的Kenwood-Roger近似值(ddf='Ken')的选项,而anova调用正好给出了我想要的结果。我还没有找到在lmer函数中请求K-R ddf近似值的方法。在lmer中也有使用K-R的方法吗?如果不是,这可能是一个很好的补充到您的包。感谢lmerTest。要使用Kenward-Roger的方差分析近似值,您可以编写方差分析(model,ddf=“kenw”),并在摘要函数中使用它,您可以类似地编写:summary(model,ddf=“kenw”)。lmerTest包使用pbkrtest包对Kenward-Roger近似的实现