Построить гистограмму и кривую функции плотности на одном графике - PullRequest
3 голосов
/ 15 января 2012

У меня есть функция плотности f, и я выполняю выборку MCMC для нее.Чтобы оценить качество выборки, мне нужно нанести hist и curve на одном графике.Проблема

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

заключается в том, что они имеют разный масштаб: частота определенного бина обычно равна сотням, а максимум функции плотности составляет около 1 или около того.То, что я хочу сделать, это настроить два графика на одной высоте, с одной осью Y слева, а другой справа.Кто-нибудь может помочь?Спасибо.

1 Ответ

7 голосов
/ 15 января 2012

Используйте аргумент prob=TRUE для hist:

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

Также см. этот вопрос SO

...