Bioequivalence and Bioavailability Forum

Phoenix: resid­u­als [General Sta­tis­tics]

posted by Helmut  – Vienna, Austria, 2014-02-21 20:25 (4079 d 21:57 ago) – Posting: # 12472
Views: 25,362

Hi Lucas,

in the following a setup for PHX6.3+ (Clayton/Leslie’s data given in Chow/Liu’s table 6.9.1; validated against table 8.2.3); sequences coded as TR|RT. It seems complicated first, but works also for unba­lanced sequences. This is not PHX’ standard model, but all effects fixed (like for EMA and in SAS Proc GLM). You have to code it only once and can use it as a template in other studies.

1. Setup the BE wizard for ln-transformed data.
   Fixed effects: subject+period+treatment+sequence+subject(sequence),
   delete Random1: subject(sequence)
   
2. Execute and navigate to Output Data > Residuals > Send To > Data > Data Wizard > Add Filter
   • Exclude everything except Dependent, Subject, Predicted, Residual (i.e., keep them mapped to None)
     
   • Built In > Add > Action: Include | Find: RT | Search Area Column: Sequence > OK
     
   • Built In > Add > Action: Include | Find: 1 | Search Area Column: Period > OK
     
   • Execute and rename to RT Filter
3. Final Results > Result > Send To > Data > Split Worksheet
   • Map Sequence to Sort
     
   • ☐ Create unique values worksheets (uncheck the box!)
   • Execute and rename to RT Count
4. Navigate back to Bioequivalence > Output Data > Final Variance Parameters > Send To > Data > Join Worksheets
   • Map Dependent to Sort, Estimate to Source Column, leave Parameter mapped to None
     
   • Map Worksheet 2 > RT Filter.Result
     
   • Map Dependent to Sort and all variables to Source Column
     
   • Execute and rename to RT MS intra
5. Final Results > Result > Send To > Data > Join Worksheets
   • Map Sequence to Sort and all other variables to Source Column
     
   • Map Worksheet 2 > RT Count.Result
     
   • Map Sequence to Sort and Count to Source Column
     
   • Execute and rename to RT data
6. Final Results > Result > Send To > Data > Data Wizard > Add Transformation
   • Transformation Type: Custom | Transformation: Custom Function | New Column Name: StudResIntra |
     Formula: Residual/sqrt(Estimate*(Count-1)/(2*Count))
     
   • Add Filter > Exclude Residual, Estimate, Count
     
   • Execute and rename to RT Stud Res intra
7. Repeat steps #2 – #6 for sequence TR. You should end up with TR Filter, TR Count, TR MS intra, TR data, TR Stud Res intra.
   
8. Navigate back to RT Stud Res intra > Final Results > Result > Send To > Data > Append Worksheets
   • Map all variables to Source Column
     
   • Map Worksheet 2 > TR Stud Res intra.Result
     
   • Map all variables to Source Column
     
   • Execute and rename to Stud Res intra

Using the standard plotting functions (see above), we get:

 

   



In Phoenix Pre-Release V1.4 (Build 6.4.0.511, 2013-12-20) still no overlays in box plots and his­to­grams. Maybe I will post some R-code later.

---

Eidt: I have to correct myself: In Phoenix (all versions) these are standardized (sometimes call internally studentized) residuals, not studentized (externally studentized) residuals.

—
Dif-tor heh smusma 🖖🏼 Довге життя Україна! 
Helmut Schütz


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

Complete thread:

• Qualitative analysis of the ANOVA residuals Lucas 2014-02-19 14:27 [General Sta­tis­tics] 
  • Qualitative analysis of the ANOVA residuals Helmut 2014-02-19 15:22
    • Qualitative analysis of the ANOVA residuals Lucas 2014-02-19 15:58
      • Qualitative analysis of the ANOVA residuals Helmut 2014-02-19 16:25
        • Qualitative analysis of the ANOVA resid­uals Lucas 2014-02-19 17:20
        • Qualitative analysis of the ANOVA resid­u­als Lucas 2014-02-20 19:07
          • First hints for Phoenix Helmut 2014-02-20 20:11
          • Phoenix: resid­u­alsHelmut 2014-02-21 19:25
            • Phoenix: resid­u­als BRN 2014-03-04 15:15
              • Phoenix: resid­u­als Helmut 2014-03-04 16:52
            • Phoenix: resid­u­als Rosy 2018-10-15 07:03
              • Software validation Helmut 2018-10-15 11:43
                • Software validation Rosy 2018-10-18 00:38
                  • Avoid raw residuals Helmut 2018-10-18 17:14
                    • Avoid raw residuals Rosy 2018-10-18 18:05
          • Phoenix → R Helmut 2014-02-23 02:20
            • Phoenix → R Lucas 2014-02-25 19:28
  • Qualitative analysis of the ANOVA residuals ElMaestro 2014-02-19 15:32
    • Qualitative analysis of the ANOVA residuals Helmut 2014-02-19 16:06
    • Qualitative analysis of the ANOVA residuals Lucas 2014-02-19 16:13
      • Qualitative analysis of the ANOVA resid­u­als ElMaestro 2014-02-19 20:40