Bioequivalence and Bioavailability Forum • lme() in bear

lme() in bear [🇷 for BE/BA]

posted by ElMaestro – Denmark, 2015-04-18 13:28 (3661 d 08:47 ago) – Posting: # 14698
Views: 29,075

Hi yung-jin,

❝ ❝ 1. Could you try and fit fixed effects as log(Cmax) ~ 0+ drug + seq + prd, and then extract your drug effects directly from that.

❝

❝ OK. add an intercept?

No, take away an intercept.

❝ Fixed effects: log(Cmax) ~ 0 + seq + prd + drug

In a sense this isn't what I suggested. Please try and fit it with drug specified as the first effect. In stead of
log(Cmax) ~ 0 + seq + prd + drug
see if you can make sure it somehow fits
log(Cmax) ~ 0 + drug + prd + seq
or equally good
log(Cmax) ~ 0 + drug + seq + prd

order of effects has lexical importance and determines how many effect estimates you get per term. Accordingly, when you fit sequence as the first term you will not see one effect estimate for each drug.

Good luck and sorry if my initial suggestion was not too well described.

—
Pass or fail!
ElMaestro

Thread locked

Complete thread:

Reproducing Bear results in Phoenix mittyri 2015-04-14 08:35
- Bug in bear? Helmut 2015-04-14 14:07
  - Bug in bear? mittyri 2015-04-14 20:43
    - Bug in bear? Helmut 2015-04-15 02:17
  - use the same dataset? yjlee168 2015-04-20 10:13
- Bugs in bear with replicated demo data set? yjlee168 2015-04-15 09:38
- Bugs fixed in bear? yjlee168 2015-04-17 01:21
  - Bugs fixed in bear? Helmut 2015-04-17 14:25
    - Bugs fixed in bear? yjlee168 2015-04-17 18:50
      - Bugs fixed in bear? ElMaestro 2015-04-17 23:14
        
        lme() in bear yjlee168 2015-04-18 11:12
        
        lme() in bearElMaestro 2015-04-18 11:28
        
        lme() in bear yjlee168 2015-04-18 12:02
        
        lme() in bear ElMaestro 2015-04-18 13:42
        
        Mean means Helmut 2015-04-18 14:08
        
        Mean means ElMaestro 2015-04-18 16:01
        
        Mean means Helmut 2015-04-19 00:59
        
        Mean means ElMaestro 2015-04-19 09:14
        
        LL and AIC Helmut 2015-04-19 11:08
        
        Confused ElMaestro 2015-04-19 11:42
        
        lsmeans for mixed model in R yjlee168 2015-04-19 23:36
        
        lsmeans() & lme() Helmut 2015-04-20 01:28
        
        lsmeans() & lme() yjlee168 2015-04-20 10:23
        lsmeans() & lme() ElMaestro 2015-04-20 10:35
        
        keep it simple! Helmut 2015-04-20 14:34
        
        keep it simple! ElMaestro 2015-04-20 15:39
        
        ML vs. REML Helmut 2015-04-20 16:35
        
        lme() with NA in R yjlee168 2015-04-19 23:44
      - Dataset Helmut 2015-04-18 13:30
        
        Dataset and lme() in bear yjlee168 2015-04-18 23:14