blob: 687af3247e06a7d380d07b328afc7ee30d098df3 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
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)))
|
