EZR– category –

Rで、一元配置分散分析を step by step で計算してみた。 lm() と Anova() を使えばあっという間だが、具体的な一つ一つの計算を自分で組み立ててみるとどうか？ 教科書の例題に沿って確認した。 R で一元配置分散分析を実行するための例題データ サンプル...

検証試験において、三群以上の平均値を比較したいときに、単純に二群比較を繰り返すと有意水準が甘くなる。 有意水準の調整によって簡単に処理する方法がボンフェローニ （Bonferroni）の方法とその進化版だ。 Rでボンフェローニ型のP値調整で多重比較を行...

ダネット検定は、比較対照群といくつかの実験群を多重比較する方法。 Rでダネット検定をするにはどうしたらよいか？ Rでダネット検定をするには？ まずmultcompパッケージをインストール。 インストールは一回だけでOK。 install.packages("multcomp") 次...

非劣性試験のサンプルサイズ計算 生存時間データの場合 非劣性試験のサンプルサイズ計算 生存時間データの場合 例 1 各記号の意味は以下のとおり。 t：観察期間。 S1とS0：それぞれ治療群とコントロール群の生存率。 dF：Freedmanの方法による各群必要な死...

非劣性試験のエンドポイントが連続量で、各群の平均値を求める場合、サンプルサイズはどのように計算するか？ RとEZRで行う方法。 非劣性試験における非劣性検定とは？ 劣っていないことを積極的に証明する検定。 臨床的に意味のある差を設定して、それよ...

非劣性検定は劣っていないことを証明する検定。 割合の非劣性検定のサンプルサイズ計算はどうやるか？ 割合の非劣性検定のサンプルサイズ計算 各変数の意味は、以下のとおりとする。 pA：試験薬の有効率とする。 pB：標準薬の有効率とする。 pB.bar：帰無...

生存時間解析のサンプルサイズ計算の方法 Cox の比例ハザードモデルを使う前提の計算 生存時間解析のサンプルサイズ計算 その 1 R でスクリプトを書いてみた。 S0がコントロール、S1が治療群、それぞれの生存率。 dが、両群合わせて合計の死亡者数。 sampl...

カプランマイヤー曲線の書き方 R と EZR の方法。 カプランマイヤー曲線の書き方説明のためのサンプルデータ カプランマイヤー曲線を描きたい場合、イベント発生のデータとイベント発生までもしくは観察終了までの時間のデータが必要。 イベント発生データ...