Changeset 718 for trunk/yat/regression/LinearWeighted.cc

Timestamp:

Dec 26, 2006, 10:56:26 AM (16 years ago)

Author:

Jari Häkkinen

Message:

Addresses #170.

File:

• trunk/yat/regression/LinearWeighted.cc (modified) (2 diffs)

trunk/yat/regression/LinearWeighted.cc

@@ -42 +42 @@
   }
 
+  double LinearWeighted::alpha(void) const
+  {
+    return alpha_;
+  }
+
+  double LinearWeighted::alpha_err(void) const
+  {
+    return sqrt(alpha_var_);
+  }
+
+  double LinearWeighted::beta(void) const
+  {
+    return beta_;
+  }
+
+  double LinearWeighted::beta_err(void) const
+  {
+    return sqrt(beta_var_);
+  }
+
   void LinearWeighted::fit(const utility::vector& x,
                            const utility::vector& y,
@@ -64 +84 @@
   }
 
+  double LinearWeighted::m_x(void) const
+  {
+    return ap_.x_averager().mean();
+  }
+
+  double LinearWeighted::m_y(void) const
+  {
+    return ap_.y_averager().mean();
+  }
+
+  double LinearWeighted::mse(void) const
+  {
+    return mse_;
+  }
+
+  double LinearWeighted::prediction_error(const double x, const double w) const
+  {
+    return sqrt(alpha_var_ + beta_var_*(x-m_x_)*(x-m_x_)+s2(w));
+  }
+
+  double LinearWeighted::s2(double w) const
+  {
+    return s2_/w;
+  }
+
+  double LinearWeighted::standard_error(const double x) const
+  {
+    return sqrt(alpha_var_ + beta_var_*(x-m_x_)*(x-m_x_) );
+  }
+
+  double LinearWeighted::sxx(void) const
+  {
+    return ap_.x_averager().sum_xx_centered();
+  }
+
+  double LinearWeighted::sxy(void) const
+  {
+    return ap_.sum_xy_centered();
+  }
+
+  double LinearWeighted::syy(void) const
+  {
+    return ap_.y_averager().sum_xx_centered();
+  }
+
 }}} // of namespaces regression, yat, and theplu