谁帮我看下这段代码，神没意思，怎么就能做出这个图的？？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

##绘制密度曲线

x <- seq(from = -3,to = 3,by = 1/20)
#标准正态分布密度函数
y <- dnorm(x)

plot(x,y,main = "Standard Normal Distribution",type = "l",
     ylab = "Density",xlab = "Quantile")

abline(h = 0)

region.x <- x[1 <= x & x <= 2]
region.y <- y[1 <= x & x <= 2]

region.x <- c(region.x[1],region.x,tail(region.x,1))
region.y <- c( 0 , region.y , 0)

polygon(region.x,region.y,dendity = 10)

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：没意思

不多说了，见代码，比较出来的结果

1. ##绘制密度曲线
   
2. 
   
3. x <- seq(from = -3,to = 3,by = 1/20)
   
4. #标准正态分布密度函数
   
5. y <- dnorm(x)
   
6. # 作一个标准正态图，type="l"是指作线图
   
7. plot(x,y,main = "Standard Normal Distribution",type = "l",
   
8.      ylab = "Density",xlab = "Quantile")
   
9. # 在plot基础上加一条y=0的水平线，h是代表horizon，是水平的意思
   
10. # 如果加垂直线，就是v=...,verticle的意思
    
11. abline(h = 0)
    
12. # 这里是确定你要选择的区域，首先选择x在1到2之间的
    
13. region.x <- x[1 <= x & x <= 2]
    
14. # 这是：当x在1到2之间时，y的范围
    
15. region.y <- y[1 <= x & x <= 2]
    
16. #我加了红线，你可以看出来，是哪一块区域
    
17. lines(region.x,region.y,col="red",lwd=2)
    
18. # 对比阶段：
    
19. polygon(region.x,region.y,density = 10)
    
20. # 如果没有这两句，就显示的是，上半部分
    
21. # 加入后，实际上，在原来的基础上，加了两个点：
    
22. # (region.x[1],0)和(tail(region.x,1),0)多了这两个点
    
23. # 其实这两个点在y=0这条线上
    
24. region.x1 <- c(region.x[1],region.x,tail(region.x,1))
    
25. region.y1 <- c( 0 , region.y , 0)
    
26. # 这个，就是从包含了红线，到y=0上的两个点，形成的一个多边形
    
27. polygon(region.x1,region.y1,density = 10)

复制代码

jiangbeilu 发表于 2015-12-9 10:27
不多说了，见代码，比较出来的结果

太谢谢你们了

ryoeng 发表于 2015-12-9 02:27
可以参考一下geom_polygon in ggplot2。

太谢谢你们了！喜欢这个论坛，我也希望能帮别人解决问题，现在技术不行

1. library(ggplot2)
   
2. x <- seq(from = -3,to = 3,by = 1/20)
   
3. y <- dnorm(x)
   
4. region.x <- c(1,x[1 <= x & x <= 2],2)
   
5. region.y <- c(0,y[1 <= x & x <= 2],0)
   
6. plot_df<-data.frame(x=x,y=y)
   
7. region_df<-data.frame(region.x=region.x,region.y=region.y)
   
8. p1<-qplot(x,y,data=plot1_df,geom="line")+geom_hline(y=0)+theme_bw()
   
9. p1+geom_polygon(aes(x=region.x,y=region.y),data=region_df)

复制代码

这一题的ggplot做法，但怎样做density那种阴影效果？请论坛大神指教。