2017-06-22 92 views
0

我正在閱讀spatstat教科書中的文本,並嘗試使用ppm學習模型擬合。 我創建了carteisan座標作爲協變量的模型。然後我想看看在模型上只有一個協變量的影響,必須提供協變量「x」的值(作爲effectfun的參數)

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」

有人可以告訴我什麼是我的錯誤。謝謝。

更新:當我只使用一個協變量,並且使用那個協變量的效應函數時,沒有錯誤。當我使用兩個協變量時,我想檢查一個協變量的效應函數,我在這個問題中得到這個錯誤。

+0

請記住讓代碼示例具有可重現性,因此可以輕鬆地將粘貼複製到R會話中。你還沒有定義'chicago_ppp',你可以直接使用'spatstat'中的一些內置數據集。 –

回答

0


爲了能夠計算的估計強度爲 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)) 

+0

謝謝。這回答了我的問題。我會用可重複的代碼改進我的下一個問題。 – BKS

相關問題