### Kuvaus

Alkuperäiskieli | englanti |
---|---|

Lehti | Annals of Applied Statistics |

Vuosikerta | 11 |

Numero | 3 |

Sivut | 1375-1402 |

Sivumäärä | 28 |

ISSN | 1932-6157 |

DOI - pysyväislinkit | |

Tila | Julkaistu - syyskuuta 2017 |

OKM-julkaisutyyppi | A1 Alkuperäisartikkeli tieteellisessä aikakauslehdessä, vertaisarvioitu |

### Tieteenalat

- 112 Tilastotiede

### Lainaa tätä

*Annals of Applied Statistics*,

*11*(3), 1375-1402. https://doi.org/10.1214/17-AOAS1029

}

*Annals of Applied Statistics*, Vuosikerta 11, Nro 3, Sivut 1375-1402. https://doi.org/10.1214/17-AOAS1029

**Gaussian processs framework for temporal dependence and discrepancy functions in Ricker-type population growth models.** / Hartmann, Marcelo; Hosack, Geoffrey R; Hillary, Richard M; Vanhatalo, Jarno Petteri.

Tutkimustuotos: Artikkelijulkaisu › Artikkeli › Tieteellinen › vertaisarvioitu

TY - JOUR

T1 - Gaussian processs framework for temporal dependence and discrepancy functions in Ricker-type population growth models

AU - Hartmann, Marcelo

AU - Hosack, Geoffrey R

AU - Hillary, Richard M

AU - Vanhatalo, Jarno Petteri

PY - 2017/9

Y1 - 2017/9

N2 - Density dependent population growth functions are of central importance to population dynamics modelling because they describe the theoretical rate of recruitment of new individuals to a natural population. Traditionally these functions are described with a fixed functional form with temporally constant parameters and without species interactions. The Ricker stock-recruitment model is one such function that is commonly used in fisheries stock assessment. In recent years, there has been increasing interest in semi-parametric and temporally varying population growth models. The former are related to the general statistical approach of using semi-parametric discrepancy functions, such as Gaussian processes (GP), to model deviations of data around the expected parametric function. In the latter, the reproductive rate, which is a key parameter describing the population growth rate, is assumed to vary in time. In this work, we introduce how these existing Ricker population growth models can be formulated under the same statistical approach of hierarchical GP models. We also show how the time invariant semi-parametric approach can be extended and combined with the time varying reproductive rate using a GP model. Then we extend these models to the multispecies setting by incorporating cross-covariances among species with a continuous time covariance structure using the linear model of coregionalization. As a case study, we examine the productivity of three Pacific salmon populations. We compare the alternative Ricker population growth functions using model posterior probabilities and leave-one-out cross validation predictive densities. Our results show substantial temporal variation in maximum reproductive rates and reveal temporal dependence among the species, which have direct management implications. However, our results do not support inclusion of semi-parametric discrepancy function and they suggest that the semi-parametric discrepancy functions may lead to challenges in parameter identifiability more generally.

AB - Density dependent population growth functions are of central importance to population dynamics modelling because they describe the theoretical rate of recruitment of new individuals to a natural population. Traditionally these functions are described with a fixed functional form with temporally constant parameters and without species interactions. The Ricker stock-recruitment model is one such function that is commonly used in fisheries stock assessment. In recent years, there has been increasing interest in semi-parametric and temporally varying population growth models. The former are related to the general statistical approach of using semi-parametric discrepancy functions, such as Gaussian processes (GP), to model deviations of data around the expected parametric function. In the latter, the reproductive rate, which is a key parameter describing the population growth rate, is assumed to vary in time. In this work, we introduce how these existing Ricker population growth models can be formulated under the same statistical approach of hierarchical GP models. We also show how the time invariant semi-parametric approach can be extended and combined with the time varying reproductive rate using a GP model. Then we extend these models to the multispecies setting by incorporating cross-covariances among species with a continuous time covariance structure using the linear model of coregionalization. As a case study, we examine the productivity of three Pacific salmon populations. We compare the alternative Ricker population growth functions using model posterior probabilities and leave-one-out cross validation predictive densities. Our results show substantial temporal variation in maximum reproductive rates and reveal temporal dependence among the species, which have direct management implications. However, our results do not support inclusion of semi-parametric discrepancy function and they suggest that the semi-parametric discrepancy functions may lead to challenges in parameter identifiability more generally.

KW - 112 Statistics and probability

KW - Marginal likelihood

KW - Model evidence

KW - Population growth

KW - Fisheries

KW - Density dependence

KW - Temporal dependence

KW - Interspecific dependence

U2 - 10.1214/17-AOAS1029

DO - 10.1214/17-AOAS1029

M3 - Article

VL - 11

SP - 1375

EP - 1402

JO - Annals of Applied Statistics

JF - Annals of Applied Statistics

SN - 1932-6157

IS - 3

ER -