Modelado geoestadístico en R, con ejemplo de Kriging Ordinario y variograma omnidireccional. Los datos se pueden bajar en el siguiente enlace mega.nz/#!Uc8S...
Muchísimas gracias Maximiliano por este vídeo, ha sido muy informativo... sin embargo, quisiera hacerte unas preguntas sobre los procedimientos que aplicaste en el mismo. 1.- ¿Por qué los datos deben seguir una distribución normal para poder aplicar el Kriging? 2.- ¿Es posible realizar una validación cruzada en R empleando para ello todos los modelos de forma automática?, es decir, que la validación me permita elegir el mejor modelo que se adapte a mis datos mediante esta validación, así como los parámetros, rango, meseta y pepita... 3.- ¿Dentro del gstat es posible realizar otros tipos de interpolación, como el IDW, o Splines? 4.- Una última pregunta... ¿Es posible que la validación me arroje no sólo los valores del RMSE, sino también, los de R2 del ajuste del modelo y el MAE, o valor medio absoluto?.. Disculpa tantas preguntas, es que estoy estudiando geoestadística como parte de mi formación como Geomático, y me llama poderosamente la atención, peeero, tengo muchas dudas....
Miguel, que dicha que te sirvió el video. Antes de responder las preguntas te aclaro que para nada soy un experto, sino mas bien un entusiasta en la materia. 1.- No es que tienen que seguir una distribución normal, se podría hacer Kriging sobre datos no normales, es solo que como base de la estadística, se asumen que los procesos generadores son normales y por lo tanto los procedimientos y métodos están basados en este supuesto. Se recomienda que sean mas o menos normales por este criterio 2.- Lamentablemente no hay forma de realizar la validación cruzada automáticamente, ya que la validación va a depender del modelo y los parámetros para cada modelo, se podría hacer una función propia o hacer un loop para generar lo que quieres, pero no hay nada pre-programado 3.- Dentro de gstat se puede hacer IDW pero no splines, ve la documentación del comando 'idw' 4.- De hecho la validación no te arroja ningún valor, el RMSE lo calculé yo por mi cuenta, asi que se podrían calcular los valores que quieras mientras tengas la formula Saludos, y espero algo de esto te ayude
Hola maximiliano muy buen video, sin embargo tengo dudas, en el ejemplo de tu base de datos , por ejemplo para unvalor de potasio le sacas el logaritmo cno sus coordenadas, en mi caso tengo un conjunto de datos de un sensor de humedad para un dia, que valor deberia ir junto a la coordenada del sensor, la media, el que mas se repite, ojala me ayudes , gracias, saludos desde Perú
Hola, disculpa la tardanza. En esos casos seria mejor usar la media, o se me ocurre escoges una hora y escoges el valor de esa hora para todos los puntos
al introducir el comando : (p2=spplot(ok,"var1.se",asp=1,col.regions=cm.colors(80), main="desviacion", scales=list(draw=T),xlab="x[m]",ylab="y[m]")) me arroja: Error in `[.data.frame`(obj@data, zcol) : undefined columns selected sabes cual podría ser el problema?
Buenas, estoy trabajando con unos datos de resistencia a la penetración (r_mpa) y humedad volumétrica del suelo (HV). Cuando voy a aplicar la función mapview, me tira un error que dice "Error: at least one of the following columns: r_mpa, log_mpa, HV, log_HV not found in x[[i]]"
Roberto Clemente Aguirre Elizondo ok, burst=T es mejor usarlo cuando se está traficando una variable cualitativa y quieres una capa por cada categoría, entonces hay que especificar el argumento zcol = “variable” que quieres graficar. Si la variable es continua no es recomendable usar burst
Hola Maximiliano! Excelente video. Siga adelante con el buen trabajo!
Hola, muchas gracias por compartir este video quizás podrías ayudarme con una tutoría (podría escribirte por favor) ??
Muchísimas gracias Maximiliano por este vídeo, ha sido muy informativo... sin embargo, quisiera hacerte unas preguntas sobre los procedimientos que aplicaste en el mismo.
1.- ¿Por qué los datos deben seguir una distribución normal para poder aplicar el Kriging?
2.- ¿Es posible realizar una validación cruzada en R empleando para ello todos los modelos de forma automática?, es decir, que la validación me permita elegir el mejor modelo que se adapte a mis datos mediante esta validación, así como los parámetros, rango, meseta y pepita...
3.- ¿Dentro del gstat es posible realizar otros tipos de interpolación, como el IDW, o Splines?
4.- Una última pregunta... ¿Es posible que la validación me arroje no sólo los valores del RMSE, sino también, los de R2 del ajuste del modelo y el MAE, o valor medio absoluto?..
Disculpa tantas preguntas, es que estoy estudiando geoestadística como parte de mi formación como Geomático, y me llama poderosamente la atención, peeero, tengo muchas dudas....
Miguel, que dicha que te sirvió el video. Antes de responder las preguntas te aclaro que para nada soy un experto, sino mas bien un entusiasta en la materia.
1.- No es que tienen que seguir una distribución normal, se podría hacer Kriging sobre datos no normales, es solo que como base de la estadística, se asumen que los procesos generadores son normales y por lo tanto los procedimientos y métodos están basados en este supuesto. Se recomienda que sean mas o menos normales por este criterio
2.- Lamentablemente no hay forma de realizar la validación cruzada automáticamente, ya que la validación va a depender del modelo y los parámetros para cada modelo, se podría hacer una función propia o hacer un loop para generar lo que quieres, pero no hay nada pre-programado
3.- Dentro de gstat se puede hacer IDW pero no splines, ve la documentación del comando 'idw'
4.- De hecho la validación no te arroja ningún valor, el RMSE lo calculé yo por mi cuenta, asi que se podrían calcular los valores que quieras mientras tengas la formula
Saludos, y espero algo de esto te ayude
Maximiliano, muy buen video, gracias. Maximiliano, me queda la duda de que es ese tres que pones al final del comando dint. Saludos y muchas gracias.
Hola maximiliano muy buen video, sin embargo tengo dudas, en el ejemplo de tu base de datos , por ejemplo para unvalor de potasio le sacas el logaritmo cno sus coordenadas, en mi caso tengo un conjunto de datos de un sensor de humedad para un dia, que valor deberia ir junto a la coordenada del sensor, la media, el que mas se repite, ojala me ayudes , gracias, saludos desde Perú
Hola, disculpa la tardanza. En esos casos seria mejor usar la media, o se me ocurre escoges una hora y escoges el valor de esa hora para todos los puntos
Excelente video, disculpa el archivo no esta disponible podrías subirlo te agradecería mucho.
Hola, no seria mas exacto utilizar un variograma direccional? que se ajuste en base a las anisotropias
Buenas, se podría resolver una función f() de variables variables, y después comparar el resultado con la respuesta analítica?
al introducir el comando : (p2=spplot(ok,"var1.se",asp=1,col.regions=cm.colors(80),
main="desviacion",
scales=list(draw=T),xlab="x[m]",ylab="y[m]"))
me arroja: Error in `[.data.frame`(obj@data, zcol) : undefined columns selected
sabes cual podría ser el problema?
Yo creo la variable “var1.se” después de realizar el kriging, fíjate de que la hayas creado antes de intentar graficarla
Buenas, estoy trabajando con unos datos de resistencia a la penetración (r_mpa) y humedad volumétrica del suelo (HV). Cuando voy a aplicar la función mapview, me tira un error que dice "Error: at least one of the following columns: r_mpa, log_mpa, HV, log_HV not found in x[[i]]"
podes compartir el comando que estas intentando correr? estas especificando 'zcol'?
@@maxgav13 Mi set de datos se llama mpa y el comando es mapView(mpa,burst=T, hide=T)
Roberto Clemente Aguirre Elizondo ok, burst=T es mejor usarlo cuando se está traficando una variable cualitativa y quieres una capa por cada categoría, entonces hay que especificar el argumento zcol = “variable” que quieres graficar. Si la variable es continua no es recomendable usar burst
Hi, is there a way How to use negative exponential covariogram
Petr Hrobař hi, I haven’t looked into this but you could check all of the variogram models available in the show.vgms() function
El enlace esta roto.. no hay archivos
Profe, podria reponer los datos del enlace? Gracias
puedes pasar la base de datos
y el scritp
hola, poderia corrigir o link? Nao esta baixando..., belo vídeo
Buen día Maximiliano,
¿Es posible tener acceso al script?
La funcion, coordinates (datos), me da error, 'no encuentra objeto x'
Villalbaruben 11 “x” tiene que ser el nombre de la columna que tiene la coordenada “x” de los datos, igual para “y”