2012-03-16から1日間の記事一覧

R R

頻度が指数分布する乱数の発生. rm(list=ls()) set.seed(1) N <- 1000 #発生させる個数 La <- 1/140 #指数分布のパラメータλ x <- round(rexp(N, rate=La),0) #Rの組み込み関数を使う x2 <- round(-(log(1-runif(N)))/La,0) #一様乱数から変換する mean(x) …