[R] Retrieve regression summary results after rq

John Kane jrkrideau at inbox.com
Wed Sep 26 17:04:31 CEST 2012


summary(fit1) perhaps/

John Kane
Kingston ON Canada


> -----Original Message-----
> From: cacamendes85 at gmail.com
> Sent: Wed, 26 Sep 2012 11:39:30 -0300
> To: r-help at r-project.org
> Subject: [R] Retrieve regression summary results after rq
> 
> 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 at 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.

____________________________________________________________
FREE 3D MARINE AQUARIUM SCREENSAVER - Watch dolphins, sharks & orcas on your desktop!




More information about the R-help mailing list