R: plm - jährliche Fixeffekte - Jahres- und Quartalsdaten
Ich habe ein Problem beim Einrichten eines Paneldatenmodells.
Hier sind einige Beispieldaten:
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ühre die folgende Regression mit 'id' als Einzelindex und 'year' als Zeitindex durch:
reg1 <- plm(y ~ x, data=data,index=c("id", "year"), model="within",effect="time")
Leider erhalte ich folgenden Fehler:
doppelte Paare (Zeit-ID) Fehler in pdim.default (Index [[1]], Index [[2]]):
Um das zu umgehen, benutze ich die kombinierte Variable 'y_q':
reg1 <- plm(y ~ x, data=data,index=c("id", "y_q"), model="within",effect="time")
Aber hier ist mein Problem - ich möchte nur festgelegte Effekte für das Jahr und nicht für das Quartal.
Gibt es eine andere Möglichkeit, die frühere Ausgabe zu umgehen, anstatt den Index 'y_q' zu erstellen?
Vielen Dank im Voraus für jede Hilfe!