[RsR] hubers m-estimator in R / SPSS

Rudolf Dutter R@Dutter @end|ng |rom tuw|en@@c@@t
Wed Feb 23 21:53:41 CET 2011


Martin,

I agree fully with you. It reminds me that some 30 years ago I puzzled
around why SPSS produces negative variances. It was not too difficult to
find out the reason by trial and error. Conclusion: We didn't recommend
the usage SPSS at a scientific level any more, and some years ago, we
cancelled all licenses from our institute. I am not surprised that
Hubers algorithm remains mysterious in that package.
Sorry to say that.

Rudi

Martin Maechler wrote:
>>>>>> "MH" == Manfred Hammerl <manfred.hammerl using edu.uni-graz.at>
>>>>>>     on Wed, 23 Feb 2011 17:14:51 +0100 writes:
>>>>>>             
>
>     MH> thanks. of course, that was my idea too, but i couldn't find out which 
>     MH> scale estimator spss uses. no manual, book, paper or whatever contains 
>     MH> information about that. i tried other programs (like R or the ACM-tool), 
>     MH> which use MAD, SMAD or an iterated standard deviation, but i always get 
>     MH> different results with spss...
>
> so why on earth are you using SPSS at all?
> ...
> actually it's hard to believe that they don't document what
> their M-estimate does ...
> but that would really belong to an SPSS help list rather than an
> R one, right?
>
> Martin
>
>     MH> manfred
>
>
>     MH> Am 22.02.2011 00:42, schrieb Matias Salibian-Barrera:
>     >> Hello Manfred,
>     >> 
>     >> I'm not familiar with SPSS, but based on my experience, I'll go out on a
>     >> limb and say that the difference may be on the scale estimator that it's
>     >> used (either a preliminary estimator like the MAD) or a simultaneously
>     >> computer M-scale. Hopefully the SPSS manual will have some details.
>     >> 
>     >> Matias
>     >> 
>     >> 
>     >> On 11-02-18 06:35 PM, Manfred Hammerl sat down at the computer and wrote:
>     >>> hello, i'm new to robust statistics but found out very quick that R
>     >>> (used huber, hubers and huberM) and SPSS (huber m-estimator) calculate
>     >>> different location estimates, given the same tuning constant k. since
>     >>> the differences a really not very small, i wanted to get some detailed
>     >>> information about this but i couldn't find out which algorithm SPSS uses
>     >>> to calculate hubers estimator so far...
>     >>> 
>     >>> does anyone know something about SPSS's huber function?
>     >>> 
>     >>> greetings,
>     >>> 
>     >>> manfred
>     >>> 
>     >>> _______________________________________________
>     >>> R-SIG-Robust using r-project.org mailing list
>     >>> https://stat.ethz.ch/mailman/listinfo/r-sig-robust
>     >> 
>
>
>     MH> -- 
>     MH> MMag. Manfred Hammerl, BBakk.
>     MH> manfred.hammerl using edu.uni-graz.at
>
>     MH> _______________________________________________
>     MH> R-SIG-Robust using r-project.org mailing list
>     MH> https://stat.ethz.ch/mailman/listinfo/r-sig-robust
>
> _______________________________________________
> R-SIG-Robust using r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-robust
>
>
>   

-- 
~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=
From: Prof. Dr. Rudolf Dutter
      Dept. of Statistics and Probability Theory
      Vienna University of Technology,
      Wiedner Hauptstr. 7 ("Goldenes Lamm")
      A-1040 Vienna, Austria
      Tel. +43 1 58801/10730
      FAX  +43 1 58801/10798
      E-Mail: R.Dutter using tuwien.ac.at
      Internet: http://www.statistik.tuwien.ac.at/public/dutt/
~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=~=




More information about the R-SIG-Robust mailing list