Thyroid cancer
Race Yes No
White 92 21
Black 15 20
与黑色相比,要获得白色的OR癌症:
(92/21) / (15/20) = 5.84127
当我使用R时,我会得到相反的情况:或黑色的5.84127与白色相比!?
贝洛是R代码:
Thyroid \<- as.table(rbind(c(92,21),c(15, 20)))
dimnames(Thyroid ) \<- list(
race=c("white","black"),
cancer=c("yes","no")
)
Thyroid \<- as.data.frame((Thyroid))`
Thyroid$race \<- relevel(Thyroid$race, ref = "black")`
glm.thy \<- glm(cancer \~ race, weights = Freq,`
family = binomial(link = 'logit'),`
data = Thyroid)`
summary(glm.thy)`
exp(coef(glm.thy))
这就是我得到的:
(Intercept) racewhite
1.3333333 0.1711957
或= 0.1711957?任何想法都是相关的。 谢谢我尝试了R中的GLM,我预计将插入0.1711957的5.84127。
我还尝试了Fisher测试,该测试提供了合理的或:
Thyroid <- as.table(rbind(c(92,21),c(15, 20)))
dimnames(Thyroid ) <- list(
race=c("white","black"),
cancer=c("yes","no"))
print(Thyroid )
fisher.test(Thyroid )
I使用此命令获得了预期的输出:
<- relevel(Thyroid$cancer, ref = "No")``