From 834bb45be7e74553f51da00bad95595666586a7c Mon Sep 17 00:00:00 2001 From: Elizabeth Hunt Date: Tue, 10 Oct 2023 14:44:33 -0600 Subject: problem one hw 4 --- src/vector,least-squares.lisp | 14 ++++++++++++++ 1 file changed, 14 insertions(+) create mode 100644 src/vector,least-squares.lisp (limited to 'src/vector,least-squares.lisp') diff --git a/src/vector,least-squares.lisp b/src/vector,least-squares.lisp new file mode 100644 index 0000000..687af32 --- /dev/null +++ b/src/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))) + -- cgit v1.3