Wright.crit(n, k, nit)
Wright.crit(n=10,k=2,nit=50) # nit is set to 50 for fast execution in the example. # nit=10000 is recommended as in Wright (2000)
Run the code above in your browser using DataCamp Workspace