[R] package for measurement error models

Prof Brian Ripley ripley at stats.ox.ac.uk
Wed Aug 11 08:19:12 CEST 2010


On Tue, 10 Aug 2010, Carrie Li wrote:

> Thank you both!
> 
> I found that the model with both x and y have measurement errors should be
> pretty common in practice. But it seems to me that there is no a simple
> solution for it...(I mean, a ready-to-use package or program handing this
> model fitting problem. )

Hmm, have you actually looked at my reference?  We provided a program 
in 1987.

> One more question, is any reference that considers weighted measurement
> error models (use the known variance of measurement errors as weights) ?

That's not really how you use known variances, which is the main case 
covered in our paper.

>  
> (I've search on google..got nothing...)
> 
> Thanks for your sharing opinions again! I appreciate!
> 
> Carrie--
> 
> On Tue, Aug 10, 2010 at 7:54 AM, Prof Brian Ripley <ripley at stats.ox.ac.uk>
> wrote:
>       On Tue, 10 Aug 2010, peter dalgaard wrote:
> 
>
>             On Aug 10, 2010, at 3:52 AM, Carrie Li wrote:
>
>                   Thanks. I found the code in the link you
>                   gave me very helpful.
>                   But, I just have few questions regarding
>                   the code.
>                   It seems to me that in (from
>                   wikipdeia)Deming regression, it assumes
>                   that
>                   the ratios of the variances of two
>                   measurement errors are constant for all
>                   pairs of (x_i, y_i). However, if the
>                   ratios are not constant, (i.e. the
>                   variances of measurement are
>                   heterogeneous) , is it still appropriate
>                   to use
>                   Deming regression ?
> 
>
>             In a word, no.
>
>             One way of looking at it is that as the ratio of
>             variances varies from 0 to infinity, the analysis
>             goes from regression of y on x to (inverse)
>             regression of x on y, and those give different
>             results, not just numerically but also
>             asymptotically. I.e., getting the ratio wrong gives
>             an inconsistent estimate; getting it wrong for some
>             of the data, as is bound to happen if you assume it
>             constant and it isn't, will also give a inconsistent
>             estimate. Unless, that is, you can find a definition
>             of "average ratio" that eliminates the bias, but I
>             don't think it is worth the paperwork.
>
>             Rather, I'd suggest direct minimization of the SSR
>             (from the Wikipedia page), noting that you can plug
>             in x_i^* as a function of beta also if the
>             _individual_ ratios are known. (I get the feeling
>             that someone must have been here before, so possibly
>             others can fill in the gaps?) For modest sample
>             sizes, it might also be possible to
> 
> 
> Yes, people have been there before. Mike Thompson and I published a
> now-much-cited-in-analytical-chemistry paper in The Analyst in 1987. A
> companion paper was rejected by a mainstream statistics journal as
> 'already known', but the journal editor was unable to get any prior
> publication out of the referee.
>
>       formulate the problem as a nonlinear model and use nls().
> 
> 
> Direct minimization is simple enough.
> 
>
>       --
>       Peter Dalgaard
>       Center for Statistics, Copenhagen Business School
>       Solbjerg Plads 3, 2000 Frederiksberg, Denmark
>       Phone: (+45)38153501
>       Email: pd.mes at cbs.dk  Priv: PDalgd at gmail.com
> 
> 
> --
> 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 272866 (PA)
> Oxford OX1 3TG, UK                Fax:  +44 1865 272595
> 
> 
> 
>

-- 
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 272866 (PA)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595


More information about the R-help mailing list