Как я могу найти модель из таблицы contigency?

Question

У меня есть таблица смежности с 3 путями, скажем, A, B, C, все двоичные данные (1 или 0) и несколько различных чисел в таблице;например:

                                  A
                            1           0
                            D           D
                          1   0       1   0

           1    C    1   21   15     11    30
     B               0   15   11      5    8
           0    C    1    5   42     3     14
                     0    4    3      2    15

D - это успех или неудача.

Я хочу использовать r, чтобы найти модель успеха, скажем, prob (D = 1) = a + bx1 + cx2+ dx3, x1, x2, x3 - переменные A, B, C.

Так что же я могу сделать?

Я пробовал model = glm (D ~ A + B + C, бином)

но не работает.

Answer 1

Вы должны преобразовать свою таблицу непредвиденных расходов (я полагаю, это действительно то, что у вас есть) в data.frame.

# Sample data
n <- 200
d <- data.frame(
  A = sample(0:1, n, replace=TRUE),
  B = sample(0:1, n, replace=TRUE),
  C = sample(0:1, n, replace=TRUE),
  D = sample(0:1, n, replace=TRUE)
)
d <- table(d)

# Put the data in a data.frame
library(reshape)
d <- melt(d)

# Compute the regression
glm(D ~ A+B+C, data=d, family=binomial(), weights=value)