Let π be a finite or countable set. Given a matrix F = (Fij)i,jβπ of distribution functions on R and a quasistochastic matrix Q = (qij)i,jβπ, i.e. an irreducible nonnegative matrix with maximal eigenvalue 1 and associated unique (modulo scaling) positive left and right eigenvectors u and v, the matrix renewal measure βnβ₯0Qn β F*n associated with Q β F := (qijFij)i,jβπ (see below for precise definitions) and a related Markov renewal equation are studied. This was done earlier by de Saporta (2003) and Sgibnev (2006, 2010) by drawing on potential theory, matrix-analytic methods, and Wiener-Hopf techniques. In this paper we describe a probabilistic approach which is quite different and starts from the observation that Q β F becomes an ordinary semi-Markov matrix after a harmonic transform. This allows us to relate Q β F to a Markov random walk {(Mn, Sn)}nβ₯0 with discrete recurrent driving chain {Mn}nβ₯0. It is then shown that renewal theorems including a Choquet-Deny-type lemma may be easily established by resorting to standard renewal theory for ordinary random walks. The paper concludes with two typical examples.