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

posted by ElMaestro  – Denmark, 2015-04-20 12:35 (3583 d 07:12 ago) – Posting: # 14714
Views: 27,832

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:

UA Flag
Activity
 Admin contact
23,376 posts in 4,912 threads, 1,662 registered users;
255 visitors (0 registered, 255 guests [including 21 identified bots]).
Forum time: 18:48 CET (Europe/Vienna)

There are sadistic scientists who hurry to hunt down errors
instead of establishing the truth.    Marie Curie

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