模拟从$X$分配中提取(R)

模拟从$X$分配中提取(R),r,R,我有一个随机变量$x$的pdf$f(x)=4x^3$,其中我需要从分布中模拟抽取 我的解决方案包括从pdf中查找cdf(第1期): 既然您已经将您的pdf标识为beta版,只需使用rbeta进行采样 s1 <- rbeta(5000,4,1) s1既然您已经将您的pdf识别为beta版,只需使用rbeta进行采样 s1 <- rbeta(5000,4,1) s1既然您已将您的pdf识别为beta版,只需使用rbeta进行采样即可。完美!所以rbeta(1,shape1=4,sha

我有一个随机变量$x$的pdf$f(x)=4x^3$,其中我需要从分布中模拟抽取

我的解决方案包括从pdf中查找cdf(第1期):


既然您已经将您的
pdf
标识为beta版,只需使用
rbeta
进行采样

s1 <- rbeta(5000,4,1)

s1既然您已经将您的
pdf
识别为beta版,只需使用
rbeta
进行采样

s1 <- rbeta(5000,4,1)

s1既然您已将您的
pdf
识别为beta版,只需使用
rbeta
进行采样即可。完美!所以
rbeta(1,shape1=4,shape2=1)
?如果给定的pdf不遵循已知的分发。如何从随机变量(在R中)模拟绘图?请参阅帮助文件
?rbeta
,并确认您使用的是相同的参数化。对于任意pdf,如果您无法解析求解,则始终可以使用拒绝采样。因为您已将您的
pdf
识别为beta版,只需使用
rbeta
进行采样。完美!所以
rbeta(1,shape1=4,shape2=1)
?如果给定的pdf不遵循已知的分发。如何从随机变量(在R中)模拟绘图?请参阅帮助文件
?rbeta
,并确认您使用的是相同的参数化。对于任意pdf,如果无法解析求解,则始终可以使用拒绝采样。
pdf <- function(x) 4*x^3 # on [0,1]
p <- runif(50000)
dp <- pdf(p)
s2 <- p[runif(50000) < dp/max(dp)]