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

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

posted by yjlee168 – Kaohsiung, Taiwan, 2015-04-20 12:23 (3714 d 12:38 ago) – Posting: # 14713
Views: 28,459

Dear Helmut,

Yes, this message is great.

❝ That’s not correct – see help(models)!

Thanks. Got it.

❝ ...

❝ TotalData[, !(names(TotalData) %in% c("AUC0t", "AUC0INF", "lnCmax",

❝ "lnAUC0t","lnAUC0INF"))]

❝ muddle <- lme(log(Cmax) ~ drug+seq+prd, random = ~drug-1|subj,

❝ data=TotalData)

I need one more line to run lsmeans here. muddle.rg1<-ref.grid(muddle,data=TotalData); otherwise, it would crash for no data. I don't know why.

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

❝ Bingo for the means! The SEs are still off (compared to PHX).

❝ Level Estimate StdError DF


❝   R   7.328141 0.039852 12

❝   T   7.391880 0.032656 12

❝ At least something to start from.

OK.

❝

Post #500 this year.

Congratulations! :ok:

—
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

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