Mixed models, did someone figure out lme? [R for BE/BA]

posted by PharmCat  – Russia, 2020-08-10 22:19 (251 d 02:48 ago) – Posting: # 21853
Views: 1,376

» Hi all,

Hi ElMaestro, it's me again.

Bates went to Julia and work on MixedModels.jl. I and many people asked how to make custom G and R structure in this package. As I know - there are no method to make some like CSH or VC for R matrix, and all issues on github to make it is closed. :angry: I think this is the position of developers.

MixedModels.jl package for Julia seems like port of lme package for R project, I think lme have equivalent limitation because PIRLS used in both cases.

From my side was made attempt to build package for bioequivalence replicate design in Julia. It work fine for all basic designs (sometimes better than SPSS), but have some limitations: no choice for variance structure, only CSH+VC, 2 formulations, strict model. You can find it here. Source is open and I try to document is clearly.

PS There are some ways to make product for general purpose (like SPSS), but it is very big job even to make good tipization for data and process. Also make good API it's a monster task. So, some more techniks I will try to implement to ReplicateBE.jl (gradient and hessian computation, some structure choice, linear algebra optimizations).

Complete thread:

 Admin contact
21,419 posts in 4,475 threads, 1,509 registered users;
online 6 (0 registered, 6 guests [including 3 identified bots]).
Forum time: Monday 01:08 CEST (Europe/Vienna)

There is one certainty in drug development
and statistics that one can depend on:
the data are always late.    Scott Patterson and Byron Jones

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