Bioequivalence and Bioavailability Forum • lsmeans() & lme()

lsmeans() & lme() [🇷 for BE/BA]

posted by ElMaestro – Denmark, 2015-04-20 12:35 (4075 d 19:11 ago) – Posting: # 14714
Views: 32,150

Gentlemen,

❝ lsmeans(muddle, "drug", cov.reduce=F, weights="equal")

❝

❝ drug lsmean SE df lower.CL upper.CL

❝ 1 7.328141 0.02504047 12 7.273582 7.382699

❝ 2 7.391880 0.03485385 12 7.315940 7.467820

❝

❝ Results are averaged over the levels of: seq, prd

❝ Confidence level used: 0.95

Well, it seems you are finally getting there, although it really isn't pretty. I think all this is about treatment effects in perception must be Least Squares Means whatever that term truly means. I don't uncritically subscribe to that view. R actually got those treatment effects right but it did not report what was expected or hoped for by those who had a love affair with SAS' invention.

Note that 'averaged over the levels' part. In a nutshell, this reminds me that one day I will add a package called 'marginal means' and it will do absolutely nothing except alias the lsmeans function to a function called marginal.means and suddenly everything will make a lot more sense. I will receive the Fields Medal for it. At least. Plus 17 Michelin Stars and the Golden Palms.

—
Pass or fail!
ElMaestro

Thread locked

Complete thread:

Reproducing Bear results in Phoenix mittyri 2015-04-14 08:35 [🇷 for BE/BA]
- 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 bear ElMaestro 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