Stable release for cran

Fixed documentation, and reverted back to using a positive constraint on the diagonal of Z. Without this, there's a very small number of cases where the parameter estimates would 'flip' within a chain, leading to identifiability and convergence issues