R– category –
-
アンケート調査には何人以上必要か R で計算する方法
アンケート調査には何人以上必要か アンケート調査には何人以上必要か 無限母集団の場合 よくある質問。 「アンケート調査は何人にとったらよいでしょうか。」 「世論調査のサンプル数は本当に意味があるのでしょうか。」 「視聴率調査のサンプル数はどう... -
R で割合の推定・検定に必要なサンプル数を計算する方法
割合に関するサンプルサイズ計算の方法 推定精度と検定の 2 種類あり 推定精度に基づく方法 相対精度deltaで，割合を推定する際のサンプルサイズ計算は、以下のような R スクリプトで計算できる。 myPsize <- function(p,delta){ n <- 4*(1-p)/(p*de... -
R で対応のある t 検定に必要なサンプル数を計算する方法
対応のある t 検定は、前後比較をするデザインの時に用いる。 対応のある t 検定のサンプルサイズ、サンプル数はどのくらいであればよいか？ 対応のある t 検定のサンプルサイズ計算に必要な数値 サンプルサイズ計算に必要なのは、 検出力 標準化された前... -
R で t 検定に必要なサンプル数を計算する方法
t 検定のサンプル数計算の方法 t 検定のサンプル数計算（1:1の場合） R の power.t.test()を用いる。 必要な情報は、検出力をどうするかと、二群の平均値の差。 標準偏差で割った標準化された差が必要だ。 検出力は慣例で80%＝0.8が多い。 90%＝0.9もよく... -
R と EZR でカプランマイヤー曲線を書く方法
カプランマイヤー曲線の書き方 R と EZR の方法。 カプランマイヤー曲線の書き方説明のためのサンプルデータ カプランマイヤー曲線を描きたい場合、イベント発生のデータとイベント発生までもしくは観察終了までの時間のデータが必要。 イベント発生データ... -
R でグラフの軸目盛・軸ラベル・色の変更と値表示をする方法
統計ソフトRのグラフ上に数値を表示させたり、軸の目盛間隔や軸ラベルを変更する方法を紹介。 軸の数値をすべて縦にする las=1 軸ラベルのスタイル（the style of axis labels）。1は縦軸も横軸も縦方向の数字にする意味。 例： library(MASS) boxplot(bwt... -
R で QQ プロットを書く方法
正規確率プロット QQプロット は、正規分布をしているかどうか、目視で確認するための方法。 図で正規分布のチェック qqnorm()を使う。 正規確率プロットを描く。 x軸が正規分布の論理的な分位数、y軸が実測値の分位数。 対角線上にまっすぐに、プロットが... -
R で連続データの 2 群比較を行う方法
実験群とコントロール群、リスク因子あり群となし群、介入群と非介入群、など二群比較 を R で実行する方法。 群ごとに平均値・標準偏差・中央値を求める 平均値を計算するなら tapply(var1, grp, mean) を使用する。 grpのグループごとに、var1の平均値を... -
非劣性マージンの決め方と R で非劣性検定を実行する方法
いままでの方法と比べて、格段にいいとか、際立っていいとか、じゃなくてもいい場合がある。 ダメじゃなければいい。 劣っていなければいい。 劣っていなければいいという検定方法が、非劣性（ひれっせい）検定だ。 じゃ、どういうときが劣っていないって... -
R で出版バイアスの検定を行う方法
Begg検定、Egger検定、Macaskill検定という出版バイアスの検定方法の解説。 出版バイアス検定のための個々の研究データ 出版バイアスを検定するためのサンプルデータは以下の通り。 a <- c(3,7,5,102,28,4,98,60,25,138,64,45,9,57,25,65,17) n1 <- ...
