GMACRO
binpower
note
note Gütefunktion bei Binomialtest (Aufgabe 52):
note   Eingabe:  K1 = Stichprobenumfang n
note             K2 = linke krit. alpha-Schranke k_l zu Bi(n,theta_0)
note             K3 = rechte krit. alpha-Schranke k_r zu Bi(n,theta_0)
note   Ausgabe:  C1: theta = 0.01 (0.01) 0.99
note             C2: beta_l(theta) = Pr{Bi(n,theta) < k_l}
note             C3: beta_r(theta) = Pr{Bi(n,theta) > k_r} 
note
note
name k1 'n'
name k2 'k_l'
name k3 'k_r'
print k1-k3
note Sind die Eingabewerte in Ordnung (Y/N)?
yesno k4
if k4=0
  exit
endif
name c1 'theta_'
name c2 'beta_l(theta)'
name c3 'beta_r(theta)'
name k5 'theta'
do k4 = 1 : 99
   let 'theta'=k4*0.01
   # Berechnung von beta_l(theta)
   let k6='k_l'-1
   CDF k6 k7;
   Binomial 'n' k5. 
   let c1(k4)='theta'
   let c2(k4)=k7
   # Berechnung von beta_r(theta)
   CDF 'k_r' k8;
   Binomial 'n' k5.    
   let c3(k4)=1-k8
enddo
ENDMACRO