3

私はRを初めて使用します。一連の変数に対して一変量ロジスティック回帰を実行するための次のコードを見つけました。私がやりたいことは、以下のロジスティック回帰コードと同様に、従属変数に対する変数のリストに対してカイ二乗検定を実行することです。変数の可能なすべての組み合わせを作成することを含むそれらのいくつかを見つけましたが、それを機能させることができません。理想的には、変数の 1 つ (X) が同じであることを望みます。

R の for ループを使用したカイ二乗分析

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)
       })
4

1 に答える 1

2

ベクター内の文字列がデータセットlapplyの列名に含まれていることを確認してください。icu

データをダウンロードするとうまくいきますicu：

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

をの列に変更statusします。いくつかの変数の例を次に示します。STAicu

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]))

または、変数のすべての組み合わせに対するカイ 2 乗検定:

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

これは機能しますか？

于 2013-09-17T13:15:28.557 に答える