Я пытаюсь провести логистический анализ c. В этом анализе данные выглядят следующим образом:
Y: двоичный (0,1)
А для независимых переменных они следующие:
X1 : DrugA: категориальная переменная. Принимали ли пациенты препарат A (да или нет)
X2: DrugA_Con c: непрерывная переменная. концентрации Drug A , который также включает ноль.
X3: DrugB: категориальная переменная. Принимали ли пациенты препарат B (да или нет).
Обратите внимание, что во входном наборе данных пациенты относятся к одной из 3 групп:
- Не сделал принимать либо препарат А, либо препарат В (плацебо)
- принимал препарат А, но НЕ препарат Б (только препарат А)
- принимал препарат А и препарат В (комбинация)
Пожалуйста, смотрите код ниже, который генерирует представление набора данных независимых переменных:
DrugA<-data.frame(DrugA=factor(1))
DrugA_Conc<-data.frame(DrugA_Conc=rep(seq(from = 0, to = 100, length.out = 3)))
DrugB<-data.frame(DrugB=factor(c(0,1)))
mergd<-merge(DrugA,DrugA_Conc)
mergd2<-merge(mergd,DrugB)
DrugA_2<-data.frame(DrugA=factor(0))
DrugA_Conc_2<-data.frame(DrugA_Conc=rep(0,length.out = 3))
DrugB_2<-data.frame(DrugB=factor(0))
mergd3<-merge(DrugA_2,DrugA_Conc_2)
mergd4<-merge(mergd3,DrugB_2)
mergd5<-rbind(mergd2,mergd4)
Мой вопрос заключается в следующем: Как я могу кодировать модель logisti c в R, чтобы предсказать эффект DrugA_Con c на Y с присутствием DrugB и без него, и прогнозируемый выходной сигнал для 4 возможных сценариев ios (3, перечисленных выше) +
НЕ принимал препарат А, но принимал препарат В (препарат В один).
Я пробовал следующие коды в R, но он не работал:
m1 <- glm(Y~ DrugB+DrugA*DrugA_Conc, data=all, family="binomial")
m2 <- glm(Y~ DrugB+I(DrugA*DrugA_Conc), data=all, family="binomial")
Обратите внимание, что я не пытаюсь оценить взаимодействие между DrugA и DrugA_Con c, но скорее, пытаясь создать прогнозы для 4 возможных сценариев комбинации лекарств ios (один препарат A, оба препарата, ни один из препаратов, один препарат B) на основе набора данных, который имеет первые 3 сценария ios. Также можно ли применить код, если я возьму журнал DrugA_Con c (я имею в виду log(DrugA_Conc))
В SAS я обнаружил, что для кодирования можно применить трюк (умножение DrugA*DrugA_Conc) модель, подобная следующей, которая позволяет прогнозировать 4 сценария ios (рисунок 1):
model Y(event='1') = DrugB DrugA*DrugA_Conc/
С уважением,