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