Title: | An implementation of the EcoTroph Ecosystem modelling approach |
---|---|
Description: | An approach and software for modelling marine and freshwater ecosystems. It is articulated entirely around trophic levels. EcoTroph's key displays are bivariate plots, with trophic levels as the abscissa, and biomass flows or related quantities as ordinates. Thus, trophic ecosystem functioning can be modelled as a continuous flow of biomass surging up the food web, from lower to higher trophic levels, due to predation and ontogenic processes. Such an approach, wherein species as such disappear, may be viewed as the ultimate stage in the use of the trophic level metric for ecosystem modelling, providing a simplified but potentially useful caricature of ecosystem functioning and impacts of fishing. This version contains catch trophic spectrum analysis (CTSA) function and corrected versions of the mf.diagnosis and create.ETmain functions. |
Authors: | J. Guitton, M. Colleter, P. Gatti, and D. Gascuel |
Maintainer: | Jerome Guitton <[email protected]> |
License: | GPL |
Version: | 1.6.1 |
Built: | 2024-11-05 02:55:55 UTC |
Source: | https://github.com/polehalieutique/ecotroph |
convert.list2tab enables to convert the list object returned by the create.ETdiagnosis function into a list of data.frames. These data.frames contain calculated variables by TL class and combinations of effort multipliers.
check.table(ecopath)
check.table(ecopath)
ecopath |
is the input table used in ET (possibly based on Ecopath data). The different variables are the group name, its trophic level, biomass, production, catches, omnivory index and accessibility (fraction of the group that can be catch assuming an infinite fishing effort). |
data(ecopath_guinee) check.table(ecopath_guinee)
data(ecopath_guinee) check.table(ecopath_guinee)
This function enables the verification of input tables based on EwE data and used in the EcoTroph routine. A template is provided in the example: data(ecopath_guinee).
convert.list2tab(diagn.list)
convert.list2tab(diagn.list)
diagn.list |
is the list object returned by the create.ETdiagnosis function. |
data(ecopath_guinee) Liste=create.ETdiagnosis(create.ETmain(ecopath_guinee)) Tab=convert.list2tab(Liste)
data(ecopath_guinee) Liste=create.ETdiagnosis(create.ETmain(ecopath_guinee)) Tab=convert.list2tab(Liste)
ET-Transpose provides a picture of an ecosystem under a given fishing mortality. ET-Diagnosis is a routine simulating how this baseline ecosystem would be impacted by increasing or decreasing fishing effort. Fishing effort can be modified per fleet and/or trophic group. Ecosystem-wide effects of altering fishing effort include potential changes in biomass, accessible biomass, production, kinetics and catch trophic spectra, as well as impacts on the mean trophic level of the catch and biomass. Additionally, ET-Diagnosis constitutes a useful exploratory tool for ecosystem-based management. It simulates how reducing or increasing fishing effort and/or preferentially targeting different trophic levels could improve yield at the ecosystem scale. Lastly, ET-Diagnosis allows to view how different assumptions on ecosystem functioning (biomass input control, top-down effect) affect both trophic level specific and ecosystem-wide properties in relation to fishing.
create.ETdiagnosis(data, Mul_eff = NULL, Group = NULL, fleet.of.interest = NULL, same.mE = NULL, B.Input=NULL, Beta = NULL, TopD = NULL, FormD = NULL, TLpred = NULL)
create.ETdiagnosis(data, Mul_eff = NULL, Group = NULL, fleet.of.interest = NULL, same.mE = NULL, B.Input=NULL, Beta = NULL, TopD = NULL, FormD = NULL, TLpred = NULL)
data |
is the list object returned by the create.ETmain function. |
Mul_eff |
is a vector of fishing effort multipliers that the user wants to test. Mul_eff must contain the value 1 (reference state). By default, the function simulates a range of fishing effort multipliers from 0 to 5 for each fleet. |
Group |
is a character vector of trophic groups that the user specifically wants to impact by changing associated fishing efforts. By default, all trophic groups are equally impacted. |
fleet.of.interest |
is a character vector of fleet(s) that the user specifically wants to impact by changing associated fishing efforts (default =NULL). This argument is of particular interest if there are more than two fleets because it limits the mE combinations to be tested, and thus the associated computation time. |
same.mE |
is a logical argument (default=F), if TRUE the same effort multipliers are simultaneously applied to all fleets. |
B.Input |
is a logical argument (default=F), if TRUE the "Biomass input control" equation is accounted for in EcoTroph equations. |
Beta |
is a coefficient expressing the extent of the biomass input control. Beta=0 refers to an ecosystem where all secondary production originates from grazing on primary producers, and Beta=1 to an ecosystem where detritus and/or recruitment contribute to a major part of the biomass input (default=0.2). |
TopD |
is a coefficient expressing the top-down control, i.e. the fraction of the natural mortality depending on predator abundance. It varies between 0 and 1. The user can specify a numeric value, which is applied to each trophic level, or a numeric vector (of the same length as TL classes), i.e. a value for each TL (default=0.4). |
FormD |
is a shape parameter varying between 0 and 1. It defines the functional relationship between prey and predators. The value 1 refers to a situation where predators abundance has a linear effect on the speed of the flow of their preys. The user can specify a numeric value, which is applied to each trophic level, or a numeric vector (of the same length as TL classes), i.e. a value for each TL (default=0.5). |
TLpred |
is the trophic level that the user considers to be the "predator" trophic classes start. The default value is 3.5. |
Fleets' names used in the argument 'fleet.of.interest' are the catch column names of the Ecopath input data.frame (e.g. 'catch.1' or 'catch.ind').
This function returns a list of elements referring to each simulated combination of fishing effort multipliers. Each element is a list of two types of results:
- Variables characterizing the state and functioning of the modeled ecosystem: biomass, flow, kinetic, catches (total and per fleet) and fishing mortality per trophic level.
- Summary statistics (contained in the ET_Main_diagnose): absolute and relative (in comparison with the reference state) total biomass, flow, catches.
plot.ETdiagnosis and plot.ETdiagnosis_isopleth to plot the principle graphics resulting from the create.ETdiagnosis function, create.ETmain to create a list of tables used as input in the create.ETdiagnosis function.
data(ecopath_guinee) #Impacts of global changes in fishing efforts multipliers (in the range 0-5) create.ETdiagnosis(create.ETmain(ecopath_guinee),same.mE=TRUE) #Test of all the combinations of fishing effort multipliers per fleet #(in the range 0-5) create.ETdiagnosis(create.ETmain(ecopath_guinee)) #With biomass input control create.ETdiagnosis(create.ETmain(ecopath_guinee),B.Input=TRUE) #Impacts of changing fishing effort against Barracudas+ and Carangids groups create.ETdiagnosis(create.ETmain(ecopath_guinee), Mul_eff=(seq(0,5,.1)),Group=c('Barracudas+','Carangids'))
data(ecopath_guinee) #Impacts of global changes in fishing efforts multipliers (in the range 0-5) create.ETdiagnosis(create.ETmain(ecopath_guinee),same.mE=TRUE) #Test of all the combinations of fishing effort multipliers per fleet #(in the range 0-5) create.ETdiagnosis(create.ETmain(ecopath_guinee)) #With biomass input control create.ETdiagnosis(create.ETmain(ecopath_guinee),B.Input=TRUE) #Impacts of changing fishing effort against Barracudas+ and Carangids groups create.ETdiagnosis(create.ETmain(ecopath_guinee), Mul_eff=(seq(0,5,.1)),Group=c('Barracudas+','Carangids'))
This function enables the creation of the ET-Main table (summarizing the principal results/variables in function of the TL classes) and other intermediate tables of the ET-Transpose routine. It provides a picture of an ecosystem under a given fishing mortality.
create.ETmain(ecopath, smooth_type=NULL, sigmaLN_cst=NULL, pas=NULL, shift=NULL, smooth_param=NULL)
create.ETmain(ecopath, smooth_type=NULL, sigmaLN_cst=NULL, pas=NULL, shift=NULL, smooth_param=NULL)
ecopath |
is the input table used in ET (possibly based on Ecopath data). The different variables are the group name, its trophic level, biomass, production, catches, omnivory index and accessibility (fraction of the group that can be catch assuming an infinite fishing effort). |
smooth_type |
is a parameter of the create.smooth function. It defines the type of sigma calculation for the lognormal distribution. The value for this parameter is 1, 2 or 3. By default smooth_type=1, this defines a constant sigma. By choosing smooth_type=2, the user has the possibility to put a sigmaLN=smooth_param*ln(TL-0.05), with smooth_param=0.07 and shift=0.95 by default. Smooth_type=3 corresponds to the use of the calculated Omnivory Index (OI) divided by the associated mean TL as sigmaLN. |
sigmaLN_cst |
is a parameter of the create.smooth function. It defines the value of the constant sigma of the lognormal distribution for smooth_type=1. By default, sigmaLN_cst=0.12. |
pas |
is a parameter of the create.smooth function. It defines the splitting of the TL classes. |
shift |
is a parameter of the create.smooth function. It defines the beginning of the smooth function and allows the substraction of 0.05 in the sigma calculation accounting for the half interval range of the trophic class. |
smooth_param |
is a parameter of the create.smooth function. It defines the slope of the log-linearly increase of the TL variability with the mean trophic level of the group for smooth_type=2. SigmaLN(TL) is thus defined as sigmaLN(TL)=smooth_param*ln(TL-0.05). |
This function returns a list containing: the ET-Main table, intermediate matrices (biomass, accessible biomass, flowP...) and a list of matrices corresponding to the different fisheries catches.
plot.ETmain to create the principle graphics resulting from the create.ETmain function, create.smooth to create the Smooth table used in this function, Transpose to convert data referring to groups into data referring to TL classes.
data(ecopath_guinee) create.ETmain(ecopath_guinee) #Use of the second smooth type create.ETmain(ecopath_guinee,smooth_type=2)
data(ecopath_guinee) create.ETmain(ecopath_guinee) #Use of the second smooth type create.ETmain(ecopath_guinee,smooth_type=2)
create.smooth is used to create a smooth function. This function enables the conversion of data pertaining to specific taxa or functional groups into data by trophic classes. The main assumption of this Smooth function is that the distribution of the biomass (or catch...) of a trophic group around its mean trophic level follows a lognormal curve. The curve is defined by a mean (the mean TL of the trophic group) and a standart deviation (sigma), which is a measure of the trophic level variability within the group. The distribution is then defined by the lognormal function LN(mean TL, sigma).
create.smooth(tab_input, smooth_type=NULL, sigmaLN_cst=NULL, pas=NULL, shift=NULL, smooth_param=NULL)
create.smooth(tab_input, smooth_type=NULL, sigmaLN_cst=NULL, pas=NULL, shift=NULL, smooth_param=NULL)
tab_input |
is the input table used in ET (possibly based on Ecopath data). The different variables are the group name, its trophic level, biomass, production, catches, omnivory index and accessibility (fraction of the group that can be catch assuming an infinite fishing effort). |
smooth_type |
defines the type of sigma calculation for the lognormal distribution. Values of this parameter are 1, 2 or 3. By default smooth_type=1, this defines a constant sigma. By choosing smooth_type=2, the user has the possibility to implement a sigmaLN=smooth_param*ln(TL-0.05), with the parameter smooth_param=0.07 and shift=0.95 by default. Smooth_type=3 corresponds to the use of the omnivory index (OI) in the sigmaLN calculation (sigmaLN=OI/TL). |
sigmaLN_cst |
defines the value of the constant sigma of the lognormal distribution in case of smooth_type=1. By default, sigmaLN_cst=0.12. |
pas |
defines the splitting of the TL classes. By default, pas=0.1. |
shift |
defines the beginning of the smooth function and allows the substraction of 0.05 in the sigma calculation accounting for the half interval range of the trophic class. By default, with a constant sigmaLN (smooth_type=1), shift=1.8; with a function defined sigmaLN (smooth_type=2), shift=0.95; and with sigmaLN=OI/TL (smooth_type=3), shift=0. |
smooth_param |
defines the slope of the log-linear increase of the TL variability with the mean trophic level of the group. SigmaLN(TL) is thus defined as sigmaLN(TL)=smooth_param*ln(TL-0.05). By default, smooth_param=0.07. |
The user has the possibility to define sigmaLN for each trophic group and also adjust the LN distribution with the smooth_type, sigmaLN_cst, smooth_param, shift and pas parameters. Different choices are available : a constant sigma, a function defined sigma (log-linear increase) , or a sigma equal to the omnivory index divided by the associated mean TL.
create.smooth returns a table of the TL distribution within a trophic class. This table enables the calculation of Trophic Spectra, it is used in the Transpose function.
plot.smooth to plot the Smooth function, Transpose to build trophic spectra, plot.Transpose to plot the trophic spectra.
data(ecopath_guinee) create.smooth(ecopath_guinee) create.smooth(ecopath_guinee,sigmaLN_cst=0.11) create.smooth(ecopath_guinee,smooth_type=2,pas=0.2)
data(ecopath_guinee) create.smooth(ecopath_guinee) create.smooth(ecopath_guinee,sigmaLN_cst=0.11) create.smooth(ecopath_guinee,smooth_type=2,pas=0.2)
CTSA.catch.input is used to create inputs for the CTSA.forward function. It is a list of data.frames referring to catches per fleet formatted with TL classes in rows and trophic groups in columns.
CTSA.catch.input(catch.group,smooth_type=NULL,sigmaLN_cst=NULL, pas=NULL,shift=NULL,smooth_param=NULL)
CTSA.catch.input(catch.group,smooth_type=NULL,sigmaLN_cst=NULL, pas=NULL,shift=NULL,smooth_param=NULL)
catch.group |
is a data.frame containing: a column group_name, column(s) referring to the catches of each fleet (named 'catch.1', 'catch.2'...), a column TL specifying the mean TL of each group, and optionally a column OI (omnivory index) used for smooth_type=3. |
smooth_type |
is a parameter of the create.smooth function. It defines the type of sigma calculation for the lognormal distribution. The value for this parameter is 1, 2 or 3. By default smooth_type=1, this defines a constant sigma. By choosing smooth_type=2, the user has the possibility to put a sigmaLN=smooth_param*ln(TL-0.05), with smooth_param=0.07 and shift=0.95 by default. Smooth_type=3 corresponds to the use of the calculated Omnivory Index (OI) divided by the associated mean TL as sigmaLN. |
sigmaLN_cst |
is a parameter of the create.smooth function. It defines the value of the constant sigma of the lognormal distribution for smooth_type=1. By default, sigmaLN_cst=0.12. |
pas |
is a parameter of the create.smooth function. It defines the splitting of the TL classes. |
shift |
is a parameter of the create.smooth function. It defines the beginning of the smooth function and allows the substraction of 0.05 in the sigma calculation accounting for the half interval range of the trophic class. |
smooth_param |
is a parameter of the create.smooth function. It defines the slope of the log-linearly increase of the TL variability with the mean trophic level of the group for smooth_type=2. SigmaLN(TL) is thus defined as sigmaLN(TL)=smooth_param*ln(TL-0.05). |
CTSA.catch.input returns a list of data.frames, referring to catches per fleet formatted with TL classes in rows and trophic groups in columns.
create.smooth, Transpose and CTSA.forward.
data(ecopath_guinee) catch.group=ecopath_guinee[,c("group_name","TL","catch.1","catch.2")] Y_test <- CTSA.catch.input(catch.group) Y_test
data(ecopath_guinee) catch.group=ecopath_guinee[,c("group_name","TL","catch.1","catch.2")] Y_test <- CTSA.catch.input(catch.group) Y_test
E_MSY_0.1 computes two indices of exploitation: Emsy or Fmsy (maximum sustainable yield), and E0.1 or F0.1 ("start" of full exploitation) per TL class.
E_MSY_0.1(data, Mul_eff=NULL, B.Input=NULL, Beta=NULL,TopD=NULL, FormD=NULL, TLpred=NULL, maxTL=NULL)
E_MSY_0.1(data, Mul_eff=NULL, B.Input=NULL, Beta=NULL,TopD=NULL, FormD=NULL, TLpred=NULL, maxTL=NULL)
data |
is the list object returned by the create.ETmain function. |
Mul_eff |
is a parameter of the create.ETdiagnosis function. It is a vector of fishing effort multipliers that the user wants to test. Mul_eff must contain the value 1 (reference state). By default, the function simulates a range of fishing effort multipliers from 0 to 5 for each fleet. |
B.Input |
is a parameter of the create.ETdiagnosis function. It is a logical argument (default=F), if TRUE the "Biomass input control" equation is accounted for in EcoTroph equations. |
Beta |
is a parameter of the create.ETdiagnosis function. It is a coefficient expressing the extent of the biomass input control. Beta=0 refers to an ecosystem where all secondary production originates from grazing on primary producers, and Beta=1 to an ecosystem where detritus and/or recruitment contribute to a major part of the biomass input (default=0.2). |
TopD |
is a parameter of the create.ETdiagnosis function. It is a coefficient expressing the top-down control, i.e. the fraction of the natural mortality depending on predator abundance. It varies between 0 and 1. The user can specify a numeric value, which is applied to each trophic level, or a numeric vector (of the same length as TL classes), i.e. a value for each TL (default=0.4). |
FormD |
is a parameter of the create.ETdiagnosis function. It is a shape parameter varying between 0 and 1. It defines the functional relationship between prey and predators. The value 1 refers to a situation where predators abundance has a linear effect on the speed of the flow of their preys. The user can specify a numeric value, which is applied to each trophic level, or a numeric vector (of the same length as TL classes), i.e. a value for each TL (default=0.5). |
TLpred |
is a parameter of the create.ETdiagnosis function. It is the trophic level that the user considers to be the "predator" trophic classes start. The default value is 3.5. |
maxTL |
is a numeric string indicating the maximum TL for which indices are computed. |
For any TL class, if E0.1 and/or Emsy value(s) is(are) equal to the maximum effort multiplier tested (max(Mul_eff)), then E/F0.1 and/or E/Fmsy are set equal to NA.
The E_MSY_0.1 function returns a data.frame containing Fmsy, Emsy, F0.1 and E0.1 per TL class.
data(ecopath_guinee) E_MSY_0.1(create.ETmain(ecopath_guinee))
data(ecopath_guinee) E_MSY_0.1(create.ETmain(ecopath_guinee))
E0.1
E0.1(TL, Y, Fish_mort)
E0.1(TL, Y, Fish_mort)
TL |
Trophic level |
Y |
Catches |
Fish_mort |
Fishing mortality |
This example dataset is extracted from the 2004 Guinean Ecopath model (Gascuel et al., 2009). It provides a template for the input table formatting, the wanted variables names and the different capabilities of this package (used in function's examples).
data(ecopath_guinee)
data(ecopath_guinee)
A data.frame with 35 observations on the following 8 variables.
group_name
a character vector corresponding to the names of the trophic groups used in the Ecopath model. Has obligatory to be written 'group_name'.
TL
a numeric vector corresponding to the trophic level of the associated trophic groups. Has obligatory to be written 'TL'.
biomass
a numeric vector corresponding to the biomass of the associated trophic groups. Has obligatory to be written 'biomass'.
prod
a numeric vector corresponding to the production on biomass ratio. For the Detritus groups (no P/B value entered in Ecopath), put 0 as a value. Has obligatory to be written 'prod'.
catch.1
a numeric vector corresponding to the catch of the artisanal fleet. A value must be entered for all groups, with a 0-value if no catch are made. Has obligatory to be written 'catch.something'.
catch.2
a numeric vector corresponding to the catch of the industrial fleet. A value must be entered for all groups, with a 0-value if no catch are made. Has obligatory to be written 'catch.something'
accessibility
a numeric vector corresponding to the fraction of the trophic group that can be catch assuming an infinite fishing effort. Has obligatory to be written 'accessibility'.
OI
a numeric vector corresponding to the omnivory index calculated by Ecopath for each trophic group. Has obligatory to be written 'OI'.
No NA are accepted in the dataset (0 for the P/B of the detritus groups, 0 for the catch...). Follow the instructions stated in the variables descriptions. Different fleets can be entered in the model using the following system: catch.1, catch.2, catch.whatyouwant... If there is only one fleet, you just have to put catch as a variable name.
Gascuel et al. (2009) Impact de la peche sur l'ecosysteme marin de Guinee - Modelisation EwE 1985/2005 -
data(ecopath_guinee) ecopath_guinee names(ecopath_guinee)
data(ecopath_guinee) ecopath_guinee names(ecopath_guinee)
E_MSY
Emsy(TL, Y, Fish_mort)
Emsy(TL, Y, Fish_mort)
TL |
Trophic level |
Y |
Catches |
Fish_mort |
Fishing mortality |
Effort multiplier diagnosis
mf.diagnosis( x, ET_Main, catch.list, TL_out, fleet, n.fleet, Fish_mort_ref, Fish_mort_acc_ref, B.Input, Beta, TopD, FormD, TLpred, n.TL, range.TLpred, lim.high.TL, range.highTL )
mf.diagnosis( x, ET_Main, catch.list, TL_out, fleet, n.fleet, Fish_mort_ref, Fish_mort_acc_ref, B.Input, Beta, TopD, FormD, TLpred, n.TL, range.TLpred, lim.high.TL, range.highTL )
x |
is the list object returned by the create.ETdiagnosis function. |
ET_Main |
is the list object returned by the create.ETmain function. |
catch.list |
catches for each fleet |
TL_out |
Maximum TL to consider |
fleet |
List of available fleet |
n.fleet |
Number of fleet |
Fish_mort_ref |
Fishing mortality |
Fish_mort_acc_ref |
Accessible Fisshing mortality |
B.Input |
Biomass |
Beta |
Beta parameters of the formula |
TopD |
Parameters of the formula |
FormD |
Parameters of the formula |
TLpred |
Trophic level for predators |
n.TL |
Number of trophic levels |
range.TLpred |
Range of trophic level for predators |
lim.high.TL |
Limit for high trophic levels |
range.highTL |
Range for High trophic levels |
This function enables to plot the mixed impacts of changes in fishing effort for two fleets (or groups of fleets).
plot_ETdiagnosis_isopleth( x, fleet1 = NULL, fleet2 = NULL, var = NULL, n.level = NULL, relative = NULL, name.fleet1 = NULL, name.fleet2 = NULL, color = NULL, ask = interactive() )
plot_ETdiagnosis_isopleth( x, fleet1 = NULL, fleet2 = NULL, var = NULL, n.level = NULL, relative = NULL, name.fleet1 = NULL, name.fleet2 = NULL, color = NULL, ask = interactive() )
x |
is the list object returned by the create.ETdiagnosis function. |
fleet1 |
is a character vector of fleets for which fishing efforts are equally changed. |
fleet2 |
is a second character vector of fleets for which fishing efforts are equally changed. Fishing efforts remain unchanged for fleets not assigned in fleet1 or fleet2. If fleet2 is NULL, all fleets not assigned in fleet1 are assigned in fleet2. If the argument 'fleet.of.interest' has been assigned in the function create.ETdiagnosis, fleet1=fleet.of.interest and thus fleet2 is composed of the remaining fleet(s) (not assigned in fleet.of.interest). |
var |
is a character vector of plotted variables (TOT_biomass,TOT_biomass_acc,Y,Y_fleet1,Y_fleet2,TL_TOT_biomass,TL_TOT_biomass_acc,TL_Y,TL_Y_fleet1,TL_Y_fleet2).All the variables are plotted by default. |
n.level |
is a numeric string, specifying the number of plotted isopleth areas (7 is the default value). |
relative |
is a logical string (by default relative=F), specifying if the variables have to be plotted in absolute or relative values (in comparison with reference state, Mul_eff=1). Note that if relative=TRUE, mean trophic level in biomass or catches (TL_TOT_biomass,TL_Y,...) are not plotted. |
name.fleet1 |
is a character string used to implement x-axis name. By default name.fleet1='fleet 1'. |
name.fleet2 |
is a character string used to implement y-axis name. By default name.fleet2='fleet 2'. If the argument fleet.of.interest has been assigned in the function create.ETdiagnosis, name.fleet1 = 'fleet of interest' and name.fleet2 = 'other fleets'. |
color |
is a vector of colors, the length of this vector should be equal to the value of levels. By default, color=rainbow(n=levels). |
ask |
default value is interactive. Parameter used to enable the user to control the display of each graph. |
Fleets' names used in the arguments 'fleet1' and 'fleet2' are the catch column names of ecopath input dataframe (e.g. 'catch.1' or 'catch.ind').
data(ecopath_guinee) diagn.list=create.ETdiagnosis(create.ETmain(ecopath_guinee)) plot_ETdiagnosis_isopleth(diagn.list,fleet1='catch.1',fleet2='catch.2') plot_ETdiagnosis_isopleth(diagn.list,fleet1='catch.1',fleet2='catch.2',relative=TRUE)
data(ecopath_guinee) diagn.list=create.ETdiagnosis(create.ETmain(ecopath_guinee)) plot_ETdiagnosis_isopleth(diagn.list,fleet1='catch.1',fleet2='catch.2') plot_ETdiagnosis_isopleth(diagn.list,fleet1='catch.1',fleet2='catch.2',relative=TRUE)
This function enables the creation of the principle graphics resulting from the create.ETdiagnosis function. The function returns the principal plots of the global ET-Diagnosis routine: the graphics of the biomass, accessible biomass...rates for the different effort multipliers, the Biomass Trophic Spectra (BTS) for the different effort multipliers, the B/Bref(mE=1) and Y/Yref graphs for the main TL classes and the Catch Trophic Spectra (CTS) (global and per fleet).
## S3 method for class 'ETdiagnosis' plot( x, scale = NULL, maxrange = NULL, legend.cex = NULL, ask = interactive(), ... )
## S3 method for class 'ETdiagnosis' plot( x, scale = NULL, maxrange = NULL, legend.cex = NULL, ask = interactive(), ... )
x |
is the list object returned by the create.ETdiagnosis function. |
scale |
is the scale parameter of the Biomass Trophic Spectra, can be log or by default the standard scale of results. |
maxrange |
is the maximum TL wanted for the x-axis. By default maxrange = 5.5. |
legend.cex |
defines the value of the cex for the legend. |
ask |
default value is interactive. Parameter used to enable the user to control the display of each graph. |
... |
plot other arguments |
data(ecopath_guinee) diagn.list<-create.ETdiagnosis(create.ETmain(ecopath_guinee),same.mE=TRUE) plot(diagn.list)
data(ecopath_guinee) diagn.list<-create.ETdiagnosis(create.ETmain(ecopath_guinee),same.mE=TRUE) plot(diagn.list)
This function enables the display of the principle plots resulting from the create.ETmain function: Biomass Trophic Spectra, Accessible Biomass Trophic Spectra, Catch by fleet Trophic spectra, Total Catch Trophic Spectra and other summary plots.
## S3 method for class 'ETmain' plot( x, scale1 = NULL, scale2 = NULL, scale3 = NULL, legend.cex = NULL, ask = interactive(), ... )
## S3 method for class 'ETmain' plot( x, scale1 = NULL, scale2 = NULL, scale3 = NULL, legend.cex = NULL, ask = interactive(), ... )
x |
is the list of tables returned by the create.ETmain function. |
scale1 |
defines the scale of the Biomass plots y-axis: can be log or not. |
scale2 |
defines the scale of the Accessible Biomass plots y-axis: can be log or not. |
scale3 |
defines the scale of the Catch by fleet plots y-axis: can be log or not. |
legend.cex |
defines the value of the cex for the legend. |
ask |
default value is interactive. Parameter used to enable the user to control the display of each graph. |
... |
plot other arguments |
The function returns the principal graphics of the global ET-Transpose routine: the Biomass Trophic Spectra, the Accessible Biomass Trophic Spectra and other graphics, notably the Catch Trophic Spectra.
create.smooth function to create the Smooth, Transpose to calculate the data transposition into trophic spectra, create.ETmain to create a list of table including the ET-Main table.
data(ecopath_guinee) plot(create.ETmain(ecopath_guinee),scale1=log) plot(create.ETmain(ecopath_guinee),scale1=log,scale3=log)
data(ecopath_guinee) plot(create.ETmain(ecopath_guinee),scale1=log) plot(create.ETmain(ecopath_guinee),scale1=log,scale3=log)
plot.smooth is used to plot the Smooth function. This function enables the user to see the TL distributions around their mean trophic levels.
## S3 method for class 'smooth' plot(x, ...)
## S3 method for class 'smooth' plot(x, ...)
x |
is the table returned by the create.smooth function. |
... |
plot other arguments |
create.smooth function to create the Smooth, Transpose to calculate the data transposition into trophic spectra.
data(ecopath_guinee) plot(create.smooth(ecopath_guinee)) plot(create.smooth(ecopath_guinee,smooth_type=2))
data(ecopath_guinee) plot(create.smooth(ecopath_guinee)) plot(create.smooth(ecopath_guinee,smooth_type=2))
This function returns the two principal plots of the Transpose function : a plot by group and the associated Trophic Spectra (CTS, BTS...).
## S3 method for class 'Transpose' plot(x, title = NULL, scale = NULL, legend.cex = NULL, ...)
## S3 method for class 'Transpose' plot(x, title = NULL, scale = NULL, legend.cex = NULL, ...)
x |
is the table returned by the Transpose function. |
title |
defines the title of the graph. |
scale |
defines the scale of the y-axis: can be log or not. |
legend.cex |
defines the value of the cex for the legend. |
... |
plot other arguments |
create.smooth function to create the Smooth, plot.smooth to plot the smooth function, Transpose to calculate the data transposition into trophic spectra.
data(ecopath_guinee) smoothed<-create.smooth(ecopath_guinee) plot(Transpose(smoothed,ecopath_guinee,"biomass"),scale=log) plot(Transpose(smoothed,ecopath_guinee,"catch.1"),title="Small Scale Fishery Catch")
data(ecopath_guinee) smoothed<-create.smooth(ecopath_guinee) plot(Transpose(smoothed,ecopath_guinee,"biomass"),scale=log) plot(Transpose(smoothed,ecopath_guinee,"catch.1"),title="Small Scale Fishery Catch")
This function loads input data from an xml file created by the user, or exported from the EwE EcoTroph plug-in, or from a web service associated to a database populated with parameters of several EwE models.
read.ecopath.model(filename)
read.ecopath.model(filename)
filename |
is the address of the file the user wants to import. |
This function returns a data.frame containing all the column needed to run the EcoTroph R package.
check.table to control the reliability of the dataset.
function used to compute pB for the higest trophic levels @param compteur @param pb.mf @param TL_out @param range.highTL
regPB(compteur, pb.mf, TL_out, range.highTL)
regPB(compteur, pb.mf, TL_out, range.highTL)
function used to compute pB for the higest trophic levels and accessible biomass @param compteur @param pb.mf @param TL_out @param range.highTL
regPB.ac(compteur, pb.mf.ac, TL_out, range.highTL)
regPB.ac(compteur, pb.mf.ac, TL_out, range.highTL)
This function enables an other calculation for the sigma of the create.smooth function. Sigma is calculated on the base of a saturation function reflecting a biological reasoning about the variability of the TL within trophic classes: the variability increases with the TL and reaches a plateau after a certain TL.
saturation(sigma_inf = NULL, coeff = NULL, pas = NULL)
saturation(sigma_inf = NULL, coeff = NULL, pas = NULL)
sigma_inf |
defines the value of the curve's plateau. |
coeff |
defines the value of the slope. |
pas |
defines the splitting of the TL classes. |
By default sigma is constant. This function enables an other user choice reflecting a different reasoning.
saturation returns a vector of values for the sigma used in the create.smooth function.
create.smooth function to create the Smooth, plot.smooth to plot the smooth function.
plot(saturation()) lines(saturation(0.2)) text(48,0.18,"sigma_inf=0.2") lines(saturation(coeff=0.5)) text(48,0.35,"coeff=0.5")
plot(saturation()) lines(saturation(0.2)) text(48,0.18,"sigma_inf=0.2") lines(saturation(coeff=0.5)) text(48,0.35,"coeff=0.5")
Transpose enables the conversion of data pertaining to specific taxa or functionnal groups into data by trophic class. Data can represent catches, biomasses or production in order to produce continuous distributions of those variables over trophic levels.
Transpose(tab_smooth, tab_input, column)
Transpose(tab_smooth, tab_input, column)
tab_smooth |
is the table returned by the create.smooth function. |
tab_input |
is the input table based on Ecopath data or on independent data. The different variables are the group name, its trophic level, biomass, production on biomass ratio, catches, omnivory index and accessibility (fraction of the group that can be catch assuming an infinite fishing effort) if the input table corresponds to an EwE model. In other case, to simply build trophic spectra, only the group names, their trophic levels and related variables are necessary. |
column |
is the tab_input table column name of the variable the user wants to transpose (for example "biomass" or "catch"). |
create.smooth function to create the Smooth, plot.smooth to plot the smooth function, plot.Transpose to plot the associated trophic spectra.
data(ecopath_guinee) Transpose(create.smooth(ecopath_guinee),ecopath_guinee,"biomass") Transpose(create.smooth(ecopath_guinee),ecopath_guinee,"catch.1")
data(ecopath_guinee) Transpose(create.smooth(ecopath_guinee),ecopath_guinee,"biomass") Transpose(create.smooth(ecopath_guinee),ecopath_guinee,"catch.1")