We construct a kind of preserving positivity numerical methods for stochastic CIR equation based on the stochastic DG method. Our proposed methods inherit the positive property of the exact solution. We also prove that the methods are convergent with the strong order 1. Finally, some numerical experiments are given to verify our theoretical results.