/
rec_effort_div_analysis.R
31 lines (23 loc) · 1.41 KB
/
rec_effort_div_analysis.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
# rec div analysis
#```{r, eval = F, include = T, echo = T}
REC_DATA <- read.csv('X:/gdepiper/ESR2018/SOE/Data/Rec_Days_Fished_2018.csv', as.is=TRUE)
REC_DATA$P_Shore <- -(as.numeric(gsub("","",'""',REC_DATA$Shore)) /
as.numeric(gsub("","",'""',REC_DATA$All_Modes))) *
log(as.numeric(gsub("","",'""',REC_DATA$Shore)) /
as.numeric(gsub("","",'""',REC_DATA$All_Modes)))
REC_DATA$P_Private <- -(as.numeric(gsub("","",'""',REC_DATA$Private_Rental)) /
as.numeric(gsub("","",'""',REC_DATA$All_Modes))) *
log(as.numeric(gsub("","",'""',REC_DATA$Private_Rental)) /
as.numeric(gsub("","",'""',REC_DATA$All_Modes)))
REC_DATA$P_Party <- -(as.numeric(gsub("","",'""',REC_DATA$Party_Charter)) /
as.numeric(gsub("","",'""',REC_DATA$All_Modes))) *
log(as.numeric(gsub("","",'""',REC_DATA$Party_Charter)) /
as.numeric(gsub("","",'""',REC_DATA$All_Modes)))
REC_DATA$Value <- exp(REC_DATA$P_Shore+REC_DATA$P_Private+REC_DATA$P_Party)
REC_DATA$Region[REC_DATA$Region=='"MID-ATLANTIC"'] <- 'MA'
REC_DATA$Region[REC_DATA$Region=='"NORTH ATLANTIC"'] <- 'NE'
E_SHANNON <- subset(REC_DATA, select=c('Time','Region','Value'))
E_SHANNON$Units <- 'Effective Shannon'
E_SHANNON$Var <- 'Recreational fleet effort diversity across modes'
E_SHANNON$Source <- 'MRIP effort time series, processed to generate diversity measure.'
#```