The 22 Ne(p,γ) 23 Na reaction in NeNa cycle plays an important role in the production of only stable sodium isotope 23 Na. This nucleus is processed by the NeNa cycle during hot bottom burning (HBB) in asymptotic giant branch (AGB) stage of low metallicity intermediate mass stats (4 MO ≤ M ≤ 6 MO). Recent measurements have addressed the uncertainty in the thermonuclear reaction rate of this reaction at relevant astrophysical energies through the identification of low lying resonances at Ep = 71,105, 156.2, 189.5 and 259.7 keV. In addition, precise measurements of low energy behaviour of the non-resonant capture has also been performed and the contribution of the sub-threshold resonance at 8664 keV excitation in 23 Na has been established. Here, in this article, we have presented a systematic R-matrix analysis of direct capture to the bound states and the decay of the sub-threshold resonance at 8664 keV to the ground state of 23 Na. A finite range distorted wave Born approximation (FRDWBA) calculation has been performed for 22 Ne( 3 He,d) 23 Na transfer reaction data to extract the asymptotic normalization coeeficients (ANC-s) required to estimate the nonresonant capture cross sections or astrophysical S-factor values in R-matrix analysis. Simultaneous R-matrix analysis constrained with ANC-s from transfer calculation reproduced the astrophysical S-factor data over a wide energy window. The S DC tot (0) = 48.8±9.5 keV.b compares well with the result of Ferraro, et al. and has a lower uncertainty. The resultant thermonuclear reaction is slightly larger in 0.1 GK ≤ T ≤ 0.2 GK temperature range but otherwise in agreeent with Ferraro, et al..