正态分布密度曲线如图1所示:x1=seq(-4,4,0.001),下面我们就来说一说关于r语言做正态分布直方图步骤?我们一起去了解并探讨一下这个问题吧!
![r语言做正态分布直方图步骤(用r语言模拟正态分布密度曲线和分位数函数)](http://img.studyofnet.com/upimg/630059608.jpg)
r语言做正态分布直方图步骤
正态分布密度曲线如图1所示:
x1=seq(-4,4,0.001)
x2=seq(-3,5,0.001)
y1=dnorm(x1,0,1)
y2=dnorm(x2,mean=1,sd=1)#y1,y2表示两条均值不同的正态密度函数曲线
plot(x1,y1,type="l")
lines(x2,y2,col="red")
abline(v=0)
abline(v=1,col='red')
#(a)下侧分位数如图二所示:
x=seq(-4,4,length=1000)
y=dnorm(x,0,1)
plot(x,y,type="l")
for(p in seq(0.01,0.95,length=50)){
p=0.95#把p的值放在里面等于没有循环
xp=qnorm(p)#xp:分位数
points(xp,0,col="red",pch=16,cex=1)
#cex:点的大小,pch=16:说明点是实心的,pch等于的数字不同,则对应的符合不同,16是一个实心点
x1=x[x<=xp]
y1=dnorm(x1)
lines(x1,y1,type="h",col="blue")
}
#(b)上侧分位数如图三所示:
x=seq(-4,4,length=1000)
y=dnorm(x,0,1)
plot(x,y,type="l")
p=0.95#如果p=0.5,则下面要改成lower.tail=T
xp=qnorm(p)#xp:分位数
points(xp,0,col="red",pch=16,cex=1)
x1=x[x>=xp]
y1=dnorm(x1)
lines(x1,y1,type="h",col="blue")
,