In this paper, we present a new hybrid high-order method for the Sobolev equation. For a given positive integer k, we adopt the Pk polynomials to approximate the discrete unknown solutions on mesh elements and faces. In our method, we use reconstruction to formulate the gradient term and add suitable stabilization term to enforce the stability and high order accuracy. We analyze the full-discrete formulation of the hybrid high-order method. The stability, energy error estimate with k+1 order and L2 estimate with k+2 order are proved. At last, numerical results show the efficiency of our method.