Использование texreg (или аналогичного) для создания таблицы регрессии для вывода кластерной стандартной функции ошибок (cl)

Question

Я использовал отличный пакет texreg производить готовые к LaTeX таблицы регрессии из подобранных объектов модели, но, похоже, она не совместима с различными функциями, которые корректируют мои стандартные ошибки для кластеризации. Некоторые поддельные данные и код приводят пример и сообщение об ошибке ниже.

Любые мысли о том, как получить желаемый результат (аналогично тому, что я получаю от texreg)?

x  = rnorm(1000)
IDs = ceiling(seq(from = .1, to = 10,length.out=1000))
s = c(rep(rnorm(1),100),rep(rnorm(1),100),rep(rnorm(1),100),rep(rnorm(1),100),rep(rnorm(1),100),rep(rnorm(1),100),rep(rnorm(1),100),rep(rnorm(1),100),rep(rnorm(1),100),rep(rnorm(1),100))
y = 3*x + 1.5^(IDs)*rnorm(n=1000,sd=s^2)*x + rnorm(1000)*.3
IDs = as.factor(IDs)
d = data.frame(x,IDs,y)
m = lm(y~IDs+x,data=d)
summary(m)

library(texreg)
texreg(m,omit.coef="IDs")

\begin{table}
\begin{center}
\begin{tabular}{l c }
\hline
            & Model 1 \\
\hline
(Intercept) & $0.12$       \\
            & $(4.50)$     \\
x           & $5.28^{***}$ \\
            & $(1.41)$     \\
\hline
R$^2$       & 0.02         \\
Adj. R$^2$  & 0.01         \\
Num. obs.   & 1000         \\
\hline
\multicolumn{2}{l}{\scriptsize{\textsuperscript{***}$p<0.001$, 
  \textsuperscript{**}$p<0.01$, 
  \textsuperscript{*}$p<0.05$}}
\end{tabular}
\caption{Statistical models}
\label{table:coefficients}
\end{center}
\end{table}

cl   <- function(dat,fm, cluster){
           cluster = as.numeric(cluster)
       attach(dat, warn.conflicts = F)
           library(sandwich)
        library(lmtest)
           M <- length(unique(cluster))
           N <- length(cluster)
           K <- fm$rank
           dfc <- (M/(M-1))*((N-1)/(N-K))
           uj  <- apply(estfun(fm),2, function(x) tapply(x, cluster, sum));
           vcovCL <- dfc*sandwich(fm, meat=crossprod(uj)/N)
           coeftest(fm, vcovCL) }

result = cl(d,m,IDs)
texreg(result,omit.coef="IDs")

Ошибка в (функция (классы, fdef, mtable): невозможно найти унаследованный метод для функции "extract" для подписи "coeftest" '

Answer 1

Раздел 5.5 виньетки пакета предлагает решение. Виньетка пакета была опубликована в виде статьи в журнале статистического программного обеспечения. Его можно найти здесь: http://www.jstatsoft.org/v55/i08/.

Более конкретно, надежные стандартные ошибки и значения p должны быть извлечены из вашего result матрица и передана texreg через override.se а также override.pval аргументы:

se <- result[, 2]
pval <- result[, 4]
screenreg(    # display the results in the R console
    m, 
    omit.coef = "IDs", 
    override.se = se, 
    override.pval = pval
)
texreg(       # for LaTeX output
    m, 
    omit.coef = "IDs", 
    override.se = se, 
    override.pval = pval
)

Еще один способ сделать это - извлечь коэффициенты из исходной модели, сохранить их в texreg объект, манипулировать этим объектом, а затем передать его texreg функция:

tr <- extract(m)
tr@pvalues <- result[, 4]
tr@se <- result[, 2]
screenreg(tr)  # or texreg including your original arguments...