newlogoPPT-大阪大学大学院基础工学研究科基础工学部.ppt

1、GLM利用反復測定解析,株式会社 SAS統計解析研究室小玉 奈津子,2,(1),始間違一応GLM、一般線形（ Linear Model）解析。 （SAS最古） 一般化線形（Generalized Linear Model) GENMOD利用。,3,GLM(2),GLM 線形回帰分析 （単回帰重回帰 =REG） 分散分析 （場合OK) 多変量分散分析 ( H-F G-G 調整 単変量型 Wilks 多変量型) 変量効果 ( =MIXED） 反復測定分散分析 (.） =MIXED） 普通多重比較,4,経時被験者1次単位、時点被験者2次単位分割実験考。基本的指定 単変量型 被験者入子（） 効果指定

2、TEST 要因検定行,GLM利用 昔指定分割型実験解析,GROUP RESP SUBJECT TIME 1 10 001 1 1 13 001 2 省略 3 23 022 1 3 36 022 2PROC GLM DATA=data1; CLASS time group subject; MODEL resp=group subject(group) time time*group subject(time group); TEST H=group E= subject(group); TEST H=time E= subject(time group); TEST H=time*group

3、E= subject(time group);RUN;,5,GLM利用 反復測定分散分析 (1),基本的指定 多変量型 REPEATED 利用分析指定,GROUP SUBJECT TIME1 TIME2 TIME3 1 001 12 24 55 1 002 34 53 47 2 003 23 32 33 省略PROC GLM DATA=data1; CLASS group; MODEL time1 time2 time3=group; /*REPATED*/ REPEATED time 3 ; RUN;,6,GLM利用 反復測定分散分析 (2),REPEATED 単変量型検定、多変量型検定行

4、構文 REPEATED 要因名 水準 変換 / 例: REPEATED time 3 POLYNOMOAL / PRINTE （ 時点3 直交多項式型対比 / 球面性検定),7,GLM利用 反復測定分散分析 (3),REPATED出力) 各時点分散分析表 時点、通常分散分析行 必要場合、次MODELNOUNI指定。( ODS可能） PROC GLM; CLASS a; MODEL t1-t3=a / NOUNI; REPEATED time; RUN;,8,GLM利用 反復測定分散分析 (4),2) 多変量分散分析（多変量型検定） Wilks Phillai Hotelling-lawley

5、Roy最大根 出力。 * 手法p値異場合。 Wilks良使答 実際？,9,GLM利用 反復測定分散分析 (5),3) 球面性検定（ Sphericity Tests) Mauchly球面性検定 2検定結果出力。 Transformed Variates REPEATED直交対比以外 指定場合利用 （例:PROFILE) Orthogonal Components 対比直交対比選択場合利用 （対比何指定場合、利用）,10,GLM利用 反復測定分散分析 (6),4) 単変量型検定 （修正） 被験者内要因検定対、球面性仮定成立場合利用2調整行。 G-G （Greenhouse-Geisser)調整p

6、値 H-F （Huynh-Feldt)調整p値 （一緒対応値出力）,11,解析例:1 宿題, 分割型解析（持方単変量型） PEPEATED利用単変量型 （持方多変量型） 多変量型 時点間多重比較 単純主効果検定,12,参考,高橋行雄大橋靖雄芳賀敏郎 (1989）SAS実験解析 東京大学出版会（内容古） Littel, Milliken, Stroup,Wolfinger (1996) SAS System for Mixed Models （MIXED入門書） 松山 裕 山口 拓洋 訳 （2001） 医学統計線形混合 社 （混合入門書） 岸本淳司 （1997) 経時解析 第回 日本SAS会論文

7、集（例貰） 愛知学院大学千野直仁先生 反復測定分散分析 基本応用 （参考） http:/www.aichi-gakuin.ac.jp/%7Echino/anova/contents.html,13,理論聞場合、狩野先生細SAS統計的話場合明日来岸本（Boss）使用法、,質疑応答,15,1Version8便利機能ODS,ODS（） SAS出力結果制御新機能。 ODS利用、今出 GLM結果（例、最小2乗平均） 出力。 出力形式 * HTML形式 *RTF形式 *SAS ,16,GLM利用 反復測定分散分析 2,REPEATED主変換 CONTRAST （基準水準） 水準対象水準場合利用。 （ RE

8、PEATED、最後水準基準。） 例 REPEATED time 3 CONTRAST(1) 指定場合、 M行列以下作成。 t1 t2 t3treatmnt_2 -1.000000000 1.000000000 0.000000000treatmnt_3 -1.000000000 0.000000000 1.000000000,17,GLM利用 反復測定分散分析 , IDENTITY 変換 各時点別変数扱検定行。 例: REPEATED time 3 IDENTITY 指定場合、対比行列以下作成。 t1 t2 t3 treatmnt_1 1.000000000 0.000000000 0.000000000 treatmnt_2 0.000000000 1.000000000 0.000000000 treatmnt_3 0.000000000 0.000000000 1.000000000,18,GLM利用 反復測定分散分析 , PROFILE 変換 時点間差対比作成。 例: REPEATED time 3 PROFILE 指定場合、 対比行列以下作成。 t1 t2 t3 treatmnt_1 1.000000000 -1.000000000 0.000000000 treatmnt_2 0.000000000 1.000000000 -1.000000000,