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 (
markbermej...@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@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es
[[alternative HTML version deleted]]
___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es