A cluster expansion of the statistical mechanical density operator for a general linear chain model with nearest-neighbor interactions is made. This expansion is then shown to lead to an expansion of a generalized transfer matrix, whose maximum eigenvalue is the per-site partition function. A number of computational features, as well as some illustrative examples, of this approach are described.