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@r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es