Оценка относительного риска в спатстате

Question

У меня возникают проблемы при вычислении оценки относительного риска (relrisk.ppp) двух точечных моделей: одна с четырьмя отметками в прямоугольной области angular, а другая с двумя отметками в круглой области.

Для первого шаблона с четырьмя отметками я могу получить относительный риск и полученный объект в виде large imlist с 4 элементами, соответствующими каждой отметке.

Однако для второго шаблона он дает список из 10 элементов, из которых первая матрица v пуста с NA записями. Я ломаю голову над тем, что может быть не так, когда созданные объекты точечного рисунка кажутся идентичными. Любая помощь будет оценена. Спасибо.

Answer 1

Для вашего первого набора данных результатом является список объектов изображения (список из четырех объектов класса im). Для вашего второго набора данных результатом relrisk.ppp является одно изображение (объект класса im). Это поведение по умолчанию, когда существует только два возможных типа точек (два возможных значения меток). См. help(relrisk.ppp).

Во всех случаях вы просто должны иметь возможность plot и print полученного объекта. Вам не нужно проверять внутренние данные изображения.

Дополнительные пояснения: , когда существует только два возможных типа точек, поведение по умолчанию relrisk.ppp состоит в их обработке. в качестве данных case-control , где точки, относящиеся к первому типу, рассматриваются как контрольные (например, неинфицированные люди), а точки второго типа - как случаи (например, инфицированные люди). Отношение интенсивностей (делений, разделенных на контроли) оценивается как изображение.

Если вы не хотите, чтобы это произошло , задайте аргумент casecontrol=FALSE, а затем relrisk.ppp всегда будет возвращать список изображений с одним изображением для каждой возможной отметки. Каждое изображение дает пространственно изменяющуюся вероятность точки этого типа.

Все это объяснено в help(relrisk.ppp) или в книге.