search for: burn_in

...# count of events N <- rep(0, K) # count of failures 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) # } } ---------------sn...

...t(rsample, method="T10") #global test, subgroup inv and receive error message as follows. "Error in m[idx1, idx2] <- 1 : subscript out of bounds" Using the same code for a different dataset does not cause this problem. Any suggestions much appreciated! (ctr spec: burn_in = 100 n_eff = 1000 step = 16 seed = 0 tfixed = FALSE ) -- View this message in context: http://r.789695.n4.nabble.com/eRm-raschsampler-error-message-tp3770265p3770265.html Sent from the R help mailing list archive at Nabble.com.