One- and two-body densities in position space have been calculated for the atomic beryllium isoelectronic series starting from explicitly correlated multideterminant wave functions. The effects of electronic correlations have been systematically studied by comparing the correlated results with the corresponding Hartree-Fock ones. Some expectation values such as , < r(n)>, , < r(12)(n)>, , and < R-n>, where r, r(12), and R stand for the electron-nucleus, interelectronic, and two electron center of mass coordinates, respectively, have been obtained. All the calculations have been carried out by using the Monte Carlo algorithm.