Как мне создать таблицы квантильной регрессии, используя Stargazer? - PullRequest
Новая
       86

Как мне создать таблицы квантильной регрессии, используя Stargazer?

1 голос
/ 23 февраля 2020

Я вычислил следующие квантильные регрессии, используя пакет quantreg 

    qr_10 = rq(inno_DELTA ~ deDomains + R_and_D_pc + Pop_dens + Agr_GDP + Manufacturing_GDP + Service_GDP + Infr_Area_Percent + Res_pc + Debt_GDP + GOV_EXP_GDP + firms_total + factor(landkreis) + factor(jahr), tau = 0.10, data = df_ip_c)
    qr_25 = rq(inno_DELTA ~ deDomains + R_and_D_pc + Pop_dens + Agr_GDP + Manufacturing_GDP + Service_GDP + Infr_Area_Percent + Res_pc + Debt_GDP + GOV_EXP_GDP + firms_total + factor(landkreis) + factor(jahr), tau = 0.25, data = df_ip_c)
    qr_50 = rq(inno_DELTA ~ deDomains + R_and_D_pc + Pop_dens + Agr_GDP + Manufacturing_GDP + Service_GDP + Infr_Area_Percent + Res_pc + Debt_GDP + GOV_EXP_GDP + firms_total + factor(landkreis) + factor(jahr), tau = 0.5, data = df_ip_c)
    qr_75 = rq(inno_DELTA ~ deDomains + R_and_D_pc + Pop_dens + Agr_GDP + Manufacturing_GDP + Service_GDP + Infr_Area_Percent + Res_pc + Debt_GDP + GOV_EXP_GDP + firms_total + factor(landkreis) + factor(jahr), tau = 0.75, data = df_ip_c)
    qr_95 = rq(inno_DELTA ~ deDomains + R_and_D_pc + Pop_dens + Agr_GDP + Manufacturing_GDP + Service_GDP + Infr_Area_Percent + Res_pc + Debt_GDP + GOV_EXP_GDP + firms_total + factor(landkreis) + factor(jahr), tau = 0.95, data = df_ip_c)

Я пытаюсь отобразить эти регрессии, используя stargazer. Код, который я пытаюсь запустить, выглядит следующим образом: 

stargazer(qr_10,qr_25,qr_50,qr_75,qr_95, rq.se = "iid", type = "text", title="Regression Results" ,initial.zero = F,single.row=TRUE, out="table_quantile_regression.html")

Однако я получаю следующее сообщение об ошибке 

Error in base::backsolve(r, x, k = k, upper.tri = upper.tri, transpose = transpose,  : 
  singular matrix in 'backsolve'. First zero in diagonal [421]

Я предполагаю, что это сообщение об ошибке связано с Стандартная функция ошибок rq.se в stargazer, как, например, summary(qr_10, se = "iid"), работает нормально.

У кого-нибудь есть решение этой проблемы? Спасибо.

1 Ответ

0 голосов
/ 09 марта 2020

Полагаю, это из-за ваших данных. «особая матрица», указывающая на некоторую проблему с регрессией. Минимальный пример прекрасно работает, и stargazer поддерживает quantreg: 

library(quantreg)
library(stargazer)

stargazer(mtcars, type="text")

qr_10 = rq(mpg ~ factor(cyl) + disp, tau = 0.10, data = mtcars)
qr_25 = rq(mpg ~ factor(cyl) + disp, tau = 0.25, data = mtcars)
qr_50 = rq(mpg ~ factor(cyl) + disp, tau = 0.50, data = mtcars)
qr_75 = rq(mpg ~ factor(cyl) + disp, tau = 0.75, data = mtcars)
qr_95 = rq(mpg ~ factor(cyl) + disp, tau = 0.95, data = mtcars)

stargazer(qr_10,qr_25,qr_50,qr_75,qr_95, rq.se = "iid", type = "text", title="Regression Results", initial.zero = F,single.row=TRUE)
stargazer(qr_10,qr_25,qr_50,qr_75,qr_95, type = "text", title="Regression Results", initial.zero = F,single.row=TRUE)

Regression Results
=======================================================================================================
                                                Dependent variable:                                    
             ------------------------------------------------------------------------------------------
                                                        mpg                                            
                   (1)               (2)               (3)               (4)                (5)        
-------------------------------------------------------------------------------------------------------
factor(cyl)6 -2.462** (.992)   -2.105 (1.324)   -4.365*** (1.453) -8.598*** (1.605)  -14.011*** (1.024)
factor(cyl)8  -2.243 (1.735)   -2.608 (2.317)    -3.885 (2.542)   -11.437*** (2.809) -18.811*** (1.792)
disp         -.026*** (.006)   -.027*** (.009)   -.033*** (.009)     -.012 (.010)       .013* (.007)   
Constant     24.630*** (.858) 25.705*** (1.145) 29.916*** (1.256) 31.575*** (1.388)   33.011*** (.885) 
-------------------------------------------------------------------------------------------------------
Observations        32               32                32                 32                 32        
=======================================================================================================

Можете ли вы дать дополнительную информацию или сводную статистику для ваших данных?

Проверка https://stats.stackexchange.com/questions/78022/cause-of-singularity-in-matrix-for-quantile-regression и связанные посты. С наилучшими пожеланиями

...