R随机数
常见分布表:
> n=10
> #生成n个均匀分布的数
> runif(n, min=0, max=1) #min默认0 max默认1
[1] 0.41915337 0.39372552 0.89188504 0.34365708 0.99997584
[6] 0.77761968 0.64008306 0.79776374 0.55633145 0.03431137
>
> hist(x,prob=T,col=gray(.9),main='unifrom on [0,1]') #对x绘制概率直方图 设置颜色 添加标题
> curve(dunif(x,0,1),add=T) #添加密度曲线
>
> x <- rnorm(100,0,1) #生成100个标准正太分布随机数
> hist(x,prob=T,col=gray(.9),main='standard normal') #绘制概率直方图
> curve(dnorm(x,0,1),add=T) #添加密度曲线
>
> par(mfrow=c(1,3))
p=0.25
for( n in c(10,20,50)){
x=rbinom(100,n,p)
hist(x,prob=T,main=paste("n =",n))
xvals=0:n
points(xvals,dbinom(xvals,n,p),type="h",lwd=3)
}
par(mfrow=c(1,1))
>x=rexp(100,1/10) # 生成100 个均值为10 的指数分布随机数
hist(x,prob=T,col=gray(0.9),main='exponential')
curve(dexp(x,1/10),add=T)