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

posted by Helmut Homepage – Vienna, Austria, 2014-02-21 19:25 (2364 d 14:31 ago) – Posting: # 12472
Views: 16,400

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 🖖
Helmut Schütz
[image]

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

Complete thread:

Activity
 Admin contact
21,009 posts in 4,379 threads, 1,460 registered users;
online 28 (0 registered, 28 guests [including 12 identified bots]).
Forum time: Thursday 10:56 CEST (Europe/Vienna)

If you obey all the rules,
you will miss all the fun.    Katharine Hepburn

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