2017-07-15 96 views
0

我有一個數據幀(「數學」)這樣的時間沒有偏差只有變量時空模型(有三種不同的方法,儘管僅示出一個) - dataframePLM包中的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,隨機效應模型如上所述。

+1

請參閱此鏈接以瞭解如何製作可重現的示例:https://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example – www

+0

嘗試創建一個「pdata .frame「,並使用它的'data'參數。看包裝的vingette如何做到這一點。 – Helix123

回答

1

這可能是您的數據的問題: 看起來,變量VerbalScoreMethod不會因人而異。因此,對於Swamy-Arora RE模型(默認),無法計算必需的方差。受影響的變量退出模型,這裏是所有RHS變量,並且您得到(不是非常具體的)錯誤消息empty model

您可以使用命令pvar()檢查每個人的變化。

如果我的假設是真實的,並且仍然想要估計隨機效應模型,則必須使用不依賴於方差內的不同隨機效應估計值,例如,嘗試Wallace-Hussain估計量(random.method="walhus")。