A systematic time-dependent perturbation scheme for classical canonical systems is developed based on a Wick's theorem for thermal averages of time-ordered products. The occurrence of the derivatives with respect to the canonical variables noted by Martin, Siggia, and Rose implies that two types of Green's functions have to be considered, the propagator and the response function. The diagrams resulting from Wick's theorem are "double graphs" analogous to those introduced by Dyson and also by Kawasaki, in which the response-function lines form a "tree structure" completed by propagator lines. The implication of a fluctuation-dissipation theorem on the self-energies is analyzed and compared with recent results by Deker and Haake.