R语言均匀分布

最新推荐文章于 2025-04-12 01:00:00 发布

眼里进沙子好疼

最新推荐文章于 2025-04-12 01:00:00 发布

阅读量8.4k

点赞数 1

分类专栏：笔记

本文链接：https://blog.youkuaiyun.com/weixin_44850594/article/details/104988625

版权

R语言均匀分布

问题：

设总体X~U(-2,2)，即总体X服从[-2,2]上的均匀分布。设Y_1=X_1，Y_2=X_1+X_2，Y_3=X_1+X_2+X_3，Y_4=∑_(i=1)^10X_i ，Y_5=∑_(i=1)^20X_i 。独立重复上述试验1000次。试画出Y_1,⋯,Y_5的直方图、分布的核密度曲线、相应的正态分布密度函数。

##代码

n<-c(1,2,3,10,20)
n
N<-1000
a<--2:2
a
data<-list()
data
y<-matrix(rep(0,N*length(n)),byrow=T,ncol=N)#产生一个N*n长度的全0矩阵
y
for(i in 1:length(n)){
   
  data[[i]]<-matrix(runif(i*N),ncol=i,byrow=T,nrow=N)
  y[i,]<-(apply(data[[i]],1,sum)-i*mean)/(sqrt(i)*sd)
}

#Y1=x1
hist(y[1,],probability=T,xlab=expression(paste(Y[1],"=",x[1])),main="",ylim=c(0,0.5),ylab="密度",cex=1.2,lwd=5,col=c("dark blue","dark red","dark orange","green","purple"