Псевдо R ^ 2 от Регрессионных Квантилей для панельных данных - PullRequest
Новая
       23

Псевдо R ^ 2 от Регрессионных Квантилей для панельных данных

0 голосов
/ 03 апреля 2019

Я использую пакет "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 здесь

Iможет получить результат через 

summary(crem.fit1)

Однако я хотел бы знать Pseudo R2 (соответствие), который, как представляется, сводная команда обеспечивает для регрессий OLS, но не в случае данных панели в квантилерегрессия.Кто-нибудь знает, как извлечь Pseudo R2?

1 Ответ

0 голосов
/ 03 апреля 2019

Я не знаю пакет 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

pR ^ 2 будет 0,268.Надеюсь, это поможет.

...