Freedom for replicate bears [🇷 for BE/BA]

posted by ElMaestro  – Denmark, 2011-01-27 22:50 (5270 d 03:11 ago) – Posting: # 6502
Views: 67,387

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
23,425 posts in 4,928 threads, 1,683 registered users;
43 visitors (0 registered, 43 guests [including 11 identified bots]).
Forum time: 03:02 CEST (Europe/Vienna)

You can do one of two things; just shut up,
which is something I don’t find easy,
or learn an awful lot very fast,
which is what I tried to do.    Jane Fonda

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