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 distribu