Since it is not a bug it /must/ be a feature! [Software]
❝ And I have a question about Studentized residuals (Outlier analysis for 4-Period study, EMA): why are we filtering the predicted values as a result of LME object 'LME Ref variability' only for periods 1 and 2 (2 sequences only for 2x2x4)?
❝ Filter->(nseq=1 and period=2) or (nseq=2 and period=1) or (nseq=3 and period=1)
Don’t get confused by the red part: The template handles also the partial replicate where we have three sequences. If you have 4-period full replicate data, only the green part will be evaluated. Navigate in the Object Browser to
Analysis
└─ ABEL
└─ Outlier Analysis
├─ Prep Steps
│ ├─ Complete reference data_Result
│ ├─ LME Ref variability_Residuals
│ ├─ LME Ref variability_Final Variance Parameters
│ ├─ Subjects by Sequence and Period
│ ├─ Join _ Residuals and Nk
│ └─ Join MSE_MSintra
└─ Studentized Residuals → Step 1 Filter
Options
below check ☑ Retain Intermediate Results
and execute.In
Results:Intermediate Results:Step 1 Filter.Results Worksheet
you see the results of all subjects who had two administrations of R (at this point incomplete data are already excluded) coming from either the first period (in sequence RTRT) or the second (in sequence TRTR). Works as designed. 
❝ What about periods 3 and 4?
The mean of residuals within particular subjects should be zero. Therefore, if one is examined the other one is not informative since it is only of the opposite sign (|Res1| = |Res3| ∨ |Res2| = |Res4|).
Check above the object
LME Ref variability_Residuals
where you have both raw residuals (before studentization: Predicted–Observed
), f.i. for the EMA’s “dataset I”:Subject Period Sequence Observed Predicted Residual
1 1 BABA 7.734541 7.472175 0.262366
1 3 BABA 7.204848 7.467214 –0.262366
2 2 ABAB 7.859143 7.721813 0.137330
2 4 ABAB 7.788792 7.926122 –0.137330
With my setup you should get two subjects tagged as potential outliers (#45 with –5.25 and #52 with +3.21). Modify the filter to
(nseq=1 and period=4) or (nseq=2 and period=3)
and see what happens: Same subjects tagged, same |studentized residuals| but with opposite signs compared to the original setup. ч.т.д.BTW, expectations {μ, σ²} of studentized residuals are {0, 1}. On my machine (Win7 64bit, both with the 64 and 32bit versions of Phoenix/WinNonlin 6.4.0.768) I get 9.577·10–16 and 1.014 for this dataset. Not that bad.
- […] it can be seen that there is a perfect negative correlation (i.e., ρ = –1) between the two intra-subject residuals within the same subject. Thus, it is sufficient to examine intra-subject residuals at the first period […]
Dif-tor heh smusma 🖖🏼 Довге життя Україна!
![[image]](https://static.bebac.at/pics/Blue_and_yellow_ribbon_UA.png)
Helmut Schütz
![[image]](https://static.bebac.at/img/CC by.png)
The quality of responses received is directly proportional to the quality of the question asked. 🚮
Science Quotes
Complete thread:
- Templates in Phoenix/WinNonlin EMA ABEL – Residuals? mittyri 2016-03-19 22:27 [Software]
- Since it is not a bug it /must/ be a feature!Helmut 2016-03-20 00:34