We used QCD factorization for the hadronic matrix elements to show that the existing data, in particular the branching ratios BR ( ?J/?K) and BR ( ?J/??), can be accounted for this approach. We analyzed the decay within the framework of QCD factorization. We have complete calculation of the relevant hard-scattering kernels for twist-2 and twist-3. We calculated this decays in a special scale ( ) and in two schemes for Wilson coefficients in NLO. We considered three functions for . The twist-3 contribution involves logarithmically divergent integral and also we considered canceling divergent. The obtained results are in agreement with available experimental data.