I have some data and a model, e.g.,
$H_0: \xi \sim \mathcal N(\mu_1 > \mu_0, \sigma_0)$, $\mu_0, \sigma_0$ are fixed,
in other words, this is the similar to the situation with "unequal means, alternative=greater".
How can I calculate the likelihood of the data under this model? Numerical integration could help, but often this integral is divergent (upper limit is infinity). I can put some prior on possible $\mu_1$ values, but it is another problem.