# Actual residuals versus simulated residuals for testing spatial autocorrelation

I have a dataset of observations at multiple sites repeated within and at multiple times. I need to test for spatial autocorrelation in the residuals of my model, which I have used time as a fixed effect for. Do I use means per site of my actual residuals extracted from the model for a Moran's Index test by the ape package, or should I use simulated residuals calculated by the recalculateResiduals function from DHARMa (since the individual points don't have unique x,y values)? Any examples I've seen have used DHARMa, but basically I'm not sure why you would use simulated residuals instead of actual residuals? Thanks.