1. (important) change prior mupv = 0.1 from 1; to make it spike at 0
2. (important) correct domega for general K
3. (important) correct gamma prior for sigma (scale instead of rate)
4. add mu, sigma, omega save and plot
1. add multinomial random draw
1. normalized prob in rcat function
2. updateg function in fortran
3. replace NA with 0 in R function to pass to fortran
4. use updateg1g2 fortran function in R, decrease 1/3 the time