我想在R中创建自己的概率密度函数来模拟论文中的某些内容。
它在某种程度上类似于指数分布,但我真正想做的是将指数分布重新定义为“修改过的”...
有这样的方法吗?感谢。
我想模拟这个:
b(x)= {(μ/ p)(e ^( - (μx-q)/ p),x> q(xbar),否则为0}
xbar是一个x,其上方有一条线,平均值
答案 0 :(得分:0)
您是否知道从下面截断的指数分布仍然是指数分布?
答案 1 :(得分:0)
我相信你可以用?sample
做你想做的事。使用已知的分布函数b(x)
生成概率向量,例如bprob
,然后
sample(x,size, replace=TRUE,prob=bprob)
有一些非常有趣的方法可以从任意分布中生成样本。例如,参见Normal RandomNumbers:UsingMachine AnalysisTo Choosethe BestAlgorithm W. H. PAYNE Washlngton州立大学(网上某处),或C中的数字食谱 科学计算第二版的艺术 William H. Press,Saul A. Teukolsky等,第7.3章