[英]Modified exponential distribution in R
我想在R中創建自己的概率密度函數,以模擬論文中的某些內容。
它在某種程度上類似於指數分布,但我真正想做的是將指數分布重新定義為“修改的”分布。
有這種方法嗎? 謝謝。
我想模擬一下:
b(x)= {(µ / p)(e ^(-(µx-q)/ p),x> q(xbar),否則為0}
xbar是一個x,上面有一條線,表示平均值
您是否知道從下面截斷的指數分布仍然是指數分布?
我相信您可以使用?sample
來做您想做的事情。 使用已知的分布函數b(x)
生成概率向量,例如bprob
,然后
sample(x,size, replace=TRUE,prob=bprob)
有一些非常有趣的方法可以從任意分布中生成樣本。 參見,例如,正常隨機數:使用機器分析選擇最佳算法WH PAYNE華盛頓州立大學(網上某處),或C中的數字食譜《科學計算的藝術》第二版William H. Press,Saul A. Teukolsky等,第7.3章
聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.