Мне нужно рассчитать меру, называемую взаимной информацией. Прежде всего, мне нужно вычислить другую меру, называемую энтропией, например, совместную энтропию х и у:Рассчитать объем по графику оценки двумерной плотности ядра
-∬p(x,y)·log p(x,y)dxdy
Таким образом, чтобы вычислить p(x,y)
, я использовал оценщик плотности ядра (в этом случае, функция kde2d
, и она вернула Z
значения (вероятность наличия х и у в этом окне).
Опять же, теперь у меня есть матрица Z
значений [1x100] x [1x100]
, это равно мой p(x,y)
. Но я должен интегрировать его, открыв объем под поверхностью (интеграл doble), но я не нашел способа сделать это. Функция quad2d
, чтобы вычислить двойной квадратурной не работает, потому что я только интегрированный числовой матрицы p(x,y)
, и это дает мне постоянный ....
Каждый знает, что-то найти, что объем/вычислить двойной интеграл?
Изображение сюжета из persp3d
:
Спасибо всем !!!!
Было бы проще, если бы я дискретировал _a priori_ и подсчитал частоту на ширине бункеров? Будет ли это производить гораздо больше ошибок? Благодарю. –
Я тестировал ваш метод против гистограммы, и что-то неправильное может произойти. Когда я тестирую ваш метод на 1000 случайных выборок (как вы объясните выше), энтропия дает мне значения около нуля. Но, по сравнению с теорией, это неверно, потому что это должно дать мне log (n). Я использовал какой-то другой пакет энтропии, и он дает правильный результат ~ log (n). Вы знаете, что происходит? Благодаря! –
Спасибо! Это справедливо для гауссовских нормальных ядер, но как насчет неопределенных плотностей? Тест на нормальную плотность должен давать хорошие результаты, но в моем случае при тестировании на финансовых временных рядах (не стохастическом) эти результаты значительно изменились, особенно из-за полосы пропускания (defaut - гауссовский). Я думаю, что у меня проблемы с хвостами дистрибутива. Проблемы, связанные с ошибками, состоят в том, что у меня есть обширная матрица для вычисления в попарно (временные ряды), поэтому я буду рассматривать большую визуализацию. Еще раз спасибо за помощь! –