У меня есть данные по выживанию (биномиальные) для людей через морфологическое пространство (формы рыбьих черепов). Я хотел бы разработать тест, чтобы увидеть, являются ли локальные области смерти в этом месте больше, чем ожидалось, случайно. Вероятность одной группы наблюдаемых смертей может быть просто вычислена из биномиального распределения (односторонний). Однако меня интересует сравнение с отверстиями одинаковой величины в любом месте пространства, а не только в наблюдаемой позиции (тест с множественными хвостами). Каков наилучший способ закодировать это в R?
Так, например: учитывая 100 выживших из 1000 особей, каждая из которых имеет уникальную позицию в двумерном пространстве, я могу сделать 1000 случайных выборок и рандомизировать, какие особи выжили в каждой, сохраняя постоянным общее число выживших их положение в двумерном пространстве в каждом образце.
В каждой из этих выборок (100 случайных выживших в фиксированном пространстве из 1000 человек) я бы хотел найти самый большой пространственный регион, в котором нет выживших. В простейшем случае эта область может быть ограничена эллипсом. На данный момент, я не уверен, как лучше закодировать эту статистику теста?
Значения размера наблюдаемых отверстий можно затем сравнить с количеством случайных образцов, содержащих отверстия равного или большего размера.
Спасибо за вашу помощь!