[RsR] lmrob: lmrob.S()-Warning is not in the result returned

[Martin Maechler]

>>>>> "ReneL" == Locher René (lor) <lor using zhwin.ch>
>>>>>     on Tue, 5 Jun 2007 18:24:47 +0200 writes:

    ReneL> Dear Matias When I run the following code I get 2
    ReneL> warnings: 

    ReneL> lmrob.S(): refinements did not converge (tol=1e-07) in 200 iterations

     library(robustbase)
     test <- data.frame(lconc=c(3.081910,3.165475,2.501436,3.650658,
				3.044522,3.660994),
                        dist=c(100,180,280,30,220,6))
     res <- lmrob(lconc~dist, dat=test)

    ReneL> res$converged gives TRUE as the IRWLS iterations after the lmrob.S() have converged.

    ReneL> It would be nice to have also the convergence status of lmrob.S() in the result object of lmrob().

I agree. 

HOWEVER, I cannot see / reproduce the convergence warnings
problems you mention, neither the above, nor the one you sent
yesterday.

Are you *REALLY* sure, you are not using an older version of
"robustbase" or mixing versions of robustbase,
or are using 'lmrob' (there was a package of that name) and
'robustbase' simultaneously?

Note: For such reports, do use sessionInfo(). 

Regards,
Martin

    ReneL> Kind regards

    ReneL> René

 
    ReneL> René Locher                                   E-Mail: rene.locher using zhwin.ch
    ReneL> Institut für Datenanalyse und Prozessdesign   Tel:    052/ 267 7810
    ReneL> Zürcher Hochschule Winterthur                 Fax:    052/ 268 7810
    ReneL> Technopark / Jägerstrasse 2                   
    ReneL> Postfach
    ReneL> CH-8400 Winterthur                            http://www.idp.zhwin.ch
 
    ReneL> P.S.

    ReneL> Here is the version of my system on which I have tested the example:

    ReneL> version              _                           
    ReneL> ## platform       i386-pc-mingw32             
    ReneL> ## arch           i386                        
    ReneL> ## os             mingw32                     
    ReneL> ## system         i386, mingw32               
    ReneL> ## status                                     
    ReneL> ## major          2                           
    ReneL> ## minor          5.0                         
    ReneL> ## year           2007                        
    ReneL> ## month          04                          
    ReneL> ## day            23                          
    ReneL> ## svn rev        41293                       
    ReneL> ## language       R                           
    ReneL> ## version.string R version 2.5.0 (2007-04-23)

    ReneL> library(help=robustbase)  
    ReneL> ## Version:       0.2-7