R/Likelihood_thetas.R

Defines functions Likelihood_thetas

Documented in Likelihood_thetas

Likelihood_thetas <-
function(thetas,covmat) {
		cholcov <- chol(covmat)
		invcholcov <- MASS::ginv(cholcov)
		logsqrtdet <- sum(log(diag(cholcov)))
		-(1/2)*(colSums((crossprod(invcholcov,thetas))^2))-logsqrtdet
	}

Try the BEDASSLE package in your browser

Any scripts or data that you put into this service are public.

BEDASSLE documentation built on April 11, 2022, 1:07 a.m.