Calculating Moran's I on large point dataframe in R

I have a dataset of about 77,000 points of which I want to test autocorrelation. I have tried to compute the test statistic using the testSpatialAutocorrelation tool in the DHARMa package, which uses “a simulation-based approach to create readily interpretable scaled (quantile) residuals for fitted (generalized) linear mixed models.” However when I try to run the test, I run into a memory error. What is a way to easily calculate Moran’s I for a large dataset?