In this paper we define the fractional Cox-Ingersoll-Ross process as Xt :=is a fractional Brownian motion with an arbitrary Hurst parameter H ∈ (0, 1). We prove that Xt satisfies the stochastic differential equation of the form dXt = (k − aXt)dt + σ √ Xt • dB H t , where the integral with respect to fractional Brownian motion is considered as the pathwise Stratonovich integral. We also show that for k > 0, H > 1/2 the process is strictly positive and never hits zero, so that actually Xt = Y 2 t . Finally, we prove that in the case of H < 1/2 the probability of not hitting zero on any fixed finite interval by the fractional Cox-Ingersoll-Ross process tends to 1 as k → ∞.