  The $^{22}$Ne(p,$\gamma$)$^{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 M$_O$ $\leq$ M
$\leq$ 6 M$_O$). 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 E$_p$ = 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
non-resonant 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$_{tot}^{DC}$(0) = 48.8$\pm$9.5 keV.b compares well with
the result of Ferraro, {\it et al.} and has a lower uncertainty. The resultant
thermonuclear reaction is slightly larger in 0.1 GK $\le$ T $\le$ 0.2 GK
temperature range but otherwise in agreeent with Ferraro, {\it et al.}.
