求助小伙伴，在用R软件算核密度图时，如何将核密度图和正态曲线图画在一张图中？谢谢

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

求助小伙伴，在用R软件算核密度图时，如何将核密度图和正态曲线图画在一张图中？谢谢

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：核密度图 核密度 曲线图 r软件 小伙伴

1. x <- rnorm(100)
   
2. plot(density(x))
   
3. lines(density(rnorm(100)))

复制代码

1. curve(dnorm,-3,3)
   
2. x <- rnorm(100)
   
3. lines(density(x),col="red")

复制代码

crystal8832 发表于 2020-6-21 09:27

你好，plot(density(x))和lines(density(rnorm(100))有什么区别，不都是对随机生成的x
画密度图了，我试了一下，但出来的两个图不一样

maths_hjxk 发表于 2020-6-21 09:49

这个01是正态密度曲线，03是数据本身的核密度曲线吧？

北京姑娘 发表于 2020-6-21 11:05
这个01是正态密度曲线，03是数据本身的核密度曲线吧？

是的。01为理论曲线，03为估计的核密度曲线

maths_hjxk 发表于 2020-6-21 17:34
是的。01为理论曲线，03为估计的核密度曲线

嗯嗯 谢谢你

1. #ggplot版本
   
2. 
   
3. library(ggplot2)
   
4. 
   
5. df = data.frame(x = rnorm(1000, 2, 1))
   
6.   
   
7. 
   
8. ggplot(df, aes(x)) +
   
9.   geom_histogram(aes(y = stat(density)), colour='white',
   
10.                  fill='skyblue', bins = 50) +
    
11.   stat_function(
    
12.     fun = dnorm,
    
13.     args = list(mean = mean(df$x), sd = sd(df$x)),
    
14.     lwd = 2,
    
15.     col = 'red'
    
16.   )+
    
17.   theme_classic()
    
18. 
    
19. #lattice 版本
    
20. 
    
21. library(lattice)
    
22. 
    
23. histogram(~x, data = df, nint=50, type="density",
    
24.           panel = function(x,...){
    
25.             panel.histogram(x,
    
26.                             col='skyblue',
    
27.                             border='white',...)
    
28.             panel.mathdensity(dmath = dnorm,
    
29.                               col = "red",
    
30.                               lwd = 4,
    
31.                               args = list(mean=mean(x), sd=sd(x)))
    
32.           })

复制代码

crystal8832 发表于 2020-6-21 09:27

你好 还想请教你一下 如何利用pch在核密度线上添加小方块  小三角等 以区别不同的核密度线

1. x1<-seq(-3,3,by=0.1)
   
2. y<-dnorm(x1)
   
3. plot(x1,y,col=3,pch=3,lty=2,ylim=c(0,0.6))
   
4. x <- rnorm(100)
   
5. points(density(x),col=4,pch=2)
   
6. legend(0, 0.6, c("normal pdf","density estimation"), col = c(3,4), lty =
   
7. c(2,-1), pch = c(3, 2), merge = TRUE, bg='gray90')

复制代码

不知道是不是你要的。