2012-03-16 ■ 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) # 144.383 mean(x2) # 142.306