Модель случайных эффектов в r

Question

Я хотел бы знать соответствующую команду для запуска модели со случайным эффектом под программным обеспечением R.

Я пробовал:

library(plm)
Model=plm(df$IndVar ~ df$densite + df$sexe + df$age_cat + df$production, data=df, index=c("village"), model="random")

Но я получил такой ответ:

Error in lm.fit(x, y, offset = offset, singular.ok = singular.ok, ...) : 
  0 (non-NA) cases
In addition: Warning message:
In Ops.factor(e1, e2) : ‘-’ not meaningful for factors

Тоже пробовал:

library(nlme)            
Model=lme(df$IndVar ~ df$densite + df$sexe + df$age_cat + df$production, random=~1|village)

И получил:

Error in eval(predvars, data, env) : object 'IndVar' not found

Что делать? Нужна ваша помощь.