h******3 发帖数: 190 | 1 MS都是用来fit linear mixed model的。用过这几个function的同学能告诉我区别么? |
p***r 发帖数: 920 | 2 developer 不一样,模型的输入规则稍有差别,lmer能计算东西的东西更多一些 |
w********3 发帖数: 1503 | 3 gls不清楚,lmer比lme好用我觉得。。。有时候lme跑不出结果来,用lmer就行了。。
。 |
y***g 发帖数: 10422 | 4 gls: Linear model using generalized least squares. The errors are allowed to
be correlated and/or have unequal variances.
不过我从来没用过。。。
lme: Linear Mixed Effect Model
lmer: Generalized Linear Mixed Effect Model
所在的 library 不同,也就是 developer 不一样。gls 和 lme 是 nlme 里面的。
lmer 是 lme4 里面的。
我去年做一个数据分析的时候发现,同样的 data 和 model,在32位和64位版本的
lmer/glmer 上竟然能出现很不同的结果。我的 model 是 logistic mixed effect
model。后来就不太敢用它。好在我需要做 non-linear mixed effect model 的时候也
不多。 |
h******3 发帖数: 190 | 5 Thanks!
OK.所以linear mixed model的时候lmer和lme是通用的
to
【在 y***g 的大作中提到】 : gls: Linear model using generalized least squares. The errors are allowed to : be correlated and/or have unequal variances. : 不过我从来没用过。。。 : lme: Linear Mixed Effect Model : lmer: Generalized Linear Mixed Effect Model : 所在的 library 不同,也就是 developer 不一样。gls 和 lme 是 nlme 里面的。 : lmer 是 lme4 里面的。 : 我去年做一个数据分析的时候发现,同样的 data 和 model,在32位和64位版本的 : lmer/glmer 上竟然能出现很不同的结果。我的 model 是 logistic mixed effect : model。后来就不太敢用它。好在我需要做 non-linear mixed effect model 的时候也
|
h******3 发帖数: 190 | 6 但是lmer的结果没有直接的p-value显示。这点很不方便。。
【在 w********3 的大作中提到】 : gls不清楚,lmer比lme好用我觉得。。。有时候lme跑不出结果来,用lmer就行了。。 : 。
|
B******5 发帖数: 4676 | 7 developer一样,一个S3,一个S4罢了
to
【在 y***g 的大作中提到】 : gls: Linear model using generalized least squares. The errors are allowed to : be correlated and/or have unequal variances. : 不过我从来没用过。。。 : lme: Linear Mixed Effect Model : lmer: Generalized Linear Mixed Effect Model : 所在的 library 不同,也就是 developer 不一样。gls 和 lme 是 nlme 里面的。 : lmer 是 lme4 里面的。 : 我去年做一个数据分析的时候发现,同样的 data 和 model,在32位和64位版本的 : lmer/glmer 上竟然能出现很不同的结果。我的 model 是 logistic mixed effect : model。后来就不太敢用它。好在我需要做 non-linear mixed effect model 的时候也
|
B******5 发帖数: 4676 | 8 developer故意不给的,他觉得p-value不合理
【在 h******3 的大作中提到】 : 但是lmer的结果没有直接的p-value显示。这点很不方便。。
|
h******3 发帖数: 190 | 9 汗。为啥呢?
【在 B******5 的大作中提到】 : developer故意不给的,他觉得p-value不合理
|
B******5 发帖数: 4676 | |
w******4 发帖数: 488 | 11 非常感谢,Doug的解释真是精彩啊。。
【在 B******5 的大作中提到】 : https://stat.ethz.ch/pipermail/r-help/2006-May/094765.html : btw 你生日是1223?
|
o****n 发帖数: 281 | 12 一直在用SAS做LMM,大家是必须用R做么?结果不同怎么处理,需要和SAS的相应PROC的
结果比对么? |
h******3 发帖数: 190 | 13 虽然不是完全懂,但是多谢了~
之前都没想过lme output里的不同的df是怎么被算出来的,确实有些奇怪哈。
And about your 无厘头 question......, your hunch is right......
【在 B******5 的大作中提到】 : https://stat.ethz.ch/pipermail/r-help/2006-May/094765.html : btw 你生日是1223?
|
B******5 发帖数: 4676 | 14 哈哈,不错,回头跟你一起庆祝生日
【在 h******3 的大作中提到】 : 虽然不是完全懂,但是多谢了~ : 之前都没想过lme output里的不同的df是怎么被算出来的,确实有些奇怪哈。 : And about your 无厘头 question......, your hunch is right......
|