<html>
<head>
<meta http-equiv="content-type" content="text/html; charset=utf-8">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<p>Dear Scilab users</p>
<p>I've tried to use histplot, but I probably just don't understand
what it does or what it is supposed to do. I hope you can help me
do this in a smarter way. Here's to code which does the correct
stuff - but as you can see, I have to "manually" count what goes
into each histogram category - and afterwards I also have to
manually "paint" the histplot. It's rather complicated for being
such a simple task.<span
style="color:rgb(100,174,100);font-style:italic;"><br>
</span></p>
<p><span style="color:rgb(100,174,100);font-style:italic;">//
sample.sce</span><span style="color:rgb(0,0,0);"><br>
</span></p>
<p><span style="color:rgb(0,0,0);">Bl</span> <span
style="color:rgb(92,92,92);">=</span> <span
style="color:rgb(74,85,219);">[</span></p>
<p><span style="color:rgb(74,85,219);"> </span><span
style="color:rgb(188,143,143);">7.007</span> <span
style="color:rgb(188,143,143);">7.065</span> <span
style="color:rgb(188,143,143);">7.047</span> <span
style="color:rgb(188,143,143);">6.999</span> <span
style="color:rgb(188,143,143);">6.999</span> <span
style="color:rgb(188,143,143);">7.025</span> <span
style="color:rgb(188,143,143);">6.998</span> <span
style="color:rgb(188,143,143);">7.075</span> <span
style="color:rgb(188,143,143);">7.035</span> <span
style="color:rgb(255,170,0);">..</span><br>
</p>
<p> <span style="color:rgb(188,143,143);">7.037</span> <span
style="color:rgb(188,143,143);">7.071</span> <span
style="color:rgb(188,143,143);">7.031</span> <span
style="color:rgb(188,143,143);">7.069</span> <span
style="color:rgb(188,143,143);">7.032</span> <span
style="color:rgb(188,143,143);">7.02</span> <span
style="color:rgb(188,143,143);">7.047</span> <span
style="color:rgb(188,143,143);">7.025</span> <span
style="color:rgb(188,143,143);">7.026</span> <span
style="color:rgb(74,85,219);">]</span><span
style="color:rgb(0,0,0);">;</span></p>
<p><span style="color:rgb(0,0,0);"></span><span
style="color:rgb(0,0,0);">avgbl</span> <span
style="color:rgb(92,92,92);">=</span> <span
style="color:rgb(174,92,176);text-decoration:underline;">mean</span><span
style="color:rgb(74,85,219);">(</span><span
style="color:rgb(0,0,0);">Bl</span><span
style="color:rgb(74,85,219);">)</span><span
style="color:rgb(0,0,0);">;</span><span
style="color:rgb(74,85,219);"><br>
</span></p>
<p><span style="color:rgb(74,85,219);">[</span><span
style="color:rgb(0,0,0);">nearest</span><span
style="color:rgb(0,0,0);">,</span><span
style="color:rgb(0,0,0);">i</span><span
style="color:rgb(74,85,219);">]</span> <span
style="color:rgb(92,92,92);">=</span> <span
style="color:rgb(50,185,185);">min</span><span
style="color:rgb(74,85,219);">(</span><span
style="color:rgb(50,185,185);">abs</span><span
style="color:rgb(74,85,219);">(</span><span
style="color:rgb(0,0,0);">Bl</span><span
style="color:rgb(92,92,92);">-</span><span
style="color:rgb(0,0,0);">avgbl</span><span
style="color:rgb(74,85,219);">)</span><span
style="color:rgb(74,85,219);">)</span><span
style="color:rgb(0,0,0);">;</span> <span
style="color:rgb(100,174,100);font-style:italic;">// </span><span
style="color:rgb(100,174,100);font-style:italic;">'</span><span
style="color:rgb(100,174,100);font-style:italic;">nearest</span><span
style="color:rgb(100,174,100);font-style:italic;">'</span><span
style="color:rgb(100,174,100);font-style:italic;"> is
unimportant, we need i</span><span style="color:rgb(0,0,0);"><br>
</span></p>
<p><span style="color:rgb(0,0,0);">nearest</span> <span
style="color:rgb(92,92,92);">=</span> <span
style="color:rgb(0,0,0);">Bl</span><span
style="color:rgb(74,85,219);">(</span><span
style="color:rgb(0,0,0);">i</span><span
style="color:rgb(74,85,219);">)</span><span
style="color:rgb(0,0,0);">;</span> <span
style="color:rgb(100,174,100);font-style:italic;">// This is our
</span><span style="color:rgb(100,174,100);font-style:italic;">'</span><span
style="color:rgb(100,174,100);font-style:italic;">Master</span><span
style="color:rgb(100,174,100);font-style:italic;">'</span><span
style="color:rgb(100,174,100);font-style:italic;"> for the BIAS
study</span><span style="color:rgb(0,0,0);"><br>
</span></p>
<p><span style="color:rgb(0,0,0);">Bl_diff</span> <span
style="color:rgb(92,92,92);">=</span> <span
style="color:rgb(0,0,0);">Bl</span> <span
style="color:rgb(92,92,92);">-</span> <span
style="color:rgb(0,0,0);">nearest</span><span
style="color:rgb(0,0,0);">;</span><span
style="color:rgb(0,0,0);"><br>
</span></p>
<p><span style="color:rgb(0,0,0);">minBl_diff</span> <span
style="color:rgb(92,92,92);">=</span> <span
style="color:rgb(50,185,185);">floor</span><span
style="color:rgb(74,85,219);">(</span><span
style="color:rgb(188,143,143);">100</span><span
style="color:rgb(92,92,92);">*</span><span
style="color:rgb(50,185,185);">min</span><span
style="color:rgb(74,85,219);">(</span><span
style="color:rgb(0,0,0);">Bl_diff</span><span
style="color:rgb(74,85,219);">)</span><span
style="color:rgb(74,85,219);">)</span><span
style="color:rgb(92,92,92);">/</span><span
style="color:rgb(188,143,143);">100</span> <span
style="color:rgb(92,92,92);">-</span> <span
style="color:rgb(188,143,143);">0.005</span><span
style="color:rgb(0,0,0);">;</span><span
style="color:rgb(0,0,0);"><br>
</span></p>
<p><span style="color:rgb(0,0,0);">maxBl_diff</span> <span
style="color:rgb(92,92,92);">=</span> <span
style="color:rgb(50,185,185);">ceil</span><span
style="color:rgb(74,85,219);">(</span><span
style="color:rgb(188,143,143);">100</span><span
style="color:rgb(92,92,92);">*</span><span
style="color:rgb(50,185,185);">max</span><span
style="color:rgb(74,85,219);">(</span><span
style="color:rgb(0,0,0);">Bl_diff</span><span
style="color:rgb(74,85,219);">)</span><span
style="color:rgb(74,85,219);">)</span><span
style="color:rgb(92,92,92);">/</span><span
style="color:rgb(188,143,143);">100</span> <span
style="color:rgb(92,92,92);">+</span> <span
style="color:rgb(188,143,143);">0.005</span><span
style="color:rgb(0,0,0);">;</span><span
style="color:rgb(0,0,0);"><br>
</span></p>
<p><span style="color:rgb(0,0,0);">bins</span> <span
style="color:rgb(92,92,92);">=</span> <span
style="color:rgb(74,85,219);">(</span><span
style="color:rgb(0,0,0);">maxBl_diff</span> <span
style="color:rgb(92,92,92);">-</span> <span
style="color:rgb(0,0,0);">minBl_diff</span><span
style="color:rgb(74,85,219);">)</span><span
style="color:rgb(92,92,92);">*</span><span
style="color:rgb(188,143,143);">100</span><span
style="color:rgb(0,0,0);">;</span> <span
style="color:rgb(100,174,100);font-style:italic;">// we create 9
bins</span><span style="color:rgb(0,0,0);"><br>
</span></p>
<p><span style="color:rgb(0,0,0);">dividers</span> <span
style="color:rgb(92,92,92);">=</span> <span
style="color:rgb(174,92,176);text-decoration:underline;">linspace</span><span
style="color:rgb(74,85,219);">(</span><span
style="color:rgb(0,0,0);">minBl_diff</span><span
style="color:rgb(0,0,0);">,</span><span
style="color:rgb(0,0,0);">maxBl_diff</span><span
style="color:rgb(0,0,0);">,</span><span
style="color:rgb(0,0,0);">bins</span><span
style="color:rgb(92,92,92);">+</span><span
style="color:rgb(188,143,143);">1</span><span
style="color:rgb(74,85,219);">)</span><span
style="color:rgb(0,0,0);">;</span> <span
style="color:rgb(100,174,100);font-style:italic;">// 10 dividers</span><span
style="color:rgb(0,0,0);"><br>
</span></p>
<p><span style="color:rgb(0,0,0);">bin_vec</span><span
style="color:rgb(74,85,219);">(</span><span
style="color:rgb(188,143,143);">1</span><span
style="color:rgb(74,85,219);">)</span> <span
style="color:rgb(92,92,92);">=</span> <span
style="color:rgb(50,185,185);">sum</span><span
style="color:rgb(74,85,219);">(</span><span
style="color:rgb(0,0,0);">Bl_diff</span> <span
style="color:rgb(92,92,92);"><</span> <span
style="color:rgb(0,0,0);">dividers</span><span
style="color:rgb(74,85,219);">(</span><span
style="color:rgb(188,143,143);">2</span><span
style="color:rgb(74,85,219);">)</span><span
style="color:rgb(74,85,219);">)</span><span
style="color:rgb(0,0,0);">;</span> <span
style="color:rgb(100,174,100);font-style:italic;">// start
initialize bin_vec</span><span style="color:rgb(160,32,240);"><br>
</span></p>
<p><span style="color:rgb(160,32,240);">for</span> <span
style="color:rgb(0,0,0);">i</span><span
style="color:rgb(92,92,92);">=</span><span
style="color:rgb(188,143,143);">2</span><span
style="color:rgb(255,170,0);">:</span><span
style="color:rgb(0,0,0);">bins</span> <span
style="color:rgb(160,32,240);">do</span><br>
</p>
<p> <span style="color:rgb(0,0,0);">bin_vec</span><span
style="color:rgb(74,85,219);">(</span><span
style="color:rgb(0,0,0);">i</span><span
style="color:rgb(74,85,219);">)</span> <span
style="color:rgb(92,92,92);">=</span> <span
style="color:rgb(50,185,185);">sum</span><span
style="color:rgb(74,85,219);">(</span><span
style="color:rgb(0,0,0);">Bl_diff</span> <span
style="color:rgb(92,92,92);"><</span> <span
style="color:rgb(0,0,0);">dividers</span><span
style="color:rgb(74,85,219);">(</span><span
style="color:rgb(0,0,0);">i</span><span
style="color:rgb(92,92,92);">+</span><span
style="color:rgb(188,143,143);">1</span><span
style="color:rgb(74,85,219);">)</span><span
style="color:rgb(74,85,219);">)</span> <span
style="color:rgb(92,92,92);">-</span> <span
style="color:rgb(50,185,185);">sum</span><span
style="color:rgb(74,85,219);">(</span><span
style="color:rgb(0,0,0);">Bl_diff</span> <span
style="color:rgb(92,92,92);"><</span> <span
style="color:rgb(0,0,0);">dividers</span><span
style="color:rgb(74,85,219);">(</span><span
style="color:rgb(0,0,0);">i</span><span
style="color:rgb(74,85,219);">)</span><span
style="color:rgb(74,85,219);">)</span><span
style="color:rgb(160,32,240);"><br>
</span></p>
<p><span style="color:rgb(160,32,240);">end</span> <span
style="color:rgb(100,174,100);font-style:italic;">// bin_vec
contains the correct frequency count</span><span
style="color:rgb(174,92,176);text-decoration:underline;"><br>
</span></p>
<p><span style="color:rgb(174,92,176);text-decoration:underline;">scf</span><span
style="color:rgb(74,85,219);">(</span><span
style="color:rgb(74,85,219);">)</span><span
style="color:rgb(0,0,0);">;</span><span
style="color:rgb(50,185,185);"><br>
</span></p>
<p><span style="color:rgb(50,185,185);">drawlater</span><span
style="color:rgb(74,85,219);">(</span><span
style="color:rgb(74,85,219);">)</span><span
style="color:rgb(0,0,0);">;</span><br>
</p>
<p> <span style="color:rgb(100,174,100);font-style:italic;">//
histplot(dividers,Bl_diff,style=2); // WRONG FREQUENCY COUNT</span><span
style="color:rgb(160,32,240);"><br>
</span></p>
<p><span style="color:rgb(160,32,240);">for</span> <span
style="color:rgb(0,0,0);">i</span><span
style="color:rgb(92,92,92);">=</span><span
style="color:rgb(188,143,143);">1</span><span
style="color:rgb(255,170,0);">:</span><span
style="color:rgb(0,0,0);">bins</span> <span
style="color:rgb(160,32,240);">do</span><br>
</p>
<p> <span style="color:rgb(0,0,0);">xi</span> <span
style="color:rgb(92,92,92);">=</span> <span
style="color:rgb(74,85,219);">(</span><span
style="color:rgb(0,0,0);">dividers</span><span
style="color:rgb(74,85,219);">(</span><span
style="color:rgb(0,0,0);">i</span><span
style="color:rgb(92,92,92);">+</span><span
style="color:rgb(188,143,143);">1</span><span
style="color:rgb(74,85,219);">)</span> <span
style="color:rgb(92,92,92);">+</span> <span
style="color:rgb(0,0,0);">dividers</span><span
style="color:rgb(74,85,219);">(</span><span
style="color:rgb(0,0,0);">i</span><span
style="color:rgb(74,85,219);">)</span><span
style="color:rgb(74,85,219);">)</span><span
style="color:rgb(92,92,92);">/</span><span
style="color:rgb(188,143,143);">2</span><span
style="color:rgb(0,0,0);">;</span> <span
style="color:rgb(100,174,100);font-style:italic;">// middle of
bin container</span><br>
</p>
<p> <span style="color:rgb(0,0,0);">xlo</span> <span
style="color:rgb(92,92,92);">=</span> <span
style="color:rgb(0,0,0);">xi</span> <span
style="color:rgb(92,92,92);">-</span> <span
style="color:rgb(188,143,143);">0.004</span><span
style="color:rgb(0,0,0);">;</span><br>
</p>
<p> <span style="color:rgb(0,0,0);">xhi</span> <span
style="color:rgb(92,92,92);">=</span> <span
style="color:rgb(0,0,0);">xi</span> <span
style="color:rgb(92,92,92);">+</span> <span
style="color:rgb(188,143,143);">0.004</span><span
style="color:rgb(0,0,0);">;</span><br>
</p>
<p> <span
style="color:rgb(174,92,176);text-decoration:underline;">plot</span><span
style="color:rgb(74,85,219);">(</span><span
style="color:rgb(74,85,219);">[</span><span
style="color:rgb(0,0,0);">xlo</span> <span
style="color:rgb(0,0,0);">xlo</span><span
style="color:rgb(74,85,219);">]</span><span
style="color:rgb(0,0,0);">,</span><span
style="color:rgb(74,85,219);">[</span><span
style="color:rgb(188,143,143);">0</span> <span
style="color:rgb(0,0,0);">bin_vec</span><span
style="color:rgb(74,85,219);">(</span><span
style="color:rgb(0,0,0);">i</span><span
style="color:rgb(74,85,219);">)</span><span
style="color:rgb(74,85,219);">]</span><span
style="color:rgb(0,0,0);">,</span><span
style="color:rgb(188,143,143);">'</span><span
style="color:rgb(188,143,143);">-k</span><span
style="color:rgb(188,143,143);">'</span><span
style="color:rgb(74,85,219);">)</span><span
style="color:rgb(0,0,0);">;</span><br>
</p>
<p> <span
style="color:rgb(174,92,176);text-decoration:underline;">plot</span><span
style="color:rgb(74,85,219);">(</span><span
style="color:rgb(74,85,219);">[</span><span
style="color:rgb(0,0,0);">xlo</span> <span
style="color:rgb(0,0,0);">xhi</span><span
style="color:rgb(74,85,219);">]</span><span
style="color:rgb(0,0,0);">,</span><span
style="color:rgb(74,85,219);">[</span><span
style="color:rgb(0,0,0);">bin_vec</span><span
style="color:rgb(74,85,219);">(</span><span
style="color:rgb(0,0,0);">i</span><span
style="color:rgb(74,85,219);">)</span> <span
style="color:rgb(0,0,0);">bin_vec</span><span
style="color:rgb(74,85,219);">(</span><span
style="color:rgb(0,0,0);">i</span><span
style="color:rgb(74,85,219);">)</span><span
style="color:rgb(74,85,219);">]</span><span
style="color:rgb(0,0,0);">,</span><span
style="color:rgb(188,143,143);">'</span><span
style="color:rgb(188,143,143);">-k</span><span
style="color:rgb(188,143,143);">'</span><span
style="color:rgb(74,85,219);">)</span><span
style="color:rgb(0,0,0);">;</span><br>
</p>
<p> <span
style="color:rgb(174,92,176);text-decoration:underline;">plot</span><span
style="color:rgb(74,85,219);">(</span><span
style="color:rgb(74,85,219);">[</span><span
style="color:rgb(0,0,0);">xhi</span> <span
style="color:rgb(0,0,0);">xhi</span><span
style="color:rgb(74,85,219);">]</span><span
style="color:rgb(0,0,0);">,</span><span
style="color:rgb(74,85,219);">[</span><span
style="color:rgb(188,143,143);">0</span> <span
style="color:rgb(0,0,0);">bin_vec</span><span
style="color:rgb(74,85,219);">(</span><span
style="color:rgb(0,0,0);">i</span><span
style="color:rgb(74,85,219);">)</span><span
style="color:rgb(74,85,219);">]</span><span
style="color:rgb(0,0,0);">,</span><span
style="color:rgb(188,143,143);">'</span><span
style="color:rgb(188,143,143);">-k</span><span
style="color:rgb(188,143,143);">'</span><span
style="color:rgb(74,85,219);">)</span><span
style="color:rgb(0,0,0);">;</span><br>
</p>
<p> <span style="color:rgb(50,185,185);">xfpoly</span><span
style="color:rgb(74,85,219);">(</span><span
style="color:rgb(74,85,219);">[</span><span
style="color:rgb(0,0,0);">xlo</span> <span
style="color:rgb(0,0,0);">xlo</span> <span
style="color:rgb(0,0,0);">xhi</span> <span
style="color:rgb(0,0,0);">xhi</span> <span
style="color:rgb(0,0,0);">xlo</span><span
style="color:rgb(74,85,219);">]</span><span
style="color:rgb(0,0,0);">,</span><span
style="color:rgb(74,85,219);">[</span><span
style="color:rgb(188,143,143);">0</span> <span
style="color:rgb(0,0,0);">bin_vec</span><span
style="color:rgb(74,85,219);">(</span><span
style="color:rgb(0,0,0);">i</span><span
style="color:rgb(74,85,219);">)</span> <span
style="color:rgb(0,0,0);">bin_vec</span><span
style="color:rgb(74,85,219);">(</span><span
style="color:rgb(0,0,0);">i</span><span
style="color:rgb(74,85,219);">)</span> <span
style="color:rgb(188,143,143);">0</span> <span
style="color:rgb(188,143,143);">0</span><span
style="color:rgb(74,85,219);">]</span><span
style="color:rgb(0,0,0);">,</span><span
style="color:rgb(50,185,185);">color</span><span
style="color:rgb(74,85,219);">(</span><span
style="color:rgb(188,143,143);">"</span><span
style="color:rgb(188,143,143);">grey</span><span
style="color:rgb(188,143,143);">"</span><span
style="color:rgb(74,85,219);">)</span><span
style="color:rgb(74,85,219);">)</span><span
style="color:rgb(0,0,0);">;</span><span
style="color:rgb(160,32,240);"><br>
</span></p>
<p><span style="color:rgb(160,32,240);">end</span><span
style="color:rgb(50,185,185);"><br>
</span></p>
<p><span style="color:rgb(50,185,185);">xtitle</span><span
style="color:rgb(74,85,219);">(</span><span
style="color:rgb(188,143,143);">"</span><span
style="color:rgb(188,143,143);">"</span><span
style="color:rgb(0,0,0);">,</span><span
style="color:rgb(188,143,143);">"</span><span
style="color:rgb(188,143,143);">Bias study, $Bl$ - mean($Bl$)</span><span
style="color:rgb(188,143,143);">"</span><span
style="color:rgb(0,0,0);">,</span><span
style="color:rgb(188,143,143);">"Count</span><span
style="color:rgb(188,143,143);">"</span><span
style="color:rgb(74,85,219);">)</span><span
style="color:rgb(0,0,0);">;</span><span
style="color:rgb(50,185,185);"><br>
</span></p>
<p><span style="color:rgb(50,185,185);">drawnow</span><span
style="color:rgb(74,85,219);">(</span><span
style="color:rgb(74,85,219);">)</span><span
style="color:rgb(0,0,0);">;</span></p>
<p>Best regards,</p>
<p>Claus<br>
</p>
</body>
</html>