randomization [🇷 for BE/BA]

posted by ElMaestro  – Denmark, 2018-02-24 18:32 (2223 d 22:57 ago) – Posting: # 18457
Views: 13,760

Hello Ana Cristina,

I am thinking you can play around with this:


Ana.Cristinas.Randomizer2=function(Seed, N, WriteFile=F)
{
  if (N %% 4 != 0) {cat("N not divisible by four.!\n");return("Error.")}
  Gender=rep(c("F", "M"), N/2)
  Subject=c(1:N)
  Seqx=c(rep("TR", N/2), rep("RT", N/2))
  SeqF=sample(Seqx, replace=F)
  SeqM=sample(Seqx, replace=F)
  Seq=NULL
  for (i in 1:(N/2)) Seq=c(Seq, SeqF[i], SeqM[i])
  TrtP1=substr(Seq,1,1)
  TrtP2=substr(Seq,2,2)
  Rslt=data.frame(Subject, Gender, Seq, TrtP1, TrtP2)
  L1=paste("\n\nDate: ", format(Sys.Date(), format="%B %d %Y"),
     " Seed=", Seed, "\n\n", sep="")
  L2="\n\n\n\n\n______________________________\n"
  L3="Signature, date\n"
  if (WriteFile)
  {
    FN=paste(format(Sys.Date(), format="%B %d %Y"), " ", Seed,".txt", sep="")
    sink(file=FN)
    cat(L1)
    print(Rslt, row.names=F)
    cat(L2)
    cat(L3)   
    sink()
  }
  cat(L1)
  print(Rslt)
  return("Success.")
}

##can be run with e.g.
Ana.Cristinas.Randomizer2(12341234, 16, TRUE)


When the last argument is TRUE the function will print the randomisation code into a text file which is dated, contains the seed and which you can print and date+sign by hand. The filename is unique to the combination of date and seed.

Pass or fail!
ElMaestro

Complete thread:

UA Flag
Activity
 Admin contact
22,957 posts in 4,819 threads, 1,638 registered users;
103 visitors (0 registered, 103 guests [including 6 identified bots]).
Forum time: 17:29 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