Skip to content
/ atsar Public

Applied time series analysis in R with Stan. Allows fast Bayesian fitting of multivariate time-series models.

Notifications You must be signed in to change notification settings

atsa-es/atsar

Repository files navigation

[![atsar status

badge](https://atsa-es.r-universe.dev/badges/atsar)](https://atsa-es.r-universe.dev)[![R-CMD-check](https://github.com/nwfsc-timeseries/atsar/workflows/R-CMD-check/badge.svg)](https://github.com/nwfsc-timeseries/atsar/actions)

[![Codecov test

coverage](https://codecov.io/gh/atsa-es/atsar/branch/master/graph/badge.svg)](https://app.codecov.io/gh/atsa-es/atsar?branch=master)

<style> .nav{ border:1px solid #ccc; border-width:1px 0; list-style:none; margin:0; padding:0; text-align:center; } .nav li{ display:inline-block; } .nav a{ display:inline-block; padding:5px; } </style>

The atsar R package implements Bayesian time series models using Stan, primarily for illustrative purposes and teaching (University of Washington’s Fish 507, Winter quarters). The Stan webpage, and appropriate citation guidelines are here.

INSTALL

Install atsar from the atsa-es r-universe repository (binaries for Windows and Mac-Intel)

install.packages('atsar', repos = c('https://atsa-es.r-universe.dev', 'https://cloud.r-project.org'))

You can build the development version of the package from the source here. Note you need to use this if you have a M1/M2 Mac.

# install.packages("remotes")
remotes::install_github("nwfsc-timeseries/atsar")

EXAMPLE

Simulate data:

library(rstan)
#> Warning: package 'rstan' was built under R version 4.3.2
#> Loading required package: StanHeaders
#> Warning: package 'StanHeaders' was built under R version 4.3.2
#> 
#> rstan version 2.32.6 (Stan version 2.32.2)
#> For execution on a local, multicore CPU with excess RAM we recommend calling
#> options(mc.cores = parallel::detectCores()).
#> To avoid recompilation of unchanged Stan programs, we recommend calling
#> rstan_options(auto_write = TRUE)
#> For within-chain threading using `reduce_sum()` or `map_rect()` Stan functions,
#> change `threads_per_chain` option:
#> rstan_options(threads_per_chain = 1)
library(atsar)
set.seed(123)
s = cumsum(rnorm(50))
plot(s)

Fit several models to this data:

# Regression, no slope
regression_model = fit_stan(y = s, x = model.matrix(lm(s~1)), model_name="regression")

# Regression, with slope
regression_model = fit_stan(y = s, x = model.matrix(lm(s~seq(1,length(s)))), model_name="regression")

# AR(1) time series model
ar1_model = fit_stan(y = s, est_drift=FALSE, P = 1, model_name = "ar")

# ARMA(1,1) time series model
arma1_model = fit_stan(y = s, model_name = "arma11")

# univariate ss model -- without drift but mean reversion estimated
ss_model = fit_stan(y = s, model_name = "ss_ar", est_drift=FALSE)

To see the Stan mode code behind each of these, look in the inst/stan folder on the GitHub repository. Note that fit_stan.R does some data preparation to deal with Stan not accepting NAs in the data.

DOCUMENTATION

  • ATSA lab book - Many applications are covered in our Applied Time Series Analysis book developed from the labs in our course.
  • ATSA course website - We have lectures and all material from our course on our course website.
  • Additional information can be found on the NWFSC time series page which includes several additional books and packages, NWFSC time series page

CITATION

Ward, E.J., M.D. Scheuerell, and E.E. Holmes. 2018. ‘atsar’: Applied Time Series Analysis in R: an introduction to time series analysis for ecological and fisheries data with Stan. DOI

NOAA Disclaimer

This repository is a scientific product and is not official communication of the National Oceanic and Atmospheric Administration, or the United States Department of Commerce. All NOAA GitHub project code is provided on an ‘as is’ basis and the user assumes responsibility for its use. Any claims against the Department of Commerce or Department of Commerce bureaus stemming from the use of this GitHub project will be governed by all applicable Federal law. Any reference to specific commercial products, processes, or services by service mark, trademark, manufacturer, or otherwise, does not constitute or imply their endorsement, recommendation or favoring by the Department of Commerce. The Department of Commerce seal and logo, or the seal and logo of a DOC bureau, shall not be used in any manner to imply endorsement of any commercial product or activity by DOC or the United States Government.

NOAA Fisheries

U.S. Department of Commerce | National Oceanographic and Atmospheric Administration | NOAA Fisheries