inst/apps/opinion_dynamics/global.R

In swarm-lab/teachR: R Package For Distributing Web Applications For Educational Use

pacman::p_load(shiny, shinythemes, simecol, ggvis, dplyr)

grid_sys <- function(time, init, parms) {
  o1 <- init
  o1[o1 < 0] <- 0
  
  o2 <- init
  o2[o2 > 0] <- 0
  
  p1 <- eightneighbors(o1) / 8
  p2 <- p1 - eightneighbors(o2) / 8
  
  r <- runif(nrow(init) * ncol(init))
  
  idx1 <- r <= p1
  idx2 <- r <= p2 & r > p1
  
  init[idx1] <- 1 * parms$w1
  init[idx2] <- -1 * parms$w2
  
  init
}