[R-es] Duda R

[Marcos Bermejo]

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]]