проп.тест, кондиционирование на другой класс - PullRequest
Новая
       13

проп.тест, кондиционирование на другой класс

0 голосов
/ 19 февраля 2019

Я запускаю prop.test на R. Вот таблица непредвиденных расходов и выходные данные prop.test.«а» - это имя фрейма данных.'cluster' и 'AD_2' - две дихотомические переменные. 

 table(a$cluster,a$AD_2)

         no  yes
   neg 1227  375
   pos  546  292

 prop.test(table(a$cluster,a$AD_2))

    2-sample test for equality of proportions with continuity 
 correction

data:  table(a$cluster, a$AD_2)
X-squared = 35.656, df = 1, p-value = 2.355e-09
alternative hypothesis: two.sided
95 percent confidence interval:
0.07510846 0.15362412
sample estimates:
   prop 1    prop 2 
0.7659176 0.6515513

Выборочные оценки обусловлены тем, что AD_2 является "нет", как видно из таблицы сопряженности, то есть 0,7659176 = 1227 / (1227 + 375)и 0,6515513 = 546 / (546 + 292).Поскольку $ cluster == pos положительное событие и AD_2 == да фактор риска, я бы хотел изменить соотношение пропорций на AD_2, равное 'yes'.

1 Ответ

0 голосов
/ 20 февраля 2019

R таблицы по сути являются матрицами.Функция prop.test может обрабатывать матрицы, поэтому используйте те же данные с переключенными столбцами: 

> prop.test( matrix(c( 375,292,1227,546), 2))

    2-sample test for equality of proportions with continuity correction

data:  matrix(c(375, 292, 1227, 546), 2)
X-squared = 35.656, df = 1, p-value = 2.355e-09
alternative hypothesis: two.sided
95 percent confidence interval:
 -0.15362412 -0.07510846
sample estimates:
   prop 1    prop 2 
0.2340824 0.3484487

Я думаю, что другой метод мог бы заключаться в обмене столбцами с: 

table(a$cluster,a$AD_2)[ , 2:1]
...