Bear to bear interval with 90% confidence [🇷 for BE/BA]
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
Balanced dataset with SAS only showing Cmax here
Unbalanced dataset using bear
Unbalanced with SAS only showing Cmax here
Yes, you're right. Sorry about my poor memory. Thank you for your suggestions.
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
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:
- Bear to bear interval with 90% confidence d_labes 2009-04-01 17:00
- Bear to bear interval with 90% confidence ElMaestro 2009-04-01 19:01
- Bear to bear interval with 90% confidence yjlee168 2009-04-01 20:38
- Bear to bear interval with 90% confidence ElMaestro 2009-04-01 21:28
- Bear to bear interval with 90% confidence yjlee168 2009-04-01 20:38
- Bear to bear interval with 90% confidence yjlee168 2009-04-01 20:15
- Bear to bear interval with 90% confidence ElMaestro 2009-04-01 21:47
- Bear to bear interval with 90% confidenceyjlee168 2009-04-04 23:11
- Bear to bear interval with 90% confidence ElMaestro 2009-04-05 13:23
- Bear to bear interval with 90% confidence yjlee168 2009-04-06 07:50
- Data set d_labes 2009-04-07 08:46
- Data set yjlee168 2009-04-07 13:44
- Bear to bear interval with 90% confidence ElMaestro 2009-04-05 13:23
- Bear to bear interval with 90% confidenceyjlee168 2009-04-04 23:11
- Models (not necessarily nice looking young woman) d_labes 2009-04-06 10:49
- Models (not necessarily nice looking young woman) yjlee168 2009-04-07 14:05
- Bear to bear interval with 90% confidence ElMaestro 2009-04-01 21:47
- Bear to bear interval with 90% confidence ElMaestro 2009-04-01 19:01