h <- matrix(c(1,2,2,1), nrow=2,ncol=2)
MultiHard(1:2, h)
# prints a sensible description of itself
data(betacells)
h <- 15.0 * matrix(c(NA,1,1,NA), nrow=2,ncol=2)
ppm(betacells, ~1, MultiHard(c("off","on"), h))
# fit the stationary multitype hardcore process to `betacells'
Run the code above in your browser using DataLab