Slightly off topic, but related :-) [Design Issues]

posted by Shuanghe  – Spain, 2019-12-09 13:00 (1624 d 22:12 ago) – Posting: # 20964
Views: 5,975

Hi ElMaestro,

❝ Note e.g. that the two fits have different residuals and residual df's, which to me means incomplete subjects are not deleted (R does not know and is not being told something is incomplete; the full rank design matrix is still invertible and so on).


What about the following?
X3 <- Xm[-7,]
M3 <- lm(logCmax ~ factor(Seq)+factor(Subj)+factor(Trt)+factor(Per), data=X3)
anova(M3)
lsmeans(M3, "Trt")
confint(pairs(lsmeans(M3, "Trt"), reverse =F), level=0.9)

Xm has missing period (1) for subject 4, X3 has no subject 4. compare anova(M2) and anova(M3), residual and df of residual are same. 90% CI also same. So wouldn't it mean that R deleted the extra period (2) of subject 4 in Xm automatically when doing BE evaluation? Lsmeans are different, so subject 4 period 2 was kept for that calculation. I woulds say that this behaviour is the same as SAS.

All the best,
Shuanghe

Complete thread:

UA Flag
Activity
 Admin contact
23,034 posts in 4,835 threads, 1,647 registered users;
34 visitors (0 registered, 34 guests [including 9 identified bots]).
Forum time: 12:12 CEST (Europe/Vienna)

Give me a fruitful error any time, full of seeds, bursting with its own corrections.
You can keep your sterile truth for yourself.    Vilfredo Pareto

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