Nothing
logPriorDensT = function(N, p, particles, priorList){
nu = particles$nu
particlesStar = reparamT(N=N, p=p, particles=particles)
#
gdls.list = priorList$gdls
# nuIndices = numeric(N)
# for(iN in 1:N){
# nuIndices[iN] = which(gdls.list == nu[iN])
# }
nuIndices = match(nu, gdls.list)
logPriorDens = # 0 + # xi
# - nSphereVolume(p, 1, LOG=T) - 0.5 * particlesStar$log.detOmega + # delta
(-(p+1)/2) * particlesStar$log.detSigma + # Sigma
# (-2) * apply(particlesStar$h,1,sum) + # |J1|
# (-1) * apply(log(particlesStar$h), 1, sum) + # |J|
priorList$nulogpriors[nuIndices] # nu
logPriorDens[which(particlesStar$pos == 0)] = rep(-Inf, sum(particlesStar$pos == 0))
return(logPriorDens)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.