目的よく教科書に載っている「相関係数が\(\rho\)のときの散布図」をRで作る。*この記事はこちらのウェブサイトに置かれています。Rのコードに関して前嶋直樹さん(東京大学大学院)よりアドバイスをいただきました。記して感謝いたします。パッケージの準備#以下のパッケージを使用する。入っていない場合はインストール。#別解のやり方をするときに使用。手順2次元正規乱数を発生させる関数を定義する(http://cse.naro.affrc.go.jp/takezawa/r-tips/r/60.html)。return(data.frame(x=x,y=y))}この関数を用いて、相関係数-1から+1まで合計21個のサンプルを生成する。ここではそれぞれサンプルサイズは5000とした。cordata_mini)}なおこの部分はpurrrというパッケージを使って書くこともできる。今回はコメントアウトしておく。magrittr::set_colnames(c("x","y","cor"))生成したサンプルの相関係数が真の値に近いかどうかを確認しておく。生成したサンプルのサイズが大きいほど、真の値に近い値を得られる。今回はそれぞれ5000ケースとたくさんのケースを抽出したので、いずれもそれなりに真の値に近い値が得られると思われる。1.0000000サンプルサイズ(size)が大きければ大きいほど生成したサンプルから得られる相関係数は真の値に近づくが、その分あとで作成する散布図がつぶれてしまうので難しいところ。相関係数の値ごとに、散布図をプロット。相関係数-1の場合は自明なので省略。(geom_point).いくつかピックアップして抜き出すと次のようになる。(geom_point).。