[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