来自回归分位数的伪R ^ 2用于面板数据

Question

我正在使用“ rqpd”包在R中运行以下分位数回归的面板数据：

library(rqpd)
data(bwd)
cre.form1 <- dbirwt ~ smoke + dmage + agesq + novisit + pretri2 + 
pretri3 | momid3 |smoke
crem.fit1 <- rqpd(cre.form1, panel(method="cre", taus <- c(0.9)),data=bwd)

rqpd软件包格式here

我可以通过

得到结果

summary(crem.fit1)

但是，我想知道Pseudo R2（拟合优度）-摘要命令似乎为OLS回归提供了此功能，但对于分位数回归中的面板数据而言却不是。 有人知道如何提取伪R2吗？

Answer 1

我不知道rqpd软件包。通常，pR ^ 2为：

  

1-偏差/null.deviance

我使用的是标准逻辑回归示例：

utils::data(anorexia, package = "MASS")
anorex.1 <- glm(Postwt ~ Prewt + Treat + offset(Prewt),
                family = gaussian, data = anorexia)

当我计算pR ^ 2时，我正在使用：

library(broom)
library(dplyr)
anorex.1 %>% glance() %>%
  mutate(pR2 = 1 - deviance/null.deviance)

输出为：

# A tibble: 1 x 8
  null.deviance df.null logLik   AIC   BIC deviance df.residual   pR2
          <dbl>   <int>  <dbl> <dbl> <dbl>    <dbl>       <int> <dbl>
1         4525.      71  -240.  490.  501.    3311.          68 0.268

p R ^ 2为0.268。希望有帮助。