Bioequivalence and Bioavailability Forum • example

example [🇷 for BE/BA]

posted by martin – Austria, 2008-10-16 22:02 (6445 d 04:10 ago) – Posting: # 2547
Views: 26,233

dear ace !

here is a rather strange (artificial) example generated in a simulation study (R version 2.6.2 for windows XP)

best regards

martin


time <- c(0.08263463, 0.52824754, 1.17358567, 1.39767757, 2.47663194,  3.19031555, 3.97925368, 6.95466074, 7.42732263, 12.08734425, 12.35692940, 13.63305248, 21.64587886, 29.33632297)

conc <- c(1.175213, 8.748377, 12.650698, 12.741200, 12.965722, 13.648985, 13.633996, 11.110552, 16.232875, 6.747145, 7.156791, 7.776978, 3.241740, 1.528261)



dat <- data.frame(time=time, conc=conc)



# without fishing for NA's

n_TTT_ARS=0

r.adj2=0

for (i in

(nrow(dat)-2):(min(seq_along(dat$time)[dat$time>=dat$time[which.max(dat$conc)]*2])))

 {

   check <- summary(lm(log(conc)~time,dat[i:nrow(dat),]))$adj.r.squared #modification

   if ((r.adj2 - check <(0.0001))) {

      n_TTT_ARS = nrow(dat)-i+1

      r.adj2 = summary(lm(log(conc)~time,dat[i:nrow(dat),]))$adj.r.squared

   }

}



# fishing for NA's

n_TTT_ARS=0

r.adj2=0

for (i in

(nrow(dat)-2):(min(seq_along(dat$time)[dat$time>=dat$time[which.max(dat$conc)]*2])))

 {

   check <- summary(lm(log(conc)~time,dat[i:nrow(dat),]))$adj.r.squared #modification

   if (!is.na(check) && (r.adj2 - check < (0.0001))) {

      n_TTT_ARS = nrow(dat)-i+1

      r.adj2 = summary(lm(log(conc)~time,dat[i:nrow(dat),]))$adj.r.squared

   }

}

Complete thread:

update: bear v1.1.4 for R (a free tool for ABE) yjlee168 2008-09-17 13:33 [🇷 for BE/BA]
- update: bear v1.1.4 for R (a free tool for ABE) Helmut 2008-09-17 16:59
  - update: bear v1.1.4 for R (a free tool for ABE) yjlee168 2008-09-17 21:03
- update: bear v1.1.4 for R (a free tool for ABE) martin 2008-09-19 13:13
  - update: bear v1.1.4 for R (a free tool for ABE) yjlee168 2008-09-22 18:19
  - update: bear v1.1.4 for R (a free tool for ABE) Helmut 2008-09-23 00:46
    - update: bear v1.1.4 for R (a free tool for ABE) martin 2008-09-27 09:37
      - Package lee for R Helmut 2008-09-27 14:05
  - update: bear v1.1.4 for R (a free tool for ABE) yjlee168 2008-09-23 10:07
    - update: bear v1.1.4 for R (a free tool for ABE) martin 2008-09-27 09:25
      - update: bear v1.1.4 for R (a free tool for ABE) yjlee168 2008-09-29 12:03
        
        update: bear v1.1.4 for R (a free tool for ABE) Helmut 2008-09-29 12:10
        
        update: bear v1.1.4 for R (a free tool for ABE) yjlee168 2008-09-29 13:07
        
        modified function for selection of time points martin 2008-10-11 15:30
        
        modified function for selection of time points yjlee168 2008-10-11 19:45
        
        last 3 time points martin 2008-10-12 00:04
        
        modified function for selection of time points yjlee168 2008-10-16 11:20
        
        modified function for selection of time points Aceto81 2008-10-16 12:09
        
        modified code martin 2008-10-16 12:29
        
        modified code Aceto81 2008-10-16 14:11
        
        abs(AIC) martin 2008-10-16 15:03
        examplemartin 2008-10-16 20:02
        
        example Aceto81 2008-10-17 12:30