2011-10-23 31 views
0

我想在R中用我自己的概率密度函數來模擬紙上的某些內容。R中的修改的指數分佈

它是某種類似於指數分佈,但我真正想要做的是重新定義指數分佈成「修飾」一個......

是否有這樣的方式來做到這一點?謝謝。


我想模擬這樣的:

B(X)= {(μ/ P)(E ^( - (μX-Q)/ P)中,x> Q(XBAR)和0否則}

XBAR是它上面,意味着線的X,平均

+2

確實存在:定義一個名爲'rmodexp'(或某些這樣的)的函數。然而,隨機偏差的產生並不總是很容易 - 對分佈進行相當小的修改會使模擬或模擬效率變得更加困難。如果你確切地告訴我們你想要的是什麼(例如,你在修改後的函數的密度函數或累積密度函數是什麼?)有人可能會幫助... –

+1

這會更適合於http://stats.stackexchange。 COM? –

+0

我沒有時間完全解決這個問題,但這看起來像是縮放和移位的指數。我不完全理解這個符號,但是......我懷疑'mu/p * rexp(n)+ q/p'會讓你走到最後。你可以參考一下紙張本身嗎? –

回答

0

您是否知道從下面截斷指數分佈保持指數分佈?

0

我相信你可以做你想要與?sample什麼我們你e已知分佈函數b(x)生成概率的向量,說bprob,然後

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

有從任意分佈產生的樣本一些非常有趣的方法。例如,請參閱:普通隨機數字:使用機器分析來選擇最佳算法 WH PAYNE WashlngtonState University(網上某處)或數字食譜C 科學計算的藝術第二版 William H. Press,Saul A. Teukolsky等人,章節7.3