Update likelihood.R