Bear to bear interval with 90% confidence [🇷 for BE/BA]

posted by yjlee168 Homepage – Kaohsiung, Taiwan, 2009-04-05 01:11 (5918 d 16:17 ago) – Posting: # 3467
Views: 13,106

Dear ElMaestro,

We have tried your method with lm() and it worked perfectly and was just exactly the same as the SAS output using balanced and unbalanced data. We also tried the method (this link) pointed by D. Labes, and it didn't work. We're still trying to find out why it didn't work. It could be due to that the dataset we used to test included a missing value. I will post the final result here later if necessary. Now I post the results obtained from your suggested method here with balanced and unbalanced dataset.

Balanced dataset with bear
      Statistical analysis (ANOVA(lm))
[...] **************** Classical (Shortest) 90% C.I. for lnCmax *****************

  CI90_lower Point_estimated CI90_upper
1     97.550         106.184    115.582 [...]  
**************** Classical (Shortest) 90% C.I. for lnAUC0t ****************

  CI90_lower Point_estimated CI90_upper
1     94.459         106.640    120.391 [...]
**************** Classical (Shortest) 90% C.I. for lnAUC0INF **************

  CI90_lower Point_estimated CI90_upper
1     94.478         106.335    119.681[...]

Balanced dataset with SAS only showing Cmax here
[...]Least Squares Means

                                                                 H0:LSMean1=
                                                       lncmax      LSMean2
                                         drug          LSMEAN       Pr > |t|

                                         1         7.33071429         0.2313
                                         2         7.39071429


                                                lncmax
                                  drug          LSMEAN      90% Confidence Limits

                                  1           7.330714        7.270747     7.390681
                                  2           7.390714        7.330747     7.450681

Least Squares Means for Effect drug

                                             Difference
                                                Between    90% Confidence Limits for
                                 i    j           Means       LSMean(i)-LSMean(j)

                                 1    2       -0.060000       -0.144806     0.024806


                                             drug coding: 1-Ref. & 2-Test
[...] Summary:
lnCmax:
point->106.184
lowerCI->97.55
upperCI->115.582

lnAUC0t:
point->106.640
lowerCI->94.459
upperCI->120.391

lnAUC0inf:
point->106.335
lowerCI->94.478
upperCI->119.681


Unbalanced dataset using bear
[...] **************** Classical (Shortest) 90% C.I. for lnCmax ****************

  CI90_lower Point_estimated CI90_upper
1     96.260         105.553    115.745[...] **************** Classical (Shortest) 90% C.I. for lnAUC0t ****************

  CI90_lower Point_estimated CI90_upper
1     92.360         105.202    119.830[...] **************** Classical (Shortest) 90% C.I. for lnAUC0INF **************

  CI90_lower Point_estimated CI90_upper
1     92.419         104.915    119.100[...]


Unbalanced with SAS only showing Cmax here
{...]
Least Squares Means

                                                                 H0:LSMean1=
                                                       lncmax      LSMean2
                                         drug          LSMEAN       Pr > |t|

                                         1         7.32916667         0.3149
                                         2         7.38321429


                                                lncmax
                                  drug          LSMEAN      90% Confidence Limits

                                  1           7.329167        7.263994     7.394339
                                  2           7.383214        7.318042     7.448387


                                         Least Squares Means for Effect drug

                                             Difference
                                                Between    90% Confidence Limits for
                                 i    j           Means       LSMean(i)-LSMean(j)

                                 1    2       -0.054048       -0.146215     0.038120

                                             drug coding: 1-Ref. & 2-Test  [...] Summary as follows:
lnCmax:
point->105.55
lowerCI->96.260
upperCI->115.745

lnAUC0t:
point->105.202
lowerCI->92.36
upperCI->119.83

lnAUC0inf:
point->104.915
lowerCI->92.419
upperCI->119.100


❝ *: Didn't dlabes some time last year explain this? I remember I tried

❝ really hard to challenge lm with an unbalanced dataset to see if I could

❝ get it "wrong" with lm due to imbalance in a 2,2,2-study, but I failed to

❝ do that.


Yes, you're right. Sorry about my poor memory. Thank you for your suggestions.

All the best,
-- Yung-jin Lee
bear v2.9.2:- created by Hsin-ya Lee & Yung-jin Lee
Kaohsiung, Taiwan https://www.pkpd168.com/bear
Download link (updated) -> here

Complete thread:

UA Flag
Activity
 Admin contact
23,424 posts in 4,927 threads, 1,671 registered users;
16 visitors (0 registered, 16 guests [including 12 identified bots]).
Forum time: 17:29 CEST (Europe/Vienna)

Truth and clarity are complementary.    Niels Bohr

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