相関係数のイメージ

使用するパッケージ(パッケージのインストール

library(mvtnorm)
library(tidyverse)

相関係数 $0, 0.2, 0.5, 0.7, 0.8, 0.9$ を作図します

matrixで分散共分散行列を作成

rmvnorm関数でそれぞれの相関をもつ乱数を各300個作成

plot関数で作図

for関数で図の繰り返し

fun <- function(k){
    matrix(c(1, k, k, 1), ncol=2) %>%
    mvtnorm::rmvnorm(n=300, mean=c(0,0), sigma=) %>%
    plot( ,xlab =k, ylab="")  
    }

par(mfrow = c(2, 3))
for(k in c(0.0, 0.2, 0.5, 0.7, 0.8, 0.9)){fun(k)}
par(mfrow = c(1, 1))

ダメ出し 間違い、分かりにくい部分などのご意見をお待ちします

タイトルとURLをコピーしました