How to pass control to previous iteration in for loop in R -


i running loop in r (as part of power analysis model ran r2jags). @ point want know if mcmc chains have converged, if not want skip iteration of loop. however, don't want skip next iteration, want loop start same iteration again. using command 'next', skipping iteration. how tell loop iteration? below whole code, little bit worried about:

  if(up1 > 1.1 | up2 > 1.1 | up3 > 1.1 | up4 > 1.1)    next

and whole code:

g.vec <- rep(na, 1000) #-----model------  model1 <- function(){      # likelihood   (c in 1:16){     d.diff[c] ~ dnorm(mu, tau)   }    # priors   mu ~ dnorm(0, 10)   tau <- 1/(sd*sd)   sd ~ dunif(0, 10) }  #---- loop--------  (i in 1:1000){   #--- data simulation --------- # paired data sim_bf <- rtnorm(16, mean = 0.4949, sd = 0.12, lower = 0.1) sim_af <- rtnorm(16, mean = 0.3959, sd = 0.12, lower = 0.1) diff = sim_bf - sim_af   #------setting jags--------------  data <- list(d.diff = diff) params <- c("mu", "tau", "sd") inits <- function(){   list(mu = rnorm(1),        sd = rlnorm(1))  }  #-------run jags----------------  output <- jags(data=data,                 # inits=inits,                 parameters.to.save=params,                n.iter=1000,                 n.burn=100,                 n.chains=2,                 n.thin=1,                model.file=model1,                progress.bar = "gui")  #-------convergence checker---------------- output.mcmc <- as.mcmc(output) x <- gelman.diag(output.mcmc) up1 <- x$psrf[1,2]        # approximate convergence diagnosed when upper limit close 1 up2 <- x$psrf[2,2]  up3 <- x$psrf[3,2]  up4 <- x$psrf[4,2]      if(up1 > 1.1 | up2 > 1.1 | up3 > 1.1 | up4 > 1.1)    next                   # 1 sided t-test   lo = output$bugsoutput$summary[2,4]   g.vec[i] <- ifelse((lo < 0), 0, 1) }  <- table(g.vec) g <- a[2]/1000 g

instead of

for(i in 1:1000){...     ...     if(...         next

do

i <- 0 while(i <=1000){...     <- i+1     ...     if(...         <- i-1         break

Comments

Post a Comment

Popular posts from this blog

java - Spring Data JPA: Why findOne(id) executing delete query internally? -

i have 1 application using @manytomany annotation link 2 different entities 1 parent entity. not sure if have implemented properly. when execute .findone() on parent entity's repository find it, @ first fetch data along executing delete queries internally delete record mapping table (a table jpa creates after adding @manytomany annotation.). following generated sql queries , code snippet: java class @entity @table(name = "inventory") public class inventory { @id long id; // many other properties @manytomany(fetch=fetchtype.lazy) private list<stream> streams; @manytomany(fetch=fetchtype.lazy) private list<subject> subjects; } there no other annotation i've put in stream , subject entity class link inventory class. when run app, create following tables: inventory stream subject inventory_stream inventory_subject (inventory_stream , inventory_subject, 2 tables created automatically, don't k...
Read more

python - Mongodb How to add addtional information when aggregating? -

i beginner , wrote line pipeline works, wanna add other information output, screen name , or number of tweets.i tried add under $group gave me syntax error everytime here pipeline: def make_pipeline(): # complete aggregation pipeline pipeline = [ { '$match': { "user.statuses_count": {"$gt":99 }, "user.time_zone": "brasilia" } }, { "$group": { "_id": "$user.id", "followers": { "$max": "$user.followers_count" } } }, { "$sort": { "followers": -1 } }, { "$limit" : 1 } ]; i using on example : { "_id" : objectid("5304e2e3cc9e684aa98bef97"), "text" : "first week of school on :p", "in_reply_to_status_id" : null, "retweet_count" : null, ...
Read more

java - Incorrect order of records in M-M relationship in hibernate -

Image
i have 4 records in bridge table in mysql database. using @manytomany in hibernate, accurate order of employees incorrect incorrect of projects. following required material figure out doing mistake. this source. @entity @table(name="project") public class projectbean { @id @generatedvalue @column(name="project_id") private int projectid; @column(name="title") private string projecttitle; @manytomany(mappedby="projects") private collection<employeebean> employees; public projectbean() { employees = new arraylist<employeebean>(); } //getters & setters } @entity @table(name="employee") public class employeebean { @id @generatedvalue @column(name="employee_id") private int employee_id; @column(name="first_name") private string ...
Read more