In this paper, we consider a class of backward doubly stochastic differential equations (BDSDE for short) with general terminal value and general random generator. Those BDSDEs do not involve any forward diffusion processes. By using the techniques of Malliavin calculus, we are able to establish the $L^p$-Hölder continuity of the solution pair. Then, an implicit numerical scheme for the BDSDE is proposed and the rate of convergence is obtained in the $L^p$-sense. As a by-product, we obtain an explicit representation of the process $Y$ in the solution pair to a linear BDSDE with random coefficients.