获取nlme或lme4中固定效果的广义最小二乘均值
-
27-10-2019 - |
题
aov
对象的最小二乘均值及其标准误差可以通过model.tables
函数获得:
通用标签
我想知道如何从nlme
对象或lme4
对象获取具有标准误差的广义最小二乘均值:
通用标签
任何评论和提示将不胜感激。谢谢
解决方案
lme和nlme适合最大似然或受限最大似然(后者是默认值),因此您的结果将基于这两种方法中的任意一种
summary(fm1Machine)
将为您提供包括均值和标准错误的输出:
通用标签
因为您已将固定效果与截距拟合,所以您会在固定效果结果中获得截距项,而不是MachineA的结果。MachineB和MachineC的结果与截距相反,因此要获得MachineB和MachineC的均值,请将每个值加到截距平均值上。但是标准错误不是您想要的。
要获取所需信息,请对模型进行拟合,以使其在固定效果中没有截距项(请参阅固定效果结尾处的-1
:
通用标签
这将为您提供所需的平均值和标准错误输出: 通用标签
其他提示
引用的道格拉斯·贝茨(Douglas Bates)
http://markmail.org/message/dqpk6ftztpbzgekm
“我强烈怀疑,对于大多数用户而言,lsmeans的定义是“我使用lsmeans语句时从SAS获得的数字”。我建议获得此类数字的方法是购买SAS许可证并使用SAS适合您的模型。“
不隶属于 StackOverflow