R: избежать summary.plm

Я использую R для запуска моделирования Монте-Карло, изучающего производительность оценщиков панельных данных. Поскольку я буду проводить большое количество испытаний, мне нужно получить от моего кода хотя бы приличную производительность.

Использование Rprof в 10 испытаниях моей симуляции показывает, что значительная часть времени тратится в звонках на summary.plm . Первые несколько строк Rprofsummary приведены ниже:

$by.total
                            total.time total.pct self.time self.pct
"trial"                          54.48     100.0      0.00      0.0
"coefs"                          53.90      98.9      0.06      0.1
"model.matrix"                   36.72      67.4      0.10      0.2
"model.matrix.pFormula"          35.98      66.0      0.06      0.1
"summary"                        33.82      62.1      0.00      0.0
"summary.plm"                    33.80      62.0      0.08      0.1
"r.squared"                      29.00      53.2      0.02      0.0
"FUN"                            24.84      45.6      7.52     13.8

Я вызываю summary в своем коде, потому что мне нужно получить стандартные ошибки оценок коэффициентов, а также коэффициентов сами (которые я мог получить только из объекта plm). Мой вызов выглядит как

regression <- plm(g ~ y0 + Xit, data=panel_data, model=model, index=c("country","period"))

coefficients_estimated <- summary(regression)$coefficients[,"Estimate"]
ses_estimated <- summary(regression)$coefficients[,"Std. Error"]

. У меня неприятное ощущение, что это огромная трата времени процессора, но я недостаточно знаю о том, как R делает что-то, чтобы избежать вызова summary. Я' Буду признателен за любую информацию о том, что здесь происходит за кулисами, или за какой-либо способ сократить время, необходимое для этого.

7
задан landroni 15 August 2014 в 07:15
поделиться