Update: S×F vari­­ance test [Regulatives / Guidelines]

posted by AngusMcLean – USA, 2014-12-08 19:18 (3424 d 19:28 ago) – Posting: # 14020
Views: 12,257

I have a recent email from Ana over at Pharsight: Please see below and let me know your thoughts. Evidently there has been confusion in the past and I am wondering whether it has been resolve (below).

"Some additional discussion about the questions that you asked about BE has taken place within Pharsight. I just wanted to give you an update.

1. We entered an enhancement on our tracking system “Provide additional output in Pop/Indiv Bioequivalence - sigmaD, sigmaWT (for subject by formulation interaction variance)”. WinNonlin is calculating this but not providing it in the output.
2. For certain designs the documentation has some information that might be useful to you:

For an RTRT/TRTR study, if the default model for replicated crossovers is used, which has a random model using type FA0(2), look at the Parameter and the Group_Level columns on the Parameter Key output workbook to determine which variance parameter is for the reference formulation and which is for the test formulation. In the following equations, it is assumed that the first output parameter is for the test formulation.

The within-subject variance of the test, sig_WT^2, and the within-subject variance of the reference, sig_WR^2, are in the output as:
sig_WT^2 = Var(Period*Formulation*Subject)_21
sig_WR^2 = Var(Period*Formulation*Subject)_22

The intra-subject CV's can be computed as:
IntraCV_T = 100%*sqrt(exp(sig_WT^2)-1)
IntraCV_R = 100%*sqrt(exp(sig_WR^2)-1)

The between-subject variance of test, sig_BT^2, the between-subject variance of reference, sig_BR^2, and the G matrix for the random model, can be computed as follows. This still assumes that the test formulation precedes the reference. The '_11' notation on the lambda parameters (which indicates that they are for the first random model and first group) has been dropped to improve readability:

G(1,1) = sig_BT^2 = lambda(1,1)^2
G(2,2) = sig_BR^2 = lambda(2,1)^2+lambda(2,2)^2
G(1,2) = rho*sig_BT*sig_BR = lambda(1,1)*lambda(2,1) (= between subject covariance for test and reference)

The subject-by-formulation interaction variance in general is sig_D^2 = sig_BT^2 + sig_BR^2 - rho*sig_BT*sig_BR, and can be computed from the above equations as:
sig_D^2=G(1,1)+G(2,2)-2*G(1,2)

The inter-subject CV's can be computed as:
InterCV_T = 100%*sqrt(exp(sig_BT^2)-1)
InterCV_R = 100%*sqrt(exp(sig_BR^2)-1)

If the CSH variance structure is used instead of FA0(2), the following equations change. This still assumes the first output parameter is for test and the second for reference.

G(1,1) = sig_BT^2 = cshSD(1)^2
G(2,2) = sig_BR^2 = cshSD(2)^2
G(1,2) = cshCorr * cshSD(1) * cshSD(2)"


The parameter name or units is not in the file I have: it is from the FDA template, it is the same file for EMEA.

Here is the output from Phoenix. Now one has to look at Ana's email

Dependent Units Parameter Estimate
Data lambda(1,1)_11 0.85299471
Data lambda(1,2)_11 0.82840678
Data lambda(2,2)_11 8.34E-07
Data Var(Period*Formulation*Subject)_21 0.20211811
Data Var(Period*Formulation*Subject)_22 0.11739421


Herr is my paramater key and the REFERENCE IS first. My file is RTRT

Dependent Units Parameter Source Type Bands Subject_Term Group_Term Group_Level
Data lambda(1,1)_11 Random No-Diag Factor Analytic 2 Subject
Data lambda(1,2)_11 Random No-Diag Factor Analytic 2 Subject
Data lambda(2,2)_11 Random No-Diag Factor Analytic 2 Subject
Data Var(Period*Formulation*Subject)_21 Repeated Variance Components Subject Formulation Formulation_R
Data Var(Period*Formulation*Subject)_22 Repeated Variance Components Subject Formulation Formulation_T


That is why the parameters are switched from yours.


Edit: Merged with two later post. Note that you can edit your original post within 24 hours. [Helmut]
Thread locked

Complete thread:

UA Flag
Activity
 Admin contact
22,993 posts in 4,828 threads, 1,658 registered users;
65 visitors (0 registered, 65 guests [including 6 identified bots]).
Forum time: 15:46 CEST (Europe/Vienna)

So far as I can remember,
there is not one word in the Gospels
in praise of intelligence.    Bertrand Russell

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