PLM包中的R - 包括对每个个体

我有一个数据帧（“数学”）这样的时间没有偏差只有变量时空模型（有三种不同的方法，尽管仅示出一个） - dataframe PLM包中的R - 包括对每个个体

我正在尝试为MathScore创建一个多级增长模型，其中VerbalScore是一个独立的，时间不变的随机效应。

我相信将R代码应与此类似 -

random <- plm(MathScore ~ VerbalScore + Method, data=math, index=c("id","Semester"), 
       model="random")

但是，运行以下错误这段代码的结果：

Error in plm.fit(object, data, model = "within", effect = effect) :
empty model

我相信这是与指数的问题，因为代码将运行，如果我使用：

random <- plm(MathScore ~ VerbalScore + Method + Semester, data=math, index="id", 
       model="random")

我将不胜感激关于如何创建一个多文件vel，随机效应模型如上所述。

2017-07-15 Eric

请参阅此链接以了解如何制作可重现的示例：https://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example – www

尝试创建一个“pdata .frame“，并使用它的'data'参数。看包装的vingette如何做到这一点。 – Helix123

这可能是您的数据的问题：看起来，变量VerbalScore和Method不会因人而异。因此，对于Swamy-Arora RE模型（默认），无法计算必需的方差。受影响的变量退出模型，这里是所有RHS变量，并且您得到（不是非常具体的）错误消息empty model。

您可以使用命令pvar()检查每个人的变化。

如果我的假设是真实的，并且仍然想要估计随机效应模型，则必须使用不依赖于方差内的不同随机效应估计值，例如，尝试Wallace-Hussain估计量（random.method="walhus"）。

2017-07-16 09:42:11 Helix123

回答