Как ksdensity вычисляет каждую точку данных xi?

Когда я использую [f,xi] = ksdensity(x) в Matlab, я получаю оценку плотности вероятности, f и точки оценки xi, в которых ksdensity вычисляет f.

У меня вопрос: как рассчитывается/определяется каждая точка xi? Есть ли формула?

В центре документации говорится: По умолчанию используется 100 равноотстоящих точек, которые охватывают диапазон данных в x. Итак, они охватывают диапазон, но это не объясняет, как рассчитываются.

Большое спасибо!

Хуан


person Juan    schedule 06.10.2013    source источник
comment
Я полагаю xi(i)=min(x)+i*(max(x)-min(x))/100   -  person pseudoDust    schedule 06.10.2013


Ответы (1)


Стандартный метод получения равноотстоящих точек в MATLAB — использование команды linspace. linspace(a,b,n) генерирует n линейно расположенных точек между a и b включительно.

Так что это, скорее всего, эквивалентно xi = linspace(min(x),max(x)) (количество точек по умолчанию 100 inlinspace).

person ramanan    schedule 06.10.2013
comment
Большое спасибо @pseudoDust и ramanan, но ни один из ваших ответов не работает. Если я использую в качестве примера: Data=[365629;258213;30986;340841;82628;92957;10329;454454;0;51643;309855;82628;51643;10329;41314;41314;330512;0;92957;227227;123942;433797;423469;464783;826280;433797;206570;433797;61971;123942;82628;0;113614;960551;557739;113614;113614]; [PDF,xi]=ksdensity(Data,'kernel','epanechnikov','npoints',10); xi_ramanan= linspace(min(Data),max(Data),10); for i=1:100; xi_pseudoDust(i)=min(Data)+i*(max(Data)-min(Data))/100; end; результаты не совпадают с xi ksdensity =/ - person Juan; 07.10.2013
comment
Ok. Теперь яснее, что вы показали свои данные. Аномалия связана с наличием 0s. Из-за них ваши данные будут частично перекрываться с отрицательными значениями, что приведет к тому, что xi не будет изначально присутствовать в Data. Более подробное объяснение приведено на stackoverflow.com/questions/ 17410155/ - person ramanan; 07.10.2013
comment
Спасибо! Я читал эту ветку раньше, но мой вопрос остается прежним: как рассчитывается каждый xi? =П - person Juan; 07.10.2013