Bioequivalence and Bioavailability Forum • Modelling Parallel bears

Modelling Parallel bears [🇷 for BE/BA]

posted by yjlee168 – Kaohsiung, Taiwan, 2010-04-23 23:14 (5474 d 12:34 ago) – Posting: # 5207
Views: 58,040

Dear d_labes,

Thank you for your message.

❝ IMHO definitely yes.

O.k.

❝ To ask it again: What is the idea behind your model, where does it come from? Do you have a reference for me? lm() or lme() is only a question of considering some effects as random or not.

No reference for this. Simply try to translate the codes of Power to know into R. When we were developing bear for parallel BE study, we tried to get exactly the same results as Power to know or WinNonlin. The codes of Power to know for parallel BE study looks like


PROC GLM DATA=EXAMPLE;

CLASS SUBJ TRT;

MODEL LAUCT LAUCI LCMAX=TRT;

ESTIMATE ‘A vs. B’ TRT 1–1;

LSMEAN TRT;

RUN;

So we used lme() (not lm()) to meet this requirement for validation purpose.
lmeCmax_ss<-lme(Cmax_ss ~ drug, random=~1|subj, data=Data, method="REML" )
Just like analyzing a replicate crossover. I can not find the results obtained from analyzing with WinNonlin using a linear mixed effects model. As far as I could remember, we got pretty similar results as with WinNonlin at that moment. Then when search the answer for your question int his thread, I find that we probably can simply use lm() as what we do in non-replicate 2x2x2 crossover BE study, such as
lmCmax_ss<- lm(Cmax_ss ~ drug , data=TotalData)
Which model is more appropriate for analyzing parallel BE study with R?

—
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:

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 bearsyjlee168 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
        Dataset Helmut 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