lmer / lme [General Sta­tis­tics]

posted by Helmut Homepage – Vienna, Austria, 2019-12-25 20:12 (2352 d 19:23 ago) – Posting: # 21030
Views: 13,021

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 🖖🏼 Довге життя Україна! [image]
Helmut Schütz
[image]

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

Complete thread:

UA Flag
Activity
 Admin contact
23,653 posts in 4,991 threads, 1,570 registered users;
413 visitors (0 registered, 413 guests [including 20 identified bots]).
Forum time: 16:35 CEST (Europe/Vienna)

I’m all in favor of the democratic principle
that one idiot is as good as one genius, but I draw the line
when someone takes the next step and concludes
that two idiots are better than one genius.    Leo Szilard

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