R: le plm année les effets fixes de l'année et le trimestre de données

Je vais avoir un problème de configuration d'un modèle de données de panel.

Voici quelques exemples de données:

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)

- Je exécuter la régression suivante à l'aide de " id "comme l'indice et "année" comme l'indice de temps:

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

Malheureusement, j'obtiens l'erreur suivante:

dupliquer les couples (id) Erreur dans la pdim.par défaut(index[[1]], index[[2]]) :

De sorte à obtenir autour de cela, j'utilise le combiné variable qui est "y_q':

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

Mais voici mon problème, je ne veux avoir les effets fixes de l'année et pas l'année-trimestre.

Est-il une autre façon de se déplacer le plus tôt question au lieu de faire de l'introduction de la gestion de l'indice de y_q'?

Merci d'avance pour toute aide!

voir aussi stackoverflow.com/questions/28789767/...

OriginalL'auteur Brad | 2013-04-11

  1. 3

    Dans un panneau de configuration, vous n'avez généralement pas à avoir une certaine valeur en double pour chaque couple de l'id de l'année.

    Dans votre trimestrielle de données, il sera difficile de calculer un an modèles à effets fixes sans agréger vos données afin de les rendre annuellement.

    Vérifier les exemples ici pour voir comment les données doivent être formatées pour panneau de modélisation de données.

    Ici est oneway pour le faire : 

    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

    OriginalL'auteur dickoa