PhregプロシジャでカスタムHRを出す方法

PHREG プロシジャでカスタムHRを出す方法

よくあること  Time to event型のエンドポイント試験において「投与群と因子の交互作用もみてくれない？」とよく聞 かれる  交互作用を考慮したハザード比はphregプロシジャのデフォルトではでてこない

CONTRAST STATEMENTを使い解決 因子がカテゴリ変数だった場合を想定（数値でも応用可） まず、SASがカテゴリ変数をどうコーディングするか確認する 例： TRTP （A=Active, P=Placebo） SEX (F=Female , M=Male) SAS内部で0，1のダミー変数（design variable）を作ってくれる Classステートメントはアルファベットなら後、数字なら大きいものを参照レベルにする(０とする) ClassステートメントでRef=オプションを使うことで参照レベルを指定することもできる

SAS内部でDESIGN MATRIXを作ってくれる  すべての水準を0，1で表すとTRTP（2水準）×SEX（2水準）の4パターンがある  Design matrixの行と行を引くことでカスタムHRを出せる  水準が3ある場合のDesign matrix statement # TRTP SEX β(trtp) β(sex) β(trtp)*β(sex) 1A M 1 0 0 2P M 0 0 0 3A F 1 1 1 4P F 0 1 0

カスタムHR CASE 1 男性被験者においてACTIVEとPLACEBOを比べる  まず比べたい行を特定する。今回の場合Statement#１と２を比べる  行１ － 行２すると（1, 0, 0）となる  カスタムHRとp値がでる statement # TRTP SEX β(trtp) β(sex) β(trtp)*β(sex) 1A M 1 0 0 2P M 0 0 0 3A F 1 1 1 4P F 0 1 0

カスタムHR CASE 2 ACTIVE群内で男女の比較をする  比べたい行を特定する。今回の場合Statement#１と３を比べる  行3 － 行1すると（0, 1, 1）となる  カスタムHRとp値がでる statement # TRTP SEX β(trtp) β(sex) β(trtp)*β(sex) 1A M 1 0 0 2P M 0 0 0 3A F 1 1 1 4P F 0 1 0

カスタムHR CASE 2 チェック ACTIVE群内で男女の比較をする