R/samplexitrue.R

In mvst: Bayesian Inference for the Multivariate Skew-t Model

samplexitrue = function(y, X, N, particles, priorList){
 p = ncol(y)
 #
 theta = dget('Output/Samples/theta.txt')
 xiTrue = theta$xi
 xi = matrix(rep(xiTrue, each=N), ncol=p, byrow=F)
 log.dxi = rep(0, N)
 #
 return(list(values=xi, log.dq=log.dxi))
}