作者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