Services & Resources / Wolfram Forums / MathGroup Archive

MathGroup Archive 2013

[Date Index] [Thread Index] [Author Index]

Search the Archive

Error on slope and intercept

  • To: mathgroup at
  • Subject: [mg131774] Error on slope and intercept
  • From: nicolasgebbo at
  • Date: Wed, 2 Oct 2013 05:54:16 -0400 (EDT)
  • Delivered-to:
  • Delivered-to:
  • Delivered-to:
  • Delivered-to:

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