Yes, the density matrix is in general a complex matrix, but it is also hermitian. This means that the imaginary parts of \rho and \rho^\dagger cancel out when summing over all i's and j's, leaving only the real part of \rho (the density matrix) multiplied by the real LCAO basis functions. The DOS is therefore real (has no imaginary components).