Simul Ants questions [Regulatives / Guidelines]

posted by d_labes  – Berlin, Germany, 2011-11-04 16:46 (4930 d 06:52 ago) – Posting: # 7606
Views: 8,939

Dear EM,

load(nlme)

 nConservative=0
 nAntiConservative=0
 
 ...
 for (iter in 1:100)
 {
 
 ## T is 1, R1 is 2 and R2 is 3
 lny1 = rnorm(n=36, mean=95, sd=0.15*95)
 for (i in 1:36) if (trt1[i]==2) lny1[i]=rnorm(n=1, mean = 90, sd = 0.10*90)
 for (i in 1:36) if (trt1[i]==3) lny1[i]=rnorm(n=1, mean = 110, sd = 0.15*110)
 ...


Really?
I would suggest:
If you consider 95, 90 and 110 as the arithmetic means on the untransformed scale use log(x)-sd2/2 as the mean of the logs.
If your numbers 0.15, 0.10 and 0.15 are the CV's user sd = sqrt(log(CV*CV+1)) as sd of the logs.

The manner you have written the simulation for the three treatments suggest you consider the log(metric) for the 3 treatments as independent.
Is this really the case for a cross-over where the 3 treatment were applied to each subject :ponder:? IMHO we have to consider some correlation between them.

Regards,

Detlew

Complete thread:

UA Flag
Activity
 Admin contact
23,424 posts in 4,927 threads, 1,667 registered users;
93 visitors (0 registered, 93 guests [including 9 identified bots]).
Forum time: 00:38 CEST (Europe/Vienna)

Patients may recover in spite of drugs or because of them.    John Gaddum

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