add nll computation in statstools/__init__