SE of ∆ [General Sta­tis­tics]

posted by Helmut Homepage – Vienna, Austria, 2020-10-03 00:29 (1272 d 17:39 ago) – Posting: # 21966
Views: 2,111

Dear Detlew,

❝ I think all the confusion comes from that sigmaw, sigmad, sigmadelta values including their estimates which are used by all the authors cited within this thread in a different meaning.


Quite possible.

❝ I'm not able to figure out who is who, what is what. Sorry.


F**ing terminology.

❝ The only thing I'm convinced of is that your formula (2) above is correct.

❝ If you write the confidence interval for the BE decision as

❝ PE(T-R) +- SD(d)*tval(0.95, df)


Exactly.

❝ The rest of your algebra is straight forward.

❝ And correct if you ask me ;-).


THX. Now three people agree. I even didn’t trust my rusty algebra and asked Maxima for help:

[image]


❝ BTW: the formula (2) is not the error term in the 2×2×2 crossover.


How would you call it? We use \((3)\) in PowerTOST’s BE_CI.R line 30:

sqrt(mse*ades$bkni*nc)

where for design = "2x2" ades$bkni is 0.5 and nc is sum(1/n[1]+1/n[2]).
If we agree that \(\small{\widehat{\sigma}_\textrm{w}=\sqrt{MSE}}\) * we end up with \((2)\):

      For the special case of a two-treatment, two-period crossover study in which \(\small{n_\textrm{1}}\) subjects receive the test formulation in period one and the reference formulation in period two, while \(\small{n_\textrm{2}}\) subjects receive the reference formulation in period one and the reference formulation in period two, the unbiased estimator is given by $$\small{Est.=\frac{\left(\bar{X}_\textrm{T1}+\bar{X}_\textrm{T2} \right)}{2}-\frac{\left(\bar{X}_\textrm{R1}+\bar{X}_\textrm{R2} \right)}{2}}$$ where
      \(\small{\bar{X}_\textrm{T1}=}\) the observed mean of the \(\small{n_\textrm{1}}\) observations of the test formulation in period one.
      \(\small{\bar{X}_\textrm{T2}=}\) the observed mean of the \(\small{n_\textrm{2}}\) observations of the test formulation in period two.
      \(\small{\bar{X}_\textrm{R1}=}\) the observed mean of the \(\small{n_\textrm{2}}\) observations of the reference formulation in period one.
      \(\small{\bar{X}_\textrm{R2}=}\) the observed mean of the \(\small{n_\textrm{1}}\) observations of the reference formulation in period two.
The standard error of this estimator is $$\small{SE=s\sqrt{\frac{1}{2}\left (\frac{1}{n_\textrm{1}}+\frac{1}{n_\textrm{2}} \right)}}$$ where […] \(\small{s}\) is the square root of the “error” mean square from the crossover analysis of variance, based on \(\small{\nu}\) degrees of freedom.




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,957 posts in 4,819 threads, 1,638 registered users;
80 visitors (0 registered, 80 guests [including 8 identified bots]).
Forum time: 17:08 CET (Europe/Vienna)

Nothing shows a lack of mathematical education more
than an overly precise calculation.    Carl Friedrich Gauß

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