R: plm — Jahres-fixed effects Jahres — und Quartal-Daten

Ich habe ein problem in der Einrichtung eines panel-Daten-Modell.

Ist hier einige Beispiel-Daten:

library(plm)

id <- c(1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2)
year <- c(1999,1999,1999,1999,2000,2000,2000,2000,1999,1999,1999,1999,2000,2000,2000,2000)
qtr <- c(1,2,3,4,1,2,3,4,1,2,3,4,1,2,3,4)
y <- rnorm(16, mean=0, sd=1)
x <- rnorm(16, mean=0, sd=1)

data <- data.frame(id=id,year=year,qtr=qtr,y_q=paste(year,qtr,sep="_"),y=y,x=x)

Ich führen Sie die folgende regression mit 'id' als dem jeweiligen index und "Jahr" als die Zeit, index:

reg1 <- plm(y ~ x, data=data,index=c("id", "year"), model="within",effect="time")

Leider bekomme ich die folgende Fehlermeldung:

doppelten Paare (Zeit-id) Fehler in pdim.default(index[[1]], index[[2]]) :

So, das zu umgehen, verwende ich die kombinierte variable, die 'y_q':

reg1 <- plm(y ~ x, data=data,index=c("id", "y_q"), model="within",effect="time")

Aber hier ist mein Problem-ich will nur Jahr fixe Effekte und nicht Jahr-Quartal.

Gibt es einen anderen Weg, um das frühere Problem, anstatt tiem index 'y_q'?

Dank im Voraus für jede Hilfe!!!

InformationsquelleAutor Brad | 2013-04-11
  1. 3

    In einer panel-Einstellung Sie in der Regel don ' T haben einige doppelte Wert für jedes paar-id-Jahr.

    In Ihrem quartalsweise Daten wird es schwierig sein, zu berechnen, jährlich fixed-Effekt-Modelle ohne Aggregation Ihrer Daten, um Sie jährlich.

    Überprüfen Sie die Beispiele hier zu sehen, wie Ihre Daten formatiert werden sollen für panel-Daten-Modellierung.

    Hier ist oneway zu tun : 

    require(plyr)
yeardata  <- ddply(data, .(year, id), summarize, y = mean(y),
                                                 x = mean(x))


require(plm)
reg1 <- plm(y ~ x, data = yeardata, index = c("id", "year"), model = "within", effect = "time")
fixef(reg1)

##      1999      2000 
## 0.2641997 0.0041193
    InformationsquelleAutor dickoa
Schreibe einen Kommentar