!set n=$teller
!if $printbaar >0
    white=magenta
!else
    white=white    
!endif    
!if $graad=0
    aantal_vragen=6
!else
    aantal_vragen=3
!endif 
bewerking=bewerking1.proc
wims_read_parm=$empty
data$n=<table border="1" cellpadding="3"><td>cijfer</td><td>frequentie</td><tr>
gemiddelde$n=$empty
som=0
deler=0
GRID=$empty
VTEXT=$empty
HTEXT=$empty
RECT=$empty
nivo_title=Het histogram 
# maxcounter=aantal leerlingen
maxcounter=!randint 15,40
D=1,2,2,3,3,4,4,4,4,5,5,5,5,6,6,6,6,6,6,6,6,6,6,6,7,7,7,7,7,7,7,7,7,7,7,7,8,8,8,8,9,9,10
!for s=1 to 10
    freq$s=0
!next s    
!for p=1 to $maxcounter
    cijfer$p=!randitem $D
    !for s=1 to 10
	!if $(cijfer$p)=$s
	    freq$s=$[$(freq$s)+1]
	!endif    
    !next s
    wims_read_parm=!append item $(cijfer$p) to $wims_read_parm
    som=$[$som + $(cijfer$p)]
!next p    

max=1
!for s=1 to 10
    !if $(freq$s) > $max
        max=$(freq$s)
    !endif
!next s 
!for p=1 to $max
    GRID=!append line line 0,$p,10,$p,black to $GRID
    VTEXT=!append line text black,-1,$[$p+0.2],huge,$p to $VTEXT
!next p
!for p=1 to 10
    RECT=!append line frect $[$p-1],0,$p,$(freq$p),red to $RECT   
    GRID=!append line line $p,0,$p,$max,black to $GRID
    HTEXT=!append line text black,$[$p-0.7],-0.3,huge,$p to $HTEXT
    data$n=!append line <td>$p</td><td>$(freq$p)</td><tr> to $(data$n)
!next p

plaatje$n=400,400\
transparent $white\
xrange -2,12\
yrange -2,$max+2\
$RECT \
$GRID\
$VTEXT \
$HTEXT \
textup black,-2,$[0.3*$max],huge,Frequentie\
text black,5,-1.2,huge,Cijfer\
arrow 5,-1.8,8,-1.8,20,black\
arrow -1.3,$[0.2*$max],-1.3,$[0.7*$max],20,black\
linewidth 4\
rect 0,0,10,$max,black \


data$n=$(data$n)</table>
# gemiddelde
GOED2$n=$[$som/$maxcounter]

# standaard wafwijking
# is s geen sigma
# !readproc slib/stat/deviation $wims_read_parm 
# GOED2$n=!item 2 of $slib_out
    
NW=0
!for p=1 to $maxcounter
    w=!item $p of $wims_read_parm
    W=$[($w-$(GOED2$n))^2]
    NW=$[$W + $NW]
!next p
    
GOED3$n=$[sqrt($NW/$maxcounter)]


# eerste & derde kwartiel
    
wims_read_parm=!sort numeric list of $wims_read_parm
!if $[gcd($maxcounter,2)]=2
    even=1
    # mediaan
    # !readproc slib/stat/median $wims_read_parm 
    # GOED3$n=$slib_out
    m1=!item $[0.5*$maxcounter] of $wims_read_parm
    m2=!item $[0.5*$maxcounter+1] of $wims_read_parm
    GOED4$n=$[($m2+$m1)/2]
    
    !if $[gcd((0.5*$maxcounter),2)]=2
        w1=$[0.25*$maxcounter]
	w11=$[0.25*$maxcounter+1]
	w3=$[0.75*$maxcounter]
	w33=$[0.75*$maxcounter+1]
	q1=!item $w1 of $wims_read_parm
	q11=!item $w11 of $wims_read_parm
	q3=!item $w3 of $wims_read_parm
	q33=!item $w33 of $wims_read_parm
    !else
	w1=$[0.25*$maxcounter+0.5]
	w3=$[0.75*$maxcounter+0.5]
	w33=$w3
	q1=!item $w1 of $wims_read_parm
	q11=$q1
    	q3=!item $w3 of $wims_read_parm
	q33=$q3    
    !endif
!else
    even=0
    # mediaan
    # !readproc slib/stat/median $wims_read_parm 
    # GOED3$n=$slib_out
    GOED4$n=!item $[($maxcounter+1)/2] of $wims_read_parm
    $[gcd((0.5*($maxcounter-1)),2)]
    !if $[gcd((0.5*($maxcounter-1)),2)]=2
	w1=$[0.25*$maxcounter - 0.25]
	w11=$[0.25*$maxcounter + 0.75]
	w3=$[0.75*$maxcounter + 0.25]
	w33=$[0.75*$maxcounter + 1.25]
	q1=!item $w1 of $wims_read_parm
	q11=!item $w11 of $wims_read_parm
	q3=!item $w3 of $wims_read_parm
	q33=!item $w33 of $wims_read_parm
    !else
	w1=$[0.25*$maxcounter+0.25]
	w3=$[0.75*$maxcounter+0.75]
	w33=$w3
	q1=!item $w1 of $wims_read_parm
	q11=$q1
	q3=!item $w3 of $wims_read_parm	
	q33=$q3
    !endif
!endif

    
GOED5$n=$[($q1+$q11)/2]
GOED6$n=$[($q3+$q33)/2]
GOED1$n=$maxcounter
goed1$n=Er zitten dus $maxcounter leerlingen in deze klas.
goed2$n=Het gemiddelde is dus $(GOED2$n)
goed3$n=De standaardafwijking $m_sigma is $(GOED3$n)
goed4$n=De mediaan $m_mu is $(GOED4$n)
goed5$n=Het eerste kwartiel is $(GOED5$n)
goed6$n=Het derde kwartiel is $(GOED6$n)
