2016/10/04

2020/05/04

R言語で回帰分析の重相関係数Rの信頼区間を出力

ライター：IMIN

今回は、R言語で関数lmを使った回帰分析をした際に、計算される重相関係数Rの95%信頼区間を出力する方法をご紹介いたします。R言語でRの信頼区間を出力……。なんともややこしいですね。

まずは何らかのデータを回帰分析して、ansという変数に代入します。

プログラム

ans <- lm(data$y~.,data=data)

この回帰分析の結果による予測値というのは”predict()“という関数を使えば、求めることができます。予測値を算出して、preという変数に代入しておきます。

プログラム

pre <-　predict(ans)

関数”cor.test()”を使う。

R言語で用意されている、関数”cor.test()”は、二つのデータの間に相関があるかどうか、検定する関数です。その結果として、相関係数の信頼区間も出してくれるという非常に便利なものです。

プログラム

cor.test(x,y)

とすることで、xとyの無相関検定（pearson）を行ってくれます。

今回の場合は、回帰分析による予測値と、実際のデータの実現値を引数に入れてあげれば、オーケーです。ということで、以下のようなプログラムになります。結果をcoという変数に代入しておきます。※ここで、”data$y”というのは、目的変数の実現値です。

プログラム

co <- cor.test(pre,data$y)

この中から95%信頼区間だけを抽出したい場合、以下のように記述してあげましょう。

プログラム

co95 <- co$conf.int

ちなみに、cor.test()はmethodで検定手法を変えることもできます。ピアソンの他に、ケンドール、スピアマンの順位相関係数の検定を指定することができます。

具体的には、

cor.test(x,y, method = “pearson”) …… ピアソン

cor.test(x,y, method = “kendall”) …… ケンドール

cor.test(x,y, method = “spearman”) …… スピアマン

というように記述すると、その通りの検定が行われます。

これで、回帰分析の重相関係数の95%信頼区間が出せましたので、今回はこの辺で失礼します！

※R言語入門のトップページはこちら

(totalcount 274 回, dailycount 34回 , overallcount 1,289,012 回)

ライター：IMIN

最新投稿記事

【独占】コインバンクの卜部代表、既存案件への影響限定的ー新型コロナで 2020年7月3日
【独占】データサイエンスコンペ参加者、３年後に１万人超目指すーNishikaの松田CTOインタビュー 2020年7月2日
【独占】純真学園大学・村中光副学長「G検定、全受講者に取得してもらいたい」ー「人工知能学」新設 2020年7月1日
OKI、９月下旬に「AIエッジ・コンテスト」ー「AE2100」で課題解決力競う 2020年6月30日
G検定に数学力はどれくらい必要？文系でも受かるのか？ 2020年6月26日
G検定に落ちた人、合格した人。勉強法の違いはどこにある？ 2020年6月25日
G検定の難易度・合格ライン・内容・評判を徹底調査！ 2020年6月19日
コニカミノルタ、パートナーと画像IoT開発推進、実践でAI人材育成ー高槻サイトで 2020年6月19日
「第１回全国高専ディープラーニングコンテスト」本選を８月22日に開催ーJDLA 2020年6月18日
G検定、94%が「誰かに勧めたい」―JDLA、1400人アンケート調査 2020年6月16日

週間ランキング

この記事を読んだ方はこちらの記事も読まれています

2016年7月22日ファイルからのデータ読み込みとアクセス【第2回】
2016年11月2日 Rでベクトルにその要素が無いなら追加する関数
2016年7月22日 Rでのベクトル・行列の作成と四則演算・要素の参照【第1回】
2016年7月22日実行するだけでR言語入門出来るプログラム
2016年7月27日 R言語で線形モデルによる回帰分析の自作関数例
2016年7月24日 Rでt検定行う方法とそのまま使える自作関数
2016年7月22日 Rでグラフをプロットして保存する【第5回】
2016年7月22日 R言語での自作関数の作り方・使い方【第4回】
2016年7月23日 Rで相関係数を求める方法とそのまま使える自作関数例
2016年7月22日 Rでロジスティック回帰分析そのまま使える自作関数例
2016年8月17日 Rでステップワイズ法による重回帰分析の自作関数例
2016年10月24日 Rで箱ひげ図を描く（プロットする）方法【boxplot】
2017年5月15日 RでYouden Indexを用いたROC曲線の最適カットオフ
2016年11月17日 Mac版のRでマルチバイトエラーが出た時は…
2016年7月22日 R言語の便利なコマンド集【第2.5回】
2017年2月6日 Rで散布図をプロットしてpngで保存する自作関数

COMMENT
コメントを残す
コメントをキャンセル