[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