summaryrefslogtreecommitdiff
path: root/deprecated-cl/vector,least-squares.lisp
diff options
context:
space:
mode:
Diffstat (limited to 'deprecated-cl/vector,least-squares.lisp')
-rw-r--r--deprecated-cl/vector,least-squares.lisp14
1 files changed, 14 insertions, 0 deletions
diff --git a/deprecated-cl/vector,least-squares.lisp b/deprecated-cl/vector,least-squares.lisp
new file mode 100644
index 0000000..687af32
--- /dev/null
+++ b/deprecated-cl/vector,least-squares.lisp
@@ -0,0 +1,14 @@
+(in-package :lizfcm.vector)
+
+(defun least-squares-reg (x y)
+ (let* ((n (length x))
+ (sum-y (reduce #'+ y))
+ (sum-x (reduce #'+ x))
+ (sum-xy (reduce #'+ (mapcar #'* x y)))
+ (sum-xsquared (reduce #'+ (mapcar #'* x x)))
+ (b (/ (- (* sum-y sum-xsquared) (* sum-x sum-xy))
+ (- (* n sum-xsquared) (* sum-x sum-x))))
+ (a (/ (- (* n sum-xy) (* sum-x sum-y))
+ (- (* n sum-xsquared) (* sum-x sum-x)))))
+ (list a b)))
+