相関係数のイメージ

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

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))

リハビリテーション研究に必要な統計学について、R(Windows, ChatGPT)を使って紹介してます。サンプルは全て架空のデータで作成しています。したがって解析結果は事実とは異なりますのでご了承ください。間違いなどのご指摘はコメント欄にご記入いただければ助かります。

統計学備忘録をフォローする

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

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