R/samplepsitrue.R

samplepsitrue = function(y, X, N, particles, priorList){
 p = ncol(y)
 #
 theta = dget('Output/Samples/theta.txt')
 psiTrue = theta$psi
 psi = matrix(rep(psiTrue, each=N), ncol=p, byrow=F)
 log.dpsi = rep(0, N)

 return(list(values=psi, log.dq=log.dpsi))
}

Try the mvst package in your browser

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

mvst documentation built on May 2, 2019, 1:46 p.m.