rwg.j.sim function

Simulate rwg(j) values from a random null distribution