summaryrefslogtreecommitdiff
path: root/regress.m
diff options
context:
space:
mode:
authorBenjamin Franzke <benjaminfranzke@googlemail.com>2012-05-12 15:59:20 +0200
committerBenjamin Franzke <benjaminfranzke@googlemail.com>2012-05-12 15:59:20 +0200
commit020f82903d960f379a3b127c0f46d6bf1106cb26 (patch)
tree8de5740b2d08d8e4063d639bbebf1a6f6c2de5ec /regress.m
parentcc8b965d52177eb45185d0b1ead65b55f113a687 (diff)
downloaddmc-020f82903d960f379a3b127c0f46d6bf1106cb26.tar.gz
dmc-020f82903d960f379a3b127c0f46d6bf1106cb26.tar.bz2
dmc-020f82903d960f379a3b127c0f46d6bf1106cb26.zip
regress_pred: Add degree param, to fold in quad regression
Diffstat (limited to 'regress.m')
-rw-r--r--regress.m7
1 files changed, 4 insertions, 3 deletions
diff --git a/regress.m b/regress.m
index ad06796..4d13c62 100644
--- a/regress.m
+++ b/regress.m
@@ -1,4 +1,5 @@
-function mn = regress(x, y)
- mn = [ x ones(length(x), 1) ] \ y;
+function polynom = regress(x, y, degree)
+ pot = @(x) x.^[degree:-1:0];
+ M = reshape(cell2mat(arrayfun(pot, x(:)', 'UniformOutput', false)), degree+1, [])';
+ polynom = M \ y(:);
endfunction
-