作者cbt (持续努力中~)
看板Statistics
标题[程式] R的指令glmmPQL
时间Tue Aug 14 08:17:52 2012
[软体程式类别]: R
[程式问题]: 找不到如何取出里面的值
资料处理、回归、叙述统计、logistic、probit...等
[软体熟悉度]: 中~高
[问题叙述]:
[程式范例]:
我使用R的library(MASS)中的glmmPQL
我可以顺利的跑出我想要的程式
举例来说
library(MASS)
mydata=data.frame(x1=rnorm(10),x2=rnorm(10),y=rgamma(10,1,1),ID=seq(1:10))
glmmPQL(y ~ x1+x2 , random = ~ 1 | ID, family = Gamma(link="log"), data = mydata)
其中一次的结果如下
iteration 1
iteration 2
iteration 3
iteration 4
iteration 5
iteration 6
iteration 7
Linear mixed-effects model fit by maximum likelihood
Data: mydata
Log-likelihood: NA
Fixed: y ~ x1 + x2
(Intercept) x1 x2
-1.578571696 0.968307141 -0.001152464
Random effects:
Formula: ~1 | ID
(Intercept) Residual
StdDev: 0.8865233 0.3324462
Variance function:
Structure: fixed weights
Formula: ~invwt
Number of Observations: 10
Number of Groups: 10
程式可以成功的跑出来 可是我有几个问题
1.在output中哪里可以找到gamma分配的scale parameter的估计值
2.因为我之後要做模拟跑一千次 所以我需要读出里面的几个值
尤其是random effect的StdDev 0.8865233
我知道如果我用"...$sigma"得到的是0.3324462
3.可以要求在跑的过程中不要print iteration 1,iteration 2,...等等吗
我卡在这几个问题好久了 尤其1&2
请板上的各位高手帮帮忙指点一下 非常感谢!!!
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 165.91.10.164
1F：推 blackalert:result = suppressMessages(glmmPQL(y~x1+x2, ...)) 08/14 17:17
2F：→ blackalert:intervals(result)$r$I[2] 08/14 17:17
3F：→ blackalert:第一个问题我不会，等高手回答 XD 08/14 17:19
4F：→ cbt:非常感谢楼上的回答,也麻烦其他高手帮帮我回答我的问题1,谢谢! 08/14 23:37
5F：→ ljta:1.attributes(objects) 08/15 10:59
6F：→ cbt:我使用attribute(objects) 可是还是没看到哪里有scale paramet 08/22 04:16
7F：→ cbt:er 可以请高手再指点一下吗? 谢谢! 08/22 04:16