Необходимо указать значение ковариаты «x» (в качестве аргумента для функции effectfun).

Я читаю текст из учебника spatstat и пытаюсь изучить соответствие модели с помощью ppm. Я создал модель с картографическими координатами в качестве ковариант. Затем я хотел увидеть влияние только одной ковариаты на модель,

model1 = ppm(chicago_ppp ~ x+y) plot(effectfun(model1, covname = "x"))

но я получаю ошибку

Error in effectfun(model1, covname = "x") : A value for the covariate “y” must be provided (as an argument to effect fun)

То же самое происходит, если я использую covname «y», он запрашивает «x»

Может кто-нибудь, пожалуйста, покажите мне, в чем моя ошибка. Спасибо.

ОБНОВЛЕНИЕ: когда я использую только одну ковариату, и я использую effectfun с этой одной ковариатой, ошибки не возникает. Когда я использую две ковариаты и хочу проверить эффект от одной ковариаты, я получаю эту ошибку в вопросе.


spatstat
person BKS    schedule 22.06.2017    source источник
comment
Не забудьте сделать образец кода воспроизводимым, чтобы его можно было легко скопировать в сеанс R. Вы не определили chicago_ppp, и вы также можете напрямую использовать некоторый встроенный набор данных из spatstat.   -  person Ege Rubak    schedule 23.06.2017

Ответы (1)


arrow_upward
0
arrow_downward


Чтобы иметь возможность рассчитать предполагаемую интенсивность для различных значений x, вам необходимо исправить значение для y следующим образом:

    library(spatstat)
    model <- ppm(cells ~ x + y)
    plot(effectfun(model, covname = "x", y = 0.1))


    plot(effectfun(model, covname = "x", y = 0.9))

person Ege Rubak    schedule 23.06.2017
comment
Спасибо. Это ответ на мой вопрос. Я улучшу свои следующие вопросы с помощью воспроизводимого кода. - person BKS; 26.06.2017