Стоит задача сгладить частоты в многомерном пространстве. Возможная постановка такова.
Есть двумерное пространство параметров

, в каждой точке которой можно вычислить значение частоты

.
На этом пространстве введем сетку из значения параметров

,

,

. Для этих значений параметров имеем матрицу значений частот

.
Необходимо найти оценки вероятностей для этих частот, удовлетворяющие условиям:

и

.
Одномерная задача решается методом isotonic regression.
Однако:
1) Решение получается ступенчатой функцией по

, что из практических соображений (выделенных значений

нет), и из соображений гладкости (или минимизации риска) не очень хорошо.
2) Для двумерного и многомерного случая обобщение и "разложение" через одномерное решение не очевидно.
3) Для многомерного случая возникает вопрос об эффективности алгоритма (кол-во точек в пространстве становится: 100 на 100 на 90 на 3).
Что делать? Какие книжки этому вопросу посоветуете?
Текущее решение:
к 2), 3) Использую одномерное решение. Существование других подпространств пока игнорирую.
к 1) Собираюсь сгладить 4-х точечные ступеньки с сильным скачком полиномом третьей степени, интерполяцией по крайним точкам и фитом по двум другим.