[R] intercept of 3D line? (Orthogonal regression)
m_olshansky at yahoo.com
Tue Sep 2 06:55:02 CEST 2008
I do not see why you can not use regression even in this case.
To make things more simple suppose that the exact model is:
y = a + b*x, i.e.
y1 = a + b*x1
yn = a + b*xn
But you can not observe y and x. Instead you observe
ui = xi + ei (i=1,...,n) and
vi = yi + di (i=1,...,n)
Now you have
vi = yi + di = a + b*xi + di = a + b*(ui - ei) + di
= a + b*ui + (di - b*ei)
and under regular assumptions about ei's end di's we get a standard regression problem (note that b is unknown to you but is constant).
--- On Tue, 2/9/08, William Simpson <william.a.simpson at gmail.com> wrote:
> From: William Simpson <william.a.simpson at gmail.com>
> Subject: [R] intercept of 3D line? (Orthogonal regression)
> To: r-help at r-project.org
> Received: Tuesday, 2 September, 2008, 4:53 AM
> I posted before recently about fitting 3D data x, y, z where
> all have
> error attached.
> I want to predict z from x and y; something like
> z = b0 + b1*x + b2*y
> But multiple regression is not suitable because all of x,
> y, and z have errors.
> I have plotted a 3D scatterplot of some data using rgl. I
> see that the
> data form a cigar-shaped cloud. I think multiple regression
> is only
> suitable when the points fall on a plane (forgetting about
> the error
> in x and y).
> I now know the right way how to find the best fitting plane
> to x,y,z
> data using princomp.
> But a new problem is how to get the best fitting *line*. I
> know how to do that too using princomp. But there is a
> problem: there's no way to specify a line in 3D space
> in the form
> z=f(x,y) or in other words with an intercept and slopes.
> Instead, one way to deal with the problem is to use a
> version of the line: you use an arbitrary starting point
> x0, y0, z0
> and the direction vector of your line (I know how to get
> the direction
> BUT how do I get the intercept??? At this point my lines
> just go
> through the origin.
> Do I just use $center from the princomp output modified in
> some way?
> Thanks for any help!
> R-help at r-project.org mailing list
> PLEASE do read the posting guide
> and provide commented, minimal, self-contained,
> reproducible code.
More information about the R-help