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. ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.