I am using the "rqpd" package to run the following panel data of quantile regression in 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)
```

the rqpd package form here

I can get the result via

```
summary(crem.fit1)
```

However, I would like to know the Pseudo R2 (goodness of fit) - which the summary command seems to deliver for OLS regressions, but not in case of panel data in quantile regressions. Does anybody know, how to extract the Pseudo R2?

I don't know the

`rqpd`

package. In General the pR^2 is:I am using a the standard logistic regression example:

When I calculate the pR^2, I am using:

Output would be:

pR^2 would be 0.268. Hope that helps.