R help - Sep 2012 - 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.

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!

Strange, I have never used it but presumably you are using the package ?survey??
The manual suggests summary(xx) should work.

What results or error messages do you get?


John Kane
Kingston ON Canada

> -----Original Message-----
> From: cacamendes85 at gmail.com
> Sent: Wed, 26 Sep 2012 12:39:46 -0300
> To: jrkrideau at inbox.com
> Subject: Re: [R] Retrieve regression summary results after rq
> 
> 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
>> 
>>
____________________________________________________________
GET FREE SMILEYS FOR YOUR IM & EMAIL - Learn more at
http://www.inbox.com/smileys
Works with AIM?, MSN? Messenger, Yahoo!? Messenger, ICQ?, Google Talk? and most
webmails