「2012-06-12 RのLevy過程シミュレーション関数」で
プログラムというのは,規則性のある繰り返しをうまく書くべきものなのだから...
library(adehabitat)
# help(simm.levy)
system.time({
set.seed(411)
w <- simm.levy(1:500, mu = 1.5, burst = "mu = 1.5")
u <- simm.levy(1:500, mu = 2, burst = "mu = 2")
v <- simm.levy(1:500, mu = 2.5, burst = "mu = 2.5")
x <- simm.levy(1:500, mu = 3, burst = "mu = 3")
par(mfrow=c(2,2))
lapply(list(w,u,v,x), plot, perani=FALSE)
})
は,もし,パラメータを変えてもっとたくさんのシミュレーションをするときのことを考えれば,以下のように書く方がよいのは自明。
system.time({
set.seed(411)
par(mfrow=c(2,2))
mu <- seq(1.5, 3, by=0.5)
sapply(mu, function(m) plot(simm.levy(1:500, mu=m, burst=sprintf("mu = %f", m)), perani=FALSE))
})
※コメント投稿者のブログIDはブログ作成者のみに通知されます