Discrete time-multidimensional renewal theory and applications

We develop a discrete-time renewal framework in which renewal events evolve along multiple time coordinates and the sojourn mechanism is described by a general distribution on the multi-index lattice.

The resulting processes, called multi-time renewal chains, are studied through multi-index convolution and the associated algebra of multivariate formal power series.

This algebraic formulation gives explicit representations for multi-time renewal equations, constructive coefficient formulas, and practical inversion schemes.

For computation, we combine FFT-based multidimensional convolution with Newton-type reciprocal iteration to evaluate renewal quantities on large grids.

For asymptotics, we prove strong laws and central limit theorems under proportional growth of the observation horizon, including a general central limit theorem for additive functionals and a Gaussian limit for the renewal counting process in directions with a unique rate-determining coordinate.

We also study fixed-horizon observations: the terminal age vector induces a genuinely multivariate right-censoring mechanism, leading to an exact nonparametric maximum likelihood estimator and its asymptotic normality.

Applications include a binomial--multiset identity, two-attribute warranty evaluation, alternating-renewal availability computation, and discretization-based approximations of continuous-time bivariate renewal and availability models.