模拟产生2组(每组1000个)独立的[0,1]上均匀分布的随机数,求和并画出数据的直方图,观察结果是否仍是均匀分布?产生12组求和的结果?
作者:互联网
本题验证均匀分布是否具有可加性和模拟中心极限定理。利用函数runif()产生随机数,其中min=0, max=1表示均匀分布的区间。其模拟1000个均匀分布随机数如图9所示,2组均匀分布随机数相加如图10所示。从图10可看出,2个均匀分布相加不为均匀分布,其为三角分布。下面给出相应的代码:
n <-1000
a <- runif(n, min=0, max=1)
b <- runif(n, min=0, max=1)
par(mfrow=c(2,1))
hist(a,prob=T,main="均匀分布")
hist(a+b,prob=T,main="均匀分布是否具有可加性")
m <- 12
n <- 1000
data <- matrix(nrow = n,ncol = m)
for(i in 1:m){
data[,i] = runif(n, min=0, max=1)
}
sumdata <- rowSums(data)
hist(sumdata,prob=T,main="12组均匀分布求和的分布")
标签:10,求和,均匀分布,直方图,随机数,所示,模拟,1000 来源: https://blog.csdn.net/xj4math/article/details/114766882