Normal probability plot

Function demo()
    Variable numPts = 100
    Make/o/n=(numPts) data 
    data = 10 + gnoise(1)
    statsNormPlot(data)
End
 
//==========================================================================================
// The following function produces a normal probability plot equivalent to plotting
// on Gaussian probability paper. based on a modification of statsProbPlot()
Function/S statsNormPlot(inWave)
    Wave inWave
    
    String outWinName=""
    Duplicate/O inWave,tmpWave,xWave
    Sort tmpWave,tmpWave
    Variable num=numpnts(inWave)
    xWave=(p-0.3175)/(num+0.365)
    xWave=StatsInvNormalCDF(xWave,0,1)
    Variable/G V_fitOptions=4
    CurveFit/Q line xWave /X=tmpWave /D
    Wave fit_xWave
    Display/K=1 /W=(5,45,470,405)  xWave vs tmpWave
    outWinName=UniqueName("WM_PROB_PLOT",6,0)
    DoWindow/C $outWinName
    AppendToGraph fit_xWave
    ModifyGraph mode(xWave)=4, lstyle(fit_xWave)=2, rgb(fit_xWave)=(0,0,0)
    TextBox/N=text0/F=0/A=MC/X=-26.90/Y=47.62 "Normal Probability plot" 
    Make/O ticVal = {0.0001, 0.001, 0.003, 0.01, 0.02, 0.05, 0.10, 0.25, 0.50, 0.75, 0.90, 0.95, 0.98, 0.99, 0.997, 0.999, 0.9999}
    Make/O/T ticStr = {"0.0001", "0.001", "0.003", "0.01", "0.02", "0.05", "0.10", "0.25", "0.50", "0.75", "0.90", "0.95", "0.98", "0.99", "0.997", "0.999", "0.9999"}
    Duplicate/O ticVal ticLoc
    ticLoc = StatsInvNormalCDF(ticVal[p],0,1)
    ModifyGraph userticks(left)={ticLoc,ticStr}
    Label bottom "Data Values"
    Label left "Probability"
 
    return outWinName
End