OK, I think this is done cooking for now. Added examples of fitting several other jolly-seber data augmentation models using Gibbs updates for latent states instead of updating them conditionally. Read to of readme and then scroll down to "alternative data augmentation approaches". So fast!
GitHub - benaug/Jolly-Seber-N-Prior-DA: Jolly-Seber MCMC samplers using N-prior data augmentation
Jolly-Seber MCMC samplers using N-prior data augmentation - benaug/Jolly-Seber-N-Prior-DA
github.com