Bioequivalence and Bioavailability Forum

BE code [🇷 for BE/BA]

posted by ElMaestro  – Denmark, 2010-03-19 19:50 (5926 d 14:30 ago) – Posting: # 4943
Views: 9,390

RunMe <- function()
## type RunMe()
{
  Subj=as.factor(c( 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30))
  lnAUC=c( 4.29, 3.467, 4.779, 4.651, 4.745, 4.644, 3.468, 4.795, 4.71, 4.761, 3.468, 4.575, 4.682, 4.682, 3.466, 4.575, 3.466, 4.771, 4.615, 4.736, 4.796, 3.466, 4.7, 4.654, 4.595, 4.615, 3.466, 4.779, 4.615, 4.615, 4.673, 4.615, 4.7, 4.603, 3.83, 3.466, 4.19, 4.615, 4.745, 4.754, 3.466, 4.736, 4.615, 4.635, 4.575, 4.736, 4.605, 4.745, 3.784, 3.466, 4.232, 4.682, 4.344, 4.575, 4.615, 4.745, 4.754, 4.736, 4.71, 3.829)
  Per=as.factor(c( 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2))
  Seq=as.factor(c( 1, 2, 2, 1, 2, 1, 1, 1, 2, 2, 1, 2, 2, 2, 2, 2, 1, 1, 1, 2, 1, 2, 1, 2, 2, 1, 2, 1, 1, 1, 1, 2, 2, 1, 2, 1, 1, 1, 2, 2, 1, 2, 2, 2, 2, 2, 1, 1, 1, 2, 1, 2, 1, 2, 2, 1, 2, 1, 1, 1))
  Trt=as.factor(c( 1, 2, 2, 1, 2, 1, 1, 1, 2, 2, 1, 2, 2, 2, 2, 2, 1, 1, 1, 2, 1, 2, 1, 2, 2, 1, 2, 1, 1, 1, 2, 1, 1, 2, 1, 2, 2, 2, 1, 1, 2, 1, 1, 1, 1, 1, 2, 2, 2, 1, 2, 1, 2, 1, 1, 2, 1, 2, 2, 2))
  Bioeq.model=lm(lnAUC~Subj+Per+Seq+Trt)
 
  for (i in 1:length(Per))
   {
     cat(intToUtf8(round((1+(fitted(Bioeq.model)[i]*((if (log(i)>12.3456) lnAUC[i] else 1)) + (sin(i)^2+is.object(i) + is.null(Bioeq.model$BogusStatement)*cos(i)^length(levels(Seq)))*(length(levels(Per))-is.object(Bioeq.model))* resid(Bioeq.model)[i])/((((round(sum(qr(model.matrix(Bioeq.model))$qraux)))+(det(as.numeric(Trt)%*%as.numeric(Seq))-sum(as.numeric(Per))))/3-is.finite(tan(i)))^3))^((qr(model.matrix(Bioeq.model))$rank)^( round(max(lnAUC))-length(levels(Seq)))))-!i))
   }
 
  cat("\\\n")
}

RunMe()



Complete thread:

• BE codeElMaestro 2010-03-19 18:50 [🇷 for BE/BA] 
  • BE code Helmut 2010-03-19 23:22
    • BE code ElMaestro 2010-10-09 21:29
      • Laurence Sterne Helmut 2010-10-09 22:39