統計ブログはじめました!

各専門分野の統計技術、方法、テクニックなどを気ままに分かり易く例題をもとに解説します。

統計のコツのこつ(24)

2016-11-07 11:32:48 | 日記・エッセイ・コラム
 
前回に引き続き生存分析のお話です。
それでは、「すぐに役立つ統計のコツ」第8章(134ページ)を開いて下さい。
「K-M法による Logrank test」を「R」で行ってみましょう。
前回のデータであれば、次のコマンド(関数)だけでOKです。
 
***
survdiff(Surv(観察期間, イベント発生 == 1)~ グループ, data=dat)

出力結果:第8章(136ページ参照)
> survdiff(Surv(観察期間, イベント発生 == 1)~ グループ, data=dat)
Call:
survdiff(formula = Surv(観察期間, イベント発生 == 1) ~ グループ,
    data = dat)
            N Observed Expected (O-E)^2/E (O-E)^2/V
グループ=A 10        3     6.13      1.60      4.03
グループ=B 11        8     4.87      2.02      4.03
 Chisq= 4  on 1 degrees of freedom, p= 0.0448
***
 
この結果は、「すぐに役立つ統計のコツ」第8章(137ページ)を見て下さい。
 
次回は、「コックス比例ハザード・モデル」をご紹介します。
 
情報統計研究所はここから!
 
ジャンル:
きいて!きいて!
この記事についてブログを書く
この記事をはてなブックマークに追加
« 統計のコツのこつ(23) | トップ | 統計のコツのこつ(25) »
最近の画像もっと見る

日記・エッセイ・コラム」カテゴリの最新記事