  The improved next-to-next-to-leading order (NNLO) QCD analysis of the
experimental data of the CCFR collaboration for the $xF_3$ structure function
is made. Theoretical ambiguities of the NNLO fits are estimated by means of the
Pad\'e resummation technique, which was applied both in the expanded and
non-expanded forms. The NNLO and the new N$^3$LO $\alpha_s(Q^2)$
$\bar{MS}$-matching conditions are used. In the process of the fits we are
taking into account the target mas corrections and the twist-4 $1/Q^2$-terms.
Our NNLO results for $\alpha_s(M_Z)$ values, extracted from the CCFR $xF_3$
data, are $\alpha_s(M_Z)=0.118 \pm 0.002 (stat)\pm 0.005 (syst) \pm
0.003(theory)$ provided the twist-4 contributions are fixed through the
infrared renormalon model and $\alpha_s(M_Z)=0.116^{+ 0.006}_{-0.007}(stat)\pm
0.005 (syst) \pm 0.004(theory)$ provided the twist-4 terms are considered as
the free parameters. It is shown that the extracted at the NNLO order $x$-shape
of the twist-4 correction is almost unchanged after application of the Pad\'e
resummation approximant.
