---
title: PhregプロシジャでカスタムHRを出す方法
tags:  #[第4回大阪sas勉強会]  
author: [森岡裕[SASユーザー総会世話人]](https://www.docswell.com/user/6484025)
site: [Docswell](https://www.docswell.com/)
thumbnail: https://bcdn.docswell.com/page/3EKD5L9GJD.jpg?width=480
description: [第4回大阪sas勉強会]
published: April 10, 24
canonical: https://www.docswell.com/s/6484025/KNRYNW-2024-04-10-112218
---
# Page. 1

![Page Image](https://bcdn.docswell.com/page/3EKD5L9GJD.jpg)

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


# Page. 2

![Page Image](https://bcdn.docswell.com/page/L73RK3W5J5.jpg)

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


# Page. 3

![Page Image](https://bcdn.docswell.com/page/87DY34KY7G.jpg)

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


# Page. 4

![Page Image](https://bcdn.docswell.com/page/VJP64MK2E8.jpg)

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


# Page. 5

![Page Image](https://bcdn.docswell.com/page/2EV3X9VX7Q.jpg)

カスタム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


# Page. 6

![Page Image](https://bcdn.docswell.com/page/57G3VZLRJL.jpg)

カスタム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


# Page. 7

![Page Image](https://bcdn.docswell.com/page/4EQ56LYYJP.jpg)

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