Rで解析：分布コマンドの紹介

300回目の記事は悩みましたが、私的によく使う分布のコマンドです。Rでは分布コマンドを利用してランダムに値を得る時にコマンド前にrを加えます。

ぜひ、紹介するコマンドをヘルプで確認し基本コマンドに収録されている「各種分布コマンド」を確認してみてください。Rの習得の一歩になるかもしれません。

なお、beeswarmプロットに「vipor」パッケージを利用しています。

実行コマンドはR version 4.2.2で確認しています。

実行コマンド

詳細はコメント、コマンドのヘルプを確認してください。

#&#12497;&#12483;&#12465;&#12540;&#12472;&#12398;&#12452;&#12531;&#12473;&#12488;&#12540;&#12523;
install.packages("vipor")

#&#12497;&#12483;&#12465;&#12540;&#12472;&#12398;&#35501;&#12415;&#36796;&#12415;
library("vipor")

###&#12487;&#12540;&#12479;&#20363;&#12398;&#20316;&#25104;#####
n <- 1500
TestDF <- data.frame(Group = sample(paste0("Group", 1:2), n, replace = TRUE),
                     "&#12505;&#12540;&#12479;&#20998;&#24067;" = rbeta(n, 1, 1, ncp = 0),
                     "&#27491;&#35215;&#20998;&#24067;" = rnorm(n),
                     "&#12459;&#12452;&#20108;&#20055;&#20998;&#24067;" = rchisq(n, 1),
                     "&#12509;&#12450;&#12477;&#12531;&#20998;&#24067;" = rpois(n, 1),
                     "&#12525;&#12472;&#12473;&#12486;&#12451;&#12483;&#12463;&#20998;&#24067;" = rlogis(n, 0, 1),
                     "&#19968;&#27096;&#20998;&#24067;" = runif(n, 0, 1))
#####

par(mfrow = c(2,3))

for(n in 1:6){
  vpPlot(TestDF[, 1], TestDF[, n + 1], pch = 21, col = "#4b61ba",
         main = colnames(TestDF[n + 1]), alpha = .1)
}

出力結果

分布のプロットです。

少しでも、あなたの解析が楽になりますように！！