R– category –
-
R で書いた散布図に相関係数を書き入れる方法
散布図内に相関係数と95%信頼区間を書き入れる方法 散布図を R の plot() で書く 以下のように Y 軸の変数をチルダ（ ~ ）の左側に書き、X 軸の変数を右側に書く。 plot(LogCD34perBW ~ WBC.entry, data=GCSF_CD34) すると、以下のように散布図が描かれる... -
R と EZR で同等性検定に必要なサンプルサイズを計算する方法
積極的に同等であることを証明していく同等性の検定。 サンプルサイズ計算はどのようにすればよいか？ 同等性検定はどのようにするのか？ 同等性の検定は、簡単に言えば、非劣性検定の群を入れ替えて、片側検定を2回行うことで実施できる。 劣っていないが... -
ANOVA Type I Type II Type III の違い
ANOVAには３つの種類がある。 Type I, II, IIIの３つ。 どんな時にどれを使えばよいか？ RにおけるANOVAの種類：Type I ANOVA Rのデフォルトで使えるANOVAは、anova()とaov()である。 これらはともにType I と呼ばれるANOVAである。 Type I は、複数の因子... -
R でブートストラップ ロジスティック回帰分析を行う方法
ブートストラップとは、小さいサンプルを復元抽出で何度もサンプリングして、疑似的に何度も試行したことにして、その結果から推定値の直接的な95％の範囲を求めたりする方法。 ブートストラップを使ったロジスティック回帰分析の方法をご紹介。 ブートス... -
多重代入法で必要な補完データセット数
欠測値があった場合、対処する方法がいくつかある。 そのうちの一つが多重代入法である。 多重代入法で作成すべき欠測値補完データの数はいくつがよいのか？ 多重代入法で作成するデータセットはいくつがよいと言われているか？ 多重代入法の提案者Rubinは... -
EZR でオッズ比の計算に必要なサンプルサイズを計算する方法
オッズ比のサンプルサイズ計算の方法 オッズ比計算のためのサンプルサイズ計算・検出力計算のための準備 EZRにepiR パッケージをインストールする。 Rスクリプト窓にinstall.packages("epiR")と書いて実行をクリック。 ダウンロード元を選ぶ（2024 年 10 ... -
R で分散の信頼区間を計算する方法
分散の信頼区間を R で計算する方法の解説。 分散の信頼区間を計算するスクリプト Rスクリプトは以下の通り。 var.interval <- function(data, conf.level=0.95){ df <- length(data) - 1 chilower <- qchisq((1-conf.level)/2, df) chiupper <... -
R で Box-Cox 変換を行う方法
連続量を何らかの方法で正規分布に近づける方法はいくつかある。 ここではBox-Cox変換の方法をまとめてみた。 R で Box-Cox変換を行う実例 まず、car パッケージを呼び出す。 library(car) car パッケージの中の、powerTransform()とbcpower()を使って変換... -
EZR のフィッシャーの正確検定で計算されるオッズ比の注意点
EZRのFisherの正確検定の関数でオッズ比を計算すると、通常の手計算で計算できるオッズ比と異なる数値が計算される。 これはなぜなのか？ EZRのFisherの正確検定を使わないオッズ比の計算 オッズ比は、通常は以下のような分割表だった場合、$ \frac{a/c}{b... -
R EZR エクセル SPSS で相関係数の計算・検定・必要なサンプル数を計算する方法
相関係数の計算に必要なサンプル数はいくつか？ 相関係数の検定との関係から紹介。 相関係数のサンプル数計算の前提となる検定の計算方法 相関係数の検定は、母集団の相関係数、母相関係数（ぼそうかんけいすう）がゼロかどうかの検定である。 通常はｔ分...
