!set n=$teller
bewerking=bewerking1.proc
wims_read_parm=$empty

!if $graad=0
    aantal_vragen=5
!else
    aantal_vragen=2
!endif     
counter=0    
data$n=<table border="1"><td>meting</td><td>frequentie</td><tr>
gemiddelde$n=$empty
som=0
deler=0
!if $soort=1
    nivo_title=
    counter=$maxcounter
    !for p=1 to $maxcounter
	meting=!randint 1,$afrondingsfactor
	meting$p=$[1000*(1+ $meting/$afrondingsfactor)]
	try$p=!randitem 1,2,3,4,5
	!for s=1 to $(try$p)
	    wims_read_parm=!append item $(meting$p) to $wims_read_parm
	!next s
	data$n=!append line <td>$(meting$p)&nbsp;seconden</td><td>$(try$p)</td><tr> to $(data$n)
	som=$[$som + $(meting$p)*$(try$p)]
	deler=$[$deler + $(try$p)]
    !next p
    klaar=1
!else
    nivo_title=
    !if $klaar=1
	!for p=1 to $maxcounter
	    som=$[$som + $(meting$p)*$(try$p)]
	    deler=$[$deler + $(try$p)]
	    !for s=1 to $(try$p)
		wims_read_parm=!append item $(meting$p) to $wims_read_parm
	    !next s
	    data$n=!append line <td>$(meting$p)&nbsp;seconden</td><td>$(try$p)</td><tr> to $(data$n)
	!next p
    !else
	 !exit
    !endif
!endif

data$n=$(data$n)</table>
# gemiddelde
GOED1$n=$[$som/$deler]

# 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 $deler
    w=!item $p of $wims_read_parm
    W=$[($w-$(GOED1$n))^2]
    NW=$[$W + $NW]
!next p
    
GOED2$n=$[sqrt($NW/$deler)]


# eerste & derde kwartiel
    
wims_read_parm=!sort numeric list of $wims_read_parm
!if $[gcd($deler,2)]=2
    even=1
    # mediaan
    # !readproc slib/stat/median $wims_read_parm 
    # GOED3$n=$slib_out
    m1=!item $[0.5*$deler] of $wims_read_parm
    m2=!item $[0.5*$deler+1] of $wims_read_parm
    GOED3$n=$[($m2+$m1)/2]
    
    !if $[gcd((0.5*$deler),2)]=2
        w1=$[0.25*$deler]
	w11=$[0.25*$deler+1]
	w3=$[0.75*$deler]
	w33=$[0.75*$deler+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*$deler+0.5]
	w3=$[0.75*$deler+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
    GOED3$n=!item $[($deler+1)/2] of $wims_read_parm
    $[gcd((0.5*($deler-1)),2)]
    !if $[gcd((0.5*($deler-1)),2)]=2
	w1=$[0.25*$deler - 0.25]
	w11=$[0.25*$deler + 0.75]
	w3=$[0.75*$deler + 0.25]
	w33=$[0.75*$deler + 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*$deler+0.25]
	w3=$[0.75*$deler+0.75]
	w33=$w3
	q1=!item $w1 of $wims_read_parm
	q11=$q1
	q3=!item $w3 of $wims_read_parm	
	q33=$q3
    !endif
!endif

    
GOED4$n=$[($q1+$q11)/2]
GOED5$n=$[($q3+$q33)/2]
goed1$n=Het gemiddelde is dus $(GOED1$n)
goed2$n=De standaardafwijking $m_sigma is $(GOED2$n)
goed3$n=De mediaan $m_mu is $(GOED3$n)
goed4$n=Het eerste kwartiel is $(GOED4$n)
goed5$n=Het derde kwartiel is $(GOED5$n)
