Я использую pmvnorm для получения тройного интегрирования для многомерной функции нормальной плотности в R. Диапазон интегрирования - от отрицательной бесконечности до конечных чисел для всех трех переменных. Однако иногда функция возвращает отрицательное значение, которое не имеет смысла.
Я смотрел другие посты, где люди предлагали заменить значение на ноль или повысить точность алгоритма. Это не работает в моем случае, однако.
После получения вероятности от pmvnorm, скажем, x, мне нужно преобразовать ее в log (x), чтобы выполнить оценку максимального правдоподобия. Поэтому я не могу просто заменить отрицательный х на ноль. Есть ли другое решение в этом вопросе? Могу ли я просто заменить -Inf некоторым значением, например -100, в случае нормальной плотности, где среднее значение не превышает 10, а дисперсия / ковариация также не превышает 10?
Заранее спасибо.