Estimados:
En forma resumida, tengo el siguiente problema, quiero validar los modelos esferico, gaussiano y exponencial para un semivariograma al cual se le ha realizado la transformacion box-cox a los datos originales. Alguien me puede dar un ejemplo de como deberia escribir el script para el xvalid del siguiente ejemplo???
variograma <- variog(D, max.dist=200, uvec=30, lambda=0.14146, estimator.type="modulus")
ini.vals_amarillo <- expand.grid(seq(0, 120, by=0.25), seq(0, 40, by=0.2))
vario.wls <- variofit(variograma, ini = ini.vals_amarillo, cov.model = "spherical", fix.nugget = FALSE, nugget = 50, weights = "cressie")
xv.wls <- xvalid(.......................................................)
################################################################################
he probado, pero ninguna ha resultado, enuncio algunas de ellas
·################################################################################
xv.wls1 <- xvalid(D, model = vario.wls)
xv.wls1 <- xvalid(D, coords = D$coords, data = D$data, model = vario.wls)
xv.wls1 <- xvalid(D, coords = D$coords, data = D$data+1, model = vario.wls)
xv.wls1 <- xvalid(D, coords = D$coords, data = D$data+1, model = vario.wls, reestimate=true, +vario.obj=variograma)
xv.wls1 <- xvalid(D, coords = D$coords, data = D$data+1, model = vario.wls, reestimate=true, +vario.obj=variograma), ini = ini.vals_amarillo, cov.model = "gaussian", fix.nugget = FALSE, nugget = 50, weights = +"cressie")
#############################################################################################Todo lo anterior me llevo a la siguiente idea:
#############################################################################################
si transformo mi matriz original D (en otras palabras, transformo por boxcox la data=D$data de mi geodata) y llamo a la nueva matriz Dt entonces:
variograma_t <- variog(Dt, max.dist=200, uvec=30, estimator.type="modulus") #### LA MATRIZ TRANSFORMADA
deberia ser lo mismo que:
variograma <- variog(D, max.dist=200, uvec=30, lambda=0.14146, estimator.type="modulus") # LA MATRIZ ORIGINAL CON TRANSFORMACION BOX-COX
pero resulta que NO.., alguien me puede explica porque??????
la idea iba a que estas dos expresiones deberain ser iguales
vario.wls <- variofit(variograma, ini = ini.vals_amarillo, cov.model = "spherical", fix.nugget = FALSE, nugget = 50, weights = "cressie")
vario.wls_t <- variofit(variograma_t, ini = ini.vals_amarillo, cov.model = "spherical", fix.nugget = FALSE, nugget = 50, weights = "cressie")
por ende DE FORMA DIRECTA:
el xvlid(Dt, vario.wls_t)
AGRADEXZCO CUALQUIER COMENTARIO Y AYUDA
SALUDOS CORDIALES