Not understood [R for BE/BA]

posted by ElMaestro  – Belgium?, 2020-08-05 06:13 (118 d 18:04 ago) – Posting: # 21824
Views: 6,727

Thanks PharmCat,

funny, I read the ref. you mention (Lindstrom & Bates) yesterday and I did not understand much of it. Will try again.

» Problem1: when you have big matrix invercing is slow as hell, but if you inverce by blocks (you can do it because it block-diagonal) it can be done faster.

Sounds neat, will check if I can somehow do such a thing.

» Problem2: In BE case you have many equal Zi matrices, so in logREML calculation (if you calc by subject) you inverce one matrix many times (because V = ZGZ'+R, Z can be different, G and R not changing subject by subject). If you cache result - you can seriously increase performance.

My approach is very simple: I make a case for saying that when T is not replicated while R is replicated, then it is not possible to make a meaningful model when V=ZGZt+R (as long as we count on R having the withins and G having the betweens). If we insest on ZGZt+R then we get irrelevant variance components associated with T whichever way we parameterise the covariance matrix. Therefore, I am working in V directly, meaning I am writing the blocks into it directly.
This may still be compatible with your proposal, though.

I could be wrong, but...

Best regards,
ElMaestro

No, of course you do not need to audit your CRO if it was inspected in 1968 by the agency of Crabongostan.

Complete thread:

Activity
 Admin contact
21,225 posts in 4,427 threads, 1,482 registered users;
online 4 (1 registered, 3 guests [including 2 identified bots]).
Forum time: Wednesday 00:18 UTC (Europe/Vienna)

A big data-analyst is an expert
in producing misleading conclusions from huge datasets.
It is much more efficient to use a statistician,
who can do the same with small ones.    Stephen Senn

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