d_labes
★★★

Berlin, Germany,
2012-07-10 15:02
(4279 d 10:42 ago)

Posting: # 8914
Views: 9,454
 

 randomizeBE - randomness control [🇷 for BE/BA]

Dear All!

In the R-package randomizeBE you can control the randomness of the list by a Wald-Wolfowitz runs test as Helmut has suggested in this post :cool:.
This suggestion was aimed to avoid some 'ordered' lists beside the fact that such lists are of course obtainable by dicing, also with low probability.

Internally randomizeBE uses the sequence numbers for the runs test. The runs test is originally formulated for dichotomous input only.
It is generalized to other cases using the dichotomizing of the values via <=median and > median. See f.i. german wikipedia. The english version suggests to omit values = median. IMHO this is not appropriate here where we are dealing with sequences.

Here some observations that the randomness control does not work that good in case of more than 2 sequences.

Example: 3 sequences numbered 1, 2, 3 (substitute a character representation whatever fits your needs. F.i. latin square "ABC", "BCA" and "CAB").
The 'alternating' random lists analogous to the 2x2 case give the following results:
    list                  runs p.value
 1 2 3 1 2 3 1 2 3 1 2 3   0.2502
 1 3 2 1 3 2 1 3 2 1 3 2   0.2502
 2 1 3 2 1 3 2 1 3 2 1 3   0.0438
 2 3 1 2 3 1 2 3 1 2 3 1   0.0438
 3 1 2 3 1 2 3 1 2 3 1 2   0.2502
 3 2 1 3 2 1 3 2 1 3 2 1   0.2502

Majority of them considered as compatible with null hypothesis: sequence is random.

The ordered list on the other hand is correctly considered as 'non-random':
    list                  runs p.value
 1 1 1 1 2 2 2 2 3 3 3 3   0.0028


BTW: You needn't dicing with function RL4() until you get one of the above lists to see the runs p-value. Use the function runs.pvalue() and a numeric representation of the sequences.
Example:
rl <- rep(c(1,2,3),4)
runs.pvalue(rl)


BTW2: The runs test is implemented via normal approximation like the implementations in the R-packages tseries or lawstat. The drawback is that it by design will give only approximate p-values for such small lists as discussed here. If these approximate values are conservative or liberal I don't know.

Regards,

Detlew
Helmut
★★★
avatar
Homepage
Vienna, Austria,
2012-07-10 16:20
(4279 d 09:25 ago)

@ d_labes
Posting: # 8915
Views: 7,289
 

 randomizeBE - randomness control

Dear Detlew,

interesting; I will chew on it time allowing. ;-)

❝ BTW2: The runs test is implemented via normal approximation like the implementations in the R-packages tseries or lawstat. The drawback is that it by design will give only approximate p-values for such small lists as discussed here. If these approximate values are conservative or liberal I don't know.


In my aged program I used exact values for N ≤200 according to

Swed FS and C Eisenhart
Tables for testing randomness of grouping in a sequence of alternatives
Ann Math Statist 14, 66–87 (1943)

These tables are reproduced in the “Wissenschaftliche Tabellen Geigy 3” p165, 8th Ed 1983. If you you don’t have them let me know and I will send you a scan.

Part of the tables at [image] books.

Dif-tor heh smusma 🖖🏼 Довге життя Україна! [image]
Helmut Schütz
[image]

The quality of responses received is directly proportional to the quality of the question asked. 🚮
Science Quotes
d_labes
★★★

Berlin, Germany,
2012-07-10 18:58
(4279 d 06:46 ago)

@ Helmut
Posting: # 8917
Views: 7,379
 

 runs test - critical values

Dear Helmut,

Tables for testing randomness of grouping in a sequence of alternatives

❝ Ann Math Statist 14, 66–87 (1943) These tables are reproduced in the “Wissenschaftliche Tabellen Geigy 3” p165, 8th Ed 1983. If you you don’t have them let me know and I will send you a scan.


Thanks for your kind offer.

But working with critical values is against the spirit of the implementation.
Here the user can choose an alpha value for the randomness control.

Regards,

Detlew
d_labes
★★★

Berlin, Germany,
2012-07-12 12:34
(4277 d 13:11 ago)

@ Helmut
Posting: # 8938
Views: 7,387
 

 randomizeBE - randomness control

Dear Helmut, dear All,

