// Test des signes en agronomie
// simulation des rendements des parcelles sous deux engrais A et A
// et réalisation du test des signes déterministe
n=10;
alpha=0.05;
// On suppose les rendements gaussiens de paramètres m et sigma
// choisis comme suit 
moy=10+3.*rand(1:n,'uniform');
sig=2+rand(1:n,'uniform');
Z=rand(2,n,'normal');
X=moy+sig.*Z(1,:);
Y=moy+sig.*Z(2,:);
// affichage des rendements
printf('%6s%13s%13s','Champ','Rendement A','Rendement B')
printf('%6d%13.2f%13.2f\n',(1:n)',X',Y')
//
// Calculs pour le test:
U=bool2s(X<Y);
S=sum(U);
salpha=invfdrbin(n,0.5,1-alpha); 
// Cette instruction suppose que l'on ait déjà chargé Binomiales.sci
// à adapter donc!
if S<=salpha then
   disp("Acceptation de l''hypothèse de non amélioration du rendement")
else
   disp("Rejet de l''hypothèse de non amélioration du rendement")
end