[R] Retrieve regression summary results after rq

Camila Mendes cacamendes85 at gmail.com
Wed Sep 26 17:39:46 CEST 2012


This does not work at all...

On Wed, Sep 26, 2012 at 12:04 PM, John Kane <jrkrideau at inbox.com> wrote:
> 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!
> Check it out at http://www.inbox.com/marineaquarium
>
>




More information about the R-help mailing list