Rosetta Code

Polynomial regression: Difference between revisions

Page Discussion

Polynomial regression (view source)

Revision as of 21:05, 7 November 2021

49 bytes added ,  2 years ago
m
→‎{{header|Phix}}: Replaced IupPlot with IupGraph, added online link
m (→‎{{header|FreeBASIC}}: fixed indent + removed empty lines)
m (→‎{{header|Phix}}: Replaced IupPlot with IupGraph, added online link)
Line 1,688:
=={{header|Phix}}==
{{trans|REXX}}
{{libheader|Phix/pGUIonline}}▼
{{libheader|Phix/pGUI}}
You can run this online [http://phix.x10.mx/p2js/Polynomial_regression.htm here].
<!--<lang Phix>(phixonline)-->
<span style="color: #000080;font-style:italic;">-- (just add ", AXS_YSCALE=LOG10" for a nice log scale)demo\rosetta\Polynomial_regression.exw</span>▼
<span style="color: #008080;">includewith</span> <span style="color: #000000008080;">pGUIjavascript_semantics</span> <span style="color: #0000FF000080;">.</span><span font-style="color: #000000italic;">e-- DEV MINSIZE not honoured</span>▼
<span style="color: #008080;">constant</span> <span style="color: #000000;">x</span> <span style="color: #0000FF;">=</span> <span style="color: #0000FF;">{</span><span style="color: #000000;">0</span><span style="color: #0000FF;">,</span><span style="color: #000000;">1</span><span style="color: #0000FF;">,</span><span style="color: #000000;">2</span><span style="color: #0000FF;">,</span><span style="color: #000000;">3</span><span style="color: #0000FF;">,</span><span style="color: #000000;">4</span><span style="color: #0000FF;">,</span><span style="color: #000000;">5</span><span style="color: #0000FF;">,</span><span style="color: #000000;">6</span><span style="color: #0000FF;">,</span><span style="color: #000000;">7</span><span style="color: #0000FF;">,</span><span style="color: #000000;">8</span><span style="color: #0000FF;">,</span><span style="color: #000000;">9</span><span style="color: #0000FF;">,</span><span style="color: #000000;">10</span><span style="color: #0000FF;">},</span>
<span style="color: #000000;">y</span> <span style="color: #0000FF;">=</span> <span style="color: #0000FF;">{</span><span style="color: #000000;">1</span><span style="color: #0000FF;">,</span><span style="color: #000000;">6</span><span style="color: #0000FF;">,</span><span style="color: #000000;">17</span><span style="color: #0000FF;">,</span><span style="color: #000000;">34</span><span style="color: #0000FF;">,</span><span style="color: #000000;">57</span><span style="color: #0000FF;">,</span><span style="color: #000000;">86</span><span style="color: #0000FF;">,</span><span style="color: #000000;">121</span><span style="color: #0000FF;">,</span><span style="color: #000000;">162</span><span style="color: #0000FF;">,</span><span style="color: #000000;">209</span><span style="color: #0000FF;">,</span><span style="color: #000000;">262</span><span style="color: #0000FF;">,</span><span style="color: #000000;">321</span><span style="color: #0000FF;">},</span>
Line 1,694 ⟶ 1,700:
<span style="color: #008080;">function</span> <span style="color: #000000;">regression</span><span style="color: #0000FF;">()</span>
<span style="color: #004080;">atom</span> <span style="color: #7060A8000000;">xm</span> <span style="color: #0000FF;">=</span> <span style="color: #000000;">0</span><span style="color: #0000FF;">,</span> <span style="color: #7060A8000000;">ym</span> <span style="color: #0000FF;">=</span> <span style="color: #000000;">0</span><span style="color: #0000FF;">,</span> <span style="color: #000000;">x2m</span> <span style="color: #0000FF;">=</span> <span style="color: #000000;">0</span><span style="color: #0000FF;">,</span> <span style="color: #000000;">x3m</span> <span style="color: #0000FF;">=</span> <span style="color: #000000;">0</span><span style="color: #0000FF;">,</span> <span style="color: #000000;">x4m</span> <span style="color: #0000FF;">=</span> <span style="color: #000000;">0</span><span style="color: #0000FF;">,</span> <span style="color: #000000;">xym</span> <span style="color: #0000FF;">=</span> <span style="color: #000000;">0</span><span style="color: #0000FF;">,</span> <span style="color: #000000;">x2ym</span> <span style="color: #0000FF;">=</span> <span style="color: #000000;">0</span>
<span style="color: #008080;">for</span> <span style="color: #000000;">i</span><span style="color: #0000FF;">=</span><span style="color: #000000;">1</span> <span style="color: #008080;">to</span> <span style="color: #000000;">n</span> <span style="color: #008080;">do</span>
<span style="color: #004080;">atom</span> <span style="color: #000000;">xi</span> <span style="color: #0000FF;">=</span> <span style="color: #000000;">x</span><span style="color: #0000FF;">[</span><span style="color: #000000;">i</span><span style="color: #0000FF;">],</span>
<span style="color: #000000;">yi</span> <span style="color: #0000FF;">=</span> <span style="color: #000000;">y</span><span style="color: #0000FF;">[</span><span style="color: #000000;">i</span><span style="color: #0000FF;">]</span>
<span style="color: #7060A8000000;">xm</span> <span style="color: #0000FF;">+=</span> <span style="color: #000000;">xi</span>
<span style="color: #7060A8000000;">ym</span> <span style="color: #0000FF;">+=</span> <span style="color: #000000;">yi</span>
<span style="color: #000000;">x2m</span> <span style="color: #0000FF;">+=</span> <span style="color: #7060A8;">power</span><span style="color: #0000FF;">(</span><span style="color: #000000;">xi</span><span style="color: #0000FF;">,</span><span style="color: #000000;">2</span><span style="color: #0000FF;">)</span>
<span style="color: #000000;">x3m</span> <span style="color: #0000FF;">+=</span> <span style="color: #7060A8;">power</span><span style="color: #0000FF;">(</span><span style="color: #000000;">xi</span><span style="color: #0000FF;">,</span><span style="color: #000000;">3</span><span style="color: #0000FF;">)</span>
Line 1,706 ⟶ 1,712:
<span style="color: #000000;">x2ym</span> <span style="color: #0000FF;">+=</span> <span style="color: #7060A8;">power</span><span style="color: #0000FF;">(</span><span style="color: #000000;">xi</span><span style="color: #0000FF;">,</span><span style="color: #000000;">2</span><span style="color: #0000FF;">)*</span><span style="color: #000000;">yi</span>
<span style="color: #008080;">end</span> <span style="color: #008080;">for</span>
<span style="color: #7060A8000000;">xm</span> <span style="color: #0000FF;">/=</span> <span style="color: #000000;">n</span>
<span style="color: #7060A8000000;">ym</span> <span style="color: #0000FF;">/=</span> <span style="color: #000000;">n</span>
<span style="color: #000000;">x2m</span> <span style="color: #0000FF;">/=</span> <span style="color: #000000;">n</span>
<span style="color: #000000;">x3m</span> <span style="color: #0000FF;">/=</span> <span style="color: #000000;">n</span>
Line 1,713 ⟶ 1,719:
<span style="color: #000000;">xym</span> <span style="color: #0000FF;">/=</span> <span style="color: #000000;">n</span>
<span style="color: #000000;">x2ym</span> <span style="color: #0000FF;">/=</span> <span style="color: #000000;">n</span>
<span style="color: #004080;">atom</span> <span style="color: #000000;">Sxx</span> <span style="color: #0000FF;">=</span> <span style="color: #000000;">x2m</span><span style="color: #0000FF;">-</span><span style="color: #7060A8;">power</span><span style="color: #0000FF;">(</span><span style="color: #7060A8000000;">xm</span><span style="color: #0000FF;">,</span><span style="color: #000000;">2</span><span style="color: #0000FF;">),</span>
<span style="color: #000000;">Sxy</span> <span style="color: #0000FF;">=</span> <span style="color: #000000;">xym</span><span style="color: #0000FF;">-</span><span style="color: #7060A8000000;">xm</span><span style="color: #0000FF;">*</span><span style="color: #7060A8000000;">ym</span><span style="color: #0000FF;">,</span>
<span style="color: #000000;">Sxx2</span> <span style="color: #0000FF;">=</span> <span style="color: #000000;">x3m</span><span style="color: #0000FF;">-</span><span style="color: #7060A8000000;">xm</span><span style="color: #0000FF;">*</span><span style="color: #000000;">x2m</span><span style="color: #0000FF;">,</span>
<span style="color: #000000;">Sx2x2</span> <span style="color: #0000FF;">=</span> <span style="color: #000000;">x4m</span><span style="color: #0000FF;">-</span><span style="color: #7060A8;">power</span><span style="color: #0000FF;">(</span><span style="color: #000000;">x2m</span><span style="color: #0000FF;">,</span><span style="color: #000000;">2</span><span style="color: #0000FF;">),</span>
<span style="color: #000000;">Sx2y</span> <span style="color: #0000FF;">=</span> <span style="color: #000000;">x2ym</span><span style="color: #0000FF;">-</span><span style="color: #000000;">x2m</span><span style="color: #0000FF;">*</span><span style="color: #7060A8000000;">ym</span><span style="color: #0000FF;">,</span>
<span style="color: #000000;">B</span> <span style="color: #0000FF;">=</span> <span style="color: #0000FF;">(</span><span style="color: #000000;">Sxy</span><span style="color: #0000FF;">*</span><span style="color: #000000;">Sx2x2</span><span style="color: #0000FF;">-</span><span style="color: #000000;">Sx2y</span><span style="color: #0000FF;">*</span><span style="color: #000000;">Sxx2</span><span style="color: #0000FF;">)/(</span><span style="color: #000000;">Sxx</span><span style="color: #0000FF;">*</span><span style="color: #000000;">Sx2x2</span><span style="color: #0000FF;">-</span><span style="color: #7060A8;">power</span><span style="color: #0000FF;">(</span><span style="color: #000000;">Sxx2</span><span style="color: #0000FF;">,</span><span style="color: #000000;">2</span><span style="color: #0000FF;">)),</span>
<span style="color: #000000;">C</span> <span style="color: #0000FF;">=</span> <span style="color: #0000FF;">(</span><span style="color: #000000;">Sx2y</span><span style="color: #0000FF;">*</span><span style="color: #000000;">Sxx</span><span style="color: #0000FF;">-</span><span style="color: #000000;">Sxy</span><span style="color: #0000FF;">*</span><span style="color: #000000;">Sxx2</span><span style="color: #0000FF;">)/(</span><span style="color: #000000;">Sxx</span><span style="color: #0000FF;">*</span><span style="color: #000000;">Sx2x2</span><span style="color: #0000FF;">-</span><span style="color: #7060A8;">power</span><span style="color: #0000FF;">(</span><span style="color: #000000;">Sxx2</span><span style="color: #0000FF;">,</span><span style="color: #000000;">2</span><span style="color: #0000FF;">)),</span>
<span style="color: #000000;">A</span> <span style="color: #0000FF;">=</span> <span style="color: #7060A8000000;">ym</span><span style="color: #0000FF;">-</span><span style="color: #000000;">B</span><span style="color: #0000FF;">*</span><span style="color: #7060A8000000;">xm</span><span style="color: #0000FF;">-</span><span style="color: #000000;">C</span><span style="color: #0000FF;">*</span><span style="color: #000000;">x2m</span>
<span style="color: #008080;">return</span> <span style="color: #0000FF;">{</span><span style="color: #000000;">C</span><span style="color: #0000FF;">,</span><span style="color: #000000;">B</span><span style="color: #0000FF;">,</span><span style="color: #000000;">A</span><span style="color: #0000FF;">}</span>
<span style="color: #008080;">end</span> <span style="color: #008080;">function</span>
Line 1,735 ⟶ 1,741:
<span style="color: #7060A8;">printf</span><span style="color: #0000FF;">(</span><span style="color: #000000;">1</span><span style="color: #0000FF;">,</span><span style="color: #008000;">" %2d %3d %3g\n"</span><span style="color: #0000FF;">,{</span><span style="color: #000000;">x</span><span style="color: #0000FF;">[</span><span style="color: #000000;">i</span><span style="color: #0000FF;">],</span><span style="color: #000000;">y</span><span style="color: #0000FF;">[</span><span style="color: #000000;">i</span><span style="color: #0000FF;">],</span><span style="color: #000000;">f</span><span style="color: #0000FF;">(</span><span style="color: #000000;">x</span><span style="color: #0000FF;">[</span><span style="color: #000000;">i</span><span style="color: #0000FF;">])})</span>
<span style="color: #008080;">end</span> <span style="color: #008080;">for</span>
<span style="color: #000080;font-style:italic;">-- And a simple plot (re-using x,y from above)</span>
<span style="color: #7060A8008080;">IupShowinclude</span> <span style="color: #0000FF000000;">(pGUI</span><span style="color: #0000000000FF;">dlg.</span><span style="color: #0000FF000000;">)e</span>▼
<span style="color: #008080;">ifinclude</span> <span style="color: #7060A8000000;">platformIupGraph</span><span style="color: #0000FF;">()!=.</span><span style="color: #004600000000;">JSe</span> <span style="color: #008080;">then</span> ▼
▼
<span style="color: #7060A8008080;">IupPlotBeginfunction</span> <span style="color: #0000FF000000;">(get_data</span><span style="color: #0000000000FF;">plot(</span><span style="color: #0000FF004080;">,Ihandle</span> <span style="color: #000000;">0graph</span><span style="color: #0000FF;">)</span>▼
<span style="color: #7060A8004080;">IupPlotAddinteger</span> <span style="color: #0000FF;">({</span><span style="color: #000000;">plotw</span><span style="color: #0000FF;">,</span> <span style="color: #000000;">xh</span><span style="color: #0000FF;">[}</span> <span style="color: #0000000000FF;">i=</span> <span style="color: #0000FF7060A8;">],IupGetIntInt</span><span style="color: #0000FF;">(</span><span style="color: #000000;">ygraph</span><span style="color: #0000FF;">[,</span><span style="color: #000000008000;">i"DRAWSIZE"</span><span style="color: #0000FF;">])</span>▼
<span style="color: #008080;">constant</span> <span style="color: #000000;">x</span> <span style="color: #0000FF7060A8;">=IupSetInt</span> <span style="color: #0000FF;">{(</span><span style="color: #000000;">0graph</span><span style="color: #0000FF;">,</span><span style="color: #000000008000;">1"YTICK"</span><span style="color: #0000FF;">,</span><span style="color: #000000008080;">2iff</span><span style="color: #0000FF;">,(</span><span style="color: #000000;">3h</span><span style="color: #0000FF;">,<</span><span style="color: #000000;">4240</span><span style="color: #0000FF;">,?</span><span style="color: #000000008080;">5iff</span><span style="color: #0000FF;">,(</span><span style="color: #000000;">6h</span><span style="color: #0000FF;">,<</span><span style="color: #000000;">7150</span><span style="color: #0000FF;">,?</span><span style="color: #000000;">880</span><span style="color: #0000FF;">,:</span><span style="color: #000000;">940</span><span style="color: #0000FF;">,):</span><span style="color: #000000;">1020</span><span style="color: #0000FF;">}))</span>▼
<span style="color: #004080;">Ihandle</span> <span style="color: #000000008080;">dlgreturn</span> <span style="color: #0000FF;">={{</span> <span style="color: #7060A8000000;">IupDialogx</span><span style="color: #0000FF;">(,</span><span style="color: #000000;">ploty</span><span style="color: #0000FF;">,</span><span style="color: #008000004600;">`TITLE="simple plot", RASTERSIZE=640x480`CD_RED</span><span style="color: #0000FF;">)}}</span>▼
<span style="color: #008080;">end</span> <span style="color: #008080;">forfunction</span>▼
<span style="color: #7060A8;">IupOpen</span><span style="color: #0000FF;">()</span>▼
<span style="color: #008080004080;">forIhandle</span> <span style="color: #000000;">igraph</span> <span style="color: #0000FF;">=</span> <span style="color: #000000;">1IupGraph</span> <span style="color: #0080800000FF;">to(</span> <span style="color: #7060A8000000;">lengthget_data</span><span style="color: #0000FF;">(,</span><span style="color: #000000008000;">x"RASTERSIZE=640x440"</span><span style="color: #0000FF;">)</span> <span style="color: #008080;">do</span>▼
<span style="color: #0000FF7060A8;">{}IupSetAttributes</span> <span style="color: #0000FF;">=(</span> <span style="color: #7060A8000000;">IupPlotEndgraph</span><span style="color: #0000FF;">(,</span><span style="color: #000000008000;">plot"XTICK=1,XMIN=0,XMAX=10"</span><span style="color: #0000FF;">)</span>▼
<span style="color: #0040807060A8;">IhandleIupSetAttributes</span> <span style="color: #000000;">plot</span> <span style="color: #0000FF;">=(</span> <span style="color: #7060A8000000;">IupPlotgraph</span><span style="color: #0000FF;">(,</span><span style="color: #008000;">"GRIDYTICK=YES20, MARGINLEFTYMIN=500, MARGINBOTTOMYMAX=40320"</span><span style="color: #0000FF;">)</span>▼
<span style="color: #004080;">Ihandle</span> <span style="color: #000000;">dlg</span> <span style="color: #0000FF;">=</span> <span style="color: #7060A8;">IupDialog</span><span style="color: #0000FF;">(</span><span style="color: #000000;">graph</span><span style="color: #0000FF;">,</span><span style="color: #008000;">`TITLE="simple plot"`</span><span style="color: #0000FF;">)</span>
<span style="color: #7060A8;">IupSetAttributes</span><span style="color: #0000FF;">(</span><span style="color: #000000;">dlg</span><span style="color: #0000FF;">,</span><span style="color: #008000;">"MINSIZE=245x150"</span><span style="color: #0000FF;">)</span>
<span style="color: #7060A8;">IupShow</span><span style="color: #0000FF;">(</span><span style="color: #7060A8000000;">IupClosedlg</span><span style="color: #0000FF;">()</span>▼
<span style="color: #008080;">if</span> <span style="color: #7060A8;">platform</span><span style="color: #0000FF;">()!=</span><span style="color: #004600;">JS</span> <span style="color: #008080;">then</span>
<span style="color: #7060A8;">IupMainLoop</span><span style="color: #0000FF;">()</span>▼
<span style="color: #0080807060A8;">endIupClose</span> <span style="color: #0080800000FF;">if()</span>▼
<span style="color: #008080;">end</span> <span style="color: #008080;">if</span>
<!--</lang>-->
{{out}}
A(plus a simple graphical plot, as per [[Polynomial_regression#Racket|Racket]]:)▼
<pre>
y=3x^2+2x+1
Line 1,753 ⟶ 1,783:
10 321 321
</pre>
===plot===
▲A simple plot, as per [[Polynomial_regression#Racket|Racket]]:
▲{{libheader|Phix/pGUI}}
<!--<lang Phix>-->
▲ <span style="color: #008080;">include</span> <span style="color: #000000;">pGUI</span><span style="color: #0000FF;">.</span><span style="color: #000000;">e</span>
▲
▲ <span style="color: #008080;">constant</span> <span style="color: #000000;">x</span> <span style="color: #0000FF;">=</span> <span style="color: #0000FF;">{</span><span style="color: #000000;">0</span><span style="color: #0000FF;">,</span><span style="color: #000000;">1</span><span style="color: #0000FF;">,</span><span style="color: #000000;">2</span><span style="color: #0000FF;">,</span><span style="color: #000000;">3</span><span style="color: #0000FF;">,</span><span style="color: #000000;">4</span><span style="color: #0000FF;">,</span><span style="color: #000000;">5</span><span style="color: #0000FF;">,</span><span style="color: #000000;">6</span><span style="color: #0000FF;">,</span><span style="color: #000000;">7</span><span style="color: #0000FF;">,</span><span style="color: #000000;">8</span><span style="color: #0000FF;">,</span><span style="color: #000000;">9</span><span style="color: #0000FF;">,</span><span style="color: #000000;">10</span><span style="color: #0000FF;">}</span>
<span style="color: #008080;">constant</span> <span style="color: #000000;">y</span> <span style="color: #0000FF;">=</span> <span style="color: #0000FF;">{</span><span style="color: #000000;">1</span><span style="color: #0000FF;">,</span><span style="color: #000000;">6</span><span style="color: #0000FF;">,</span><span style="color: #000000;">17</span><span style="color: #0000FF;">,</span><span style="color: #000000;">34</span><span style="color: #0000FF;">,</span><span style="color: #000000;">57</span><span style="color: #0000FF;">,</span><span style="color: #000000;">86</span><span style="color: #0000FF;">,</span><span style="color: #000000;">121</span><span style="color: #0000FF;">,</span><span style="color: #000000;">162</span><span style="color: #0000FF;">,</span><span style="color: #000000;">209</span><span style="color: #0000FF;">,</span><span style="color: #000000;">262</span><span style="color: #0000FF;">,</span><span style="color: #000000;">321</span><span style="color: #0000FF;">}</span>
▲ <span style="color: #7060A8;">IupOpen</span><span style="color: #0000FF;">()</span>
▲ <span style="color: #004080;">Ihandle</span> <span style="color: #000000;">plot</span> <span style="color: #0000FF;">=</span> <span style="color: #7060A8;">IupPlot</span><span style="color: #0000FF;">(</span><span style="color: #008000;">"GRID=YES, MARGINLEFT=50, MARGINBOTTOM=40"</span><span style="color: #0000FF;">)</span>
▲ <span style="color: #000080;font-style:italic;">-- (just add ", AXS_YSCALE=LOG10" for a nice log scale)</span>
▲ <span style="color: #7060A8;">IupPlotBegin</span><span style="color: #0000FF;">(</span><span style="color: #000000;">plot</span><span style="color: #0000FF;">,</span> <span style="color: #000000;">0</span><span style="color: #0000FF;">)</span>
▲ <span style="color: #008080;">for</span> <span style="color: #000000;">i</span><span style="color: #0000FF;">=</span><span style="color: #000000;">1</span> <span style="color: #008080;">to</span> <span style="color: #7060A8;">length</span><span style="color: #0000FF;">(</span><span style="color: #000000;">x</span><span style="color: #0000FF;">)</span> <span style="color: #008080;">do</span>
▲ <span style="color: #7060A8;">IupPlotAdd</span><span style="color: #0000FF;">(</span><span style="color: #000000;">plot</span><span style="color: #0000FF;">,</span> <span style="color: #000000;">x</span><span style="color: #0000FF;">[</span><span style="color: #000000;">i</span><span style="color: #0000FF;">],</span> <span style="color: #000000;">y</span><span style="color: #0000FF;">[</span><span style="color: #000000;">i</span><span style="color: #0000FF;">])</span>
▲ <span style="color: #008080;">end</span> <span style="color: #008080;">for</span>
▲ <span style="color: #0000FF;">{}</span> <span style="color: #0000FF;">=</span> <span style="color: #7060A8;">IupPlotEnd</span><span style="color: #0000FF;">(</span><span style="color: #000000;">plot</span><span style="color: #0000FF;">)</span>
▲ <span style="color: #004080;">Ihandle</span> <span style="color: #000000;">dlg</span> <span style="color: #0000FF;">=</span> <span style="color: #7060A8;">IupDialog</span><span style="color: #0000FF;">(</span><span style="color: #000000;">plot</span><span style="color: #0000FF;">,</span><span style="color: #008000;">`TITLE="simple plot", RASTERSIZE=640x480`</span><span style="color: #0000FF;">)</span>
▲ <span style="color: #7060A8;">IupShow</span><span style="color: #0000FF;">(</span><span style="color: #000000;">dlg</span><span style="color: #0000FF;">)</span>
▲ <span style="color: #008080;">if</span> <span style="color: #7060A8;">platform</span><span style="color: #0000FF;">()!=</span><span style="color: #004600;">JS</span> <span style="color: #008080;">then</span>
▲ <span style="color: #7060A8;">IupMainLoop</span><span style="color: #0000FF;">()</span>
▲ <span style="color: #7060A8;">IupClose</span><span style="color: #0000FF;">()</span>
▲ <span style="color: #008080;">end</span> <span style="color: #008080;">if</span>
<!--</lang>-->
 
=={{header|PowerShell}}==
Petelomax
7,794

edits

Retrieved from "https://rosettacode.org/wiki/Special:MobileDiff/53620"