replicateBE solution with interactions [General Sta­tis­tics]

posted by mittyri – Russia, 2019-12-23 15:30 (1188 d 01:46 ago) – Posting: # 21020
Views: 6,023

Dear Nastya,

whenever possible use simple solutions, at least till the moment you must dive into mixed models galaxy ;-)

And more: interaction terms are not the same in SAS and R
take a look at method.B implementation inside replicateBE package
I just modified it a bit

library(replicateBE)
library(nlme)
alpha <- 0.05
options(contrasts=c("contr.treatment","contr.poly"))
Dataset <- rds01
Dataset$group<- factor(ifelse(as.numeric(levels(Dataset$subject))[Dataset$subject]<31, 1,2))
modB <- lme(log(PK) ~ sequence + group + sequence:group + period + period%in%group + treatment, random = ~1 | subject, na.action = na.omit, data = Dataset)

EMA.B <- summary(modB)
PE <- EMA.B$tTable["treatmentT", "Value"] # exp(PE) = 1.157275
CI <- exp(PE + c(-1, +1) * qt(1 - alpha, EMA.B$tTable["treatmentT", "DF"]) * EMA.B$tTable["treatmentT", "Std.Error"])
# CI = 1.071136 1.250340

Kind regards,
Mittyri

Complete thread:

UA Flag
Activity
 Admin contact
22,550 posts in 4,724 threads, 1,606 registered users;
18 visitors (0 registered, 18 guests [including 13 identified bots]).
Forum time: 17:16 CET (Europe/Vienna)

If there is an exception to any rule,
and if it can be proved by observation,
that rule is wrong.    Richard Feynman

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