package randomizeBE [R for BE/BA]

posted by Helmut Homepage – Vienna, Austria, 2019-08-26 10:35  – Posting: # 20521
Views: 1,016

Hi sury,

» […] when we generate the randomization schedule in SAS using PROC PLAN Procedure, we generally take for example if subjects are 12, blocks will be 6…

As specified in your code.
» factors blocks=6 random sequence=2 random/ noprint;

» Comparison: when we run in SAS, within the block same sequence is not repeating, but when we run in R, we are getting the sequence repeated within the blocks. even though we are getting the P-value of 1.000 in R

RTFM

If blocksize is missing it defaults to 2*number of sequences.


I guess you didn’t specify that you want a blocksize of 6. Hence, with two sequences you will get a blocksize of four:

library(randomizeBE)
print(RL4(nsubj = 12), sumry = TRUE) # blocksize missing

Randomization table          created: 2019-08-26 10:28:28
(seed: 6284375
blocksize: 4 )

 subject seqno sequence
       1     1       TR
       2     2       RT
       3     2       RT
       4     1       TR
       5     1       TR
       6     1       TR
       7     2       RT
       8     2       RT
       9     2       RT
      10     2       RT
      11     1       TR
      12     1       TR


Summary of randomisation

12 subjects randomized into 2 sequence groups.
Number of subjects in sequence groups:
RT TR
 6  6
Runs test of randomness: p.value=0.2259


print(RL4(nsubj = 12, blocksize = 6, seed = 6284375), sumry = TRUE)

Randomization table          created: 2019-08-26 10:28:29
(seed: 6284375
blocksize: 6 )

 subject seqno sequence
       1     2       RT
       2     2       RT
       3     2       RT
       4     1       TR
       5     1       TR
       6     1       TR
       7     2       RT
       8     2       RT
       9     2       RT
      10     1       TR
      11     1       TR
      12     1       TR


Summary of randomisation

12 subjects randomized into 2 sequence groups.
Number of subjects in sequence groups:
RT TR
 6  6
Runs test of randomness: p.value=0.0693



» but is it fine to get that.

So what’s the problem?

Cheers,
Helmut Schütz
[image]

The quality of responses received is directly proportional to the quality of the question asked. ☼
Science Quotes

Complete thread:

Activity
 Mix view
Bioequivalence and Bioavailability Forum |  Admin contact
19,992 posts in 4,230 threads, 1,375 registered users;
online 9 (1 registered, 8 guests [including 7 identified bots]).
Forum time (Europe/Vienna): 22:38 CET

I’m not a pessimist,
I’m just a well informed optimist.    José Saramago

The BIOEQUIVALENCE / BIOAVAILABILITY FORUM is hosted by
BEBAC Ing. Helmut Schütz
HTML5