24小时热门版块排行榜    

查看: 6362  |  回复: 25
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

大傻meier

新虫 (正式写手)

[求助] 折磨我很久的 地统计学分析,望大神指教 已有3人参与

在约52×100m的范围内按地形等差异布置了100个点,取0-10cm土样,做土壤CEC的空间变异性。
经纬度坐标在ARCGIS里投影成平面坐标,数据踢除异常值,做正态性转换。
在GS+里,参数Active lag distance(样点间最大距离的一半)和Lag class distance interval(约等于样点间最小距离)采用软件默认值,未做更改。得到的半变异函数模型为高斯模型,具体参数见表1,半变异函数图1。
在Gis里,做趋势分析,数据南北方向呈倒U型趋势(图2)。半变异云图,系统默认的为图3,资料里说“步长*步长组的值应该在最大的样点对距离的1/2左右,Lag size可由平均最近邻点距离求得(Spatial Statistics toolbox—Analyzing Patterns---Average Nearest Neighbor tool)”,故将Lag Size和number of lags更改为图4中的数值。
在地统计学向导里,插值方法选择Ordinary Kring—Prediction,移除2次趋势。更改Lag Size 和number of lags后,软件生成的半方差模型为图5,线性模型,和GS+里得到的结果很不一样。
我把GS+得到的模型参数输入,得到图6,但是Prediction error里标准均方根预测误差达到12了(图7),明显错了,

折磨我很久了,论文卡在这里,干着急上火,不知道可以问谁。没有学过GIS,不知道我这操作过程有没有什么错误。地统计学分析到底该怎么做嘞,求各位指教,我如果有说的不明白的,我再解释。谢谢,所有金币奉上

折磨我很久的 地统计学分析,望大神指教
表1.PNG


折磨我很久的 地统计学分析,望大神指教-1
图1.PNG


折磨我很久的 地统计学分析,望大神指教-2
图2.PNG


折磨我很久的 地统计学分析,望大神指教-3
图3.PNG


折磨我很久的 地统计学分析,望大神指教-4
图4


折磨我很久的 地统计学分析,望大神指教-5
图5.PNG


折磨我很久的 地统计学分析,望大神指教-6
图6.PNG


折磨我很久的 地统计学分析,望大神指教-7
图7.PNG
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

truebelief

木虫 (著名写手)

都是去除2次趋势再对残差建模的吗?lag size大了,导致没有一个明显上升的semivariogram。还有不要把sill和nugget都定死了,这样根本就没有办法fit semivariogram。

» 本帖已获得的红花(最新10朵)

5楼2014-10-26 09:19:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 26 个回答

keguigiser

新虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
大傻meier: 金币+20, 有帮助, 谢谢 2014-11-09 13:21:13
不要修改Lag Size,默认的试一试。

[ 发自手机版 http://muchong.com/3g ]
肯定自己,超越自己,创造自己
2楼2014-10-25 13:14:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

大傻meier

新虫 (正式写手)

引用回帖:
2楼: Originally posted by keguigiser at 2014-10-25 13:14:38
不要修改Lag Size,默认的试一试。

默认的是这个,也是线性模型。
explore data里的半变异云图和Geostatistical Wizard里的这一步,怎么得到的semivariogram图和Lag size也不一样呢
折磨我很久的 地统计学分析,望大神指教-8
0.PNG

3楼2014-10-25 16:50:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

keguigiser

新虫 (著名写手)

【答案】应助回帖

引用回帖:
3楼: Originally posted by 大傻meier at 2014-10-25 16:50:58
默认的是这个,也是线性模型。
explore data里的半变异云图和Geostatistical Wizard里的这一步,怎么得到的semivariogram图和Lag size也不一样呢

0.PNG
...

是一样的啊,4.17*12就是整个半变异的横坐标啊

[ 发自手机版 http://muchong.com/3g ]
肯定自己,超越自己,创造自己
4楼2014-10-25 18:38:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +37 医学老男孩 2026-08-13 87/4350 2026-08-18 06:35 by 医学老男孩
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[基金申请] 时间戳变了,能看出什么问题? +4 基诺咪客 2026-08-17 4/200 2026-08-17 16:10 by Vivilian
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +12 Ldrop2023 2026-08-13 15/750 2026-08-17 15:02 by 小豌豆_发芽
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 6/300 2026-08-16 20:07 by 奶牛小黑
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +10 majunge000 2026-08-11 12/600 2026-08-16 08:18 by xli1984
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
信息提示
请填处理意见