Mauro,
Não sei o que pode estar ocorrendo. Estou rodando sem problemas na v3.02.
Em quinta-feira, 26 de dezembro de 2013, Mauro Sznelwar<sznelwar@uol.com.br> escreveu:
>
> Eu não consegui baixar a bibliotéca sptsatat, diz que não está disponível na versão R 3.0.1 do Windows.
>
> Alexandre, bom dia!
> Não tenho experiência com a {spatstat}, mas tentei trabalhar seu código e acredito que possa lhe ser útil.
>
> ### <BEGIN>
> require(raster); require(spatstat); require(sp)
> ## Raster
> -----------------------------------------------------------
> set.seed(20130622)
> r <- raster(nc=100, nr=100)
> r <- setValues(r, round(runif(ncell(r))* 255))
> plot(r)
> r.spdf <- data.frame(coordinates(r), px=getValues(r)); names(r.spdf)
> gridded(r.spdf) <- ~x+y
> spplot(r.spdf)
> ## Amostras
> -------------------------------------------------
> parc.p <- spsample(r.spdf, 15, type="random")
> parc.c <- vector("list",length(parc.p))
> for(i in 1:length(parc.p)) {
> parc.c[[i]] <- disc(radius=5, centre=c(parc.p$x[i], parc.p$y[i]))}
> W <- do.call(union.owin,parc.c)
> W.df <- as.data.frame(as.im(W))
> W.sp <- W.df; coordinates(W.sp) <- ~x+y
> plot(W.sp)
> W.df$px <- over(W.sp, r.spdf)
> head(W.df)
> #
> # x y value px
> #1 -179.9856 -39.83543 1 11
> #2 -179.9856 -39.29907 1 76
> #3 -179.9856 -38.76271 1 76
> #4 -179.9856 -38.22635 1 76
> #5 -179.9856 -37.68998 1 34
> #6 -178.7263 -41.98087 1 157
> #
> ### <END>
> .
--