相关文章推荐
苦恼的白开水  ·  mvnrepository.com无法访问咋 ...·  3 周前    · 
苦恼的白开水  ·  Hive之窗口函数 - Lucas_zhao ·  1 月前    · 
苦恼的白开水  ·  windows关闭java进程·  8 月前    · 
苦恼的白开水  ·  [麦肯锡]：行至山顶：看见中国职场女性领导者 ...·  9 月前    · 
苦恼的白开水  ·  单细胞工具箱|singleR-单细胞类型自动 ...·  9 月前    · 
苦恼的白开水  ·  SQLSERVER ...·  9 月前    · 
苦恼的白开水  ·  DataGrip PGSQL ...·  9 月前    · 
苦恼的白开水  ·  Android开发：不会ANR？这里有ANR ...·  10 月前    · 
苦恼的白开水  ·  Torch not compiled ...·  10 月前    · 
苦恼的白开水  ·  CUDA C 和 C++ 的简单介绍 - ...·  10 月前    · 
憨厚的冲锋衣  ·  好玩的网站“彩蛋”们（程序员视角）...图形 ...·  5 小时前    · 
傻傻的毛豆  ·  程序员鸡翅的个人空间-程序员鸡翅个人主页-视频·  5 小时前    · 
茫然的春卷  ·  程序员Sunday的个人空间-程序员Sund ...·  5 小时前    · 
满身肌肉的炒饭  ·  强烈推荐程序员常用网站 - ·  5 小时前    · 
豁达的单车  ·  “错误：在关键字'CONVERT'附近有语法 ...·  5 小时前    · 
威武的香瓜  ·  [转]Entity Framework ...·  5 小时前    · 
爱笑的大脸猫  ·  在线直播-中华人民共和国民政部·  6 小时前    · 
博学的火锅  ·  武威市人民政府 各区县信息公开工作年度报告 凉州区·  6 小时前    · 
礼貌的土豆  ·  政府工作报告 厦门市思明区人民政府·  6 小时前    · 
面冷心慈的马克杯  ·  两会&上外｜黄绮：为“星星的孩子”谋教育福利·  6 小时前    · 
Code  ›  [R-sig-Geo] How to perform a variogram model in gstat?
https://stat.ethz.ch/pipermail/r-sig-geo/2012-September/016153.html
苦恼的白开水
11 月前
Tue Sep 25 14:13:37 CEST 2012
  • Previous message: [R-sig-Geo] How to perform a variogram model in gstat?
  • Next message: [R-sig-Geo] NaN / NA values in SpatialPixelDataFrame conversion
  • Messages sorted by: [ date ] [ thread ] [ subject ] [ author ] About your question, the variogram fit depend on what you want!!!, if you just want to fit a model to describe structures you can choose many options like eyefit, Index GOF (Goodness of fit), AIC (Akaike) or RSS to decide which parameters, models and values select, If you want an interpolate data with kriging then you would use a cross validation to compares values. Regards, Carolina Lang A. Investigador ----- Mensaje original ----- De: "Julia Chacon Labella" < juliachacon at gmail.com > Para: r-sig-geo at r-project.org Enviados: Martes, 25 de Septiembre 2012 4:57:08 GMT -04:00 Georgetown Asunto: [R-sig-Geo] How to perform a variogram model in gstat? Please excuse these basic questions, but I am a bit lost in geostatistics. I am trying to make a kriging of the soil properties of a plot. I am already using gstat with R and I don't have any problems with the scripts and so on...but, I don't really know how to choose the nugget, psill and range for the variogram model. Moreover, I am not sure about how to include the anisotropy in the model, I have already include a linar model in some cases, to remove any trend. I am following the steps recomended in the gstat help, but I am a bit befuddled. Let me show you my scripts (the psill, range and nugget are choosen just to make the test). coordinates(corg)= c("x", "y") v = variogram(SOC~1, corg) SSErr.a <- numeric(4) names(SSErr.a) <- c("Exponencial","Esférico"," Gausiano","Lineal") v.fit = fit.variogram(v, vgm(psill=0.18, model="Exp", range=5, nugget=0.1)) SSErr.a[1] = attr(v.fit, "SSErr") v.fit = fit.variogram(v, vgm(psill=0.18, model="Sph", range=5, nugget=0.1)) SSErr.a[2] = attr(v.fit, "SSErr") v.fit = fit.variogram(v, vgm(psill=0.18, model="Gau", range=5, nugget=0.1)) SSErr.a[3] = attr(v.fit, "SSErr") v.fit = fit.variogram(v, vgm(psill=0.18, model="Lin", range=5, nugget=0.1)) SSErr.a[4] = attr(v.fit, "SSErr") SSErr.a v = variogram(SOC~x+y, corg) SSErr.b <- numeric(4) names(SSErr.b) <- c("Exponencial","Esférico","Gausiano","Lineal") v.fit = fit.variogram(v, vgm(psill=0.18, model="Exp", range=5, nugget=0.1)) SSErr.b[1] = attr(v.fit, "SSErr") v.fit = fit.variogram(v, vgm(psill=0.18, model="Sph", range=5, nugget=0.1)) SSErr.b[2] = attr(v.fit, "SSErr") v.fit = fit.variogram(v, vgm(psill=0.18, model="Gau", range=5, nugget=0.1)) SSErr.b[3] = attr(v.fit, "SSErr") v.fit = fit.variogram(v, vgm(psill=0.18, model="Lin", range=5, nugget=0.1)) SSErr.b[4] = attr(v.fit, "SSErr") SSErr.b Thanks a lot, Cheers. Julia Chacon Labella [PhD candidate] Área de Biodiversidad y Conservación Escuela Superior de Ciencias Experimentales y Tecnología Universidad Rey Juan Carlos c/ Tulipán, s/n. E-28933 Móstoles (Madrid) España Julia [[alternative HTML version deleted]] _______________________________________________ R-sig-Geo mailing list R-sig-Geo at r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo ------------------------------------------------------------------------------------------------ Certificacion ISO 9001/2000: Sistema de Datos Bio-Pesqueros (Iquique, Valparaiso y Talcahuano) Acreditacion INN - ISO 17025: Determinacion de la edad de jurel y merluza comun
  • Previous message: [R-sig-Geo] How to perform a variogram model in gstat?
  • Next message: [R-sig-Geo] NaN / NA values in SpatialPixelDataFrame conversion
  • Messages sorted by: [ date ] [ thread ] [ subject ] [ author ]
    •  
    推荐文章
    苦恼的白开水  ·  mvnrepository.com无法访问咋回事？ -
    3 周前
    苦恼的白开水  ·  Hive之窗口函数 - Lucas_zhao
    1 月前
    苦恼的白开水  ·  windows关闭java进程
    8 月前
    苦恼的白开水  ·  [麦肯锡]：行至山顶：看见中国职场女性领导者 - 发现报告
    9 月前
    苦恼的白开水  ·  单细胞工具箱|singleR-单细胞类型自动注释（含数据版） -
    9 月前
    苦恼的白开水  ·  SQLSERVER where条件加判断语句_sqlserverwhere条件中的判断_难自拔于世界之大的博客
    9 月前
    苦恼的白开水  ·  DataGrip PGSQL 执行计划_datagrip查看执行计划_LaLaLa_OvO的博客
    9 月前
    苦恼的白开水  ·  Android开发：不会ANR？这里有ANR解析和案例！开发者社区
    10 月前
    苦恼的白开水  ·  Torch not compiled with CUDA enabled - PyTorch Forums
    10 月前
    苦恼的白开水  ·  CUDA C 和 C++ 的简单介绍 - NVIDIA 技术博客
    10 月前
    憨厚的冲锋衣  ·  好玩的网站“彩蛋”们（程序员视角）...图形类彩蛋文字类彩蛋开发者社区
    5 小时前
    傻傻的毛豆  ·  程序员鸡翅的个人空间-程序员鸡翅个人主页-视频
    5 小时前
    茫然的春卷  ·  程序员Sunday的个人空间-程序员Sunday个人主页-视频
    5 小时前
    满身肌肉的炒饭  ·  强烈推荐程序员常用网站 -
    5 小时前
    豁达的单车  ·  “错误：在关键字'CONVERT'附近有语法错误。”
    5 小时前
    威武的香瓜  ·  [转]Entity Framework 异常： 'OFFSET' 附近有语法错误。\r\n在 FETCH 语句中选项 NEXT 的用法无效 - Ace001
    5 小时前
    爱笑的大脸猫  ·  在线直播-中华人民共和国民政部
    6 小时前
    博学的火锅  ·  武威市人民政府 各区县信息公开工作年度报告 凉州区
    6 小时前
    礼貌的土豆  ·  政府工作报告 厦门市思明区人民政府
    6 小时前
    面冷心慈的马克杯  ·  两会&上外｜黄绮：为“星星的孩子”谋教育福利
    6 小时前
    今天看啥   ·   Py中国   ·   codingpro   ·   藏经阁   ·   小百科   ·   link之家   ·   卧龙AI搜索
    Code - 代码工具平台
    删除内容请联系邮箱 2879853325@qq.com
    © 2024 ~ 沪ICP备11025650号