关于matlab和r语言做线性混合效应模型的问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

matlab里面做线性混合效应模型可以用fitlme函数，r语言可以用nlme包中的lme函数，用这两个函数构建相同的混合效应模型，发现自由度不同，最后结果也相差很大，请问有同学知道怎么回事吗？
matlab调用函数：
lme=fitlme(table,'Y~g+t+g:t+(1|subject)','FitMethod','REML')
g t 和subject都设置为dummy变量
得到residual的自由度为95

r语言调用函数：
model<-lme(Y~g*t,random=~1|subject,data=table)
g t和subject都设置为factor
得到residual的自由度为76

自由度的差值看起来是和subject有关，与matlab相比，r似乎将随机效应(subject的截距)的自由度从residual中排除了...

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：MATLAB 混合效应模型 atlab matla Mat

大多数情况下，R更具有统计专业水准。

llb_321 发表于 2020-7-19 01:06
大多数情况下，R更具有统计专业水准。

我也感觉是，主要是matlab还涉及到其它功能，现在就是不得不现在matlab里面做一下统计....

llb_321 发表于 2020-7-19 01:06
大多数情况下，R更具有统计专业水准。

我也感觉是，主要是matlab还涉及到其它功能，现在就是不得不先在matlab里面做一下统计....