[R-sig-Geo] spatstat predict rhohat question

[Rolf Turner]

Thank you for the *excellently* reproducible example.  We are looking 
into it and will get back to you.

cheers,

Rolf Turner

-- 
Technical Editor ANZJS
Department of Statistics
University of Auckland
Phone: +64-9-373-7599 ext. 88276
Home phone: +64-9-480-4619


On 22/02/15 17:06, js_wvu15 wrote:
> 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!