Displaying 1 result from an estimated 1 matches for "min_sets_index".
Did you mean:
min_sets_indexes
2011 Feb 08
3
Parallel processing question ...
...e
for (i in 1:n) {
nodes <- as.numeric(sets[indexes[i]:(indexes[i+1]-1)])
node <- resample(nodes, 1, prob=lam[nodes]/sum(lam[nodes]))
N[node] = N[node] + 1
}
N
}
lambdas<- foreach(j=1:(2*burn_in), .combine=cbind) %dopar% {
N <- update_N(min_sets, min_sets_indexes, lambda)
lambda <- rgamma(K, shape=a+N, rate=bT)
lambda
if (j%%100==0) { print(j); print(lambda); print(N)}
# if (j > burn_in) {
# lambdas <- cbind(lambdas, lambda)
# }
}
---------------snip