# Quantitative Ecotoxicology, Page 35, Robust Regression on Order Statistics

Get the data from here and read it into R:

First we need to convert the vector indicating if an observation is censored to TRUE/FALSE: I store it in a new colum called ‘rem2’ (you could also overwrite rem):

Then we can run the Robust Regression on Order Statistics with the ros() function from the NADA package:

Which gives the same mean and standard deviation as the SAS-Makro (5.16 and 2.07).

Code and data are available at my github-repo under filename ‘p35’.

Written on December 1, 2012