R– category –

三群以上の平均値を多重比較したい。 でも各群の背景因子がそろっていない。 背景因子を調整しながら三群以上の平均値を多重比較するにはどうすればいいか？ R でのやり方を解説する。 共変量を調整して多重比較する方法 群ごとの背景因子は、群分け変数と...

ノンパラメトリックの多重比較をRで実施する方法。 ノンパラメトリックとは何か？ ノンパラメトリックとはパラメトリックではないという意味。 パラメトリックとは、パラメーターを使うという意味だ。 パラメーターとは、日本語では母数（ぼすう）と言われ...

フィッシャーの正確確率検定、カイ二乗検定の3群以上の比較をRで実施する方法の解説。 カイ二乗検定の3群以上の比較 三群以上の割合の比較はどうやればいいのか？ Bonferroni型のp値調整を使う方法がある。 Rで行う場合、pairwise.prop.test()という関数を...

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

Tukey HSD検定をRで行う方法の解説。 Tukey HSD検定をRで行う方法 aov()とTukeyHSD()という二つの関数を使う。 ダネット検定のときと同じように、例としてwarpbreaksというデータを使う。 ダネット検定は以下を参照。 機織りにおいて、tension（緊張、張り...

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

クロスオーバー試験のサンプルサイズ計算を R で行う方法 クロスオーバー試験とは クロスオーバー試験とは、二群どちらにも二つの介入を受けてもらう試験。 従来の薬と新薬とか、方法Aと方法Bとか、二つの方法を順番を変えて受けてもらう。 利点は、小さい...

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

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

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