R: plm - rok stałych efektów - dane za rok i kwartał
Mam problem z ustawieniem modelu danych panelu.
Oto kilka przykładowych danych:
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)
Używam następującej regresji, używając „id” jako indeksu indywidualnego i „roku” jako indeksu czasu:
reg1 <- plm(y ~ x, data=data,index=c("id", "year"), model="within",effect="time")
Niestety pojawia się następujący błąd:
duplikat par (identyfikator czasu) Błąd w pdim.default (indeks [[1]], indeks [[2]]):
Aby to obejść, używam połączonej zmiennej „y_q”:
reg1 <- plm(y ~ x, data=data,index=c("id", "y_q"), model="within",effect="time")
Ale oto mój problem - chcę mieć tylko stałe efekty roku, a nie kwartał.
Czy istnieje inny sposób obejścia wcześniejszego problemu niż tworzenie indeksu tiemu „y_q”?
Dzięki z góry za każdą pomoc!