script for computing hte likelihood by internally marginalizing over tub