initializeSMARTchain[inputsnpls, nsq, nbp, initheta, inirho, inieps] generates initial values for one MCMC chain, and returns {{theta, epsilon, rho}, treels, adpaptive parameters, snpls}. In the last column of snpls, the log likelihood given the initialized treels over-writes the values in inputsnpls. |