[R-sig-Geo] spatstat predict rhohat question

js_wvu15 jsheehan06 at gmail.com
Sun Feb 22 05:06:21 CET 2015


Hello,

I'm using spatstat 1.40-0 in R version 3.1.2 and encountered this warning
for predict.rhohat when the rhohat object is a fitted point process model
and the window is irregular:

Warning message:
In Y * lambda :
  longer object length is not a multiple of shorter object length

I managed to duplicate my warning:

library(spatstat)
X <- rpoispp(function(x,y){exp(3+3*x)})
win1 <- owin(poly=list(list(x=c(0,1,1,0), y=c(0,0,1,1)),
	list(x=c(0.6,0.4,0.4,0.6), y=c(0.2,0.2,0.4,0.4))))
X1 <- X[win1]
fit1 <- ppm(X1, ~x)
rho1 <- rhohat(fit1, "y")
plot(predict(rho1)) ###get the warning

I researched the warning but couldn't figure it out. Any help is greatly
appreciated!

Jim Sheehan
West Virginia, USA



--
View this message in context: http://r-sig-geo.2731867.n2.nabble.com/spatstat-predict-rhohat-question-tp7587812.html
Sent from the R-sig-geo mailing list archive at Nabble.com.



More information about the R-sig-Geo mailing list