ffun(y,Rho) init.all(nval,K,rand.start,mixture=F) recurse(fy,tpm,nc,epsilon) revise.ispd(tpm) revise.rho(y,gamma,m) revise.tpm(xi,mixture)