Gibbs sampling for MRFs

[davharris]

[davharris]

Partly implemented in the simulation script that makes fakedata

[davharris]

Example code:

# Gibbs sampling.  Update one species (column) at a time, conditional 
# on everyone else.
for(i in 1:iter){
  for(j in 1:n.species){
    state[ , j] = rbinom(
      n.sites,
      prob = sigmoid(inputs[, j] + state %*% lateral[, j]),
      size = 1
    )
  }
}