s*********t 发帖数: 3 | 1 有个关于Generalized Linear Mixed Models(GLMMs)的问题:
请问在R或SAS中能否得到 covariance matrix between fixed effect estimators (\
beta) and variance component estimator for random effect(\sigma)
在R或SAS中可以得到fixed effect estimators的variance matrix,也可以得到
variance component estimator(\sigma)的variance。但是能否输出它们两者的
covariance呢?
多谢~ |
w*****n 发帖数: 375 | 2 fixed effect和 random effect之间没有covariance. 如果你想知道的是intreaction
的话, 把interacation term 放到model 里 |
s*********t 发帖数: 3 | 3 在linear mixed model里,fixed effect和 random effect之间的确没有covariance。
(因为Fisher Information Matrix的off-diagonal block matrix=0) 但是在
Generalized Linear Mixed Models里,fixed effect和 random effect之间是有
covariance。(因为Fisher Information Matrix的off-diagonal block matrix !=0) |
s*r 发帖数: 2757 | 4 你有没有看过proc glimmix的manual
【在 s*********t 的大作中提到】 : 在linear mixed model里,fixed effect和 random effect之间的确没有covariance。 : (因为Fisher Information Matrix的off-diagonal block matrix=0) 但是在 : Generalized Linear Mixed Models里,fixed effect和 random effect之间是有 : covariance。(因为Fisher Information Matrix的off-diagonal block matrix !=0)
|
s*********t 发帖数: 3 | 5 看了proc mixed的manual, 没有相关输出项。proc glimmix倒没看。
我目前的做法是在R中写出loglikelihood function,然后用nlm,可以得到Hessian
matrix,再求逆;
多谢大家的回复~
【在 s*r 的大作中提到】 : 你有没有看过proc glimmix的manual
|