The initial data associated with mathematical models or equations that describe physical phenomena, may have errors. It is important to know the effect of these errors on the desired behaviour of the solutions of initial value problems. In this paper, we discuss the continuous dependence of solutions on the initial conditions for nabla fractional difference equations. We also obtain the linear variation of parameters formula for nabla fractional difference equations involving Riemann-Liouville type fractional differences.