Блог

Cluster.robust.se кажется, изменяет тип данных iv регрессии и сокращает количество наблюдений с помощью stargazer, есть какие-нибудь советы по их возвращению?

Пример:

 library(ivpack)
library(stargazer)
df<- tibble(
  indiv_id = sample(10, 100, rep = TRUE),
  group_id = sample(10, 100, rep = TRUE),
  cluster= sample(2,100,rep=TRUE)
)
iv<-cluster.robust.se(ivreg(df$indiv_id~df$group_id),df$cluster)
stargazer(iv,
          type="text")
  

Что происходит, так это cluster.robust.se() изменение класса вашей модели с ivreg на coeftest . Проверьте это:

 > class(ivreg(df$indiv_id~df$group_id))
[1] "ivreg"

> class(iv)
[1] "coeftest"
  

Это два разных объекта, поэтому stargazer выдает разные результаты, исключая количество наблюдений из coeftest выходных данных.

Одним из возможных решений является включение обоих объектов в одну таблицу, например:

 iv<-cluster.robust.se(ivreg(df$indiv_id~df$group_id),df$cluster)
myivreg <- ivreg(df$indiv_id~df$group_id)

stargazer(myivreg, iv, type="text")