我在R中使用rqpd软件包进行具有固定效果的分位数回归(quantreg软件包不支持具有固定效果(,如下所示:
reg_q1 <- rqpd::rqpd(NoneHousingExpense ~ totalIncome + treatmentMonthly +
as.factor(year) + HHsize_eq| Address,
panel(method = "pfe", taus = seq(0.02,0.98,0.02), tauw = rep(1/49,49)),
data = adj_df1, weights = adj_df1$rWeight)
并且我使用CCD_ 1来获得系数和标准误差。一年前它运行得很好,但一年后我又回到了我的代码中,summary函数只返回这样一个字符:
Length Class Mode
coefficients 246 -none- numeric
residuals 1313054 -none- numeric
coef 246 -none- numeric
X 1 matrix.csr S4
Z 1 matrix.csr S4
y 26797 -none- numeric
ids 26797 -none- numeric
panel 6 -none- list
control 7 -none- list
formula 3 formula call
call 5 -none- call
rank 1 -none- numeric
你知道我怎样才能从中得到标准错误吗?该方案最近在产出方面有变化吗?
我遇到了同样的问题。与其使用summary(reg_q1(,不如尝试使用rqpd::summary.rqpd(reg_q 1(。这为我解决了问题。
答案在这里:
summary.rqpd(reg_q1)