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

posted by PharmCat  – Russia, 2020-08-08 01:08 (161 d 18:10 ago) – Posting: # 21838
Views: 6,681

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,303 posts in 4,441 threads, 1,487 registered users;
online 7 (0 registered, 7 guests [including 5 identified bots]).
Forum time: Saturday 18:18 CET (Europe/Vienna)

If you can’t solve a problem, then there is an easier problem
you can solve: find it.    George Pólya

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