Instructions/testing.r

In vpourfaraj/lobsterCatch0.1: Models capture processes in American lobster fishery

#testing

require(lobsterCatch)
require(devtools)
load_all('~/git/lobsterCatch')
require(geosphere)
nr = 18
nc = 18
arena = matrix(0,nr,nc)
y=x=seq(1,17,1)
require(geosphere)
load_all('D:/git/lobsterCatch')

arena = matrix(0,12,12)
y=x=seq(1,11,2)
traps = expand.grid(x,y)

p = list()
p$nrowgrids = nr
p$ncolgrids = nc
p$ngrids=p$nrowgrids * p$ncolgrids
p$initlambda=1
p$initD = 3
p$smult = 0.993
p$currentZoIInit = 1

p$trapEastStart = traps[,1]
p$trapNorthStart = traps[,2]
p$ntrapsstart = length(p$trapEastStart)

p$saturationThresholdStart = 5
p$how_closeStart = .1
p$dstepstart = 5
p$trapSaturationStart = T

p$tSteps = 50

p$realizations = 1

a = SimulateLobsterMovement(p=p)

#nLobsters
uL = unique(a[[1]]$lobsters$I)
dims = length(uL)
dims = ceiling(sqrt(dims))
par( mfrow=c(dims,dims) ) # create a plot with 1 row and 2 columns to show plots side by side
for(i in uL){
with(subset(a[[1]]$lobsters, I==i),plot(EASTING, NORTHING, type='l',ylim=c(-30,30),xlim=c(-30,30)))
points(traps, pch=16, col='red')
#with(subset(outmove, I==1), text(x=EASTING[c(1,100)], y=NORTHING[c(1,100)], c('Start','End')))
}

#testing trap sat

smm = c(1,3,5,10)
gg = list()
p$initlambda = .5
for(i in 1:length(smm)){
p$saturationThresholdStart=smm[i]
a = SimulateLobsterMovement(p=p)
gg[[i]]=c(mean(apply(a[[1]]$traps,1,max)),dispersion(apply(a[[1]]$traps,2,max)))
}