Bioequivalence and Bioavailability Forum • Dataset

Dataset [🇷 for BE/BA]

posted by Helmut – Vienna, Austria, 2010-04-26 00:38 (5113 d 21:47 ago) – Posting: # 5216
Views: 55,881

Dear bears,

❝ When I revised bear today, […] (with built-in dataset in bear v2.4.2) for CIs.

Are you trying to confuse me? The built-in dataset consists of 20 (10:10) observations. For comparisons can we please stick to this dataset, lnAUC0t only?

   subj drug Cmax AUC0t AUC0INF lnCmax lnAUC0t lnAUC0INF

1     1    1 88.6  1510    1530 4.4841  7.3199    7.3330

2     2    1 52.5   883     890 3.9608  6.7833    6.7912

3     3    1 92.0  1650    1670 4.5218  7.4085    7.2226

4     4    1 56.0  1015    1050 4.0254  6.9226    6.9565

5     5    1 84.0  1556    1570 4.4308  7.3499    7.3588

6     6    1 84.8  1412    1432 4.4403  7.2528    7.2668

7     7    1 83.0  1353    1356 4.4188  7.2101    7.2123

8     8    1 96.4  1443    1450 4.5685  7.2745    7.2793

9     9    1 68.1  1299    1305 4.2210  7.1694    7.1740

10   10    1 33.5   560     570 3.5115  6.3279    6.3456

11   11    2 70.3  1284    1290 4.2528  7.1577    7.1624

12   12    2 73.5  1391    1340 4.2973  7.2378    7.2004

13   13    2 50.2   873     890 3.9160  6.7719    6.7920

14   14    2 62.2  1211    1230 4.1304  7.0992    7.1148

15   15    2 74.1  1233    1255 4.3054  7.1172    7.1349

16   16    2 60.4  1172    1182 4.1010  7.0665    7.0750

17   17    2 60.4  1172    1185 4.1010  7.0665    7.0775

18   18    2 75.3  1336    1355 4.3215  7.1974    7.2116

19   19    2 76.8  1348    1355 4.3412  7.2064    7.2116

20   20    2 82.9  1419    1425 4.4176  7.2577    7.2619

For non-bears, here is the dataset (balanced, 20 subjects) for download.

With this code

# Change to folder containg data files in RGui!

TotalData <- read.csv("PAR-BEAR.txt", header=T, sep="\t", dec=".")

attach(TotalData)

subj <- as.factor(subj)

drug <- as.factor(drug)

boxplot(lnAUC0t ~ drug, col="lightgray", yl)

points(jitter(as.numeric(drug), factor=0.25), lnAUC0t, col="blue")

result <- t.test(lnAUC0t ~ drug, var.equal = FALSE, conf.level = 0.90)

result

cat("T/R [%] with alpha 0.05 (90% CI)", sep=" ", fill=TRUE)

tbldiff <- matrix(

  c(as.numeric(exp(diff(result$estimate))),

  sort(as.numeric(exp(-result$conf.int)))),

  byrow = TRUE, nrow = 1)

dimnames(tbldiff) <- list("Ratio",

                     c("  Point Estimate",

                       "  CL90 lower",

                       "  CL90 upper" ))

round(tbldiff*100,3)

I get

        Welch Two Sample t-test



data:  lnAUC0t by drug 

t = -0.1392, df = 12.043, p-value = 0.8916

alternative hypothesis: true difference in means is not equal to 0 

90 percent confidence interval:

 -0.2199223  0.1880423 

sample estimates:

mean in group 1 mean in group 2 

        7.10189         7.11783



T/R [%] with alpha 0.05 (90% CI)

        Point Estimate   CL90 lower   CL90 upper

Ratio          101.607       82.858      124.598

As expected the CI is slightly wider as the ones obtained by simple lm or the strange lme (83.312-123.920). Have a look at the plot. ;-)

Edit: Remove the ~~red part~~ of the code. The line should read
boxplot(lnAUC0t ~ drug, col="lightgray")
THX to D. Labes for beta-testing. [Helmut]

—
Dif-tor heh smusma 🖖🏼 Довге життя Україна!
Helmut Schütz

The quality of responses received is directly proportional to the quality of the question asked. 🚮
Science Quotes

Complete thread:

Parallel bears meeting at random in infinity d_labes 2010-04-22 11:43 [🇷 for BE/BA]
- Parallel bears meeting at random in infinity ElMaestro 2010-04-22 12:53
  - Parallel groups in bear - CIs d_labes 2010-04-22 14:00
    - Parallel groups in bear - CIs ElMaestro 2010-04-22 21:47
      - Parallel groups in bear - CIs d_labes 2010-04-23 09:09
    - Parallel groups in bear - CIs yjlee168 2010-04-25 23:29
- Parallel bears meeting at random in infinity yjlee168 2010-04-22 23:09
  - Modelling Parallel bears d_labes 2010-04-23 09:12
    - Modelling Parallel bears yjlee168 2010-04-23 21:14
      - Validating vs. WinNonlin... Helmut 2010-04-24 00:28
        
        Validating vs. WinNonlin... yjlee168 2010-04-24 19:36
        Validating vs. WinNonlin... yjlee168 2010-04-26 00:09
        
        Validating vs. WinNonlin... Helmut 2010-04-26 01:29
        
        WNL in replicate BE yjlee168 2010-04-26 08:59
        
        WNL in replicate BE Helmut 2010-04-26 16:15
      - Modelling Parallel bears yjlee168 2010-04-25 19:34
        
        Modelling Parallel bears ElMaestro 2010-04-25 20:40
        DatasetHelmut 2010-04-25 22:38
        
        Dataset yjlee168 2010-04-25 22:44
        
        Dataset Helmut 2010-04-26 01:13
        
        Dataset yjlee168 2010-04-26 08:16
        
        NCA → Statistical analysis for parallel study Helmut 2010-04-26 13:12
        
        NCA → Statistical analysis for parallel study yjlee168 2010-04-26 18:43
        
        dilemma yjlee168 2010-04-26 08:41
        Equal variances d_labes 2010-04-26 09:04
        
        Equal variances yjlee168 2010-04-26 09:22
        
        GLM = Equal variances d_labes 2010-04-26 13:29
        
        GLM = Equal variances Helmut 2010-04-26 14:45
        
        I'm a believer d_labes 2010-04-26 15:58
        
        I'm a believer Helmut 2010-04-26 16:31
        
        Equal variances Helmut 2010-04-26 12:55
        
        gls() for unequal variances? d_labes 2010-04-26 16:36
        
        gls() for unequal variances? Helmut 2010-04-26 17:00
        Sims Helmut 2010-04-27 01:36
        
        Sandwich - Simsalabim d_labes 2010-04-28 10:58
        
        Sandwich - Simsalabim Helmut 2010-04-28 14:19
        parametrization of R function rlnorm martin 2010-05-02 18:22
        
        Mean of log-normal d_labes 2010-05-03 16:22
        parametrization of R function rlnorm ElMaestro 2013-07-26 21:42
        
        Martin‽ Helmut 2013-07-28 02:01