R在Mathematica中有类似于TransformedDistribution的东西吗？

Question

我有一个随机变量X和一个转换f，我想知道f(X)的概率分布函数，至少近似。在Mathematica中有TransformedDistribution，但我在R中找不到类似的东西。正如我所说，某种近似解也可以。

Answer 1

您可以查看distr包裹。例如，假设为y = x^2+2x+1，其中x通常以均值2和标准差5分布。您可以：

require(distr)
x<-Norm(2,5)
y<-x^2+2*x+1
#y@r gives random samples. We make an histogram.
hist(y@r(10000))
#y@d and y@p are the density and the cumulative functions
y@d(80)
#[1] 0.002452403
y@p(80)
#[1] 0.8891796