在R语言中使用循环来进行卡方检验

Question

我是R的新手。我找到了以下用于对一组变量进行单变量逻辑回归的代码。我想做的是运行卡方检验，将一系列变量与因变量进行比较，类似于下面的逻辑回归代码。我发现了几个涉及创建所有可能变量组合的方法，但我无法使其正常工作。理想情况下，我希望其中一个变量（X）保持不变。 使用for循环进行R中的卡方分析

lapply(c("age","sex","race","service","cancer",
         "renal","inf","cpr","sys","heart","prevad",
         "type","frac","po2","ph","pco2","bic","cre","loc"),

       function(var) {

         formula    <- as.formula(paste("status ~", var))
         res.logist <- glm(formula, data = icu, family = binomial)

         summary(res.logist)
       })

- dj_stat

---

这里的信息不足以回答这个问题。你能提供一些数据样本吗？或者制作一个可重现的例子（最好使用 set.seed()）？status，age等是你的变量吗？具体哪里出了问题？ - gung - Reinstate Monica

可能是使用R中的for循环进行卡方分析的重复问题。 - C8H10N4O2

Answer 1

您确定lapply的向量中的字符串是否在icu数据集的列名中吗？

当我下载icu数据时，它可以正常工作：

system("wget http://course1.winona.edu/bdeppa/Biostatistics/Data%20Sets/ICU.TXT")
icu <- read.table('ICU.TXT', header=TRUE)

并将 status 更改为 STA，这是 icu 中的一列。以下是您某些变量的示例：

my.list <- lapply(c("Age","Sex","Race","Ser","Can"),         
       function(var) {
         formula    <- as.formula(paste("STA ~", var))
         res.logist <- glm(formula, data = icu, family = binomial)
         summary(res.logist)
       })

这将给我一个summary.glm对象列表。例如：

lapply(my.list, coefficients)
[[1]]
               Estimate Std. Error   z value     Pr(>|z|)
(Intercept) -3.05851323 0.69608124 -4.393903 1.113337e-05
Age          0.02754261 0.01056416  2.607174 9.129303e-03

[[2]]
              Estimate Std. Error    z value     Pr(>|z|)
(Intercept) -1.4271164  0.2273030 -6.2784758 3.419081e-10
Sex          0.1053605  0.3617088  0.2912855 7.708330e-01

[[3]]
              Estimate Std. Error    z value   Pr(>|z|)
(Intercept) -1.0500583  0.4983146 -2.1072198 0.03509853
Race        -0.2913384  0.4108026 -0.7091933 0.47820450

[[4]]
              Estimate Std. Error   z value     Pr(>|z|)
(Intercept) -0.9465961  0.2310559 -4.096827 0.0000418852
Ser         -0.9469461  0.3681954 -2.571858 0.0101154495

[[5]]
                 Estimate Std. Error       z value     Pr(>|z|)
(Intercept) -1.386294e+00  0.1863390 -7.439638e+00 1.009615e-13
Can          7.523358e-16  0.5892555  1.276756e-15 1.000000e+00

如果你想进行卡方检验：

my.list <- lapply(c("Age","Sex","Race","Ser","Can"),function(var)chisq.test(icu$STA, icu[,var]))

或对所有变量组合进行卡方检验:

my.list.all <- apply(combn(colnames(icu), 2), 2, function(x)chisq.test(icu[,x[1]], icu[,x[2]]))

这个可以用吗？

- user1981275

---

嘿，谢谢，回归代码对我有用。我需要类似的东西来运行卡方检验。 - dj_stat