lmer / lme [General Sta­tis­tics]

posted by Helmut Homepage – Vienna, Austria, 2019-12-25 19:12 (279 d 14:57 ago) – Posting: # 21030
Views: 3,208

Hi Nastia,

» I still struggle to understand what is the difference between lmer and lme in this case...

The syntax of the formula for random effects is different.

mod <- lmer(logPK  ~ sequence + period + treatment +
                     (1|subject),
                     data = data)
mod <- lme(log(PK) ~ sequence + period + treatment,
                     random = ~1|subject,
                     data = data, na.action = na.omit)

Important: If you have incomplete data, in lme() you have to use na.action = na.omit because the default na.action = na.fail would stop and throw an error.

The method to extract the PE and calculate the CI is also different. Hint: Lines 53–56 and 80–84 of the source-code of method.B.R.

Dif-tor heh smusma 🖖
Helmut Schütz
[image]

The quality of responses received is directly proportional to the quality of the question asked. 🚮
Science Quotes

Complete thread:

Activity
 Admin contact
21,087 posts in 4,398 threads, 1,468 registered users;
online 7 (1 registered, 6 guests [including 3 identified bots]).
Forum time: Wednesday 11:09 CEST (Europe/Vienna)

When puzzled, it never hurts to read the primary documents 
a rather simple and self-evident principle that has, nonetheless,
completely disappeared from large sectors
of the American experience.    Stephen Jay Gould

The Bioequivalence and Bioavailability Forum is hosted by
BEBAC Ing. Helmut Schütz
HTML5