-
Notifications
You must be signed in to change notification settings - Fork 2
/
Cephalopholis_argus_cytb_Hawaii_parameters.R
49 lines (37 loc) · 1.83 KB
/
Cephalopholis_argus_cytb_Hawaii_parameters.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
#### Parameters for space_invaders.R model #####
NP <- 20 #Number of processors to use
MONTHS <- 12 #Number of "months" or portions to divide a year into
MODEL_DURATION <- 2 #Number of years to run the model
#### Demographic Parameters for Model ####
np <- NP-1 #Number of processors to use
NUM_BOOTSTRAPS <- 10 #Number of simulations to run with each source
##Native Range Demographics
#from arlequin
source.theta<-7.63797
source.theta.sd<-2.65535
## Haplotype information
source.hap<-c(194,172,24,19,14,19,15,1,1)
destination.hap<-c(159,405,3,34,0,0,0,0,0)
## Demographic Parameters
BIN <- 12 # Number of different age-classes
JUVI.MORT <- 1-0.165 # Juvenile mortality
ADULT.MORT <- 1-0.052 # Adult mortality
ADULT.FRAC <- 1-0.96 # Fraction of starting population that are adults*
# * Based on the empirical estimates of Belize 2014 lionfish sample
# on the forereef.
## Calculating Recruit per Individual
# This section takes egg and larval demographic parameters and calculates
# the monthly number of recruits per individual adult lionfish.
ADULT.FEM.FRAC <- 0.49 # Proportion of adults that are females
ADULT.FEM.MAT <- 0.79 # Proportion of females that are mature
FE <- 194577 # Fecundity - number of eggs per female per month
ME <- 0.31 # Egg mortality (days)
DE <- 3 # Egg duration (days)
ML <- 0.35 # Larval mortality (days)
DL <- 27 # Larval duration (days)
## Vectors and List for above ***** DO NOT CHANGE ****
#Mortalities and time in month for recovery
Demo.param <- c(ADULT.MORT,JUVI.MORT,ADULT.FRAC,BIN)
#Demographic parameters used to calculate recruit per individual in monthly
#time steps
RPR <- c(ADULT.FEM.FRAC,ADULT.FEM.MAT,FE,ME,DE,ML,DL)