In the present work, we discuss the differentiability properties of solutions of nabla fractional difference equations of order $\alpha$ ($0<\alpha<1$) with respect to the initial conditions. Further, we develop a nonlinear variation of parameters formula to obtain the solution of a perturbed nonlinear nabla fractional difference equation.