randomization [R for BE/BA]

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

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,408 posts in 4,123 threads, 1,325 registered users;
online 11 (0 registered, 11 guests [including 7 identified bots]).
Forum time (Europe/Vienna): 11:55 CEST

We must be careful not to confuse data with the abstractions
we use to analyze them.    William James

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