Anonymous user
Multiple regression: Difference between revisions
→{{header|zkl}}: added GSL solution
m (remove hard-coded whitespace & big tags) |
(→{{header|zkl}}: added GSL solution) |
||
Line 1,237:
=={{header|zkl}}==
Using the GNU Scientific Library:
<lang zkl>var [const] GSL=Import("zklGSL"); // libGSL (GNU Scientific Library)
height:=GSL.VectorFromData(1.47, 1.50, 1.52, 1.55, 1.57, 1.60, 1.63,
1.65, 1.68, 1.70, 1.73, 1.75, 1.78, 1.80, 1.83);
weight:=GSL.VectorFromData(52.21, 53.12, 54.48, 55.84, 57.20, 58.57, 59.93,
61.29, 63.11, 64.47, 66.28, 68.10, 69.92, 72.19, 74.46);
v:=GSL.polyFit(height,weight,2);
v.format().println();
GSL.Helpers.polyString(v).println();
GSL.Helpers.polyEval(v,height).format().println();</lang>
{{out}}
<pre>
128.81,-143.16,61.96
128.813 - 143.162x + 61.9603x^2
52.25,53.48,54.36,55.77,56.77,58.37,60.08,61.28,63.18,64.50,66.58,68.03,70.30,71.87,74.33
</pre>
Or, using Lists:
{{trans|Common Lisp}}
<lang zkl>// Solve a linear system AX=B where A is symmetric and positive definite, so it can be Cholesky decomposed.
| |||