We construct Lyapunov functionals for delay differential equation models of infectious diseases in vivo to analyze the stability of the equilibria. The Lyapunov functionals contain the terms that integrate over all previous states. An appropriate evaluation of the logarithm functions in those terms guarantees the existence of the integrals. We apply the rigorous analysis for the one-strain models to multistrain models by using mathematical induction.