pnme: A Gibbs sampler for a projected normal mixed-effects model
In bpnreg: Bayesian Projected Normal Regression Models for Circular Data

A Gibbs sampler for a projected normal mixed-effects model

pnme(
  theta_cos,
  theta_sin,
  X1,
  X2,
  Z1,
  Z2,
  ZtZ1,
  ZtZ2,
  R,
  pred,
  its,
  lag,
  burn,
  N
)

`theta_cos`	A List with the cosine of the circular dependent variable.
`theta_sin`	A List with the sine of the circular dependent variable.
`X1`	A list of fixed effect model matrices for component I.
`X2`	A list of fixed effect model matrices for component II.
`Z1`	A list of random effect model matrices for component I.
`Z2`	A list of random effect model matrices for component II.
`ZtZ1`	A list of transformed random effect model matrices for component I.
`ZtZ2`	A list of transformed random effect model matrices for component II.
`R`	A list of starting values for R.
`pred`	An empty list for likelihood computation.
`its`	An integer specifying the number of iterations
`lag`	An integer specifying the amount of lag.
`burn`	An integer specifying the number of burn-in iterations.
`N`	An integer specifying the number of burn-in iterations.