请问如何实现在正态分布密度曲线下面涂对应置信区间的颜色

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

#区间估计与假设检验之间的联系
i=0
while (i<1.5) {
  x <- seq(-3,3,length = 5000)
  y <- dnorm(x,0,1)
  a<-plot(x,y,col="red",xlim = c(-3,3),ylim = c(0,0.5),type = 'l', ylab =
            'denisity',main='区间估计与假设检验之间的联系',sub="蓝色区域为置信区间，两侧空白区域为假设检验中双侧检验的拒绝域")
  if (i<1.5){
  q=seq(0,i,length=floor(80*i))
  abline(v=c(q,-q),col="blue")
  text(-i-0.2,0.15,"→"); text(i+0.2,0.15,"←")
  text(0,0.1,"置信区间")
  text(-2.5,0.2,"P=α");text(-1.7,0.05,"α/2");text(1.7,0.05,"α/2")
  i=i+0.05
  Sys.sleep(.1)}
}

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏2 回帖

关键词：之间的联系 正态分布 区间估计 假设检验

1. x=seq(-3,3,0.001)
   
2. y1=dnorm(x,0,1)
   
3. x1=seq(-2,2,0.001)
   
4. y2=dnorm(x1,0,1)
   
5. plot(x,y1,type='l')
   
6. polygon(x=c(x,rev(x)),y=c(rep(0,length(x)),rev(y1)),col='skyblue')
   
7. polygon(x=c(x1,rev(x1)),y=c(rep(0,length(x1)),rev(y2)),col='yellow')

复制代码

1. #区间估计与假设检验之间的联系
   
2. i=0.01
   
3. while (i<1.5) {
   
4.   x <- seq(-3,3,length = 5000)
   
5.   y <- dnorm(x,0,1)
   
6.   a<-plot(x,y,col="red",xlim = c(-3,3),ylim = c(0,0.5),type = 'l', ylab =
   
7.             'denisity',main='区间估计与假设检验之间的联系',sub="蓝色区域为置信区间，两侧空白区域为假设检验中双侧检验的拒绝域")
   
8.   if (i<1.5){
   
9.   q=seq(0,i,length=floor(i*40)+1)
   
10.   yq=dnorm(q,0,1)
    
11.   for(j in 1:length(q)){
    
12.   segments(q[j],0,q[j],yq[j],col='skyblue')
    
13.   segments(-q[j],0,-q[j],yq[j],col='skyblue')
    
14.   }
    
15.   #abline(v=c(q,-q),col="skyblue")
    
16.   text(-i-0.2,0.15,"→"); text(i+0.2,0.15,"←")
    
17.   text(0,0.1,"置信区间")
    
18.   text(-2.5,0.2,"P=α");text(-1.7,0.05,"α/2");text(1.7,0.05,"α/2")
    
19.   i=i+0.05
    
20.   Sys.sleep(.9)}
    
21. }

复制代码

zhou1_20 发表于 2018-12-27 15:14

能详细说说吗？

huapeng66 发表于 2018-12-28 07:16
能详细说说吗？

看我发的第一代码块，体会下R是如何绘制区域图的就好了

zhou1_20 发表于 2018-12-27 15:14

相当不错的代码！