plm 固定效應錯誤 - R 中的空模型 - 不平衡數據

Question

我正在嘗試使用 R 進行固定效應模型

fixed1<-plm(ITEMS_STARPU ~ TREAT_SUM, data=dataHyp1.i, index=c("PRACTICE"), model = "within")

我不斷收到以下錯誤：

Error in plm.fit(formula, data, model, effect, random.method, random.models,  : empty model

結果變量是每次全科醫生診所（即診所）配發的抗生素數量，因此該指數是在實踐水平上（與患者相反）

當我使用“池化”時數據有效

我讀過之前的帖子，指出問題可能出在時間段上，但數據已經包含一個日期元素（數據是每月收集和顯示的）

Hyp1.i <- read.csv("Hypothesis 1.i.csv",  header = TRUE, na.strings = c("NA"), stringsAsFactors = FALSE)

myvarsHyp1.i<- c("PERIOD_NO", "PRACTICE", "REGIONAL_TEAM", "ITEMS_STARPU", "TREAT_SUM")
dataHyp1.i <- Hyp1.i[,myvarsHyp1.i]

table(dataHyp1.i$PERIOD_NO)
class(dataHyp1.i$PERIOD_NO)
dataHyp1.i$PERIOD_NO<-factor(dataHyp1.i$PERIOD_NO)

table(dataHyp1.i$TREAT_SUM)
class(dataHyp1.i$TREAT_SUM)

dataHyp1.i$TREAT_SUM <- factor(dataHyp1.i$TREAT_SUM, levels=c(0:2), labels=c("No Letter", "Single Letter", "Repeat Letter"))

table(dataHyp1.i$TREAT_SUM)

fixed1<-plm(ITEMS_STARPU ~ TREAT_SUM, data=dataHyp1.i, index=c("PRACTICE"), model = "within")
summary(fixed)

Answer 1

您使用的是面板數據，因此您應該同時擁有橫截面數據和時間序列。 所以我認為問題在於index=c("PRACTICE")那部分應該是您的橫截面 ID 和時間 ID 的標題，在您的情況下它是 PERIOD_NO。 所以我認為你應該首先寫index=c("PRACTICE","PERIOD_NO")並且你的面板數據可能不再是空的了