我有一个密度函数f，我为它进行MCMC采样。为了评估采样的优度，我需要在同一图表中绘制hist和curve。的问题
hist(samples);
curve(dfun,add=TRUE);
是因为它们在不同的尺度上：某个仓的频率通常为数百，而密度函数的最大值约为1左右。我想做的是在同一高度配置两个图，一个y轴在左边，另一个在右边。有人能帮忙吗？非常感谢。
将prob=TRUE参数用于hist:
hist(samples, prob=TRUE)
curve(dfun,add=TRUE)
另请参阅此SO问题