Re: [R] summary.rlm (was p-values)

From: Prof Brian Ripley <>
Date: Wed 05 Jul 2006 - 20:11:22 EST

On Wed, 5 Jul 2006, Celso Barros wrote:

> Dear All,
> When I run rlm to obtain robust standard errors, my output does not include
> p-values. Is there any reason p-values should not be used in this case? Is
> there an argument I could use in rlm so that the output does
> include p-values?

First you would have to derive a reliable theory for the tests you want p-values for. The t ratios are approximately normally distributed, but to quote realistic p-values you would need much more accurate distribtution theory.

For summary.lm we can use Student's t distribution and say the results are exact only under normality. There is no analogue for summary.rlm. Note that S does not quote p-values for summmary.glm for similar reasons: R does but I do not consider it to be a good idea.

> Thanks in advance,
> Celso
> [[alternative HTML version deleted]]
> ______________________________________________
> mailing list
> PLEASE do read the posting guide!

Brian D. Ripley,        
Professor of Applied Statistics,
University of Oxford,             Tel:  +44 1865 272861 (self)
1 South Parks Road,                     +44 1865 272866 (PA)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595

______________________________________________ mailing list
PLEASE do read the posting guide!
Received on Wed Jul 05 20:15:31 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Wed 05 Jul 2006 - 22:16:16 EST.

Mailing list information is available at Please read the posting guide before posting to the list.