data(betacells)
   # unmarked 
   X <- unmark(betacells)
   model <- ppm(X, ~1, PairPiece(seq(10,100,by=10)))
   f <- fitin(model)
   f
   plot(f)
   # marked
   # fit the stationary multitype Strauss process to `betacells'
   r <- 30.0 * matrix(c(1,2,2,1), nrow=2,ncol=2)
   model <- ppm(betacells, ~1, MultiStrauss(c("off","on"), r))
   f <- fitin(model)
   f
   plot(f)Run the code above in your browser using DataLab