Come calcolare gli errori standard di boot bootsrapped con R?

Question

Recentemente ho cambiato da STATA a R e in qualche modo fatica a trovare alcuni comandi corrispondenti. Vorrei ottenere pannello bootsrapped errori standard da un modello con effetti fissi utilizzando la libreria plm come descritto qui here per gli utenti STATA:

1. Le mie domande riguardano l'approccio in generale (sia di avvio è la libreria appropriata o la libreria (meboot) )

2. come risolvere per quel particolare errore utilizzando boot:

prima ottenere alcuni dati panel:

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")) 

Secondo:

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

Per qualche ragione, boot passerà un indice (originale qui) per plm con i valori duplicati. È necessario rimuovere tutti i valori duplicati e affermare che l'indice è univoco prima di passarlo a plm.

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