R中密度估计的等步长

时间:2019-04-03 02:19:13

标签: r probability-density

使两个分布不同,但要控制密度图的x值上的“步长”。目前,R可以提供最佳或随机的步长

尝试弄乱带宽并查看不同密度的帖子,但没有一个解决这个问题。 R帮助文档同样含糊不清。

这是我尝试过的一些代码。

data2 <- rlnormAlt(n = 100000, mean = 5.0, cv = 3.5/5.0)
data1 <- rlnormAlt(n = 100000, mean = 12.3, cv = 19.6/12.3)

x = density(data1, bw = 0.01, kernel = "gaussian", from= range(c(data1, data2))[1], to=range(c(data1, data2))[2], n = 10000)$x
y = density(data1, bw = 0.01, kernel = "gaussian", from= range(c(data1, data2))[1], to=range(c(data1, data2))[2], n = 10000)$y -
density(data2, bw = 0.01, kernel = "gaussian", from= range(c(data1, data2))[1], to=range(c(data1, data2))[2], n = 10000)$y

z <- cbind(x,y) %>% as.data.frame() %>% dplyr::filter(., y > 0, x > 4)
sapply(seq(1,10,1), function(i) z$x[i+1] - z$x[i])

您可以看到步长“有时”是一致的,但是我无法设置。

例如,我希望步长为每1或0.5等(即seq(0,50,1),我可以选择此步长)。

任何帮助将不胜感激。

0 个答案:

没有答案