inst/Assessments/LFA34Assessment/3.ginifisheryFootprint.r

In LobsterScience/bio.lobster: Lobster Stock Assessment Tools for DFO Maritimes

# Lobster Fishery Footprint
require(bio.survey)
require(bio.lobster)
p = bio.lobster::load.environment()
p$libs = NULL

require(bio.utilities)
la()

a = lobster.db('process.logs') 

#grid areas by LFA
		LFAs<-read.csv(file.path( project.datadirectory("bio.lobster"), "data","maps","LFAPolys.csv"))
		LFAgrid<-read.csv(file.path( project.datadirectory("bio.lobster"), "data","maps","GridPolys.csv"))
		attr(LFAgrid,'projection') <- "LL"
		ar = calcArea(LFAgrid)
		names(ar)[1:2] <- c('LFA','GRID_NUM')
		ar7 = subset(ar,LFA==37)
		ar6 = subset(ar,LFA==36)
		ar8 = subset(ar,LFA==38)

		ar6b = ar7
		ar6b$LFA=36
		ar6 = as.data.frame(rbind(ar6,ar6b))
		ar6 = aggregate(area~LFA+GRID_NUM,data=ar6,FUN=sum)

		ar8b = ar7
		ar8b$LFA=38
		ar8 = as.data.frame(rbind(ar8,ar8b))
		ar8 = aggregate(area~LFA+GRID_NUM,data=ar8,FUN=sum)

		ar = subset(ar, !LFA %in% c(36,38))
		ar = as.data.frame(rbind(rbind(ar,ar6),ar8))


aA = aggregate(cbind(WEIGHT_KG,NUM_OF_TRAPS)~GRID_NUM+LFA+SYEAR,data=a,FUN=sum)
aA = merge(aA, ar, by=c('LFA','GRID_NUM'),all.x=T)
aA$CPUE = aA$WEIGHT_KG / aA$NUM_OF_TRAPS
lfa = c(34)
out = list()
m= 0 
for(i in 1:length(lfa)){
		aB = subset(aA,LFA==lfa[i])
		y = unique(aB$SYEAR)

		for(j in 1:length(y)){
			print(paste(lfa[i],y[j]))
			aC = subset(aB,SYEAR==y[j])
			if(nrow(aC)>=5){
			m=m+1
			out[[m]] = c(lfa[i],y[j],giniFootprint(aC$CPUE,aC$area),giniFootprint(aC$WEIGHT_KG,aC$area),giniFootprint(aC$NUM_OF_TRAPS,aC$area))
		}
	}
}
	out = as.data.frame(do.call(rbind,out))
	names(out) = c('LFA','YEAR','CPUE','LANDINGS','EFFORT')

out = out[order(out$LFA, out$YEAR),]
write.csv(subset(out,LFA==34),file=file.path(project.figuredirectory('bio.lobster'),"LFA34Assessment",'GiniLandings34.csv'))
png(file=file.path(project.figuredirectory('bio.lobster'),"LFA34Assessment",'GiniLandings34.png'),units='in',width=10,height=8,pointsize=18, res=300,type='cairo')
with(subset(out,LFA==34),plot(YEAR,LANDINGS, pch=16,xlab='Year',ylab='Gini Index'))
with(subset(out,LFA==34),lines(YEAR,runmed(LANDINGS,k=3), lwd=3,col='salmon'))
dev.off()