Yurra-Bob: a yet not paid debt

All this News Readers are (eternally) waiting for the justification how they posted on May 19, 2006 3:51 PM about my program concerning Mann-Whitney critical values by Monte Carlo, namely

***Another incorrect program: ( ***

If not based, dishonest and defamatory statements are allowed to be posted here then any serious and performing discussion is completely impossible. The cited Readers should be compelled to prove what they stated.

______licas (Luis A. Afonso)


REM "MANN"
CLS
INPUT " size shorter "; na
INPUT " size larger "; nb: n = na + nb
nna = (1 + n) * n / 2
DIM x(nb), y(nb), join(na + nb), has(nna)
vv = 10000
FOR rpt = 1 TO vv
LOCATE 6, 50
PRINT USING "###########"; vv - rpt
RANDOMIZE TIMER
FOR i = 1 TO na
x(i) = RND: join(i) = x(i)
NEXT i
FOR ii = 1 TO nb
y(ii) = RND: join(na + ii) = y(ii)
NEXT ii
FOR jo = 1 TO na
w = join(jo): order = 1
FOR j1 = 1 TO na + nb
IF join(j1) < w THEN order = order + 1
NEXT j1
xorder = xorder + order
NEXT jo
has(xorder) = has(xorder) + 1
xorder = 0
aa = INT(rpt / 10000): a = rpt / 10000
IF aa <> a THEN GOTO 100
u(1) = .999 * rpt: u(2) = .995 * rpt
u(3) = .99 * rpt: u(4) = .975 * rpt
u(5) = .95 * rpt: u(6) = .9 * rpt
FOR e = 1 TO 6: sum = 0
FOR ij = 0 TO nna
sum = sum + has(ij)
IF sum > u(e) THEN GOTO 111
NEXT ij
LOCATE 7, 12
111 wy(e) = sum / rpt: iab = na * (n + 1) - ij
PRINT USING "### #.#### ###";
ij; wy(e); iab
NEXT e
100 NEXT rpt : END

This program provide exactly the same critical values that can be found on statistical Tables.
.

Relevant Pages