Abstract In this paper, we propose a new type of univariate and bivariate Integer-valued autoregressive model order one (INAR(1)) to approximate linear birth death process with constant rates. Under specific parametric setting, the dynamic transition probabilities probability generating function INAR(1) will converge that as length subintervals goes 0. Due simplicity Markov structure, maximum l...