統計ER

R, EZR, SPSS, KH Coder を使ったデータ分析方法を紹介するブログ。ニッチな内容が多め

EZR でスピアマンの順位相関係数の 95%信頼区間を計算する方法

EZR で Spearman の順位相関係数の 95%信頼区間を計算する方法

デフォルトでは計算できないが、パッケージをインストールすると計算できる

>>もう統計で悩むのを終わりにしませんか?


↑1万人以上の医療従事者が購読中

EZR で Spearman の順位相関係数を計算する方法

統計解析 → ノンパラメトリック検定 → 相関係数の検定(Spearman の順位相関係数の検定)を選択する

Ctrl キーを押しながら、二つの変数を選択する

OK をクリックすると、以下のように結果が表示される

Spearman の順位相関係数の点推定値及び P 値は計算されるが、95 %信頼区間は計算されない

95 % 信頼区間を計算するにはどうしたらよいか

EZR で Spearman の順位相関係数の 95 %信頼区間を計算する方法

RVAideMemoire パッケージの spearman.ci 関数を使うと計算できる

まず、RVAideMemoire パッケージをインストールする

以下の一行を R スクリプト枠に手入力、もしくは、コピペして、実行ボタンをクリックする

install.packages("RVAideMemoire")

出てきた窓内から、Japan を探して選択して、OK をクリック

インストールが終了する

次に、パッケージを呼び出し、乱数シードを設定し、信頼区間を計算する

library(RVAideMemoire) # パッケージ呼び出し

set.seed(1) # 乱数シードの設定:再現性の確保のため

spearman.ci(bfi$A1, bfi$A2) # 順位相関係数を求める二つの変数を選択

bfi$A1, bfi$A2 と書いてある、変数名の前の bfi$ は、データセット名と半角ドルマークである

これは、bfi というデータセット内の A1 とか A2 とかを指定するための記述である

ご自分のデータセットで計算するときは、「ご自分のデータセット名$変数名」 に変更してほしい

上記を 1 行ずつ、カーソルをおいて、実行すると、以下の結果が出てくる

Spearman の順位相関係数は、EZR のメニューから計算した結果と同じく、約 -0.371 である

95 % 信頼区間は、小数点以下 4 桁目で四捨五入すると -0.403 ~ -0.335 と計算された

>>もう統計で悩むのを終わりにしませんか?


↑1万人以上の医療従事者が購読中

まとめ

EZR で Spearman の順位相関係数の 95 % 信頼区間を計算する方法を紹介した

RVAideMemoire パッケージの spearman.ci 関数を使うと計算できる

参考になれば