meanwhile I have tried to implement the exact runs distribution.
If I had not made a a big mistake the results are as following for the above presented 'alternating' lists:
(Also a little bug corrected if value==median):crying:
                         asymptotic asymptotic    exact
    list                   normal    normal cc.  runs distri.
 1 2 3 1 2 3 1 2 3 1 2 3   0.2502     0.4209      0.4242
 1 3 2 1 3 2 1 3 2 1 3 2   0.2502     0.4209      0.4242
 2 1 3 2 1 3 2 1 3 2 1 3   0.0658     0.1350      0.1414
 2 3 1 2 3 1 2 3 1 2 3 1   0.2502     0.4209      0.4242
 3 1 2 3 1 2 3 1 2 3 1 2   0.0658     0.1350      0.1414
 3 2 1 3 2 1 3 2 1 3 2 1   0.2502     0.4209      0.4242

cc. denotes a continuity correction to the asymptotic normal distributed statistic like SPSS uses.

For these lists it had to be concluded that the use of the exact distribution worsens the situation. It is more unlikely that the Null "list is a random sequence" is rejected. The use of cc. approximates the exact distribution well but also leads to a conservative test.

The fact that the randomness control doesn't work in all cases seems known.
Just to cite from:
Ying Wang
ECE 461 PROJECT REPORT, MAY 2003
"Nonparametric Tests for Randomness"
online resource

"... From the above analysis, most of the tests are vulnerable to a certain set of sequences, which are deterministic but accepted as random processes."

Regards,

Detlew
Helmut
★★★
avatar
Homepage
Vienna, Austria,
2021-07-15 18:07
(987 d 07:37 ago)

@ d_labes
Posting: # 22469
Views: 2,102
 

 randomness control

Dear Detlew,

excavating an old thread… Recently I saw a deficiency letter of a European assessor.

The Applicant is asked to provide block length used in randomization as clear pattern cannot be deduced from randomization scheme.


Single site, one group. The randomization was:

Subj Seq
  1   TR
  2   TR
  3   RT
  4   TR
  5   RT
  6   TR
  7   RT
  8   RT
  9   RT
 10   TR
 11   RT
 12   TR
 13   TR
 14   RT
 15   RT
 16   TR
 17   TR
 18   TR
 19   RT
 20   RT
 21   TR
 22   TR
 23   RT
 24   RT

What was the question?

Dif-tor heh smusma 🖖🏼 Довге життя Україна! [image]
Helmut Schütz
[image]

The quality of responses received is directly proportional to the quality of the question asked. 🚮
Science Quotes
Ohlbe
★★★

France,
2021-07-15 19:59
(987 d 05:45 ago)

@ Helmut
Posting: # 22470
Views: 2,212
 

 randomness control

Dear Helmut,

❝ The Applicant is asked to provide block length used in randomization as clear pattern cannot be deduced from randomization scheme


Well, can't the assessor count up to 8 ?

SCNR :-D

Regards
Ohlbe
Helmut
★★★
avatar
Homepage
Vienna, Austria,
2021-07-15 20:32
(987 d 05:12 ago)

@ Ohlbe
Posting: # 22471
Views: 2,100
 

 randomness control

Dear Ohlbe,

❝ ❝ The Applicant is asked to provide block length used in randomization as clear pattern cannot be deduced from randomization scheme


❝ Well, can't the assessor count up to 8 ?


Given from the remainder of the deficiency letter, I wouldn’t place a bet on it. :-D

Dif-tor heh smusma 🖖🏼 Довге життя Україна! [image]
Helmut Schütz
[image]

The quality of responses received is directly proportional to the quality of the question asked. 🚮
Science Quotes
d_labes
★★★

Berlin, Germany,
2021-07-18 13:15
(984 d 12:29 ago)

@ Helmut
Posting: # 22477
Views: 2,098
 

 randomness control

Dear Helmut,

as numerous times noticed: Assessors are a strange bunch of subjects :lol3:

Regards,

Detlew
UA Flag
Activity
 Admin contact
22,957 posts in 4,819 threads, 1,636 registered users;
70 visitors (0 registered, 70 guests [including 8 identified bots]).
Forum time: 00:45 CET (Europe/Vienna)

Nothing shows a lack of mathematical education more
than an overly precise calculation.    Carl Friedrich Gauß

The Bioequivalence and Bioavailability Forum is hosted by
BEBAC Ing. Helmut Schütz
HTML5