dualtraj_constrained: dualtraj_consrained
In jtm508/bayestraj:

Estimate the dual trajectory model using a Gibbs Sampler

1 2	dualtraj_constrained(X1, X2, y1, y2, K, z, iterations, thin = 1, dispIter = 10, ll = FALSE)

`X1:`	Matrix, design matrix for series 1. 1st column should be the id.
`X2:`	Matrix, design matrix for series 2. 1st column should be the id.
`y1:`	Vector, outcomes for series 1
`y2:`	Vector, outcomes for series 2
`K:`	Integer, number of latent classes
`z:`	Matrix, K x dim(X)[2] indicator matrix indicating which variables to inlcude in each group.
`iterations:`	Integer, number of MCMC iterations
`thin:`	Integer, store every 'thin' iteration
`dispIter:`	Integer, frequency of printing the iteration number
`ll:`	Boolean, Set to TRUE to display the maximum log-likelihood over all the draws.