knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)

What is Back-Calculation

Francis (1990) defined back-calculation as "... a technique that uses a set of measurements made on a fish at one time to infer its length at an earlier time or times. Specifically, the dimensions of one or more marks in some hard part of the fish, together with its current body length, are used to estimate its length at the time of formation of each of the marks. ... The marks are generally annual rings associated with growth checks, ... ." Thus, backcalculation is the reconstruction of the length of a fish at previous ages from measurements made on calcified structures.

More thorough introductions to back-calculation are in Vigliola and Meekan (2009) and Section 11.2.3.1 in Shoup and Michaeletz (2017).

\

\


Terminology

Calcified structures are scales, spines, fin rays, otoliths, or other bones from fish that can be processed in such a way as to show marks that represent annuli, or yearly benchmarks. Hereafter, "structure" will refer to "calcified structure."

Two types of measurements can be made on structures. A radial measurement is the distance from the center of the structure (e.g., focus of scale or nucleus of otolith) to the edge of an annulus. An incremental measurement is the distance between two successive annuli. Radial measurements are required for back-calculation of fish length. Incremental measurements are used in some forms of modeling growth [e.g., Weisberg (1993) and Weisberg et al. (2010)].

Back-calculation models estimate length at previous age $i$ (i.e., $L_{i}$) from known values of length at time of capture ($L_{Cap}$), scale radius to the $i$th annulus ($R_{i}$), and scale radius at time of capture ($R_{Cap}$). Several back-calculation models rely on the relationship between $R_{Cap}$ and $L_{Cap}$. Depending on the model, a function of mean $R_{Cap}$ for a given $L_{Cap}$ (i.e., $E(R_{Cap}|L_{Cap})$ ) or mean $L_{Cap}$ for a given $R_{Cap}$ (i.e., $E(L_{Cap}|R_{Cap})$) is used. These functions are not required to be linear, but often are, and in their linear form are represented as

$$ E(L_{Cap}|R_{Cap})=a+bR_{Cap} \quad \quad \text{(1)} $$

$$ E(R_{Cap}|L_{Cap})=A+BL_{Cap} \quad \quad \text{(2)} $$

\

\


Back-Calculation Models

Four Most Common Models

The first back-calculation model was jointly developed by Knut Dahl and Einar Lea and appeared in Lea (1910). The underlying concept of the Dahl-Lea model is that growth of the calcified structure is in exact proportion to growth in length of the fish. With this, the ratio of $R_{i}$ to $R_{Cap}$ is the same as the ratio of $L_{i}$ to $L_{Cap}$. Rearrangement of this equality yields the Dahl-Lea back-calculation model

$$ L_{i}=\frac{R_{i}}{R_{Cap}}L_{Cap} \quad \quad \text{(3)} $$

The Dahl-Lea model describes a family of straight lines that pass through the origin and each observed ($R_{Cap},L_{Cap}$) point. Visually (Figure 1), the estimated $L_{i}$ for a particular fish is found by locating $R_{i}$ along the x-axis, moving vertically until the straight line for that fish is met, and then moving horizontally to the point on the y-axis.

\

data(SMBassWB,package="FSA")
wb90 <- subset(SMBassWB,yearcap=1990)
par(mfrow=c(2,2),mar=c(3.05,3.05,1.15,0.65),mgp=c(1.5,0.3,0),tcl=-0.2)
# ------------------------------------------------------------
# Helper functions
# ------------------------------------------------------------
makePlot <- function(name) {
  plot(lencap~radcap,data=wb90,
       xlim=c(0,max(radcap)),ylim=c(0,max(lencap)),
       xlab=expression(R[Cap]),ylab=expression(L[Cap]),
       xaxt="n",yaxt="n",cex=0.75)
  axis(2,at=seq(0,350,50),labels=c(0,50,NA,NA,seq(200,350,50)))
  xs <- c(0,2,6,8,10,12,14,16)
  axis(1,xs)
  mtext(name,line=0.1)
}
addCalc <- function(Ri,Li,Rc,Lc,int) {
  points(Rc,Lc,pch=19,col="red")
  lines(c(0,Rc),c(int,Lc),lwd=2,col="red")
  text(Ri,-10,paste(round(Ri,3)),pos=1,xpd=TRUE)
  arrows(Ri,-5,Ri,Li,lwd=1,length=0.1,angle=20,xpd=TRUE,col="red")
  arrows(Ri,Li,-0.1,Li,lwd=1,length=0.1,angle=20,xpd=TRUE,col="red")
  text(-0.1,Li,paste(round(Li,2)),pos=2,xpd=TRUE)
}

