[R] How to get around heteroscedasticity with non-linear leas t squares in R?

Peter Dalgaard p.dalgaard at biostat.ku.dk
Tue Feb 21 21:28:14 CET 2006


"Liaw, Andy" <andy_liaw at merck.com> writes:

> Your understanding isn't similar to mine.  Mine says robust/resistant
> methods are for data with heavy tails, not heteroscedasticity.  The common
> ways to approach heteroscedasticity are transformation and weighting.  The
> first is easy and usually quite effective for dose-response data.  The
> second is not much harder.  Both can be done in R with nls().

And there is gnls() which allows direct modelling of the variance.

        -p
 
> Andy
> 
> From: Quin Wills
> > 
> > I am using "nls" to fit dose-response curves but am not sure 
> > how to approach
> > more robust regression in R to get around the problem of the my error
> > showing increased variance with increasing dose.  
> > 
> >  
> > 
> > My understanding is that "rlm" or "lqs" would not be a good idea here.
> > 'Fairly new to regression work, so apologies if I'm missing something
> > obvious.
> > 
> >  
> > 
> > 
> > 	[[alternative HTML version deleted]]
> > 
> > ______________________________________________
> > R-help at stat.math.ethz.ch mailing list
> > https://stat.ethz.ch/mailman/listinfo/r-help
> > PLEASE do read the posting guide! 
> > http://www.R-project.org/posting-guide.html
> > 
> >
> 
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
> 

-- 
   O__  ---- Peter Dalgaard             Øster Farimagsgade 5, Entr.B
  c/ /'_ --- Dept. of Biostatistics     PO Box 2099, 1014 Cph. K
 (*) \(*) -- University of Copenhagen   Denmark          Ph:  (+45) 35327918
~~~~~~~~~~ - (p.dalgaard at biostat.ku.dk)                  FAX: (+45) 35327907




More information about the R-help mailing list