Create a new variable **normal_z**$=\Phi^{-1}(\frac{i-0.5}{2000}), i=1,2,\cdots,2000$, where $\Phi^{-1}$ is the quantile/inverse function of a standard normal distribution