Microsoft PowerPoint - SAS2012_ZHANG_0629.ppt [互換モード]

Size: px
Start display at page:

Download "Microsoft PowerPoint - SAS2012_ZHANG_0629.ppt [互換モード]"

Transcription

1 SAS による生存時間解析の実務 張方紅グラクソ スミスクライン ( 株 バイオメディカルデータサイエンス部 Practice of Survival Analysis sing SAS Fanghong Zhang Biomedical Data Science Department, GlaxoSmithKline K.K. 要旨 : SASによる生存時間解析の実務経験を共有する. データの要約 検定 推定の前提に注意し, 群間差の検定とハザード比の信頼区間に焦点を与える. 複数に存在するプログラミング方法の間の関係と使い勝手を紹介する. キワド生存時間解析 LIFETEST プロシジャロキーワード : 生存時間解析,LIFETEST プロシジャ, ログランク検定, ハザード比

2 動機づけ : 実務経験を共有 生存時間解析に報告が多い 概念 本質 浜田 (20, 生存時間解析再入門 生存時間解析のミステリーをひも解く 学術 平井 吉田 田崎 (2009,22 つの生存時間の差に対する同時信頼区間 実務レベルの経験 仮想な臨床試験 SASプログラミング 感想 まとめ 発表内容

3 仮想な臨床試験 主要目的 : 治験薬群であるT+ 併用療法と対照薬群である 単独療法とで全生存期間 (OS Overall Survival を比較する 主要評価項目 : 全生存期間 解析方法. 累積生存曲線, 生存期間中央値の推定は Kaplan-Meier 法を用いて行い Greenwood の公式を用いて95% 信頼区間を求める 2. 群間比較には施設以外の割付調整因子を層とした層別ログランク検定を用いる 3. 治療効果の推定値として ox の比例ハザードモデルを用いて群間の治療効果のハザード比とその95% 信頼区間を求める JOG プロトコールマニュアル version 2.0

4 ASO (American Society of linical Oncology R T+ Overall survival T+ N=77 N=75 Median (95% I 2.8 (9.6, (6.8, 0. ox HR (95% I 0.57 (0.40, 0.83 Two-sided p-value (Log-rank test T+

5 レポート : 結果 主要評価項目であるOSの中央値はT+ 併用群で 2.8 ヵ月 単独群で8 8. ヵ月 (HR=0.57, 95%I: , p= と有意な延長を示した. HR: Hazard Ratio I: onfidence Interval 解析段階と前提 解析段階内容前提 分布の記述要約 検定 推定 生存関数を Kaplan-Meier 法で推定し, プロットする. 推定された生存曲線から中央値を求める ノンパラメトリック法であり, 前提が必要なし ノンパラメトリック法, 群間で生存関数の違いが生存関数が群間であるかを検定する交差していない ハザード比を推定する セミパラメトリック法, 比例ハザード性

6 ログランク検定の帰無仮説 両群のOSの中央値が等しいという仮説に対する検定 両群のハザード比 =という仮説に対する検定 両群の生存関数が等しいという仮説に対する検定 中央値 ハザード比 P 値 T+ N=58 N=57 Median (95% I 5.6 ( (2.8, 5.7 HR (95% I 0.6 (0.37,.00 (ox Two-sided p-value (Log-rank test I of HR と P 値の結果と一致しない 前提が異なる 手法が対応していない :Wald vs Score

7 ここまでの結論 ログランク検定は両群の生存関数が等しいという仮説に対する検定 中央値 ハザード比 ログランク検定はそれぞれ前提が異なる, 方法が対応していない場合が多い ログランク検定 : 概念 記号 ログランク (log-rank 検定は, 生存曲線の差を ある時点ではなく観察期間全体を通じて比較し検定できる方法です. O E V 2 O-E: スコア統計量 V: O-E の分散 O: 治験薬群における観察死亡数, E: 治験薬群における期待死亡数 2

8 ログランク検定 :SAS プロシジャ 2 3 PRO LIFETEST; STRATA 治療群変数 ; PRO LIFETEST; TEST 治療群変数 ; PRO PHREG; MODEL 時間変数 * 打ち切り変数 = 治療群変数 /TIES=; 3つの方法で本質的に同等な結果違いはタイの取り扱い 大橋 浜田 (995, p.54 ログランク検定 :SAS プロシジャ PRO LIFETEST; STRATA 治療群変数 ; PRO LIFETEST; TEST 治療群変数 ; PRO PHREG; MODEL 時間変数 * 打ち切り変数 = 治療群変数 / TIES=DISRETE; PRO PHREG; MODEL 時間変数 * 打ち切り変数 = 治療群変数 / TIES=BRESLOW;/* デフォルト */ 上の表の左 右の方法はそれぞれ結果が一致ただし,PHREGプロシジャのScore 検定を使用 ollett (994, p.332

9 2 3 ログランク検定 :STRATA 文を選ぶ PRO LIFETEST; STRATA 治療群変数 ; PRO LIFETEST; TEST 治療群変数 ; PRO PHREG; MODEL 時間変数 * 打ち切り変数 = 治療群変数 /TIES=; STRATA 文を選ぶ理由. 生存関数を計算 2. スコア統計量 (O-E を出力 層別ログランク検定 :SAS プロシジャ 2 3 PRO LIFETEST; STRATA 層別変数 /GROUP= 治療群変数 ; PRO LIFETEST; STRATA 層別変数 ; TEST 治療群変数 ; PRO PHREG; MODEL 時間変数 * 打ち切り変数 = 治療群変数 /TIES=; STRATA 層別変数 ; 3 つの方法で本質的に同等な結果 SAS 9. マニュアル

10 層別ログランク検定 :GROUP オプション PRO LIFETEST; STRATA 層別変数 /GROUP= 治療群変数 ; 2 PRO LIFETEST; STRATA 治療群変数 ; BY 層別変数 ; 2つの方法で同じ結果 大橋 浜田 (995, p99-00 SAS 9.GROUP オプションが便利になった ログランク検定 : 片側. ログランク検定は両側しかない 2. 解決方法 : ログランク検定に対応する漸近正規分布に従う Z 統計量を利用 Z O-E: スコア統計量 V: O-Eの分散 O E V

11 ログランク検定 : 片側 ods output HomStats=stats /* O-E スコア統計量 */ LogHomov=cov; /*(V: O-E の分散 */ proc lifetest; strata 治療群変数 ; run; ログランク検定 : まとめ 両側 片側 層なし PRO LIFETEST; STRATA 治療群変数 ; PRO LIFETEST; 層あり STRATA 層別変数 /GROUP= 治療群変数 ; スコアO-Eとその分散 Vをデータセットに落とし,Zを計算 ods output HomStats=stats LogHomov=cov; PRO LIFETEST; 層なし STRATA 治療群変数 ; PRO LIFETEST; 層あり STRATA 層別変数 /GROUP= 治療群変数 ;

12 . ox 法 ハザード比 2. Peto 法 3. Pike 法 4. Mantel-Haenszel 法 ハザード比 : 比例ハザード性を前提 ( t exp( ( t T λ T : 治験薬群におけるハザード λ : 対照薬群におけるハザード 比例ハザード性が成り立っていれば, 群間のリスクの違いを つのパラメータに要約することが可能となる.

13 ox 法 : 部分尤度 r x L ( exp( ( ollett p.07, (3.28, x ( = 治験薬群, 0 対照薬群 O: 治験薬群における総死亡数,n と n 2 はリスク集合の大きさ n l l x L exp( ( r r n l l r n n O x x L 2 ( exp( log exp( log ( log ox 法 2 2 exp( exp( ( log ( H E O n n n O L u r E H : 治験薬群における対立仮説の下での期待死亡数

14 ˆ ML u( ˆ I of ˆ ox 法 :SAS : 最尤推定値 ML 0, O E H ML.96SE( ˆ ML I of HR exp( ˆ.96 SE ( ˆ ML PRO PHREG; MODEL 時間変数 * 打ち切り変数 = 治療群変数 /RL; ML Peto 法 β :β の最尤推定の one step 近似 O E V O-E: スコア統計量 V: O-Eの分散 O: 治験薬群における観察死亡数, E: 治験薬群における期待死亡数 大橋 浜田 (995, p.9-95

15 Peto 法 :SAS O E V SE( V I of O-E HR exp V.96 V ods output HomStats=stats LogHomov=cov; PRO LIFETEST; Pike 法 ~ O log O T / E / E T ~ SE( E T E I of OT / E HR exp log O / E T.96 E T E

16 Pike 法 :SAS I of O T / E T HR exp log.96 O / E ET E 計算対象 SAS プログラミング ods output censoredsummary=events; O PRO LIFETEST; ods output ods output HomStats=stats HomStats=stats; LogHomov=cov; O E PRO PRO LIFETEST; LIFETEST; E E=O (O E 感想 ログランク検定 ( 層別あり / なし :SAS 9. 以降は直接に出力でき, 便利になった. ハザード比 :ox 法しか直接に出力できなく, 不便のまま ハザード比に対して, 方法論の研究と比較も不十分

17 まとめ 中央値 ハザード比 ログランク検定のP 値はそれぞれ前提が異なる, 対応していない LIFETESTプロシジャSTRATA 文は使い勝手がよい ハザード比に対する方法論の比較研究が必要 参考文献 (/2 大橋靖雄 浜田知久馬 (995, 生存時間解析 SAS による生物統計. 東京大学出版会. 浜田知久馬 (20, 生存時間解析再入門 生存時間解析のミステリーをひも解く. 日本 SAS ユーザー会論文集, 平井健太 吉田祐樹 田崎武信 (2009, 2 つの生存時間の差に対する同時信頼区間. 日本 SAS ユーザー会論文集, JOG (2008 JOG プロトコールマニュアル version (202.6 アクセス ASO 202 演題速報 (202.6 アクセス

18 参考文献 (2/2 Berry G., Kitchin R. M. and Mock P. A. (99: A comparison of two simple hazard ratio estimators based on the logrank test. Statistics in Medicine, Vol. 0, ollett D. (994, Modelling Survival Data in Medical Research. hapman & Hall: London. Sato T.(992 Letters to the editor, Statistics in Medicine, i p SAS/STAT 9. User s Guide The LIFETEST Procedure

Microsoft PowerPoint - 【配布・WEB公開用】SAS発表資料.pptx

Microsoft PowerPoint - 【配布・WEB公開用】SAS発表資料.pptx 生存関数における信頼区間算出法の比較 佐藤聖士, 浜田知久馬東京理科大学工学研究科 Comparison of confidence intervals for survival rate Masashi Sato, Chikuma Hamada Graduate school of Engineering, Tokyo University of Science 要旨 : 生存割合の信頼区間算出の際に用いられる各変換関数の性能について被覆確率を評価指標として比較した.

More information

日本製薬工業協会シンポジウム 生存時間解析の評価指標に関する最近の展開ー RMST (restricted mean survival time) を理解するー 2. RMST の定義と統計的推測 2018 年 6 月 13 日医薬品評価委員会データサイエンス部会タスクフォース 4 生存時間解析チー

日本製薬工業協会シンポジウム 生存時間解析の評価指標に関する最近の展開ー RMST (restricted mean survival time) を理解するー 2. RMST の定義と統計的推測 2018 年 6 月 13 日医薬品評価委員会データサイエンス部会タスクフォース 4 生存時間解析チー 日本製薬工業協会シンポジウム 生存時間解析の評価指標に関する最近の展開ー RMST (restricted mean survival time) を理解するー 2. RMST の定義と統計的推測 2018 年 6 月 13 日医薬品評価委員会データサイエンス部会タスクフォース 4 生存時間解析チーム 日本新薬 ( 株 ) 田中慎一 留意点 本発表は, 先日公開された 生存時間型応答の評価指標 -RMST(restricted

More information

NLMIXED プロシジャを用いた生存時間解析 伊藤要二アストラゼネカ株式会社臨床統計 プログラミング グループグルプ Survival analysis using PROC NLMIXED Yohji Itoh Clinical Statistics & Programming Group, A

NLMIXED プロシジャを用いた生存時間解析 伊藤要二アストラゼネカ株式会社臨床統計 プログラミング グループグルプ Survival analysis using PROC NLMIXED Yohji Itoh Clinical Statistics & Programming Group, A NLMIXED プロシジャを用いた生存時間解析 伊藤要二アストラゼネカ株式会社臨床統計 プログラミング グループグルプ Survival analysis using PROC NLMIXED Yohji Itoh Clinical Statistics & Programming Group, AstraZeneca KK 要旨 : NLMIXEDプロシジャの最尤推定の機能を用いて 指数分布 Weibull

More information

スライド 1

スライド 1 生存時間解析における Lakatos の症例数設計法の有用性の評価 魚住龍史, * 水澤純基 浜田知久馬 日本化薬株式会社医薬データセンター 東京理科大学工学部経営工学科 Evaluation of availability about sample size formula by Lakatos on survival analysis Ryuji Uozumi,, * Junki Mizusawa,

More information

SAS_2014_zhang_3

SAS_2014_zhang_3 生物学的同等性試験における例数設計 : 正確 近似と漸近 張方紅 安藤英一グラクソ スミスクライン 株 バイオメディカルデータサイエンス部 Sample size for bioequivalence rials: Exac approximae and asympoic mehods Fanghong Zhang and Hidekazu Ando Biomedical Daa Science Deparmen

More information

Chapter 1 Epidemiological Terminology

Chapter 1 Epidemiological Terminology Appendix Real examples of statistical analysis 検定 偶然を超えた差なら有意差という P

More information

PHREG プロシジャにおける 共変量調整解析に関連したオプション機能 魚住龍史 1 * 矢田真城 2 浜田知久馬 3 1 京都大学大学院医学研究科医学統計生物情報学 2 エイツーヘルスケア株式会社 3 東京理科大学 Investigating fascinating aspects associa

PHREG プロシジャにおける 共変量調整解析に関連したオプション機能 魚住龍史 1 * 矢田真城 2 浜田知久馬 3 1 京都大学大学院医学研究科医学統計生物情報学 2 エイツーヘルスケア株式会社 3 東京理科大学 Investigating fascinating aspects associa PHREG プロシジャにおける 共変量調整解析に関連したオプション機能 魚住龍史 1 * 矢田真城 2 浜田知久馬 3 1 京都大学大学院医学研究科医学統計生物情報学 2 エイツーヘルスケア株式会社 3 東京理科大学 Investigating fascinating aspects associated with covariate-adjusted analysis using PHREG procedure

More information

Microsoft PowerPoint - R-stat-intro_12.ppt [互換モード]

Microsoft PowerPoint - R-stat-intro_12.ppt [互換モード] R で統計解析入門 (12) 生存時間解析 中篇 準備 : データ DEP の読み込み 1. データ DEP を以下からダウンロードする http://www.cwk.zaq.ne.jp/fkhud708/files/dep.csv /fkh /d 2. ダウンロードした場所を把握する ここでは c:/temp とする 3. R を起動し,2. 2 の場所に移動し, データを読み込む 4. データ

More information

<4D F736F F F696E74202D204D C982E682E892B290AE82B582BD838A E8DB782CC904D978A8BE68AD482C98AD682B782E988EA8D6C8E402E >

<4D F736F F F696E74202D204D C982E682E892B290AE82B582BD838A E8DB782CC904D978A8BE68AD482C98AD682B782E988EA8D6C8E402E > SAS ユーザー総会 2017 Mantel-Haenszel 法により調整したリスク差の信頼区間に関する一考察 武田薬品工業株式会社日本開発センター生物統計室佐々木英麿 舟尾暢男 要旨 Mantel-Haenszel 法により調整したリスク差に関する以下の信頼区間の算出方法を紹介し 各信頼区間の被覆確率をシミュレーションにより確認することで性能評価を行う Greenland 信頼区間 Sato 信頼区間

More information

Kaplan-Meierプロットに付加情報を追加するマクロの作成

Kaplan-Meierプロットに付加情報を追加するマクロの作成 Kaplan-Meier 1, 2,3 1 2 3 A SAS macro for extended Kaplan-Meier plots Kengo Nagashima 1, Yasunori Sato 2,3 1 Department of Parmaceutical Technochemistry, Josai University 2 School of Medicine, Chiba University

More information

統計的データ解析

統計的データ解析 統計的データ解析 011 011.11.9 林田清 ( 大阪大学大学院理学研究科 ) 連続確率分布の平均値 分散 比較のため P(c ) c 分布 自由度 の ( カイ c 平均値 0, 標準偏差 1の正規分布 に従う変数 xの自乗和 c x =1 が従う分布を自由度 の分布と呼ぶ 一般に自由度の分布は f /1 c / / ( c ) {( c ) e }/ ( / ) 期待値 二乗 ) 分布 c

More information

基礎統計

基礎統計 基礎統計 第 11 回講義資料 6.4.2 標本平均の差の標本分布 母平均の差 標本平均の差をみれば良い ただし, 母分散に依存するため場合分けをする 1 2 3 分散が既知分散が未知であるが等しい分散が未知であり等しいとは限らない 1 母分散が既知のとき が既知 標準化変量 2 母分散が未知であり, 等しいとき 分散が未知であるが, 等しいということは分かっているとき 標準化変量 自由度 の t

More information

Medical3

Medical3 Chapter 1 1.4.1 1 元配置分散分析と多重比較の実行 3つの治療法による測定値に有意な差が認められるかどうかを分散分析で調べます この例では 因子が1つだけ含まれるため1 元配置分散分析 one-way ANOVA の適用になります また 多重比較法 multiple comparison procedure を用いて 具体的のどの治療法の間に有意差が認められるかを検定します 1. 分析メニュー

More information

2 H23 BioS (i) data d1; input group patno t sex censor; cards;

2 H23 BioS (i) data d1; input group patno t sex censor; cards; H BioS (i) data d1; input group patno t sex censor; cards; 0 1 0 0 0 0 1 0 1 1 0 4 4 0 1 0 5 5 1 1 0 6 5 1 1 0 7 10 1 0 0 8 15 0 1 0 9 15 0 1 0 10 4 1 0 0 11 4 1 0 1 1 5 1 0 1 1 7 0 1 1 14 8 1 0 1 15 8

More information

JMP による 2 群間の比較 SAS Institute Japan 株式会社 JMP ジャパン事業部 2008 年 3 月 JMP で t 検定や Wilcoxon 検定はどのメニューで実行できるのか または検定を行う際の前提条件の評価 ( 正規性 等分散性 ) はどのメニューで実行できるのかと

JMP による 2 群間の比較 SAS Institute Japan 株式会社 JMP ジャパン事業部 2008 年 3 月 JMP で t 検定や Wilcoxon 検定はどのメニューで実行できるのか または検定を行う際の前提条件の評価 ( 正規性 等分散性 ) はどのメニューで実行できるのかと JMP による 2 群間の比較 SAS Institute Japan 株式会社 JMP ジャパン事業部 2008 年 3 月 JMP で t 検定や Wilcoxon 検定はどのメニューで実行できるのか または検定を行う際の前提条件の評価 ( 正規性 等分散性 ) はどのメニューで実行できるのかというお問い合わせがよくあります そこで本文書では これらについて の回答を 例題を用いて説明します 1.

More information

RSS Higher Certificate in Statistics, Specimen A Module 3: Basic Statistical Methods Solutions Question 1 (i) 帰無仮説 : 200C と 250C において鉄鋼の破壊応力の母平均には違いはな

RSS Higher Certificate in Statistics, Specimen A Module 3: Basic Statistical Methods Solutions Question 1 (i) 帰無仮説 : 200C と 250C において鉄鋼の破壊応力の母平均には違いはな RSS Higher Certiicate in Statistics, Specimen A Module 3: Basic Statistical Methods Solutions Question (i) 帰無仮説 : 00C と 50C において鉄鋼の破壊応力の母平均には違いはない. 対立仮説 : 破壊応力の母平均には違いがあり, 50C の方ときの方が大きい. n 8, n 7, x 59.6,

More information

PowerPoint プレゼンテーション

PowerPoint プレゼンテーション 27//26 第 4 回 医学統計勉強会 東北大学病院循環器内科 東北大学病院臨床研究推進センター 共催 東北大学大学院医学系研究科 EBM 開発学寄附講座 宮田 敏 生存時間解析生存曲線,Cox 比例ハザードモデル 生存時間解析 (survival time analysis) では, 基準となるある時点から, 目的となるイベントの発生までの時間を解析する. 例えば, ある疾患の登録研究において,

More information

スライド 1

スライド 1 SAS による二項比率の差の非劣性検定の正確な方法について 武藤彬正宮島育哉榊原伊織株式会社タクミインフォメーションテクノロジー Eact method of non-inferiority test for two binomial proportions using SAS Akimasa Muto Ikuya Miyajima Iori Sakakibara Takumi Information

More information

分析のステップ Step 1: Y( 目的変数 ) に対する値の順序を確認 Step 2: モデルのあてはめ を実行 適切なモデルの指定 Step 3: オプションを指定し オッズ比とその信頼区間を表示 以下 このステップに沿って JMP の操作をご説明します Step 1: Y( 目的変数 ) の

分析のステップ Step 1: Y( 目的変数 ) に対する値の順序を確認 Step 2: モデルのあてはめ を実行 適切なモデルの指定 Step 3: オプションを指定し オッズ比とその信頼区間を表示 以下 このステップに沿って JMP の操作をご説明します Step 1: Y( 目的変数 ) の JMP によるオッズ比 リスク比 ( ハザード比 ) の算出と注意点 SAS Institute Japan 株式会社 JMP ジャパン事業部 2011 年 10 月改定 1. はじめに 本文書は JMP でロジスティック回帰モデルによるオッズ比 比例ハザードモデルによるリスク比 それぞれに対する信頼区間を求める操作方法と注意点を述べたものです 本文書は JMP 7 以降のバージョンに対応しております

More information

青焼 1章[15-52].indd

青焼 1章[15-52].indd 1 第 1 章統計の基礎知識 1 1 なぜ統計解析が必要なのか? 人間は自分自身の経験にもとづいて 感覚的にものごとを判断しがちである 例えばある疾患に対する標準治療薬の有効率が 50% であったとする そこに新薬が登場し ある医師がその新薬を 5 人の患者に使ったところ 4 人が有効と判定されたとしたら 多くの医師はこれまでの標準治療薬よりも新薬のほうが有効性が高そうだと感じることだろう しかし

More information

Kumamoto University Center for Multimedia and Information Technologies Lab. 熊本大学アプリケーション実験 ~ 実環境における無線 LAN 受信電波強度を用いた位置推定手法の検討 ~ InKIAI 宮崎県美郷

Kumamoto University Center for Multimedia and Information Technologies Lab. 熊本大学アプリケーション実験 ~ 実環境における無線 LAN 受信電波強度を用いた位置推定手法の検討 ~ InKIAI 宮崎県美郷 熊本大学アプリケーション実験 ~ 実環境における無線 LAN 受信電波強度を用いた位置推定手法の検討 ~ InKIAI プロジェクト @ 宮崎県美郷町 熊本大学副島慶人川村諒 1 実験の目的 従来 信号の受信電波強度 (RSSI:RecevedSgnal StrengthIndcator) により 対象の位置を推定する手法として 無線 LAN の AP(AccessPont) から受信する信号の減衰量をもとに位置を推定する手法が多く検討されている

More information

講義「○○○○」

講義「○○○○」 講義 信頼度の推定と立証 内容. 点推定と区間推定. 指数分布の点推定 区間推定 3. 指数分布 正規分布の信頼度推定 担当 : 倉敷哲生 ( ビジネスエンジニアリング専攻 ) 統計的推測 標本から得られる情報を基に 母集団に関する結論の導出が目的 測定値 x x x 3 : x 母集団 (populaio) 母集団の特性値 統計的推測 標本 (sample) 標本の特性値 分布のパラメータ ( 母数

More information

JMP によるオッズ比 リスク比 ( ハザード比 ) の算出方法と注意点 SAS Institute Japan 株式会社 JMP ジャパン事業部 2008 年 3 月改定 1. はじめに本文書は JMP でオッズ比 リスク比 それぞれに対する信頼区間を求める算出方法と注意点を述べたものです この後

JMP によるオッズ比 リスク比 ( ハザード比 ) の算出方法と注意点 SAS Institute Japan 株式会社 JMP ジャパン事業部 2008 年 3 月改定 1. はじめに本文書は JMP でオッズ比 リスク比 それぞれに対する信頼区間を求める算出方法と注意点を述べたものです この後 JMP によるオッズ比 リスク比 ( ハザード比 ) の算出方法と注意点 SAS Institute Japan 株式会社 JMP ジャパン事業部 2008 年 3 月改定 1. はじめに本文書は JMP でオッズ比 リスク比 それぞれに対する信頼区間を求める算出方法と注意点を述べたものです この後の 2 章では JMP でのオッズ比 オッズ比の信頼区間の算出方法について サンプルデータを用いて解説しております

More information

<4D F736F F F696E74202D F95618A7789EF B836A F838C834E B88E38A77939D8C76322E >

<4D F736F F F696E74202D F95618A7789EF B836A F838C834E B88E38A77939D8C76322E > 204 年 9 月 26 日第 62 回日本心臓病学会学術集会モーニングレクチャー 医学統計の基礎 於 : 仙台国際センター第 9 会場 医学統計の基礎 東北大学大学院医学系研究科循環器内科学分野 宮田敏 miyata@cardio.med.tohou.ac.jp 日本心臓病学会 COI 開示 東北大学大学院医学系研究科循環器内科学宮田敏 演題発表に関連し 開示すべき CO I 関係にある企業などはありません

More information

<4D F736F F D2091E63489F190B691B68E9E8AD489F090CD2E646F6378>

<4D F736F F D2091E63489F190B691B68E9E8AD489F090CD2E646F6378> 医学統計勉強会 東北大学病院循環器内科 東北大学臨床研究推進センター共催 東北大学大学院医学系研究科 EBM 開発学寄附講座 宮田敏 Daa! daa! daa! he cried impaienly. I can' mae brics wihou clay. From The Advenure of he Copper Beeches, The Advenure of Sherloc Holmes.

More information

Statistical inference for one-sample proportion

Statistical inference for one-sample proportion RAND 関数による擬似乱数の生成 魚住龍史 * 浜田知久馬東京理科大学大学院工学研究科経営工学専攻 Generating pseudo-random numbers using RAND function Ryuji Uozumi * and Chikuma Hamada Department of Management Science, Graduate School of Engineering,

More information

連載 第 5 回 医学データの統計解析の基本 生存時間データの解析 Fundamentals of statistical analysis in biomedical research:analysis of survival data 朝倉こう子 * * 濱﨑俊光 1 生物学や医学分野において

連載 第 5 回 医学データの統計解析の基本 生存時間データの解析 Fundamentals of statistical analysis in biomedical research:analysis of survival data 朝倉こう子 * * 濱﨑俊光 1 生物学や医学分野において 連載 第 5 回 医学データの統計解析の基本 生存時間データの解析 Fundamentals of statistical analysis in biomedical research:analysis of survival data 朝倉こう子 * * 濱﨑俊光 1 生物学や医学分野において 個体 ( 例えば 培養細胞 動物あるいは被験者 患者 ) に対して化合物や医薬品といった処理を与えることにより惹き起こされる反応は

More information

Microsoft PowerPoint - SASユーザ総会2016_MRCT_送付用.pptx

Microsoft PowerPoint - SASユーザ総会2016_MRCT_送付用.pptx Multi Regional Clinical Trial の生存時間解析における地域ごとのイベント発現数の予測 淀康秀 1 浜田知久馬 2 ( 1 大日本住友製薬株式会社 2 東京理科大学 ) The prediction of the number of survival event occurrence by region in Multi-Regional Clinical Trial Yasuhide

More information

解析センターを知っていただく キャンペーン

解析センターを知っていただく キャンペーン 005..5 SAS 問題設定 目的 PKパラメータ (AUC,Cmax,Tmaxなど) の推定 PKパラメータの群間比較 PKパラメータのバラツキの評価! データの特徴 非反復測定値 個体につき 個の測定値しか得られない plasma concentration 非反復測定値のイメージ図 測定時点間で個体の対応がない 着目する状況 plasma concentration 経時反復測定値のイメージ図

More information

Microsoft PowerPoint - e-stat(OLS).pptx

Microsoft PowerPoint - e-stat(OLS).pptx 経済統計学 ( 補足 ) 最小二乗法について 担当 : 小塚匡文 2015 年 11 月 19 日 ( 改訂版 ) 神戸大学経済学部 2015 年度後期開講授業 補足 : 最小二乗法 ( 単回帰分析 ) 1.( 単純 ) 回帰分析とは? 標本サイズTの2 変数 ( ここではXとY) のデータが存在 YをXで説明する回帰方程式を推定するための方法 Y: 被説明変数 ( または従属変数 ) X: 説明変数

More information

Microsoft PowerPoint - sc7.ppt [互換モード]

Microsoft PowerPoint - sc7.ppt [互換モード] / 社会調査論 本章の概要 本章では クロス集計表を用いた独立性の検定を中心に方法を学ぶ 1) 立命館大学経済学部 寺脇 拓 2 11 1.1 比率の推定 ベルヌーイ分布 (Bernoulli distribution) 浄水器の所有率を推定したいとする 浄水器の所有の有無を表す変数をxで表し 浄水器をもっている を 1 浄水器をもっていない を 0 で表す 母集団の浄水器を持っている人の割合をpで表すとすると

More information

スライド 1

スライド 1 SASによる二項比率における正確な信頼区間の比較 原茂恵美子 1) 武藤彬正 1) 宮島育哉 2) 榊原伊織 2) 1) 株式会社タクミインフォメーションテクノロジーシステム開発推進部 2) 株式会社タクミインフォメーションテクノロジービジネスソリューション部 Comparison of Five Exact Confidence Intervals for the Binomial Proportion

More information

Microsoft PowerPoint - H17-5時限(パターン認識).ppt

Microsoft PowerPoint - H17-5時限(パターン認識).ppt パターン認識早稲田大学講義 平成 7 年度 独 産業技術総合研究所栗田多喜夫 赤穂昭太郎 統計的特徴抽出 パターン認識過程 特徴抽出 認識対象から何らかの特徴量を計測 抽出 する必要がある 認識に有効な情報 特徴 を抽出し 次元を縮小した効率の良い空間を構成する過程 文字認識 : スキャナ等で取り込んだ画像から文字の識別に必要な本質的な特徴のみを抽出 例 文字線の傾き 曲率 面積など 識別 与えられた未知の対象を

More information

(3) 検定統計量の有意確率にもとづく仮説の採否データから有意確率 (significant probability, p 値 ) を求め 有意水準と照合する 有意確率とは データの分析によって得られた統計値が偶然おこる確率のこと あらかじめ設定した有意確率より低い場合は 帰無仮説を棄却して対立仮説

(3) 検定統計量の有意確率にもとづく仮説の採否データから有意確率 (significant probability, p 値 ) を求め 有意水準と照合する 有意確率とは データの分析によって得られた統計値が偶然おこる確率のこと あらかじめ設定した有意確率より低い場合は 帰無仮説を棄却して対立仮説 第 3 章 t 検定 (pp. 33-42) 3-1 統計的検定 統計的検定とは 設定した仮説を検証する場合に 仮説に基づいて集めた標本を 確率論の観点から分析 検証すること 使用する標本は 母集団から無作為抽出されたものでなければならない パラメトリック検定とノンパラメトリック検定 パラメトリック検定は母集団が正規分布に従う間隔尺度あるいは比率尺度の連続データを対象とする ノンパラメトリック検定は母集団に特定の分布を仮定しない

More information

Microsoft PowerPoint - R-survival.ppt

Microsoft PowerPoint - R-survival.ppt 統計解析フリーソフト R 入門 R による生存時間解析 本日のメニュー R のインストール R による生存時間解析 イントロ 生存関数の推定と群間比較 競合リスクについて その他 2 R のインストール 実行ファイル R-2.6.0pat-win32.exe をダブルクリック http://cran.md.tsukuba.ac.jp/bin/windows/base/r-2.6.0pat-win32.exe

More information

骨髄移植データに関する イベントヒストリー解析 大阪電気通信大学大学院情報工学専攻 辻谷研究室 中井崇人 1

骨髄移植データに関する イベントヒストリー解析 大阪電気通信大学大学院情報工学専攻 辻谷研究室 中井崇人 1 骨髄移植データに関する イベントヒストリー解析 大阪電気通信大学大学院情報工学専攻 辻谷研究室 中井崇人 1 目次 1. はじめに 2. Multi-stateモデル 3. 平滑化スプライン ( 一般化加法モデル ) 4. 生存率の予測 2 1. はじめに 白血病は骨髄のガン化が原因 他人の骨髄を体内にいれるのでさまざまな問題がある 患者が知りたいはどのくらい生きれるか 1 年後の生存確率 3 2.

More information

<4D F736F F D208EC08CB18C7689E68A E F AA957A82C682948C9F92E82E646F63>

<4D F736F F D208EC08CB18C7689E68A E F AA957A82C682948C9F92E82E646F63> 第 7 回 t 分布と t 検定 実験計画学 A.t 分布 ( 小標本に関する平均の推定と検定 ) 前々回と前回の授業では, 標本が十分に大きいあるいは母分散が既知であることを条件に正規分布を用いて推定 検定した. しかし, 母集団が正規分布し, 標本が小さい場合には, 標本分散から母分散を推定するときの不確実さを加味したt 分布を用いて推定 検定しなければならない. t 分布は標本分散の自由度 f(

More information

EBNと疫学

EBNと疫学 推定と検定 57 ( 復習 ) 記述統計と推測統計 統計解析は大きく 2 つに分けられる 記述統計 推測統計 記述統計 観察集団の特性を示すもの 代表値 ( 平均値や中央値 ) や ばらつきの指標 ( 標準偏差など ) 図表を効果的に使う 推測統計 観察集団のデータから母集団の特性を 推定 する 平均 / 分散 / 係数値などの推定 ( 点推定 ) 点推定値のばらつきを調べる ( 区間推定 ) 検定統計量を用いた検定

More information

Microsoft PowerPoint - 【魚住】 発表原稿【Web公開用】.ppt [互換モード]

Microsoft PowerPoint - 【魚住】 発表原稿【Web公開用】.ppt [互換モード] SG (Statistical Graphics) Procedures による Kaplan-Meier プロットの作成 魚住龍史 1, * 浜田知久馬 2 1 日本化薬株式会社医薬データセンター 2 東京理科大学工学部経営工学科 Kaplan-Meier plots using Statistical Graphics Procedures Ryuji Uozumi 1, * and Chikuma

More information

生存時間データに対するベイズ流例数設計 矢田真城 1 魚住龍史 2 浜田知久馬 1 エイツーヘルスケア株式会社開発戦略本部生物統計部 2 京都大学大学院医学研究科医学統計生物情報学 3 東京理科大学工学部情報工学科 3 Bayesian sample size calculation for sur

生存時間データに対するベイズ流例数設計 矢田真城 1 魚住龍史 2 浜田知久馬 1 エイツーヘルスケア株式会社開発戦略本部生物統計部 2 京都大学大学院医学研究科医学統計生物情報学 3 東京理科大学工学部情報工学科 3 Bayesian sample size calculation for sur 生存時間データに対するベイズ流例数設計 矢田真城 魚住龍史 浜田知久馬 エイツーヘルスケア株式会社開発戦略本部生物統計部 京都大学大学院医学研究科医学統計生物情報学 3 東京理科大学工学部情報工学科 3 Bayesian sample size calculation for survival analysis Shinjo Yada, Ryuji Uozumi, and Chikuma Hamada

More information

<4D F736F F D204B208C5182CC94E497A682CC8DB782CC8C9F92E BD8F6494E48A722E646F6378>

<4D F736F F D204B208C5182CC94E497A682CC8DB782CC8C9F92E BD8F6494E48A722E646F6378> 3 群以上の比率の差の多重検定法 013 年 1 月 15 日 017 年 3 月 14 日修正 3 群以上の比率の差の多重検定法 ( 対比較 ) 分割表で表記される計数データについて群間で比率の差の検定を行う場合 全体としての統計的有意性の有無は χ 検定により判断することができるが 個々の群間の差の有意性を判定するためには多重検定法が必要となる 3 群以上の比率の差を対比較で検定する方法としては

More information

Microsoft Word - Stattext12.doc

Microsoft Word - Stattext12.doc 章対応のない 群間の量的データの検定. 検定手順 この章ではデータ間に 対 の対応のないつの標本から推定される母集団間の平均値や中央値の比較を行ないます 検定手法は 図. のようにまず正規に従うかどうかを調べます 但し この場合はつの群が共に正規に従うことを調べる必要があります 次に 群とも正規ならば F 検定を用いて等分散であるかどうかを調べます 等分散の場合は t 検定 等分散でない場合はウェルチ

More information

MedicalStatisticsForAll.indd

MedicalStatisticsForAll.indd みんなの 医療統計 12 基礎理論と EZR を完全マスター! Ayumi SHINTANI はじめに EZR EZR iii EZR 2016 2 iv CONTENTS はじめに... ⅲ EZR をインストールしよう... 1 EZR 1...1 EZR 2...3...8 R Console...10 1 日目 記述統計量...11 平均値と中央値... 11...12...15...18

More information

スライド 1

スライド 1 SAS による生存時間分布の予測 Death Note の統計学 東京理科大学 浜田知久馬, 魚住龍史 Death Note Predcton of Survval Dstrbuton usng SAS Statstcs for Death Note Chkuma Hamada Toko Unverst of Scence 要旨 分位点回帰を行うプロシジャ QUANTREG プロシジャを打ち切りがある場合に拡張した

More information

Medical3

Medical3 1.4.1 クロス集計表の作成 -l m 分割表 - 3つ以上のカテゴリを含む変数を用いて l mのクロス集計表による分析を行います この例では race( 人種 ) によってlow( 低体重出生 ) に差が認められるかどうかを分析します 人種には3つのカテゴリ 低体重出生には2つのカテゴリが含まれています 2つの変数はともにカテゴリ変数であるため クロス集計表によって分析します 1. 分析メニュー

More information

Presentation Title Goes Here

Presentation  Title Goes Here SAS 9: (reprise) SAS Institute Japan Copyright 2004, SAS Institute Inc. All rights reserved. Greetings, SAS 9 SAS 9.1.3 Copyright 2004, SAS Institute Inc. All rights reserved. 2 Informations of SAS 9 SAS

More information

目次 はじめに P.01 適用分野

目次 はじめに P.01 適用分野 Cox 比例ハザードモデル 米銀行の倒産確率の推定 2016/11/1 目次 はじめに --------------------------------------------------------------------------------------------------------------------------- P.01 適用分野 ---------------------------------------------------------------------------------------------------------------------------

More information

<4D F736F F F696E74202D2088E38A77939D8C7695D78BAD89EF313791E63589F194E497A682C695AA8A84955C2E >

<4D F736F F F696E74202D2088E38A77939D8C7695D78BAD89EF313791E63589F194E497A682C695AA8A84955C2E > 27//9 第 5 回 医学統計勉強会 東北大学病院循環器内科 東北大学病院臨床研究推進センター 共催 東北大学大学院医学系研究科 EBM 開発学寄附講座 宮田 敏 比率と分割表 疾患の発症率など, 物事の頻度 (frequency) を議論する際, 以下の三つの概念を使い分ける. 比 (ratio):a, B ( ) が存在するとき,A/B を比という. A と B は互いを含まない. 例 : 性比.BMI=

More information

カイ二乗フィット検定、パラメータの誤差

カイ二乗フィット検定、パラメータの誤差 統計的データ解析 008 008.. 林田清 ( 大阪大学大学院理学研究科 ) 問題 C (, ) ( x xˆ) ( y yˆ) σ x πσ σ y y Pabx (, ;,,, ) ˆ y σx σ y = dx exp exp πσx ただし xy ˆ ˆ はyˆ = axˆ+ bであらわされる直線モデル上の点 ( ˆ) ( ˆ ) ( ) x x y ax b y ax b Pabx (,

More information

際 正規分布に従わない観測値に対して通常の t 検定を適用した場合 どのような不都合が生じるかを考える 一般に通常の t 検定や Wilcoxon 検定などの仮説検定を行う場合 2つの処理の間に差がないことが真実であるにもかかわらず差があると主張する過誤確率 ( 第 1 種の過誤確率 ) 2つの処理

際 正規分布に従わない観測値に対して通常の t 検定を適用した場合 どのような不都合が生じるかを考える 一般に通常の t 検定や Wilcoxon 検定などの仮説検定を行う場合 2つの処理の間に差がないことが真実であるにもかかわらず差があると主張する過誤確率 ( 第 1 種の過誤確率 ) 2つの処理 連載 第 2 回 医学データの統計解析の基本 2 つの平均の比較 * 朝倉こう子 濱﨑俊光 Fundamentals of statistical analysis in biomedical research:two-sample tests for comparing means 1 基礎研究や臨床研究を問わず医学研究において 新しい化合物や治療法を発見し その性能を特徴づける場合 何らかの対照

More information

ビジネス統計 統計基礎とエクセル分析 正誤表

ビジネス統計 統計基礎とエクセル分析 正誤表 ビジネス統計統計基礎とエクセル分析 ビジネス統計スペシャリスト エクセル分析スペシャリスト 公式テキスト正誤表と学習用データ更新履歴 平成 30 年 5 月 14 日現在 公式テキスト正誤表 頁場所誤正修正 6 知識編第 章 -3-3 最頻値の解説内容 たとえば, 表.1 のデータであれば, 最頻値は 167.5cm というたとえば, 表.1 のデータであれば, 最頻値は 165.0cm ということになります

More information

<4D F736F F F696E74202D2088E38A77939D8C7695D78BAD89EF313691E63589F194E497A682C695AA8A84955C2E >

<4D F736F F F696E74202D2088E38A77939D8C7695D78BAD89EF313691E63589F194E497A682C695AA8A84955C2E > 26// 第 5 回 医学統計勉強会 東北大学病院循環器内科 東北大学病院臨床研究推進センター 共催 東北大学大学院医学系研究科 EBM 開発学寄附講座 宮田 敏 比率と分割表 疾患の発症率など, 物事の頻度 (frequency) を議論する際, 以下の三つの概念を使い分ける. 比 (ratio):a, B ( ) が存在するとき,A/B を比という. A と B は互いを含まない. 例 : 性比.BMI=

More information

Excelによる統計分析検定_知識編_小塚明_5_9章.indd

Excelによる統計分析検定_知識編_小塚明_5_9章.indd 第7章57766 検定と推定 サンプリングによって得られた標本から, 母集団の統計的性質に対して推測を行うことを統計的推測といいます 本章では, 推測統計の根幹をなす仮説検定と推定の基本的な考え方について説明します 前章までの知識を用いて, 具体的な分析を行います 本章以降の知識は操作編での操作に直接関連していますので, 少し聞きなれない言葉ですが, 帰無仮説 有意水準 棄却域 などの意味を理解して,

More information

スライド 1

スライド 1 データ解析特論重回帰分析編 2017 年 7 月 10 日 ( 月 )~ 情報エレクトロニクスコース横田孝義 1 ( 単 ) 回帰分析 単回帰分析では一つの従属変数 ( 目的変数 ) を 一つの独立変数 ( 説明変数 ) で予測する事を考える 具体的には y = a + bx という回帰直線 ( モデル ) でデータを代表させる このためにデータからこの回帰直線の切片 (a) と傾き (b) を最小

More information

Microsoft Word - å“Ÿåłžå¸°173.docx

Microsoft Word - å“Ÿåłžå¸°173.docx 回帰分析 ( その 3) 経済情報処理 価格弾力性の推定ある商品について その購入量を w 単価を p とし それぞれの変化量を w p で表 w w すことにする この時 この商品の価格弾力性 は により定義される これ p p は p が 1 パーセント変化した場合に w が何パーセント変化するかを示したものである ここで p を 0 に近づけていった極限を考えると d ln w 1 dw dw

More information

抗悪性腫瘍薬の開発における臨床試験エンドポイント 一般的には, 全生存期間 (OS) が真のエンドポイントとして考えられている. OS の定義 : ランダム化からあらゆる原因による死亡までの期間 OS を主要評価項目とした臨床試験を実施する場合, がん腫によっては, 非常に試験期間が長くなってしまう

抗悪性腫瘍薬の開発における臨床試験エンドポイント 一般的には, 全生存期間 (OS) が真のエンドポイントとして考えられている. OS の定義 : ランダム化からあらゆる原因による死亡までの期間 OS を主要評価項目とした臨床試験を実施する場合, がん腫によっては, 非常に試験期間が長くなってしまう 抗悪性腫瘍薬の開発における臨床試験エンドポイント 一般的には, 全生存期間 (OS) が真のエンドポイントとして考えられている. OS の定義 : ランダム化からあらゆる原因による死亡までの期間 OS を主要評価項目とした臨床試験を実施する場合, がん腫によっては, 非常に試験期間が長くなってしまう. 無増悪生存期間 (PFS) や無病生存期間 (DFS), 客観的奏効率など,OS よりも短期間で評価可能なエンドポイントを主要評価項目とした臨床試験の実施が検討されることが多い.

More information

Jupyter Notebook を活用したプログラムライブラリ構築の検討 吹谷芳博 1, 藤澤正樹 1 ( 1 あすか製薬株式会社 ) Examination of the program library construction using Jupyter Notebook ASKA Pharm

Jupyter Notebook を活用したプログラムライブラリ構築の検討 吹谷芳博 1, 藤澤正樹 1 ( 1 あすか製薬株式会社 ) Examination of the program library construction using Jupyter Notebook ASKA Pharm Jupyter Notebook を活用したプログラムライブラリ構築の検討 吹谷芳博 1, 藤澤正樹 1 ( 1 あすか製薬株式会社 ) Examination of the program library construction using Jupyter Notebook ASKA Pharmaceutical Co., Ltd. 要旨 : PC SAS での Jupyter Notebook

More information

SAS Enterprise Guideによるデータ解析入門

SAS Enterprise Guideによるデータ解析入門 ........ 1 / 70.... SAS Enterprise Guide Kengo NAGASHIMA Laboratory of Biostatistics, Department of Parmaceutical Technochemistry, Josai University 2010 11 16 ........ 2 / 70 (SAS / SAS Enterprise Guide

More information

情報工学概論

情報工学概論 確率と統計 中山クラス 第 11 週 0 本日の内容 第 3 回レポート解説 第 5 章 5.6 独立性の検定 ( カイ二乗検定 ) 5.7 サンプルサイズの検定結果への影響練習問題 (4),(5) 第 4 回レポート課題の説明 1 演習問題 ( 前回 ) の解説 勉強時間と定期試験の得点の関係を無相関検定により調べる. データ入力 > aa

More information

Mantel-Haenszelの方法

Mantel-Haenszelの方法 Mantel-Haenszel 2008 6 12 ) 2008 6 12 1 / 39 Mantel & Haenzel 1959) Mantel N, Haenszel W. Statistical aspects of the analysis of data from retrospective studies of disease. J. Nat. Cancer Inst. 1959; 224):

More information

ANOVA

ANOVA 3 つ z のグループの平均を比べる ( 分散分析 : ANOVA: analysis of variance) 分散分析は 全体として 3 つ以上のグループの平均に差があるか ということしかわからないために, どのグループの間に差があったかを確かめるには 多重比較 という方法を用います これは Excel だと自分で計算しなければならないので, 分散分析には統計ソフトを使った方がよいでしょう 1.

More information

untitled

untitled Data cleaning Original datan=8479 NCC/ptkgN=7958 NCC10x10E628 8139 NCC/kgPage 1-2) Covariate cleaningcovariate Page Page5-8 Kaplan-Meier method, Log-rank, Cox hazard model Overall survival Cumulative incidence

More information

スライド 1

スライド 1 . 無情報事前分布 前回 前回の復習 データの分布 ( 統計モデル を設定 ( θ モデルごとに相性のよい事前分布 ( 共役事前分布 を紹介 事後分布の計算 π (θ π ( θ, π ( θ dθ, 昔 共役事前分布を利用して 解析的に事後分布を導出簡単な形で書けない場合は ラプラス近似を利用 現在 パラメータの次元が高い, もしくは複雑な入り方をする統計モデル 数値的な方法 ( 第六回 で近似するため共役事前分布にこだわらない

More information

第7章

第7章 5. 推定と検定母集団分布の母数を推定する方法と仮説検定の方法を解説する まず 母数を一つの値で推定する点推定について 推定精度としての標準誤差を説明する また 母数が区間に存在することを推定する信頼区間も取り扱う 後半は統計的仮説検定について述べる 検定法の基本的な考え方と正規分布および二項確率についての検定法を解説する 5.1. 点推定先に述べた統計量は対応する母数の推定値である このように母数を一つの値およびベクトルで推定する場合を点推定

More information

Microsoft Word - 補論3.2

Microsoft Word - 補論3.2 補論 3. 多変量 GARC モデル 07//6 新谷元嗣 藪友良 対数尤度関数 3 章 7 節では 変量の対数尤度を求めた ここでは多変量の場合 とくに 変量について対数尤度を求める 誤差項 は平均 0 で 次元の正規分布に従うとする 単純化のため 分散と共分散は時間を通じて一定としよう ( この仮定は後で変更される ) したがって ij から添え字 を除くことができる このとき と の尤度関数は

More information

様々なミクロ計量モデル†

様々なミクロ計量モデル† 担当 : 長倉大輔 ( ながくらだいすけ ) この資料は私の講義において使用するために作成した資料です WEB ページ上で公開しており 自由に参照して頂いて構いません ただし 内容について 一応検証してありますが もし間違いがあった場合でもそれによって生じるいかなる損害 不利益について責任を負いかねますのでご了承ください 間違いは発見次第 継続的に直していますが まだ存在する可能性があります 1 カウントデータモデル

More information

スライド 1

スライド 1 線形モデルにおける CLASS ステートメントの機能 吉田早織 1 魚住龍史 2 1 日本化薬株式会社医薬データセンター 2 京都大学大学院医学研究科 The fascinating features for the CLASS in the context of linear models Saori Yoshida 1 and Ryuji Uozumi 2 1 Clinical Data Management

More information

Microsoft Word - tmp_input.doc

Microsoft Word - tmp_input.doc アバスチン 2.7.3 臨床的有効性の概要 Page 19 2.7.3.2.4 NO16966 試験 (FOLFOX-4 療法及び XELOX 療法との併用による無作為化比較試験, 一次治療 ) 本試験は転移性結腸 直腸癌の一次治療例を対象とした第 Ⅲ 相臨床試験であり, 当初,XELOX 療法 (L-OHP/ カペシタビン ) と FOLFOX-4 療法 (L-OHP/5-FU/LV) の非盲検比較試験として開始されたが

More information

Microsoft PowerPoint - åł€æžœæ”¨è«Œã†®å�ºç¤”(çfl°æ€Š)final

Microsoft PowerPoint - åł€æžœæ”¨è«Œã†®å�ºç¤”(çfl°æ€Š)final 日本製薬工業協会シンポジウム 因果推論の基礎 横浜市立大学医学部 臨床統計学 田栗正隆 1 検証的臨床試験の目的 治療法の効果を適切に評価すること 因果効果の推定に興味 薬剤の承認申請を行うための第 III 相治験では ランダム化比較試験が行われる 2 ランダム化比較試験 治療 A を対象者にランダムに割付 A を受けたかどうかのみが違う 2 群が作れる 比較可能性 (comparability)

More information

ii 3.,. 4. F. (), ,,. 8.,. 1. (75% ) (25% ) =9 7, =9 8 (. ). 1.,, (). 3.,. 1. ( ).,.,.,.,.,. ( ) (1 2 )., ( ), 0. 2., 1., 0,.

ii 3.,. 4. F. (), ,,. 8.,. 1. (75% ) (25% ) =9 7, =9 8 (. ). 1.,, (). 3.,. 1. ( ).,.,.,.,.,. ( ) (1 2 )., ( ), 0. 2., 1., 0,. 23(2011) (1 C104) 5 11 (2 C206) 5 12 http://www.math.is.tohoku.ac.jp/~obata,.,,,.. 1. 2. 3. 4. 5. 6. 7.,,. 1., 2007 ( ). 2. P. G. Hoel, 1995. 3... 1... 2.,,. ii 3.,. 4. F. (),.. 5.. 6.. 7.,,. 8.,. 1. (75%

More information

 

  早稲田大学大学院理工学研究科 博士論文概要 論文題目 Various statistical methods in time series analysis 時系列解析における種々の統計手法 申請者 天野友之 Tomoyuki AMANO 数理科学専攻数理統計学研究 007 年 月 時とともに変動する偶然量の観測値の系列である時系列の解析は近年 様々な統計手法が導入され自然工学 医学 経済学 など多方面で急速に発展している

More information

自動車感性評価学 1. 二項検定 内容 2 3. 質的データの解析方法 1 ( 名義尺度 ) 2.χ 2 検定 タイプ 1. 二項検定 官能検査における分類データの解析法 識別できるかを調べる 嗜好に差があるかを調べる 2 点比較法 2 点識別法 2 点嗜好法 3 点比較法 3 点識別法 3 点嗜好

自動車感性評価学 1. 二項検定 内容 2 3. 質的データの解析方法 1 ( 名義尺度 ) 2.χ 2 検定 タイプ 1. 二項検定 官能検査における分類データの解析法 識別できるかを調べる 嗜好に差があるかを調べる 2 点比較法 2 点識別法 2 点嗜好法 3 点比較法 3 点識別法 3 点嗜好 . 内容 3. 質的データの解析方法 ( 名義尺度 ).χ 検定 タイプ. 官能検査における分類データの解析法 識別できるかを調べる 嗜好に差があるかを調べる 点比較法 点識別法 点嗜好法 3 点比較法 3 点識別法 3 点嗜好法 : 点比較法 : 点識別法 配偶法 配偶法 ( 官能評価の基礎と応用 ) 3 A か B かの判定において 回の判定でAが選ばれる回数 kは p の二項分布に従う H :

More information

統計学 - 社会統計の基礎 - 正規分布 標準正規分布累積分布関数の逆関数 t 分布正規分布に従うサンプルの平均の信頼区間 担当 : 岸 康人 資料ページ :

統計学 - 社会統計の基礎 - 正規分布 標準正規分布累積分布関数の逆関数 t 分布正規分布に従うサンプルの平均の信頼区間 担当 : 岸 康人 資料ページ : 統計学 - 社会統計の基礎 - 正規分布 標準正規分布累積分布関数の逆関数 t 分布正規分布に従うサンプルの平均の信頼区間 担当 : 岸 康人 資料ページ : https://goo.gl/qw1djw 正規分布 ( 復習 ) 正規分布 (Normal Distribution)N (μ, σ 2 ) 別名 : ガウス分布 (Gaussian Distribution) 密度関数 Excel:= NORM.DIST

More information

PowerPoint プレゼンテーション

PowerPoint プレゼンテーション 1/X Chapter 9: Linear correlation Cohen, B. H. (2007). In B. H. Cohen (Ed.), Explaining Psychological Statistics (3rd ed.) (pp. 255-285). NJ: Wiley. 概要 2/X 相関係数とは何か 相関係数の数式 検定 注意点 フィッシャーのZ 変換 信頼区間 相関係数の差の検定

More information

数値計算法

数値計算法 数値計算法 011/5/5 林田清 ( 大阪大学大学院理学研究科 ) レポート課題 1( 締め切りは 5/5) 平均値と標準偏差を求めるプログラム 入力 : データの数 データ データは以下の 10 個 ( 例えばある月の最高気温 ( )10 日分 ) 34.3,5.0,3.,34.6,.9,7.7,30.6,5.8,3.0,31.3 出力 :( 標本 ) 平均値 標準偏差 ソースプログラムと出力結果をメイルの本文にして

More information

Slide 1

Slide 1 ODS 統 計 グラフ 機 能 を 用 いたグラフの 作 成 SAS Institute Japan 株 式 会 社 プロフェッショナルサービス 本 部 テクニカルサポート 部 深 澤 武 志 Copyright 2010 SAS Institute Inc. All rights reserved. 目 次 ODSとは ODS 統 計 グラフのご 紹 介 ODSテンプレートに 関 して ODS 統

More information

日本製薬工業協会シンポジウム 我々は Time to event data をどのように要約し, 評価するのか 生存時間データの要約指標と適用事例の紹介 2018 年 12 月 6 日医薬品評価委員会データサイエンス部会タスクフォース 4 生存時間解析チーム 塩野義製薬 ( 株 ) 長谷川貴大大鵬薬

日本製薬工業協会シンポジウム 我々は Time to event data をどのように要約し, 評価するのか 生存時間データの要約指標と適用事例の紹介 2018 年 12 月 6 日医薬品評価委員会データサイエンス部会タスクフォース 4 生存時間解析チーム 塩野義製薬 ( 株 ) 長谷川貴大大鵬薬 日本製薬工業協会シンポジウム 我々は Time to event data をどのように要約し, 評価するのか 生存時間データの要約指標と適用事例の紹介 2018 年 12 月 6 日医薬品評価委員会データサイエンス部会タスクフォース 4 生存時間解析チーム 塩野義製薬 ( 株 ) 長谷川貴大大鵬薬品工業 ( 株 ) 棚瀬貴紀 留意点 本発表は, 先日公開された 生存時間型応答の評価指標 -RMST(restricted

More information

スライド 1

スライド 1 データ解析特論第 10 回 ( 全 15 回 ) 2012 年 12 月 11 日 ( 火 ) 情報エレクトロニクス専攻横田孝義 1 終了 11/13 11/20 重回帰分析をしばらくやります 12/4 12/11 12/18 2 前回から回帰分析について学習しています 3 ( 単 ) 回帰分析 単回帰分析では一つの従属変数 ( 目的変数 ) を 一つの独立変数 ( 説明変数 ) で予測する事を考える

More information

Microsoft PowerPoint - A1.ppt [互換モード]

Microsoft PowerPoint - A1.ppt [互換モード] 011/4/13 付録 A1( 推測統計学の基礎 ) 付録 A1 推測統計学の基礎 1. 統計学. カイ 乗検定 3. 分散分析 4. 相関係数 5. 多変量解析 1. 統計学 3 統計ソフト 4 記述統計学 推測統計学 検定 ノンパラメトリック検定名義 / 分類尺度順序 / 順位尺度パラメトリック検定間隔 / 距離尺度比例 / 比率尺度 SAS SPSS R R-Tps (http://cse.aro.affrc.go.jp/takezawa/r-tps/r.html)

More information

Microsoft Word - Stattext11.doc

Microsoft Word - Stattext11.doc 章母集団と指定値との量的データの検定. 検定手順 前章で質的データの検定手法について説明しましたので ここからは量的データの検定について話します 量的データの検定は少し分量が多くなりますので 母集団と指定値との検定 対応のない 群間の検定 対応のある 群間の検定 と 3つに章を分けて話を進めることにします ここでは 母集団と指定値との検定について説明します 例えば全国平均が分かっている場合で ある地域の標本と全国平均を比較するような場合や

More information

untitled

untitled WinLD R (16) WinLD https://www.biostat.wisc.edu/content/lan-demets-method-statistical-programs-clinical-trials WinLD.zip 2 2 1 α = 5% Type I error rate 1 5.0 % 2 9.8 % 3 14.3 % 5 22.6 % 10 40.1 % 3 Type

More information

ペムブロリズマブ ( 遺伝子組換え ) 注射剤 2.7 臨床概要 臨床的有効性 ORR 海外 001 試験 パート D パート D では 治験担当医師がベースライン時点で測定可能病変ありとし 独立中央判定

ペムブロリズマブ ( 遺伝子組換え ) 注射剤 2.7 臨床概要 臨床的有効性 ORR 海外 001 試験 パート D パート D では 治験担当医師がベースライン時点で測定可能病変ありとし 独立中央判定 2.7.3.3.2.3 ORR 2.7.3.3.2.3.1 海外 001 試験 2.7.3.3.2.3.1.1 パート D パート D では 治験担当医師がベースライン時点で測定可能病変ありとし 独立中央判定機関ではなしと判定された患者が11 例 (10.7%) であったため FAS は APaT 集団とは異なる FAS BOR の解析結果を治験総括報告書第 2 版 [ 資料 5.3.5.2.1.2:

More information

<4D F736F F D2090B695A8939D8C768A E F AA957A82C682948C9F92E8>

<4D F736F F D2090B695A8939D8C768A E F AA957A82C682948C9F92E8> 第 8 回 t 分布と t 検定 生物統計学 A.t 分布 ( 小標本に関する平均の推定と検定 ) 前々回と前回の授業では, 標本が十分に大きいあるいは母分散が既知であることを条件に正規分布を用いて推定 検定した. しかし, 母集団が正規分布し, 標本が小さい場合には, 標本分散から母分散を推定するときの不確実さを加味したt 分布を用いて推定 検定しなければならない. t 分布は標本分散の自由度 f(

More information

「スウェーデン企業におけるワーク・ライフ・バランス調査 」報告書

「スウェーデン企業におけるワーク・ライフ・バランス調査 」報告書 1 2004 12 2005 4 5 100 25 3 1 76 2 Demoskop 2 2004 11 24 30 7 2 10 1 2005 1 31 2 4 5 2 3-1-1 3-1-1 Micromediabanken 2005 1 507 1000 55.0 2 77 50 50 /CEO 36.3 37.4 18.1 3-2-1 43.0 34.4 / 17.6 3-2-2 78 79.4

More information

Microsoft Word - StatsDirectMA Web ver. 2.0.doc

Microsoft Word - StatsDirectMA Web ver. 2.0.doc Web version. 2.0 15 May 2006 StatsDirect ver. 2.0 15 May 2006 2 2 2 Meta-Analysis for Beginners by using the StatsDirect ver. 2.0 15 May 2006 Yukari KAMIJIMA 1), Ataru IGARASHI 2), Kiichiro TSUTANI 2)

More information

Microsoft PowerPoint - R-stat-intro_13.ppt [互換モード]

Microsoft PowerPoint - R-stat-intro_13.ppt [互換モード] R で統計解析入門 (13) 生存時間解析 後篇 本日のメニュー 1. 競合リスクに関する解析 2. 再発事象の解析 2 復習 カプランマイヤー推定量 5 人のがん患者さんに薬物療法を行い, ガンの再発 をイベントと してカプランマイヤー推定量によりイベント発生割合を計算する 何らかの理由でイベントを発生せずに観察を終了した場合は打ち切り 時間リスクイベントイベントイベント打ち切り ( 日 ) 集合無発生割合累積発生割合

More information

PowerPoint プレゼンテーション

PowerPoint プレゼンテーション 学位論文作成のための疫学 統計解析の実際 徳島大学大学院 医歯薬学研究部 社会医学系 予防医学分野 有澤孝吉 (e-mail: karisawa@tokushima-u.ac.jp) 本日の講義の内容 (SPSS を用いて ) 記述統計 ( データのまとめ方 ) 代表値 ばらつき正規確率プロット 正規性の検定標準偏差 不偏標準偏差 標準誤差の区別中心極限定理母平均の区間推定 ( 母集団の標準偏差が既知の場合

More information

Microsoft PowerPoint - ch04j

Microsoft PowerPoint - ch04j Ch.4 重回帰分析 : 推論 重回帰分析 y = 0 + 1 x 1 + 2 x 2 +... + k x k + u 2. 推論 1. OLS 推定量の標本分布 2. 1 係数の仮説検定 : t 検定 3. 信頼区間 4. 係数の線形結合への仮説検定 5. 複数線形制約の検定 : F 検定 6. 回帰結果の報告 入門計量経済学 1 入門計量経済学 2 OLS 推定量の標本分布について OLS 推定量は確率変数

More information

*1 * Wilcoxon 2 2 t t t t d t M t N t M t n t N t n t N t d t N t t at ri

*1 * Wilcoxon 2 2 t t t t d t M t N t M t n t N t n t N t d t N t t at ri Wilcoxon H23 BioS 1 Wilcoxon 2 2.1 1 2 1 0 1 1 5 0 1 2 7 0 1 3 8 1 1 4 12 0 2 5 2 0 2 6 3 1 2 7 4 1 2 8 10 0 Wilcoxon 2.2 S 1 t S 2 t Wilcoxon H 0 H 1 H 0 : S 1 t S 2 t H 1 : S 1 t S 2 t 1 *1 *2 2.3 2.3.1

More information

Microsoft Word _肺がん統計解析資料.docx

Microsoft Word _肺がん統計解析資料.docx 治療症例数第 2 位 : (2015/1-2017/9) 統計解析資料 A) はじめに免疫治療効果の成否に大きく関与するT 細胞を中心とした免疫機構は 細胞内に進入した外来生物の排除ならびに対移植片拒絶や自己免疫疾患 悪性腫瘍の発生進展に深く関与している これら細胞性免疫機構は担癌者においてその機能の低下が明らかとなり 近年では腫瘍免疫基礎研究において各種免疫学的パラメータ解析によるエビデンスに基づいた治療手法が大きく注目されるようになった

More information

統計学 Ⅱ( 章 ( 区間推定のシミュレーション 母平均 μ の区間推定 X ~ N, のとき X T ~ 自由度 1の t分布 1 自由度 -1のt 分布の97.5% 点 :t.975 P t T t この式に T を代入する t.975 母集団

統計学 Ⅱ( 章 ( 区間推定のシミュレーション 母平均 μ の区間推定 X ~ N, のとき X T ~ 自由度 1の t分布 1 自由度 -1のt 分布の97.5% 点 :t.975 P t T t この式に T を代入する t.975 母集団 統計学 Ⅱ(16 11-1 章 11 章母集団パラメータの推定 1. 信頼区間 (1 点推定と区間推定 ( 区間推定のシミュレーション (3 母平均 μの信頼区間 (4 母比率 pの信頼区間 (5 母比率 pのより厳密な信頼区間. 点推定量の特性 (1 標本平均 X の持つ望ましい性質 ( 不偏性 (3 推定量の分散と有効性 (4 平均 乗誤差 MEと最小分散性 (5 一致性 (6 チェビシェフの不等式

More information

仮説検定を伴う方法では 検定の仮定が満たされ 検定に適切な検出力があり データの分析に使用される近似で有効な結果が得られることを確認することを推奨します カイ二乗検定の場合 仮定はデータ収集に固有であるためデータチェックでは対応しません Minitab は近似法の検出力と妥当性に焦点を絞っています

仮説検定を伴う方法では 検定の仮定が満たされ 検定に適切な検出力があり データの分析に使用される近似で有効な結果が得られることを確認することを推奨します カイ二乗検定の場合 仮定はデータ収集に固有であるためデータチェックでは対応しません Minitab は近似法の検出力と妥当性に焦点を絞っています MINITAB アシスタントホワイトペーパー本書は Minitab 統計ソフトウェアのアシスタントで使用される方法およびデータチェックを開発するため Minitab の統計専門家によって行われた調査に関する一連の文書の 1 つです カイ二乗検定 概要 実際には 連続データの収集が不可能な場合や難しい場合 品質の専門家は工程を評価するためのカテゴリデータの収集が必要となることがあります たとえば 製品は不良

More information

サーバに関するヘドニック回帰式(再推計結果)

サーバに関するヘドニック回帰式(再推計結果) 2012 年 3 月 日本銀行調査統計局 企業物価指数 サーバ に関するヘドニック回帰式 ( 再推計結果 ) 企業物価指数 サーバ の品質調整に適用するヘドニック回帰式について 1 最新のデータを用いて再推計しましたので その結果をお知らせします 1. サーバのヘドニック推計に関する基本方針 留意事項推計頻度 年 1 回 (2 月 ) 適用範囲 国内品 輸出品 輸入品に対し 同一の推計式を適用 2

More information

Microsoft PowerPoint 古川杉本SASWEB用プレゼン.ppt

Microsoft PowerPoint 古川杉本SASWEB用プレゼン.ppt ロジスティックモデルと ROC AUC 分析を 組み合わせた検査性能の評価と 疫学基本モデル評価方法 古川敏仁 杉本典子株式会社バイオスタティスティカルリサーチ Test Perforance Evaluation in Epideiological Basic Model Using ROC AUC with logistic regression Toshihito Furukawa, Noriko

More information

kubo2015ngt6 p.2 ( ( (MLE 8 y i L(q q log L(q q 0 ˆq log L(q / q = 0 q ˆq = = = * ˆq = 0.46 ( 8 y 0.46 y y y i kubo (ht

kubo2015ngt6 p.2 ( ( (MLE 8 y i L(q q log L(q q 0 ˆq log L(q / q = 0 q ˆq = = = * ˆq = 0.46 ( 8 y 0.46 y y y i kubo (ht kubo2015ngt6 p.1 2015 (6 MCMC kubo@ees.hokudai.ac.jp, @KuboBook http://goo.gl/m8hsbm 1 ( 2 3 4 5 JAGS : 2015 05 18 16:48 kubo (http://goo.gl/m8hsbm 2015 (6 1 / 70 kubo (http://goo.gl/m8hsbm 2015 (6 2 /

More information

7 / 70 8 / 70 SAS Enterprise Guide, Windows SAS, SAS,,,, SAS SAS Enterprise Guide SAS: SAS Enterprise Guide:, SAS SAS Enterprise Guide... 5 / 70 6 / 7

7 / 70 8 / 70 SAS Enterprise Guide, Windows SAS, SAS,,,, SAS SAS Enterprise Guide SAS: SAS Enterprise Guide:, SAS SAS Enterprise Guide... 5 / 70 6 / 7 3 / 70 4 / 70 SAS Enterprise Guide (SAS / SAS Enterprise Guide ) Kengo NAGASHIMA Laboratory of Biostatistics, Department of Parmaceutical Technochemistry, Josai University 2010 11 16 1 / 70 2 / 70 SAS

More information

目次 1 章 SPSS の基礎 基本 はじめに 基本操作方法 章データの編集 はじめに 値ラベルの利用 計算結果に基づく新変数の作成 値のグループ化 値の昇順

目次 1 章 SPSS の基礎 基本 はじめに 基本操作方法 章データの編集 はじめに 値ラベルの利用 計算結果に基づく新変数の作成 値のグループ化 値の昇順 SPSS 講習会テキスト 明治大学教育の情報化推進本部 IZM20140527 目次 1 章 SPSS の基礎 基本... 3 1.1 はじめに... 3 1.2 基本操作方法... 3 2 章データの編集... 6 2.1 はじめに... 6 2.2 値ラベルの利用... 6 2.3 計算結果に基づく新変数の作成... 7 2.4 値のグループ化... 8 2.5 値の昇順 降順... 10 3

More information

生命情報学

生命情報学 生命情報学 5 隠れマルコフモデル 阿久津達也 京都大学化学研究所 バイオインフォマティクスセンター 内容 配列モチーフ 最尤推定 ベイズ推定 M 推定 隠れマルコフモデル HMM Verアルゴリズム EMアルゴリズム Baum-Welchアルゴリズム 前向きアルゴリズム 後向きアルゴリズム プロファイル HMM 配列モチーフ モチーフ発見 配列モチーフ : 同じ機能を持つ遺伝子配列などに見られる共通の文字列パターン

More information