Statistics for Cosmology with HI intensity mapping: effect of higher order corrections