Стандартные ошибки начальной загрузки модели кластерных точечных процессов (kppm)

Я хотел бы сообщить о стандартной ошибке параметров кластеризации (каппа, сигма) неоднородной модели точечного процесса Томаса, которую я использовал в spatstat. Юэ и Ло (2015) сообщили, что делали это с помощью параметрической начальной загрузки. Я не очень разбираюсь в этой концепции и не применяю ее к моделям точечных процессов. Как мне это сделать?

Мое первое предположение - смоделировать свое kppm несколько раз и заново подогнать полученные смоделированные точки с теми же ковариатами. Затем вычислите стандартные ошибки из параметров кластеризации каждой последующей подгонки. Это верно? Если да, то сколько симуляций будет считаться приемлемым в этом случае? Заранее благодарим за любые указатели!


r statistics-bootstrap spatial spatstat
person KK Li    schedule 18.11.2015    source источник

Ответы (1)


arrow_upward
0
arrow_downward

По сути, ваше собственное описание полностью верно.

Мое первое предположение - смоделировать свое kppm несколько раз и заново подогнать полученные смоделированные точки с теми же ковариатами. Затем вычислите стандартные ошибки из параметров кластеризации каждой последующей подгонки.

Остается только один вопрос - сколько симуляций нужно сделать. В основном ответ такой: «Столько, сколько успеешь сделать!». Часто можно увидеть, как люди выполняют 1000 симуляций, так почему бы вам не начать с этого?

person Ege Rubak    schedule 18.11.2015
comment
Спасибо за помощь! Есть ли функция в spatstat или r, которая могла бы помочь в этом? Я подумал, что может быть функция конверта, но после прочтения документации кажется, что пользователь может предоставлять только функции, которые вычисляются по r, как Kest ... - person KK Li; 18.11.2015
comment
В настоящее время в spatstat нет функции для этого. Я не знаю, есть ли в R какие-либо пакеты, которые упростят это, чем просто написать для этого свой собственный скрипт. - person Ege Rubak; 21.11.2015