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

posted by Helmut Homepage – Vienna, Austria, 2020-07-22 17:37 (857 d 18:27 ago) – Posting: # 21779
Views: 2,517

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 🖖🏼 Довге життя Україна! [image]
Helmut Schütz
[image]

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

Complete thread:

UA Flag
Activity
 Admin contact
22,417 posts in 4,692 threads, 1,596 registered users;
17 visitors (0 registered, 17 guests [including 12 identified bots]).
Forum time: 11:05 CET (Europe/Vienna)

To know that we know what we know,
and to know that we do not know what we do not know,
that is true knowledge.    Nicolaus Copernicus

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