tappek

library("plm")
data("Grunfeld", package = "plm")

## set some observations of the dependent variable to missing
Grunfeld$inv[Grunfeld$year<=1936] <- NA
fe <- plm(inv ~ value + capital, data = Grunfeld, model = "within")

library(sandwich)
sandwich::vcovBS(fe, cluster=~firm, R = 999) # errors
# Error in vcovBS.default(fe, cluster = ~firm, R = 999) : 

tappek

# speed test matrix splitting by individual subsetting, esp. in 
# plm::pcce and plm::pmg

library("plm")
data(Produc)
pProduc <- pdata.frame(Produc)

X <- as.matrix(pProduc[ , -c(1:3)])
ind <- index(pProduc)[[1L]]