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中排除了...