Chaotic Communications have drawn a great deal of attention to the wireless communication industry and research due to its limitless meritorious features, including excellent anti-fading and anti-intercept capabilities and jamming resistance exempli gratia. Differential Chaos Shift Keying (DCSK) is of particular interest due to its low-complexity and low-power and many attractive properties. However, most of the DCSK studies reported in the literature considered the additive white Gaussian noise environment in non-cooperative scenarios. Moreover, the analytical derivations and evaluation of the error rates and other performance metrics are generally left in an integral form and evaluated using numerical techniques. To circumvent on these issues, this work is dedicated to present a new approximate error rates analysis of multi-access multiple-input multiple-output dual-hop relaying DCSK cooperative diversity (DCSK-CD) in Nakagami-m fading channels (enclosing the Rayleigh fading as a particular case). Based on this approximation, closed-form expressions for the average error rates are derived for multiple relaying protocols, namely the error-free and the decode-and-forward relaying. Testing results validate the accuracy of the derived analytical expressions.