       Error on slope and intercept

• To: mathgroup at smc.vnet.net
• Subject: [mg131774] Error on slope and intercept
• From: nicolasgebbo at googlemail.com
• Date: Wed, 2 Oct 2013 05:54:16 -0400 (EDT)
• Delivered-to: l-mathgroup@mail-archive0.wolfram.com
• Delivered-to: l-mathgroup@wolfram.com
• Delivered-to: mathgroup-outx@smc.vnet.net
• Delivered-to: mathgroup-newsendx@smc.vnet.net

I was informing myself about linear regression with errors in x and y.
I dont really understand the theory behind the following algorithm but it would be nice if someone could tell me how i get the error of the intercept and slope.

worals[x_, y_, sx_, sy_] :=
Block[{a, b, f, z, u = 1/sx, v = 1/sy, w = (sy/sx)^2},

{a, b} = (y*v).PseudoInverse@{v, x*v}; f = #.# &[(a + b*x - y) v];
While[f > (z = (x*w + (y - a) b)/(b^2 + w);
{a, b} = (y*v).PseudoInverse@{v, z*v};
f = #.# &@Join[(z - x) u, (a + b*z - y) v])];

{f, {a, b}}]

Thanks for any help

• Prev by Date: Error in slope and intercept
• Next by Date: Re: Output display by slide
• Previous by thread: Re: Error in slope and intercept
• Next by thread: small string issue