## Needed regressions and coefficients
lm.sl <- lm(radcap~lencap,data=wb90)
a <- coef(lm.sl)[1]; b <- coef(lm.sl)[2]
lm.ls <- lm(lencap~radcap,data=wb90)
c <- coef(lm.ls)[1]; d <- coef(lm.ls)[2]

## Info for fish #701
Rc701 <- 9.2219; Lc701 <- 312
## Info for fish #704
Rc <- 7.44389; Lc <- 218; Ri <- 3.49804

#-----------------------------------------------------------------------
# Plots the Dahl-Lea
#-----------------------------------------------------------------------
makePlot("Dahl-Lea")
# Fish #701 as an example -- just show line
points(Rc701,Lc701,pch=19,col="green3",cex=1.25)
lines(c(0,Rc701),c(0,Lc701),lwd=2,col="green3")
# Fish #704 as an example -- show calculation
Li <- (Ri/Rc)*Lc
addCalc(Ri,Li,Rc,Lc,0)

#-----------------------------------------------------------------------
# Plots Fraser-Lee method
#-----------------------------------------------------------------------
makePlot("Fraser-Lee")
# Fish #701 as an example -- just show line
points(Rc701,Lc701,pch=19,col="green3",cex=1.25)
lines(c(0,Rc701),c(c,Lc701),lwd=2,col="green3")
# Fish #704 as an example -- show calculation
Li <- (Ri/Rc)*(Lc-c)+c
addCalc(Ri,Li,Rc,Lc,c)

#-----------------------------------------------------------------------
# Plots the SPH method
#-----------------------------------------------------------------------
makePlot("SPH")
# Fish #701 as an example -- just show line
points(Rc701,Lc701,pch=19,col="green3",cex=1.25)
lines(c(0,Rc701),c(-a/b,Lc701),lwd=2,col="green3")
# Fish #704 as an example -- show calculation
Li <- (-a/b)+(Lc+a/b)*(Ri/Rc)
addCalc(Ri,Li,Rc,Lc,-a/b)

#-----------------------------------------------------------------------
# Plots the BPH method
#-----------------------------------------------------------------------
makePlot("BPH")
# Fish #701 as an example -- just show line
int <- (c*Lc701)/(c+d*Rc701)
points(Rc701,Lc701,pch=19,col="green3",cex=1.25)
lines(c(0,Rc701),c(int,Lc701),lwd=2,col="green3")
# Fish #704 as an example -- show calculation
int <- (c*Lc)/(c+d*Rc)
Li <- Lc*(c+d*Ri)/(c+d*Rc)
addCalc(Ri,Li,Rc,Lc,int)

Figure 1. Plot of length-at-capture versus scale radius for West Bearskin Lake Smallmouth Bass in 1990. All four methods of backcalculation are shown for fish 704 ($R_{2}$=3.498, $L_{Cap}$=218, and $R_{Cap}$=7.44389; red point and line) with calculational steps shown with the arrows. Fish 701 is shown as the green point and line for comparative purposes.

\

Fraser (1916) was the first to describe, but Lee (1920) was the first to formally derive, the back-calculation model from the concept that "the growth increment of the scale is, on the average ..., a constant proportion of the growth increment of the fish" (Francis 1990). In practice, the Fraser-Lee model modified the Dahl-Lea model by adjusting for the length of the fish when the structure forms (i.e., $L=c$ when $R=0$), that is,

