R/upPsi2T.R

upPsi2T<-function(psi2T, xiT, SI4SigmaT, cont, psi2){
	
	1 / (1 / psi2 + 2 * sum(1 / (1 + exp(-psi2T / 2) / (sqrt(det(SI4SigmaT)) * exp(xiT) * exp(-cont)))))

	}
igollini/lvm4net documentation built on June 20, 2019, 4:48 p.m.