Hamiltonian Monte Carlo in R