$$ L_{i}=\frac{R_{i}}{R_{Cap}}(L_{Cap}-a)+a \quad \quad \text{(4)} $$

\noindent where $a$ comes from the length of the fish at the time of structure formation, the intercept of the length-structure relationship regression (e.g., from Equation 1), or, when using scales, from published "standards" for a species (Carlander 1982). The Fraser-Lee model describes a family of lines with an intercept of $a$ that pass through the ($R_{Cap},L_{Cap}$) point (Francis 1990; Figure 1).

The scale proportional hypothesis (SPH) model was named by Francis (1990), but was first recognized by Whitney and Carlander (1956) when they said "{i}f the scale was 10 per cent larger when the fish was caught than the average scale for that size of fish, [then] the scale would be 10 per cent larger than normal throughout the life." If "average" and "normal" are considered to be expected values, then this hypothesis can be written as

$$ \frac{R_{i}}{E[R|L_{i}]}=\frac{R_{Cap}}{E[R|L_{Cap}]} $$

If it is assumed that the scale-length relationship is linear, then the two expected values in these ratios are computed by plugging $L_{i}$ and $L_{Cap}$, respectively, into the scale-length relationship (i.e., Equation 2) to produce

$$ \frac{R_{i}}{A+BL_{i}}=\frac{R_{Cap}}{A+BL_{Cap}} $$

which can be solved for $L_{i}$ to yield the general SPH back-calculation model

$$ L_{i} = \frac{R_{i}}{R_{Cap}}\left(L_{Cap}+\frac{A}{B}\right) - \frac{A}{B} \quad \quad \text{(5)} $$

The linear SPH model produces a family of lines that all have an intercept of $-\frac{A}{B}$ and pass through each observed ($R_{Cap},L_{Cap}$) point (Figure 1). The SPH model is the same as the Fraser-Lee model except that the intercept from Equation 1 is replaced with $-\frac{A}{B}$. The SPH model is the same as the Dahl-Lea model if $A=0$.

The body proportional hypothesis (BPH) model was also named by Francis (1990) and was also first recognized by Whitney and Carlander (1956) when they said "{i}f a fish at time of capture were 10 per cent smaller than the average fish with that size of scale, [then] the fish would be 10 per cent smaller than the expected length for the size of that scale throughout life." This hypothesis can be written as

$$ \frac{L_{i}}{E[L|R_{i}]}=\frac{L_{Cap}}{E[L|R_{Cap}]} $$

If the length-scale relationship is linear then the expected values can be found by plugging $R_{i}$ and $R_{Cap}$ into Equation 1 to get

$$ \frac{L_{i}}{a+bR_{i}}=\frac{L_{Cap}}{a+bR_{Cap}} $$

which can be solved for $L_{i}$ to yield the general BPH back-calculation model

$$ L_{i}=L_{Cap}\frac{a+bR_{i}}{a+bR_{Cap}} \quad \quad \text{(6)} $$

The linear BPH model produces a family of lines that have an intercept of $\frac{aL_{Cap}}{a+bR_{Cap}}$ and pass through each observed ($R_{Cap},L_{Cap}$) point (Figure 1). In contrast to the other back-calculation models, the BPH model uses lines with a different intercept for each fish. The linear BPH model is the same as the Dahl-Lea model if $a=0$.

\

Many Other Models

Vigliola and Meekan (2009) described 18 other models for the back-calculation of fish length. Many of these models are variations of the BPH and SPH models described above, but assuming an other than (geometrically) linear model (e.g., quadratic or other polynomials, or exponential models that are either linearized or fit with non-linear regression). Still other models model the effect of age or time on the the relationship between $L_{Cap}$ and $R_{Cap}$ or $R_{Cap}$ and $L_{Cap}$. Some of these other models have been shown to be more useful when using otoliths, modeling daily increments, or for specific fishes. See Vigliola and Meekan (2009) and Section 11.2.3.1 in Shoup and Michaeletz (2017) for more details.

\

\


References




fishR-Core-Team/RFishBC documentation built on Jan. 30, 2024, 4:42 a.m.