R肯定可以用来实现模拟,只要知道随即模拟的原理就可以了,配合相应的R语法,编出程序应该没问题。举个简单例子,用R模拟中心极限定理:
<br />
for (i in 1:200) {<br />
x = replicate(200, mean(runif(500))) <br />
plot(density((x)), main = "Demonstration of CLT", <br />
ylab = expression("Density Estimation of" ~ ~italic(bar(x))), <br />
col = "blue")<br />
}<br />
意思就是重复200次产生500个均匀分布的独立随机变量值均值,看这200个均值的分布。