removed parameter draw loop inside irf
WIP: new function for posterior marginal and prior densities