fitNull function

Fit the null distribution