Lagged auto-correlations of random_normal not nul about ncl

[yruprich]

Actually, I am facing the same problem with Python (v2.7.9 and v3.7.4):

import numpy as np

N     = 10000000
T     = 100
invT  = -1./(T-1.)
sd    = 1
av    = 0
mxlag = 10

X     = np.random.normal(av, sd, size=(N, T))
acf   = X[:,0:mxlag+1]
for i in range(N):
    acf[i,:] = [1. if l==0 else np.corrcoef(X[i,l:],X[i,:-l])[0][1] for l in range(mxlag+1)]

acf_mean=np.average(acf, axis=0)

print('mean auto-correlation of random_normal vector of length T=',T,' : ',acf_mean)
print('to be compared with -1/(T-1) = ',invT)

[yruprich]

Link to the related Python issue I opened: https://stackoverflow.com/questions/73392044/lagged-auto-correlations-of-numpy-random-normal-not-nul

[yruprich]

Actually this is not a shortcoming of the NCL function. My problem is coming from the bias in the estimate of the auto-correlation. This has been already documented back in 1954...

Reference: Marriott, F. H. C., and J. A. Pope. "Bias in the estimation of autocorrelations." Biometrika 41.3/4 (1954): 390-402 (https://www.jstor.org/stable/2332719)