我在 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)
我summary(reg_q1)$coefficients
用来获取系数和标准误差。一年前它工作得很好,但一年后我回到了我的代码,摘要函数只返回一个像这样的字符:
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
你知道我怎么能从中得到标准错误吗?软件包最近在输出方面是否发生了变化?