We study a correspondence between f(R) model of gravity in the Jordan frame and a phenomenological kind of dark energy (DE), which is known as QCD ghost DE. Since this kind of DE is not stable in the context of Einsteinian theory of gravity and Brans–Dicke model of gravity, we consider two kinds of correspondence between modified gravity and DE. By studding the dynamical evolution of model and finding relevant quantities such as, equation of state parameter, deceleration parameter, dimensionless density parameter, we show that the model can describe the present Universe and also the EoS parameter can cross the phantom divide line without needs to any kinetic energy with negative sign. Furthermore, by obtaining the adiabatic squared sound speed of the model for different cases of interaction, we show that this model is stable. Finally, we fit this model with supernova observational data in a noninteraction case and we find the best values of parameter at 1σ confidence interval as; , and . These best-fit values show that DE equation of state parameter, ωd0, can cross the phantom divide line at the present time.