In this paper, we investigate the rate of convergence of the solution $u_\varepsilon$ of the random elliptic partial difference equation $(\nabla^{\varepsilon *} a(x/\varepsilon,\omega)\nabla^\varepsilon+1)u_\varepsilon(x,\omega)=f(x)$ to the corresponding homogenized solution. Here $x\in\varepsilon Z^d$, and $\omega\in\Omega$ represents the randomness. Assuming that $a(x)$'s are…