Wie berechnet Panel bootsrapped Standardfehler mit R?

Question

Ich habe vor kurzem von STATA zu R gewechselt und habe irgendwie Mühe, entsprechende Befehle zu finden. Ich mag Panel Standardfehler von einem Fest Effect Modell wie hier beschrieben here für STATA Benutzer mit der plm Bibliothek bootsrapped erhalten:

1. Meine Fragen betreffen den Ansatz im Allgemeinen (ob das Boot ist die entsprechende Bibliothek oder die Bibliothek (meboot) )

2. Wie für diesen speziellen Fehler lösen Boot mit:

zunächst einige Panel-Daten erhalten:

library(plm) 
data(EmplUK) # from plm library 

test<-function(data, i) coef(plm(wage~emp+sector,data = data[i,], 
             index=c("firm","year"),model="within")) 

Zweitens:

library(boot) 
boot<-boot(EmplUK, test, R = 100) 

> boot<-boot(EmplUK, test, R = 100) 
    duplicate couples (time-id) 
    Error in pdim.default(index[[1]], index[[2]]) : 
    Called from: top level 
    Browse[1]> 

Answer 1

Aus irgendeinem Grund boot einen Index (Original hier) plm mit duplizierten Werte übergeben. Sie sollten alle duplizierten Werte entfernen und bestätigen, dass der Index einmalig ist, bevor Sie ihn an plm übergeben.

test <- function(data,original) { 
    coef(plm(wage~emp+sector,data = data[unique(original),], 
     index=c("firm","year"),model="within")) 
} 


boot(EmplUK, test, R = 100) 

## ORDINARY NONPARAMETRIC BOOTSTRAP 
## Call: 
## boot(data = EmplUK, statistic = test, R = 100) 
## Bootstrap Statistics : 
##  original  bias std. error 
## t1* -0.1198127 -0.01255009 0.05269375