Explicit computation of stats if not provided
Fixed readme to reflect current code
Added LN option
Added outside nonlinearity fitting routine