This paper deals with error bounds for numerical solutions of linear ordinary differential equations by global or piecewise polynomial collocation methods which are based on consideration of the involved differential operator, related matrices and the residual. It is shown that significant improvement may be obtained if direct bounds for the error in the solution are considered. The practical implementation of the theory is illustrated by a selection of numerical examples.