Huge gap in my understanding [General Sta­tis­tics]

posted by martin  – Austria, 2018-03-06 15:40 (2620 d 14:14 ago) – Posting: # 18498
Views: 5,784

Dear ElMaestro,

You may have a look at R package lsmeans for more details and this R code may help (credits to Alex) in understanding.

library(lsmeans)

### Covariance example (from Montgomery Design (8th ed.), p.656)
print(fiber)

### model
fiber.lm <- lm(strength ~ diameter + machine, data = fiber)
summary(fiber.lm)

# means versus ls-means
fiber$pred <- predict(fiber.lm, list(machine=fiber$machine, diameter=rep(mean(fiber$diameter), 15)))
aggregate(fiber$strength, by=list(fiber$machine), mean) # mean
aggregate(fiber$pred, by=list(fiber$machine), mean)     # lsmean

# lsmeans via R package lsmeans
fiber.lsm <- lsmeans(fiber.lm, "machine")
fiber.lsm

Complete thread:

UA Flag
Activity
 Admin contact
23,424 posts in 4,927 threads, 1,671 registered users;
110 visitors (0 registered, 110 guests [including 9 identified bots]).
Forum time: 06:54 CEST (Europe/Vienna)

Freedom is always and exclusively
freedom for the one
who thinks differently.    Rosa Luxemburg

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