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

该怎么办？需要您的帮助。