Metropolis Hastings in R with specified multivariate prior and likelihood

Question

I need to run a metropolis hastings algorithm in R where the prior contains four variables and the likelihood function is somewhat complex. What is the best package/function for