R/nextA.R

Defines functions nextA

nextA <- function(b,d,env)
{
     SQNP <- sqrt(env$gNP)
     
     # draws from multivariate normals can be taken using 
     # draws = mu + l * eta where l*l' = d (l is the cholesky decomposition)      
     
     #accounting for non-diffuse priors
     return(colMeans(b) + t(chol(d))%*%matrix(rnorm(env$gNIV),nrow=env$gNIV,ncol=1)/SQNP)
}

Try the RSGHB package in your browser

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

RSGHB documentation built on May 29, 2017, 6:23 p.m.