Yeah but, no but, yeah but… [General Sta­tis­tics]

posted by Helmut Homepage – Vienna, Austria, 2020-07-22 15:37 (100 d 00:07 ago) – Posting: # 21779
Views: 1,511

Hi ElMaestro,

I must confess that I don’t have the slightest idea what you have done here. Not surprising. Hazelnut-sized brain < walnut-sized brain. However:

» I was happy to look up how Chow & Liu see it; see e.g. formula 2.5.1 and 9.1.1 in the 3rd edition
$$Y_{ijk}=\mu+S_{ik}+P_j+F_{(j,k)}+{\color{Red}{C_{(j-1,k)}}}+e_{ijk} \tag{2.5.1=9.1.1}$$Did you implement the first-order carryover as well? If yes, try it without. Might explain the differences below.

» Here's my result with EMA's dataset II, if it is of interest:
»  Var.Component  Ini.value      Value
»          varWR 0.01240137 0.01211072

OK, you start with the EMA’s approach.
library(replicateBE)
CV.wR <- 0.01*method.A(data = rds02, print = FALSE, details = TRUE)[["CVwR(%)"]]
cat("varwR", signif(log(CV.wR^2+1), 7), "\n")
varwR 0.01240137

Anything goes. (© Paul Feyerabend)

With your previous REML-code I got 0.01324648 and following the FDA’s approach (intra-subject contrasts) I got 0.01298984 (in Phoenix, SAS, and by my [image]-code). Hence, I see two problems:
  1. Your result does not match the FDA’s approach. Your previous result is even closer (+2.0%) than the new one (–6.8%).
  2. Even if it would match, how would you code the FDA’s mixed model for ABE in [image]?
    That’s the most important question.

Dif-tor heh smusma 🖖
Helmut Schütz
[image]

The quality of responses received is directly proportional to the quality of the question asked. 🚮
Science Quotes

Complete thread:

Activity
 Admin contact
21,180 posts in 4,415 threads, 1,476 registered users;
online 9 (0 registered, 9 guests [including 6 identified bots]).
Forum time: Friday 14:45 CET (Europe/Vienna)

The idea is to try and give all the information to help others
to judge the value of your contribution;
not just the information that leads to judgment
in one particular direction or another.    Richard Feynman

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