package randomizeBE [R for BE/BA]

posted by Helmut Homepage – Vienna, Austria, 2019-08-26 10:35 (390 d 02:58 ago) – Posting: # 20521
Views: 2,920

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?

Dif-tor heh smusma 🖖
Helmut Schütz
[image]

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

Complete thread:

Activity
 Admin contact
21,068 posts in 4,392 threads, 1,465 registered users;
online 15 (0 registered, 15 guests [including 7 identified bots]).
Forum time: Saturday 13:33 CEST (Europe/Vienna)

The purpose of models is not to fit the data,
but to sharpen the questions.    Samuel Karlin

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