randomization [R for BE/BA]

posted by ElMaestro  – Denmark, 2018-02-24 17:32  – Posting: # 18457
Views: 5,224

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.

if (3) 4

x=c("Foo", "Bar")
b=data.frame(x)
typeof(b[,1]) ##aha, integer?
b[,1]+1 ##then let me add 1


Best regards,
ElMaestro

“(...) targeted cancer therapies will benefit fewer than 2 percent of the cancer patients they’re aimed at. That reality is often lost on consumers, who are being fed a steady diet of winning anecdotes about miracle cures.” New York Times (ed.), June 9, 2018.

Complete thread:

Activity
 Mix view
Bioequivalence and Bioavailability Forum |  Admin contact
19,478 posts in 4,133 threads, 1,333 registered users;
online 4 (0 registered, 4 guests [including 3 identified bots]).
Forum time (Europe/Vienna): 05:46 CEST

If you don’t know anything about computers,
just remember that they are machines that do exactly what you tell them
but often surprise you in the result.    Richard Dawkins

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