[R] Retrieve regression summary results after rq

Camila Mendes cacamendes85 at gmail.com
Wed Sep 26 16:39:30 CEST 2012


Hi all,

I am using quantile regression with svy design. I want to retrieve
summary regression statistics (std error, p-value), since I don't have
any in my output:

Commands:

clus1_d<- svydesign(id=~cd002_co, weights=~wtper, strata=~str, data=data)
bclus1<-as.svrepdesign(clus1_d,type="bootstrap",replicates=100)
fit1<- withReplicates(bclus1,quote(coef(rq(newm428b~sch_new_2+sch_new_3+sch_new_4,
tau=c(0.05,0.25,0.5,0.75,0.95),weights=wtper))))

Output:

coefficients(fit1)
            tau= 0.05 tau= 0.25     tau= 0.50 tau= 0.75 tau= 0.95
(Intercept)      2340      2980  3.300000e+03      3750 4360.0000
sch_new_2        185       -30 -5.500000e+01      -250 -360.0000
sch_new_3        155        20  5.032317e-16      -200 -245.0000
sch_new_4         65        20 -1.500000e+01      -200 -309.9998

Any help will be appreciated.

Best,

C.



More information about the R-help mailing list