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

posted by Helmut Homepage – Vienna, Austria, 2014-02-21 20:25 (3688 d 19:20 ago) – Posting: # 12472
Views: 20,334

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:

[image] [image]

[image]   [image]

[image]

In Phoenix Pre-Release V1.4 (Build 6.4.0.511, 2013-12-20) still no overlays in box plots and his­to­grams. :lookaround: 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 🖖🏼 Довге життя Україна! [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,636 registered users;
117 visitors (0 registered, 117 guests [including 6 identified bots]).
Forum time: 15:45 CET (Europe/Vienna)

With four parameters I can fit an elephant,
and with five I can make him wiggle his trunk.    John von Neumann

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