[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