相関係数のイメージ

2023.06.04

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

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をコピーしました