R中的修正指数分布

Question

我想在R中创建自己的概率密度函数，以模拟论文中的某些内容。

它在某种程度上类似于指数分布，但我真正想做的是将指数分布重新定义为“修改的”分布。

有这种方法吗？ 谢谢。

---

我想模拟一下：

b（x）= {（µ / p）（e ^（-（µx-q）/ p），x> q（xbar），否则为0}

xbar是一个x，上面有一条线，表示平均值

Answer 1

您是否知道从下面截断的指数分布仍然是指数分布？

Answer 2

我相信您可以使用?sample来做您想做的事情。 使用已知的分布函数b(x)生成概率向量，例如bprob ，然后

sample(x,size, replace=TRUE,prob=bprob)

有一些非常有趣的方法可以从任意分布中生成样本。 参见，例如，正常随机数：使用机器分析选择最佳算法WH PAYNE华盛顿州立大学（网上某处），或C中的数字食谱《科学计算的艺术》第二版William H. Press，Saul A. Teukolsky等，第7.3章