Как использовать lmer для анализа линейной регрессии в R - PullRequest
Новая
       16

Как использовать lmer для анализа линейной регрессии в R

0 голосов
/ 27 февраля 2020

У меня фрейм данных выглядит следующим образом: 

order <- c('1','2','3','4')  #integer
sample <- c('A01','A02','A03','A04')  #factor
Dose <- c('100','200','500','1000')   #factor
Dose.rate <- c('chronic','acute','acute','chronic')   #factor
Time.point <- c('3h','3w','3h','3w')   #factor
ctrl.irr <- c('ctrl','ctrl','irr','irr')   #factor
treatment <- c('0','0','1','1') #factor
CD45 <- c('23.45','45.13','39.93','21.35')   #numeric
Replicate <- c('1','4','3','7')   #factor
df <- data.frame(order,sample,Dose,Dose.rate,Time.point,ctrl.irr,treatment,CD45,Replicate)

Я хочу провести множественный линейный регрессионный анализ среди разных столбцов с учетом столбца CD45 в качестве основного столбца для сравнения. Как и между CD45 column и Dose, а также другими столбцами (парными).

У меня была командная строка для использования, но, к сожалению, неверный вывод. 

model9.1 <- lmer( of.CD45 ~ treatment+ Dose+ Dose.rate+ Time.point+
                   (treatment+ Dose.rate+ Time.point)*Dose+
                   (Dose+ Dose.rate+ Time.point)*treatment+
                  (Dose+ treatment+ Time.point)*Dose.rate+
                   (Dose+ treatment+ Dose.rate)* Time.point+
                   (Dose*treatment*Dose.rate*Time.point)+
                   (1| Replicate), data = df)

Есть мой код достаточно правильный или как я могу изменить его, чтобы получить правильный вывод?

...