作者andrew43 (Myrmarachne)
看板Statistics
标题[讨论] 广义线性模型的多重比较
时间Thu Apr 21 20:52:36 2011
各位先进好.
前些日子向各位求教於广义线性模型之回归分析,
各位的帮助对我的工作确实有很大的帮助,
非常感谢大家的热心.
今日小弟遇到另一个问题,
是关於广义线模的 "多重比较",
请见下例.
# R code
# 定义依变数
y <- rpois(40,5)
# 定义自变数为 4 个 level 之类别变数
x <- gl(4 , 10 , labels=c('A','B','C','D'))
# 定义自变数为 4 个 level 之类别变数
# 进行 Poisson regression
m <- glm(y~x , family=poisson) ; summary(m)
可能的结果如下:
Call:
glm(formula = y ~ x, family = poisson)
Deviance Residuals:
Min 1Q Median 3Q Max
-1.8988 -0.4599 0.0000 0.4450 1.4888
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) 1.79176 0.12910 13.879 <2e-16 ***
xB -0.18232 0.19149 -0.952 0.341
xC 0.01653 0.18182 0.091 0.928
xD -0.20252 0.19255 -1.052 0.293
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for poisson family taken to be 1)
Null deviance: 26.017 on 39 degrees of freedom
Residual deviance: 23.795 on 36 degrees of freedom
AIC: 172.30
Number of Fisher Scoring iterations: 4
可得 xB, xC 及 xD 之参数检定皆以 xA 为基准点,
例如 xB 之发生率为 xA 之 e^(-0.18232) 倍,
但未达显着差异 (p = 0.341);
若欲求得 xD 对比 xC 之发生率, 则取 e^[(-0.20252)-(0.01653)],
但接下来的检验就不知如何进行了.
请问我可以怎麽做才能得到 xD vs xC 之检定结果?
甚至我如何求得这 4 组资料所有配对之比较?
重新强制地设定基准点为 xB, xC 及 xD 再重做几次的方式是洽当的吗?
以上敬请前辈们指点.
谢谢.
--
http://apansharing.blogspot.com/
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.128.117.73