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’.