获取nlme或lme4中固定效果的广义最小二乘均值

Question

aov对象的最小二乘均值及其标准误差可以通过model.tables函数获得： 通用标签

我想知道如何从nlme对象或lme4对象获取具有标准误差的广义最小二乘均值： 通用标签

任何评论和提示将不胜感激。谢谢

Answer 1

lme和nlme适合最大似然或受限最大似然（后者是默认值），因此您的结果将基于这两种方法中的任意一种

summary(fm1Machine)将为您提供包括均值和标准错误的输出： 通用标签

因为您已将固定效果与截距拟合，所以您会在固定效果结果中获得截距项，而不是MachineA的结果。MachineB和MachineC的结果与截距相反，因此要获得MachineB和MachineC的均值，请将每个值加到截距平均值上。但是标准错误不是您想要的。

要获取所需信息，请对模型进行拟合，以使其在固定效果中没有截距项（请参阅固定效果结尾处的-1： 通用标签

这将为您提供所需的平均值和标准错误输出： 通用标签

Answer 2

引用的道格拉斯·贝茨（Douglas Bates）

http://markmail.org/message/dqpk6ftztpbzgekm

“我强烈怀疑，对于大多数用户而言，lsmeans的定义是“我使用lsmeans语句时从SAS获得的数字”。我建议获得此类数字的方法是购买SAS许可证并使用SAS适合您的模型。“