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")
可以正常工作。
有人有解决这个问题的方法吗?谢谢你。