如何绘制函数的概率密度函数？

Question

假设A遵循指数分布; B遵循Gamma分布如何绘制0.5*(A + B)的PDF

Answer 1

使用"distr"包非常简单:

library(distr)

A <- Exp(rate=3)
B <- Gammad(shape=2, scale=3)

conv <- 0.5*(A+B)

plot(conv)
plot(conv, to.draw.arg=1)

由JD Long编辑

结果图如下所示:

Answer 2

如果您只是在寻找快速图表,我通常会采用快速而肮脏的模拟方法.我做了一些绘画,在绘制时猛击高斯密度并绘制那个坏男孩的情节:

numDraws   <- 1e6
gammaDraws <- rgamma(numDraws, 2)
expDraws   <- rexp(numDraws)
combined   <- .5 * (gammaDraws + expDraws)
plot(density(combined))

输出应该看起来像这样: