Commit 0f8ba986 authored by matbuoro's avatar matbuoro
Browse files

Mise à jour CIEM

parent f482da25
rm(list=ls()) # Clear memory
year <- "2016"
######## Table 7 - Exploitation rate in the rivers Scorff ######
## SCORFF
site <- "Scorff"
stade <- "adult"
# load dataset
load(paste('~/Documents/RESEARCH/PROJECTS/ORE/Abundance/Scorff/adult/data/data_',stade,"_",year,'.Rdata',sep="")) # chargement des données
## Cm_F[t,a]: Annual number of marked fish caught by fishing per sea age category and showed at Moulin des Princes. 1:1SW, 2:MSW
## Cum_F[t,a]: Annual number of unmarked fish caught by fishing per sea age category and showed at Moulin des Princes. 1:1SW, 2:MSW
C_F_1SW <- data$Cm_F[,1] + data$Cum_F[,1] # marked + unmarked 1SW fish caugth by fishing
C_F_MSW <- data$Cm_F[,2] + data$Cum_F[,2] # marked + unmarked MSW fish caugth by fishing
# load estimations of size popualtions
load("~/Documents/RESEARCH/PROJECTS/ORE/Abundance/Scorff/adult/results/Results_adult_2016.RData")
n_1SW <- fit$median$n_1SW # medians
n_MSW <- fit$median$n_MSW # medians
Expl_rate <- cbind(
Expl_rate_1SW = (C_F_1SW / n_1SW)*100,
Expl_rate_MSW = (C_F_MSW / n_MSW)*100
)
Expl_rate <- rbind(Expl_rate,colMeans(Expl_rate))
rownames(Expl_rate) <- c(seq(1994,2016,1), "Average")
colnames(Expl_rate) <- c("1SW (%)", "MSW (%)")
write.csv(round(Expl_rate,1), file=paste('~/Documents/RESEARCH/PROJECTS/ORE/Abundance/CIEM/Table7_Scorff_',year,'.csv',sep=""))
######## Table 8 - Index rivers :spawning stock, egg deposition and attainment of CLs ######
## NIVELLE
site <- "Nivelle"
stade <- "adult"
load(paste("~/Documents/RESEARCH/PROJECTS/ORE/Abundance/Nivelle/adult/results/Results_adult_",year,".RData",sep=""))
years <- seq(1987, 2016, 1)
table <- array(, dim=c(length(years), 4))
colnames(table) <- c("1SW", "MSW", "eggs (million)", "eggs/CL")
rownames(table) <- years
#Conservation Limit:
CL = 1.44
######## Table 9 - juvenile and adult salmon numbers (estim.), in-river return rate in the monitored rivers ######
"","1SW (%)","MSW (%)"
"1994",3.7,21.3
"1995",6.4,5.7
"1996",4.9,6.9
"1997",5.3,8.2
"1998",9.8,18.5
"1999",7.8,3.3
"2000",12.7,11.9
"2001",5.8,10
"2002",2.1,4.3
"2003",0,11.9
"2004",7.7,15.4
"2005",6,20.2
"2006",7.9,30.3
"2007",10.3,12.6
"2008",4.2,16.3
"2009",1.9,11.3
"2010",7.6,4.5
"2011",7.3,5.5
"2012",9.1,18.4
"2013",4.9,6.8
"2014",8.2,9
"2015",8.7,17.1
"2016",5.4,12.2
"Average",6.4,12.2
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment