# HG changeset patch
# User Frederic Bastien <bastienf@iro.umontreal.ca>
# Date 1215448026 14400
# Node ID 74b402b5a81b91227e62c1f62994874ca3a75050
# Parent  835830e52b42202a2bb197877777a256a8422a26
small modif by yoshue

diff -r 835830e52b42 -r 74b402b5a81b linear_regression.py
--- a/linear_regression.py	Mon Jul 07 12:02:34 2008 -0400
+++ b/linear_regression.py	Mon Jul 07 12:27:06 2008 -0400
@@ -88,10 +88,11 @@
     P = LinearPredictorEquations
     XtX = T.matrix() # (n_inputs+1) x (n_inputs+1)
     XtY = T.matrix() # (n_inputs+1) x n_outputs
-    extended_input = T.prepend_scalar_to_each_row(1,P.inputs)
+    extended_input = T.prepend_scalar_to_each_row(1.,P.inputs)
     new_XtX = add_inplace(XtX,T.dot(extended_input.T,extended_input))
     new_XtY = add_inplace(XtY,T.dot(extended_input.T,P.targets))
-    
+    new_theta = T.Cholesky_solve_inplace(P.theta,XtX,XtY)  # solve linear system XtX theta = XtY 
+
 class LinearPredictor(object):
     """
     A linear predictor has parameters theta (a bias vector and a weight matrix)