[RsR] robustbase V 0.1-5

[Locher René (lor)]

Dear Martin

As you think, that this discussion (and the dataset) is interesting enough for other people working on robustness, I followed your suggestion and posted this mail for the R-SIG-robust-Community. You might recapitulate later our discussion leading to this mail for the rest of the R-SIG-robust-Community.

I have mentioned this morning during our phonecall, that the dataset included cannot be published at the moment officially with all the information statisticians would like. I am sure that a will finally get the permission to do so and you will send additional variables. For the moment I can give you only the following information:

LNOx:   log of hourly mean of NOx concentration in ambient air [ppb] next to a highly frequented motorway
LNOxEm: log of hourly sum of NOx emission of cars on this motorway in arbitrary units
sqrtWS: Square root of wind speed [m/s]

The model I used for this discussion is not the full model I am using to describe the ambient concentration at the moment. Furthermore, I am still working on the improvement of the model. Be aware that this dataset contains, contrary to the former extract, still some NAs!

Here is again the (slightly modified) code I have (finally) used to test lmrob() V 0.1-5 discussed in our mail thread:

library(robustbase)
dat <- na.omit(read.table(file="air.txt"))
dat$julday <- factor(dat$julday)
dim(dat)  ## 8197    5

r.MM01 <- lmrob(LNOx ~ LNOxEm + sqrtWS + julday,
                dat=dat,  compute.rd = FALSE)  

I got the repeated error message: 
Too many singular resamples
Aborting

and could not abort the function by pressing Control-C. A part of this problem has been fixed meanwhile by you.

Kind regards

René



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

-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: air.txt
URL: <https://stat.ethz.ch/pipermail/r-sig-robust/attachments/20060420/53b5c0d1/attachment.txt>