来自回归分位数的伪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)

但是，我想知道Poryo R2（拟合优度） - 摘要命令似乎为OLS回归提供了这一点，但在分位数回归中的面板数据中却没有。有谁知道，如何提取Pseudo R2？

Answer 1

我不知道rqpd包。一般来说，pR ^ 2是：

1 - deviance / 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

pR ^ 2将是0.268。希望有所帮助。