Площадь перекрытия двух круговых гауссовских функций - PullRequest
2 голосов
/ 14 марта 2012

Я пытаюсь написать код, который находит перекрытие между трехмерными фигурами.

Каждая фигура определяется двумя пересекающимися нормальными распределениями (одно в направлении x, другое в направлении y).

Есть ли у вас какие-либо предложения существующего кода, который решает этот вопрос или функцииЯ могу использовать для создания этого кода?Большая часть моего опыта программирования была на языке R, но я открыт для решений и на других языках.

Заранее благодарю за любые предложения и помощь!


Более длительный контекст исследования по этому вопросу: я изучаю использование акустического пространства насекомыми.Я хочу знать, будут ли случайно собранные группы насекомых иметь вызовы, которые более или менее похожи, чем мы наблюдаем в естественных сообществах (тест рандомизации).Для этого мне нужно случайным образом выбрать виды насекомых и рассчитать сходство между их призывами.

Для каждого вида у меня есть среднее значение и дисперсия для двух характеристик вызова, которые приблизительно нормально распределены.Я хотел бы использовать эти две характеристики вызова для построения трехмерного распределения вероятностей для вида.Затем я хотел бы рассчитать сумму, на которую PDF для одного вида перекрывается с другим.

Пожалуйста, примите мои извинения, если вопрос не ясен или не подходит для данного форума.

1 Ответ

1 голос
/ 13 апреля 2012

Я работаю над открытием низкомолекулярных лекарств и часто использую программу (ROCS от OpenEye Scientific Software), основанную на алгоритмах, которые представляют молекулы в виде наборов сферических гауссовых функций и вычисляют объемы пересечений.Вы можете посмотреть следующие ссылки, а также документацию ROCS:

(1) Grant and Pickup, J. Phys.Химреагент1995, 99, 3503-3510

(2) Грант, Галлардо и Пикап, J. Comp.Химреагент1996, 17, 1653-1666

...