R/a0_gibbs_rate.R

Defines functions a0_gibbs_rate

Documented in a0_gibbs_rate

a0_gibbs_rate <-
function(thetas,covmat,a0){
		cholcov <- chol(covmat)
		invcholcov <- MASS::ginv(cholcov)
		tmp <- (1/2)*(colSums((crossprod(invcholcov,thetas))^2))
		return(sum(tmp)/a0)
	}

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.