======================= T DE STUDENT ========================= 'st-test.bas 'Traduit par www.FreeNRG.info d'un programme PASCAL (ttest: Student's t-test) 'de chez: 'http://www.lecb.ncifcrf.gov/~toms/delila/delilaprograms.html DECLARE FUNCTION zip! (q AS SINGLE, i AS SINGLE, j AS SINGLE, b AS SINGLE) DECLARE FUNCTION buzz! (t AS SINGLE, n AS SINGLE) DECLARE FUNCTION ttestproba! (t AS SINGLE, ddl AS SINGLE) DECLARE SUB try (t AS SINGLE, ddl AS SINGLE) CONST ecran = 5, imprimante = 6 DIM sortie AS INTEGER DIM SHARED e AS SINGLE DIM SHARED pi AS SINGLE DIM SHARED pj2 AS SINGLE DIM SHARED pj4 AS SINGLE DIM SHARED pi2 AS SINGLE DIM SHARED dgr AS SINGLE e = EXP(1!) pi = 4! * ATN(1!) pj2 = pi / 2! pj4 = pi / 4! pi2 = pi * 2! dgr = 180! / pi CLS OPEN "scrn:" FOR RANDOM AS #ecran OPEN "lpt1:" FOR OUTPUT AS #imprimante sortie = ecran '''sortie = imprimante PRINT #sortie, "les valeurs devraient approcher 0.90:" try 3.078, 1 try 1.886, 2 try 1.638, 3 try 1.345, 14 try 1.3128, 28 try 1.3103, 30 PRINT #sortie, "" IF (sortie = ecran) THEN a$ = INPUT$(1) PRINT #sortie, "les valeurs devraient approcher 0.95:" try 6.314, 1 try 2.92, 2 try 2.353, 3 try 1.761, 14 try 1.701, 28 try 1.697, 30 PRINT #sortie, "" IF (sortie = ecran) THEN a$ = INPUT$(1) PRINT #sortie, "les valeurs devraient approcher 0.975:" try 12.706, 1 try 4.303, 2 try 3.182, 3 try 2.145, 14 try 2.048, 28 try 2.042, 30 PRINT #sortie, "" IF (sortie = ecran) THEN a$ = INPUT$(1) PRINT #sortie, "les valeurs devraient approcher 0.99:" try 31.821, 1 try 6.965, 2 try 4.541, 3 try 2.624, 14 try 2.467, 28 try 2.457, 30 PRINT #sortie, "" IF (sortie = ecran) THEN a$ = INPUT$(1) PRINT #sortie, "les valeurs devraient approcher 0.995:" try 63.7, 1 try 9.93, 2 try 5.84, 3 try 2.98, 14 try 2.76, 28 try 2.75, 30 PRINT #sortie, "" IF (sortie = ecran) THEN a$ = INPUT$(1) PRINT #sortie, "les valeurs devraient approcher 0.9995:" try 636.619, 1 try 31.598, 2 try 12.924, 3 try 4.14, 14 try 3.694, 28 try 3.646, 30 PRINT #sortie, "" END ' ' ' ' ' ' FUNCTION buzz! (t AS SINGLE, n AS SINGLE) 'var DIM rt AS SINGLE, fk AS SINGLE, ek AS SINGLE, dk AS SINGLE t = ABS(t) rt = t / SQR(n) fk = ATN(rt) IF (n = 1) THEN buzz! = 1 - fk / pj2 ELSE ek = SIN(fk) dk = COS(fk) ''IF ((round(n) MOD 2) = 1) THEN IF ((INT(n) MOD 2) = 1) THEN buzz! = 1 - (fk + ek * dk * zip!(dk * dk, 2, n - 3, -1)) / pj2 ELSE buzz! = 1 - ek * zip!(dk * dk, 1, n - 3, -1) END IF END IF END FUNCTION ' ' ' ' ' ' SUB try (t AS SINGLE, ddl AS SINGLE) DIM prob AS SINGLE SHARED sortie AS INTEGER prob = ttestproba!(t, ddl) PRINT #sortie, USING " t = ###.### "; t; PRINT #sortie, USING " ddl = ### "; ddl; PRINT #sortie, USING " P = #.#### (#.####) "; prob; 1! - prob END SUB ' ' ' ' ' ' FUNCTION ttestproba! (t AS SINGLE, ddl AS SINGLE) ttestproba! = 1! - (buzz(t, ddl) / 2!) END FUNCTION ' ' ' ' ' ' FUNCTION zip! (q AS SINGLE, i AS SINGLE, j AS SINGLE, b AS SINGLE) DIM k AS SINGLE, zz AS SINGLE, z AS SINGLE zz = 1 z = zz k = i WHILE (k <= j)'do zz = zz * q * k / (k - b) z = z + zz k = k + 2 WEND zip! = z END FUNCTION ======================= SORTIES DU PROGRAMME ========================= les valeurs devraient approcher 0.90: t = 3.078 ddl = 1 P = 0.9000 (0.1000) t = 1.886 ddl = 2 P = 0.9000 (0.1000) t = 1.638 ddl = 3 P = 0.9000 (0.1000) t = 1.345 ddl = 14 P = 0.9000 (0.1000) t = 1.313 ddl = 28 P = 0.9000 (0.1000) t = 1.310 ddl = 30 P = 0.9000 (0.1000) les valeurs devraient approcher 0.95: t = 6.314 ddl = 1 P = 0.9500 (0.0500) t = 2.920 ddl = 2 P = 0.9500 (0.0500) t = 2.353 ddl = 3 P = 0.9500 (0.0500) t = 1.761 ddl = 14 P = 0.9500 (0.0500) t = 1.701 ddl = 28 P = 0.9500 (0.0500) t = 1.697 ddl = 30 P = 0.9500 (0.0500) les valeurs devraient approcher 0.975: t = 12.706 ddl = 1 P = 0.9750 (0.0250) t = 4.303 ddl = 2 P = 0.9750 (0.0250) t = 3.182 ddl = 3 P = 0.9750 (0.0250) t = 2.145 ddl = 14 P = 0.9750 (0.0250) t = 2.048 ddl = 28 P = 0.9750 (0.0250) t = 2.042 ddl = 30 P = 0.9750 (0.0250) les valeurs devraient approcher 0.99: t = 31.821 ddl = 1 P = 0.9900 (0.0100) t = 6.965 ddl = 2 P = 0.9900 (0.0100) t = 4.541 ddl = 3 P = 0.9900 (0.0100) t = 2.624 ddl = 14 P = 0.9900 (0.0100) t = 2.467 ddl = 28 P = 0.9900 (0.0100) t = 2.457 ddl = 30 P = 0.9900 (0.0100) les valeurs devraient approcher 0.995: t = 63.700 ddl = 1 P = 0.9950 (0.0050) t = 9.930 ddl = 2 P = 0.9950 (0.0050) t = 5.840 ddl = 3 P = 0.9950 (0.0050) t = 2.980 ddl = 14 P = 0.9950 (0.0050) t = 2.760 ddl = 28 P = 0.9950 (0.0050) t = 2.750 ddl = 30 P = 0.9950 (0.0050) les valeurs devraient approcher 0.9995: t = 636.619 ddl = 1 P = 0.9995 (0.0005) t = 31.598 ddl = 2 P = 0.9995 (0.0005) t = 12.924 ddl = 3 P = 0.9995 (0.0005) t = 4.140 ddl = 14 P = 0.9995 (0.0005) t = 3.694 ddl = 28 P = 0.9995 (0.0005) t = 3.646 ddl = 30 P = 0.9995 (0.0005)
septix