dualtraj: dualtraj
In jtm508/bayestraj:

Estimate the dual trajectory model using a Gibbs Sampler

1 2	dualtraj(X1, X2, y1, y2, K1, K2, z1, z2, 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
`K1:`	Integer, number of latent classes in series 1
`K2:`	Integer, number of latent classes in series 2
`z1:`	Matrix, K1 x dim(X1)[2] indicator matrix indicating which variables to include in each group.
`z2:`	Matrix, K2 x dim(X2)[2] indicator matrix indicating which variables to include 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.