2011-07-01から1日間の記事一覧

Rの確率分布関数の"rxxxx()"を使ってみる

R

# 一様分布 help(runif) args(runif) n<-10000 R<-runif(N) hist(R) summary(R) 結果をWikiの定義と比較してみる # ポアッソン分布 help(rpois) args(rpois) n lambda R hist(R) summary(R) |

ベクトル→行列→アレイ

R

ベクトルは1次元 行列は2次元 アレイはさらに高次元 dims<-c(3,4,5,6) A<-array(1:prod(dims),dims) apply()関数は次元に応じて使える apply(A,1,sum) apply(A,2,sum) apply(A,3,sum) apply(A,c(1,2),sum) apply(A,c(2,3),sum) apply(A,c(3,1),sum) apply(…

行列データを図示する

R

2つの場合 列ベクトルの束として見る場合 2次元座標に値が入っていると見る場合 列ベクトルの束として見る場合 対応のある値の列ベクトルのペア N1<-100 V1<-runif(N1) V2<-V1+rnorm(N1)*0.1 M<-cbind(V1,V2) plot(M) 対応のある値の列ベクトルのトリオ N1…