作者andrew43 (Myrmarachne)
看板Statistics
標題Re: [程式] 多種陽性率比較
時間Sun Apr 8 14:57:57 2012
※ 引述《ProfessorX (Xavier)》之銘言：
: ------------------------------------------------------------------------
: [軟體程式類別]:
: SAS 9.2
: [程式問題]:
: 資料處理
: [軟體熟悉度]:
: 低(1~3個月)
: [問題敘述]:
: 各位先進好，小弟最近遇到四筆資料(陽性 陰性)，製做成2x4表格後，
: 用fisher檢定可以推翻H0，但是只能證明有差異存在
: 請問我該如何讓四組的陽性率差兩兩檢定(Duncan)出有差異呢?
: [程式範例]:
: data Salmonella;
: input Result $ Space $ Head;
: datalines;
: (1)Donse (+) 66
: (1)Donse (-) 30
: (2)Budai (+) 4
: (2)Budai (-) 16
: (3)Yugon (+) 1
: (3)Yugon (-) 10
: (4)Other (+) 9
: (4)Other (-) 1
: ;
: proc FREQ DATA= Salmonella;
: tables Result * Space / norow nocol nopercent fisher;
: weight Head;
: 後面我就寫不下去了.....還請高手指點!!
考慮 saturated logistic models.
我不用 sas, 所以就用 R 代替了.
但是你可以參考
http://www.ats.ucla.edu/stat/sas/library/logistic.pdf
中的 Output. 5 那段例子.
group <- c("Donse", "Budai", "Yugon", "Other")
re.positive <- c(66, 4, 1, 9)
re.negative <- c(30, 16, 10, 1)
data <- data.frame(group, re.positive, re.negative, stringsAsFactors=T)
model <- glm(
cbind(re.positive,re.negative) ~ group, family=binomial, data=data
)
summary(model)
上述結果中會以 Budai 組為基準, 得到三個 ln(勝算比) 並檢驗是否為 0.
我猜想這是你要的.
要改變基準組, 把該基準組名的第一個字前加個 0 就好.
--
http://apansharing.blogspot.com/
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 122.117.37.172
※ 編輯: andrew43 來自: 122.117.37.172 (04/08 15:09)
1F：推 ProfessorX:謝謝高手回答! 不過R是什麼@@? 04/08 15:06
2F：→ andrew43:我有幫你找到sas寫法, 請再看一次文章. 04/08 15:15
3F：推 ProfessorX:有的!我正在看! 大感激~~~!! 04/08 15:19
4F：→ andrew43:不對, 你是x教授, 不用看用感應的就知道了. XD 04/08 23:07