MTB > name c1 'weights' MTB > set c1 DATA> 32.6 33.1 32.1 32.7 31.9 31.8 32.4 32.3 32.8 33.3 32.7 32.6 DATA> end MTB > let k1=mean(c1) MTB > let k2=std(c1) MTB > let k3=count(c1) MTB > name c2 'chisq' c3 's2' c4 'mu' MTB > let k4=count(c1)-1 MTB > rand 1000 'chisq'; SUBC> chisq k4. MTB > let 's2'=k4*k2**2/'chisq' MTB > rand 1000 'mu'; SUBC> normal 0 1. MTB > let 'mu'=mean(c1)+sqrt('s2')/sqrt(k3)*'mu' MTB > plot 's2' 'mu' - * 1.05+ * - * s2 - * - * * * * - * 0.70+ ** * - 2 *** *** * - * * * * * 2* 2 *** * * * - * 2 * * 2***2 ***22* *2** * * * - 2 222 2*22**2*53322*** 2 * 0.35+ * 2***2* 244*54*5333447627336*55* 3* * *2 - ** **232*3*3567++868776++4656542 23 * - * * * *447678+++++++++++++8+66743*2* 2* - *42429+9++++++++9+++8 3622 * - * *25**324*543 2 0.00+ ----+---------+---------+---------+---------+---------+--mu 32.00 32.20 32.40 32.60 32.80 33.00 MTB > dotplot 'mu' 's2' Each dot represents 4 points .: . ::: . ::. ::::.: .:::.::::::. ::::::::::::: :::::::::::::::: .. :::::::::::::::: :: :::::::::::::::::::::::.. . .........::::::::::::::::::::::::::::::.... .. ---+---------+---------+---------+---------+---------+---mu 32.00 32.20 32.40 32.60 32.80 33.00 Each dot represents 7 points .. ::: .::::: :::::::. ::::::::. .:::::::::::. ::::::::::::::. .:::::::::::::::::::::.:...... ... ..... .. . . +---------+---------+---------+---------+---------+-------s2 0.00 0.20 0.40 0.60 0.80 1.00