search for: min_set

...es per node 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