Programming an HMM with R