From f16a58979fb59afd84f4b4dd046d50276f891cdf Mon Sep 17 00:00:00 2001
From: AJ Rossini <blindglobe@gmail.com>
Date: Thu, 22 Jan 2009 17:02:57 +0100
Subject: [PATCH] move lin regr example away until we get it solved for
 lisp-matrix.

Signed-off-by: AJ Rossini <blindglobe@gmail.com>
---
 TODO.lisp | 115 ++------------------------------------------------------------
 1 file changed, 2 insertions(+), 113 deletions(-)

diff --git a/TODO.lisp b/TODO.lisp
index c1db818..2a7616c 100644
--- a/TODO.lisp
+++ b/TODO.lisp
@@ -1,6 +1,6 @@
 ;;; -*- mode: lisp -*-
 
-;;; Time-stamp: <2009-01-20 08:19:58 tony>
+;;; Time-stamp: <2009-01-22 17:02:09 tony>
 ;;; Creation:   <2008-09-08 08:06:30 tony>
 ;;; File:       TODO.lisp
 ;;; Author:     AJ Rossini <blindglobe@gmail.com>
@@ -461,115 +461,4 @@
 			    :top-k-of-n-select))
 		:parametric
 		:partially-parametric))
-    "start of ontology")
-
-
-
-  )
-
-
-
-
-
-#+nil
-(progn  
-  ;; FIXME: Numerical linear algebra, least squares, and QR
-  ;; factorization
-
-  ;; Need to incorporate the xGEQRF routines, to support linear
-  ;; regression work.   
-
-  ;; Some issues exist in the LAPACK vs. LINPACK variants, hence R
-  ;; uses LINPACK primarily, rather than LAPACK.  See comments in R
-  ;; source for issues.  
-  
-  ;; LAPACK suggests to use the xGELSY driver (GE general matrix, LS
-  ;; least squares, need to lookup Y intent (used to be an X alg, see
-  ;; release notes).
-
-  ;; Goal is to start from X, Y and then realize that if
-  ;; Y = X \beta, then,   i.e. 8x1 = 8xp px1  + 8x1
-  ;;      XtX \hat\beta = Xt Y
-  ;; so that we can solve the equation  W \beta = Z   where W and Z
-  ;; are known, to estimate \beta.
-  (defparameter *xv*
-    (make-vector
-     8
-     :initial-contents '((1d0 3d0 2d0 4d0 3d0 5d0 4d0 6d0))))
-
-  (defparameter *xv+1*
-    (make-matrix
-     8 2
-     :initial-contents '((1d0 1d0)
-			 (1d0 3d0)
-			 (1d0 2d0)
-			 (1d0 4d0)
-			 (1d0 3d0)
-			 (1d0 5d0)
-			 (1d0 4d0)
-			 (1d0 6d0))))
-
-  (defparameter *xm*
-    (make-matrix
-     2 8
-     :initial-contents '((1d0 3d0 2d0 4d0 3d0 5d0 4d0 6d0)
-			 (1d0 2d0 3d0 4d0 5d0 6d0 7d0 8d0))))
-
-  (defparameter *y*
-    (make-vector
-     8
-     :initial-contents '((1d0 2d0 3d0 4d0 5d0 6d0 7d0 8d0))))
-
-  ;; so something like (NOTE: matrices are transposed to begin with,
-  ;; hence the incongruety)
-  (defparameter *xtx* (m* *xv* (transpose *xv*)))
-  (defparameter *xty* (m* *xv* (transpose  *y*)))
-  (defparameter *rcond* 1d0) ; condition number bound -- examples
-			   ; suggest should be zero on input? 
-
-  (defparameter *betahat*  (gelsy *xtx* *xty* *rcond*))
-  *betahat*
-
-#|
-  (#<LA-SIMPLE-VECTOR-DOUBLE (1 x 1)
-     1.293103448275862>
-     1)
-
-   ## Test case in R:
-   x <- c( 1.0, 3.0, 2.0, 4.0, 3.0, 5.0, 4.0, 6.0)
-   y <- c( 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0)
-   lm (y ~ x  -1)
-   ## => 
-   Call:
-   lm(formula = y ~ x - 1)
-
-   Coefficients:
-      x  
-   1.293  
-|#
-
-
-  ;; so something like (NOTE: matrices are transposed to begin with, hence the incongruety)
-  (defparameter *xtx* (m* *xv+1* (transpose *xv+1*)))
-  (defparameter *xty* (m* (transpose *xv+1*) (transpose  *y*)))
-  (defparameter *rcond* 1d0)
-  (defparameter *jpvt* (make-matrix 1 8))
-  (defparameter *betahat*  (gelsy *xtx* *xty* *rcond* *jpvt*))
-  ;; Kills SBCL, with:
-  ;; * Parameter 7 to routine DGELSYSPM was incorrect
-  ;; clearly a bit of error checking needed!
-
-  *betahat*
-
-
-
-  ;; which suggests one might do (modulo ensuring correct orientations)
-  (defun lm (x y)
-    (let ((betahat (gelsy (m* x (transpose x))
-			  (m* x (transpose y)))))
-      
-      (values betahat (sebetahat betahat x y))))
-  ;; to get a results list containing betahat and SEs
-
-  (values-list '(1 3 4))
-  )
\ No newline at end of file
+    "start of ontology"))
-- 
2.11.4.GIT