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