A Gibbs sampler for a projected normal regression model
pnr(theta, X1, X2, its, lag, burn)
A vector with the circular dependent variable.
A model matrix for component I.
A model matrix for component II.
An integer specifying the number of iterations
An integer specifying the amount of lag.
An integer specifying the number of burn-in iterations.