Still can't make it work [R for BE/BA]

posted by PharmCat  – Russia, 2020-08-08 01:08 (54 d 04:34 ago) – Posting: # 21838
Views: 5,898

Hi ElMaestro!

» That is why you have five elements in theta (T is not replicated, but you are trying to fit both a within and a between variance for it). I am saying we do not work on them both, in stead we replace their sum by a single variance estimate. This is the whole point of the thread here.

I understand your point. And generaly your variance component function is better because is not redundant. From other side if we have some different functions V = ƒ11) and V = ƒ22) it could be any and from any θ if we get equal V (so it could be CSH, FA0(2), manual construct).

» If you are looking at a final log likehood around 15.337 then I think we are at least talking about the same model.

I have 15.337299690629834, model is the same ;-) it is good proof: different V construction methods lead to same results.

» I guess it must be a little typo in Gurka's paper.

From time to time I find inaccuracies even in good articles and sometime it give strong headake when you try directly reproduce some methods. log(Σ|X|) instead log(|ΣX|) can give hours of frustration.:crying:

» Thanks PharmCat, pointing me towards Gurka and insisting on Xi was a great help:ok:.

I am very glad if my thoughts were useful :clap:

Complete thread:

Activity
 Admin contact
21,090 posts in 4,398 threads, 1,469 registered users;
online 16 (0 registered, 16 guests [including 7 identified bots]).
Forum time: Thursday 05:42 CEST (Europe/Vienna)

In these days, a man who says a thing cannot be done
is quite apt to be interrupted by some idiot doing it.    Elbert Green Hubbard

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