Freedom for replicate bears [🇷 for BE/BA]

posted by ElMaestro  – Denmark, 2011-01-27 22:50 (4836 d 01:22 ago) – Posting: # 6502
Views: 64,096

Dear dlabes and others,

❝ After reading again tons of Websites (seems nobody evaluates x-over studies with R) I came up with the following attempt in R:

❝   model2 <- lme(y~ drug + prd + seq,

❝        random= ~ tmt-1|subj,

❝        #different within variabilities

❝        weights=varIdent(form = ~ 1 | drug),

❝        data=BlaBla, method="REML")


I am starting to believe lme does in some cases not really reach the goal, so I switched to lmer (package lme4). Syntactically something like:
ElModel=lmer(logCmax ~ per+seq+trt + (trt-1|subj))

-which achieves the same.
It will occasionally -depending on the dataset- actually fit the data better (that is, to a higher likelihood and thus with better sigma estimates) than lme can.

And also:
Who's got a 3-period dataset (ref. replicated, but not test-replicated) with a known sWR that I can look at?

Mucho gracias
—
Pass or fail!
ElMaestro

Complete thread:

UA Flag
Activity
 Admin contact
22,993 posts in 4,828 threads, 1,655 registered users;
86 visitors (0 registered, 86 guests [including 3 identified bots]).
Forum time: 01:12 CEST (Europe/Vienna)

So far as I can remember,
there is not one word in the Gospels
in praise of intelligence.    Bertrand Russell

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