Hola Marcos.
Utilicé hace mucho tiempo la función krige, pero en este caso me temo que
el error está en los datos que estás suministrando.
Comprueba que en "plantilla.sp" está la columna llamada RadPcp, porque
el
error lo que dice es que no encuentra ese objeto (primero mira a ver si es
una variable de los datos que suministras, y al no encontrarla, lo intenta
como objeto y tampoco lo encuentra). Creo, aunque sin unos datos para
comprobarlo no estoy seguro, que el problema viene en este fragmento de
código:
plantilla.sp <- as.data.frame(radar.mx)[,1:2]
coordinates(plantilla.sp) <- ~x+y
class(plantilla.sp)
donde creo que te estás dejando la columna RadPcp fuera.
Espero que te sirva. ¡Un saludo!
Víctor Granda García
PhD. Plant Biology
Dpto BOS, Área de Fisiología Vegetal
Universidad de Oviedo
El mar., 21 jul. 2015 a las 11:31, Marcos Bermejo (<
markbermejo90 en hotmail.com>) escribió:
> Hola,
>
> Estoy intentando hacer un kriging con deriva externa con la función
> "krige()", pero me da un error:
> Error in eval(expr, envir, enclos) : object 'radPcp' not found
> El código principal es:
>
> radar.mx <- rasterToPoints(radar.cu)
> colnames(radar.mx) <-
c("x","y","radPcp")
> radar.spdf <- as.data.frame(radar.mx)
>
> lluvia.pluv <- as.data.frame(t(rbind(pluvio.h[1,log.NA],
> coord.aux[,log.NA])))
> colnames(lluvia.pluv) <- c("pluvPcp", "x",
"y")
> lluvia.pluv.spdf <- lluvia.pluv
> coordinates(lluvia.pluv.spdf) <- ~x+y
>
> #radarpx <- as(radar.spdf, "SpatialPixelsDataFrame")
>
>
> plantilla.sp <- as.data.frame(radar.mx)[,1:2]
> coordinates(plantilla.sp) <- ~x+y
> class(plantilla.sp)
>
>
> for (v in 1:length(pluvio.h[1,log.NA])) {
> lluvia.pluv.df <- data.frame(lluvia.pluv.spdf) [-v,]
> lluvia.pluv.spdf <- lluvia.pluv.df
> coordinates(lluvia.pluv.spdf) <- ~x+y
>
> # Hacemos el KED. Ver función "krige()":
> KED.rad <- krige(
> formula=pluvPcp~radPcp, # covariable -> radar
> locations=lluvia.pluv.spdf,
> newdata=plantilla.sp, # podría ser cualquier objeto
> Spatial
> model=v.fit, # modelo de semivariograma.
> maxdist=5000
> )
>
>
> ¿Me podríais ayudar?
>
> Un saludo,
>
> [[alternative HTML version deleted]]
>
> _______________________________________________
> R-help-es mailing list
> R-help-es en r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
[[alternative HTML version deleted]]