或使用r与预期的不同

问题描述 投票:0回答:1
我的2x2表:

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使用此命令获得了预期的输出:
logistic-regression glm
1个回答
0
投票
ria no中的默认值,但它在这里没有工作,我不知道为什么。我必须使用上述COCE指定参考。关于为什么R不使用“否”的任何答案,因为默认值将被赞赏。 thanks

<- relevel(Thyroid$cancer, ref = "No")``

最新问题
© www.soinside.com 2019 - 2025. All rights reserved.