\( \DeclareMathOperator{\abs}{abs} \newcommand{\ensuremath}[1]{\mbox{$#1$}} \)

_j_msc_job1_ex.wxmx

(%i1) kill(all);
\[\tag{\% o0} \mathit{done}\]
(%i1) load("stats")$
\[\mbox{}\\\mbox{0 errors, 0 warnings}\]
(%i2) mm:[4.97, 14.9, 6.69, 21.0, 7.94, 11.9, 2.43, 4.43, 0.726, 25.1, 10.5, 15.7, 5.61, 21.2, 18.5, 18.7, 4.66, 28.1, 21.0, 5.06, 15.7, 0.55, 13.9, 1.72, 9.59, 8.92, 16.7, 15.3, 20.5, 12.4];
\[\tag{mm}[4.97,14.9,6.69,21.0,7.94,11.9,2.43,4.43,0.726,25.1,10.5,15.7,5.61,21.2,18.5,18.7,4.66,28.1,21.0,5.06,15.7,0.55,13.9,1.72,9.59,8.92,16.7,15.3,20.5,12.4]\]
(%i3) test_mean(mm,'conflevel=0.95);
\[\tag{\% o3} \begin{array}{c}MEAN TEST\\ \mathit{mean\_ estimate}=12.14653333333333\\ \mathit{conf\_ level}=0.95\\ \mathit{conf\_ interval}=[9.307794616017414,14.98527205064925]\\ \mathit{method}=Exact t-test. Unknown variance.\\ \mathit{hypotheses}=H0: mean = 0 , H1: mean \neq 0\\ \mathit{statistic}=8.751228096812154\\ \mathit{distribution}=[{{\mathit{student}}_t},29]\\ {p_{\mathit{value}}}=1.242122849021143 {{10}^{-9}}\end{array}\]
(%i4) mean(mm),numer;
\[\tag{\% o4} 12.14653333333333\]
(%i5) A:lmax(mm);
\[\tag{A}28.1\]
(%i6) B:lmin(mm);
\[\tag{B}0.55\]
(%i7) Rang:A-B;
\[\tag{Rang}27.55\]
(%i8) St:Rang/6;
\[\tag{St}4.591666666666667\]
(%i10) s[1]:0; for i:1 step 1 while i<30 do ( if mm[i]<B+St then s[1]:s[1]+1 else s[1]:s[1]);
\[\tag{\% o9} 0\] \[\tag{\% o10} \mathit{done}\]
(%i11) for j:1 step 1 while j<7 do (s[j]:0, for i:1 step 1 while i<30 do ( if mm[i]<B+j*St then s[j]:s[j]+1 else s[j]:s[j]));
\[\tag{\% o11} \mathit{done}\]
(%i13) n[1]:s[1]; for k:2 step 1 while k<7 do (n[k]:s[k]-s[k-1]);
\[\tag{\% o12} 8\] \[\tag{\% o13} \mathit{done}\]
(%i14) ss:n[5];
\[\tag{ss}4\]
(%i16) xy: [[1, n[1]], [2, n[2]], [3, n[4]], [4, n[4]], [5, n[5]], [6, n[6]]];
plot2d([discrete, xy]);
\[\tag{xy}[[1,8],[2,5],[3,7],[4,7],[5,4],[6,1]]\] \[\tag{\% o16} [C:/Users/y/AppData/Local/Temp/maxout4640.gnuplot]\]
(%i22) B1:B+St;  B2:B+2*St; B3:B+3*St; B4:B+4*St; B5:B+5*St; B6:B+6*St;
\[\tag{B1}5.141666666666667\] \[\tag{B2}9.733333333333334\] \[\tag{B3}14.325\] \[\tag{B4}18.91666666666667\] \[\tag{B5}23.50833333333334\] \[\tag{B6}28.1\]
(%i23) draw2d(
   key          = "Group A",
 fill_color   = blue,
 fill_density = 0.2,
 bars([B1,8,2],[B2,5,2],[B3,7,2],[B4,7,2],[B5,4,2],[B6,1,2]),
 xaxis = true ) $
Гістограма

Риc. 1.

(%i24) smin(mm);
\[\tag{\% o24} 0.55\]
(%i25) smax(mm);
\[\tag{\% o25} 28.1\]
(%i26) range(mm);
\[\tag{\% o26} 27.55\]
(%i27) var1(mm);
\[\tag{\% o27} 57.79469018850576\]
(%i28) std1(mm);
\[\tag{\% o28} 7.60228190667156\]
(%i29) histogram(mm, nclasses = 6, fill_color= blue, fill_density = 0.5 , line_width   = 4);
\[\tag{\% o29} [\operatorname{gr2d}\left( \mathit{bars}\right) ]\]
Гістограма

Риc. 2.

Created with wxMaxima.