Свертка плотностей

Question

Если x - это данные и плотность $ hist дает мне эмпирическую плотность, как получить свертку x с самим собой. Convolve() дает мне много возможностей, но я не уверен, что использовать. Мне нужна функция для этого: http://en.wikipedia.org/wiki/Convolution#Discrete_convolution (на самом деле просто реальный случай будет достаточно)

Простой тест: свертка Bernoullies дает мне биномиальное

Pr = с (0,7, 0,3)

Правильный ответ должен быть биномиальными параметры п = 2 р = 0,3

Answer 1

Хорошо правильный ответ:

> convolve(Pr,rev(Pr),type="o") 
[1] 0.49 0.42 0.09 

бинома параметров п = 2 и р = 0,3. Таким образом, чтобы свернутые плотности можно использовать:

convolve(his$density, rev(his$density), type="o") 

Это хорошо работает с дискретными распределениями, но может работать (очень) плохо для непрерывных распределений.

Примечания: если и хочет использовать FDR cumsum() в результате свертки

Предложение: для непрерывной таблицы распределения (х)/суммы (таблица (х)) дает более точный вход для свертки