Under certain assumptions on the dependence structure of the residual lives of the insured (i.e. their independence, positive association, or negative association), in this paper we establish some laws of large numbers for the convex upper bounds, derived by the technique of comonotonicity, of the present value function of a homogeneous portfolio composed of the whole-life insurance policies.

Two methods for approximating the limiting distribution of the present value of the benefits of a portfolio of identical endowment insurance contracts are suggested. The model used assumes that both future lifetimes and interest rates are random. The first method is similar to the one presented in Parker (1994b). The second method is based on the relationship between temporary and endowment insurance contracts.

Two approaches used to model interest randomness are presented. They are the modeling of the force of interest accumulation function and the modeling of the force of interest. The expected value, standard deviation and coefficient of skewness of the present value of annuities-immediate are presented as illustrations. The implicit behavior of the force of interest under the two approaches is investigated by looking at a particular conditional expectation of the force of interest accumulation function.

An approximation of the distribution of the present value of the benefits of a portfolio of temporary insurance contracts is suggested for the case where the size of the portfolio tends to infinity. The model used is the one presented in Parker (1922b) and involves random interest rates and future lifetimes. Some justifications of the approximation are given. Illustrations for limiting portfolios of temporary insurance contracts are presented for an assumed Ornstein-Uhlenbeck process for the force of interest.