Compute the Likelihood of the PN distribution (regression)
lik_reg(X1, X2, theta, beta1, beta2, n)
the model matrix of the first component
the model matrix of the second component
a circular outcome value
estimated linear coefficients of the first component
estimated linear coefficients of the second component
sample size