R/d12_Q.R

Defines functions `d12_Q`

`d12_Q` <-
function(n,X,d1_X_theta,d1_X_logb,sigma2omega,B){   
	d12_Q <- -n*sum(diag(solve(X) %*% d1_X_theta %*% solve(X) %*% d1_X_logb)) + (2/sigma2omega)*sum(diag(solve(X) %*% d1_X_theta
             %*% solve(X) %*% d1_X_logb %*%  solve(X) %*% B))
	return(d12_Q)
}

Try the Stem package in your browser

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

Stem documentation built on May 2, 2019, 8:56 a.m.