Work on SI likelihood function