[R] nonlinear fitting when both x and y having measurement error?

Prof Brian Ripley ripley at stats.ox.ac.uk
Wed Sep 12 09:03:19 CEST 2001


On Wed, 12 Sep 2001, Etsushi Kato wrote:

> Dear r-help,
>
> I want to conduct nonlinear fitting to a data frame having x and y
> variables.  Because both x and y have measurement error, I want to
> include error term of x variable in the model.  I'm not sure but I
> think ordinary nls model only consider error term of y variable.
>
> How can I do this kind of nonlinear fitting in R.  Is there any
> examples in nls package?

That is not a least-squares problem.  Even in the simple linear case (one
x one y) it's a hard problem, one that cannot be solved without more
information (for example on the ratio of the error variances, or knowing
one of them). As far as I know there is no software available in R for
that case (although it's not hard to write).

I think you need to write down a suitable likelihood and optimize it
numerically (with optim).

This topic is not well covered in regression texts, although it is in that
by Sprent, for example.  There is a specialist text by W. A. Fuller (1987)
Measurement Error Models.


-- 
Brian D. Ripley,                  ripley at stats.ox.ac.uk
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford,             Tel:  +44 1865 272861 (self)
1 South Parks Road,                     +44 1865 272860 (secr)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._



More information about the R-help mailing list