We present autoregressive (AR) and autoregressive moving average (ARMA) processes with bivariate exponential (BE) and bivariate geometric (BG) distributions. The theory of positive dependence is used to show that in various cases, the BEAR, BGAR, BEARMA, and BGARMA models consist of associated random variables. We discuss special cases of the BEAR and BGAR processes in which the bivariate processes are stationary and have well-known bivariate exponential and geometric distributions. Finally, we fit a BEAR model to a real data set.
