On Sep 26, 2012, at 8:50 AM, John Kane wrote:
> 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?
She didn't return a model object as fit1. She return only the results of the
coef() on the rq(model) object.
Perhaps:
summ1 <- withReplicates(bclus1,
quote( summary( 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))))
Untested, and I've never used 'survey' and 'rq' packages together so this is
based entirely on analogous reasoning.)
--
David.
>
>
> John Kane
> Kingston ON Canada
>
>
>> -----Original Message-----
>> From: cacamende...@gmail.com
>> Sent: Wed, 26 Sep 2012 12:39:46 -0300
>> To: jrkrid...@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 <jrkrid...@inbox.com> wrote:
>>> summary(fit1) perhaps/
>>>
>>> John Kane
>>> Kingston ON Canada
>>>
>>>
>>>> -----Original Message-----
>>>> From: cacamende...@gmail.com
>>>> Sent: Wed, 26 Sep 2012 11:39:30 -0300
>>>> To: r-help@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@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
>
> ______________________________________________
> 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.
David Winsemius, MD
Alameda, CA, USA
______________________________________________
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.
