Исправлена ​​модель эффекта на основе заданной формулы с использованием lme4 из R - PullRequest
Новая
       11

Исправлена ​​модель эффекта на основе заданной формулы с использованием lme4 из R

1 голос
/ 01 апреля 2020

С учетом следующей формулы: 

y = β0 + β1·(x1 - 30) + β2·(x2 - 3) + β3·(x3 - 10)+ β4·x4 + β5·x6 + β6·x6·(x2-3) + e

Как я могу построить модель в R с использованием пакета, такого как lme4, et c? Заранее спасибо.

Моя пробная версия пока, но это явно не модель, как требуется формула. 

library(lme4)
# model <- lmer((y ~ (x1 -30) + (x2 - 3) + (x3-10) + x4 + x5 + (x6 | (x2-3))), data = df, REML = FALSE) this line not work
model <- lmer(y ~ x1 + (1|x2),  data = df, REML = FALSE)
model

Out: 

Linear mixed model fit by maximum likelihood [lmerMod
]
Formula: y ~ x1 + (1 | x2)
   Data: df
      AIC       BIC    logLik  deviance  df.resid 
 22427.02  22452.94 -11209.51  22419.02      4821 
Random effects:
 Groups   Name        Std.Dev.
 cal      (Intercept) 0.3649  
 Residual             2.4430  
Number of obs: 4825, groups:  cal, 4785
Fixed Effects:
(Intercept)          x1  
   19.68828      0.05818
...