biQ_function function

Q-function to replace log-likelihood function

biQ_function(par, prob_log, omega_log, X, k, n1, n2, times1, times2)

Arguments

• par: numeric vector for parameters need to be estimated
• prob_log: mixture component weights(log)
• omega_log: latent variables(log)
• X: matrix for cluster
• k: vector for the cluster number
• n1: scalar for number of column contain first trait/location etc
• n2: scalar for number of column contain second trait/location etc
• times1: vector for the x values or time points
• times2: vector for the x values or time points

Returns

the Loglikelihood value