2012-01-15 49 views
3

我有一個密度函數f,並且我對它進行MCMC採樣。爲了評估抽樣的好處,我需要在同一個圖表中繪製histcurve。的在一個圖表上繪製直方圖和密度函數曲線

hist(samples); 
curve(dfun,add=TRUE); 

的問題是,它們在不同的規模:在一定區間的頻率通常是幾百個,而最大密度函數的爲約1左右。我想要做的是在同一高度上配置兩個圖,一個y軸在左邊,另一個在右邊。誰能幫忙?謝謝。

回答

7

使用prob=TRUE參數hist

hist(samples, prob=TRUE) 
curve(dfun,add=TRUE) 

另見this SO question

+0

謝謝@Dirk Eddelbuettel,它是完美的。還有一個問題。如果我的密度函數沒有正常化呢?選項'prob'沒有幫助。 (除了手動標準化) – ziyuang 2012-01-15 19:40:38

+1

我會說你必須手動標準化。 – 2012-01-15 20:00:52