For a fixed $a \in \{1, 2, 3, \ldots\},$ the radius of starlikeness of positive order is obtained for each of the normalized analytic functions \begin{align*} \mathtt{f}_{a, ν}(z)&:= \bigg(2^{a ν-a+1} a^{-\frac{a(aν-a+1)}{2}} Γ(a ν+1) {}_a\mathtt{B}_{2a-1, a ν-a+1, 1}(a^{a/2} z)\bigg)^{\tfrac{1}{a ν-a+1}},\\ \mathtt{g}_{a, ν}(z)&:= 2^{a ν-a+1} a^{-\frac{a}{2}(aν-a+1)} Γ(a ν+1) z^{a-aν} {}_a\mathtt{B}_{2a-1, a ν-a+1, 1}(a^{a/2} z),\\ \mathtt{h}_{a, ν}(z)&:= 2^{a ν-a+1} a^{-\frac{a}{2}(aν-a+1)} Γ(a ν+1) z^{\frac{1}{2}(1+a-aν)} {}_a\mathtt{B}_{2a-1, a ν-a+1, 1}(a^{a/2} \sqrt{z}) \end{align*} in the unit disk, where ${}_a\mathtt{B}_{b, p, c}$ is the generalized Bessel function \begin{align*} {}_a\mathtt{B}_{b, p, c}(z):= \sum_{k=0}^\infty \frac{(-c)^k}{k! \; \mathrmΓ{\left( a k +p+\frac{b+1}{2}\right)} } \left(\frac{z}{2}\right)^{2k+p}. \end{align*} The best range on $ν$ is also obtained for a fixed $a$ to ensure the functions $\mathtt{f}_{a, ν}$ and $\mathtt{g}_{a, ν}$ are starlike of positive order in the unit disk. When $a=1,$ the results obtained reduced to earlier known results.

15 pages