Skip to content

Instantly share code, notes, and snippets.

Keybase proof

I hereby claim:

  • I am ozjimbob on github.
  • I am jimbob (https://keybase.io/jimbob) on keybase.
  • I have a public key whose fingerprint is 7ABC 13B3 C2B3 2583 C3A2 C03C EC2D AA9F 5376 6533

To claim this, I am signing this object:

Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
@ozjimbob
ozjimbob / engine.r
Created December 11, 2012 23:55
SPAM Model
# Run model simulation with given values for:
# initial cull, maintenance cull, cell target density, allocated budget,
# area to cull, and "control area" (area for which proportional population is calculated)
# (last two are kept separate for the spatial optimisation, where we may choose to cull a subset
# of the target control area for efficiency)
run.sim = function(curr.init.cull, curr.maint.cull, curr.target.density, budget, C.mask, ctrl.mask)
{
Nx = ncol(K); Ny = nrow(K); # gauge size of map
IC = K # note - initial population is set to carrying capacity
### Plot length summary
data=read.csv("plotBAincr.csv",stringsAsFactors=F)
data$Date.start=as.Date(data$Date.start)
data$Date.end=as.Date(data$Date.end)
plotlist=unique(data$Plot)
out_frame=data.frame(Plot="",Date.start="",Date.end="",length=0,stringsAsFactors=F)