空间和时间变异性
考虑简单模型 。 通过误差项 ei 对空间或时间结构建模。通常,空间相关性模型可以定义为
si 表示 yi 的位置,其中通过反映空间或时间的坐标指定 si。通常通过假定协方差是欧几里得距离 dijsi sj 之间)的函数,限制空间或时间结构。该协方差可以表示为 ,其中 表示观测值 yiyj 之间的相关性。
在两个或更多位置坐标的情况下,若 f(dij) 不依赖于方向,则协方差结构为各向同性。若依赖于方向,则该结构为各向异性
空间相关性结构
下面显示 JMP 中可用的空间模型的相关性结构。这些结构由参数 ρ 表示,除非另行约束,否则该参数值为正数。
当空间过程为二阶平稳时,空间相关性结构中列出的结构定义变差图。从地质统计学借鉴而来,变差图是用于描述和估计空间变异性的标准工具。它使用半方差测量作为两个观测值之间的距离 dij 的函数的空间变异性。
Z(s) 表示位置 s 处的响应值。按以下方式给出 sisj 处的两个观测值之间的半方差
定义为截距。它表示 h = 0 处的跳跃不连续性。
定义为半变差图达到 sill 处的距离。在小于极差的距离处,观测值在空间上是相关的。对于大于等于极差的距离,空间相关性实际上为零。在球形模型中,ρ 为极差。在指数模型中,3ρ 为实际极差。在高斯模型中, 为实际极差。实际极差定义为协方差下降到 95% 的 sill 处的距离。
“拟合混合模型”报表 - 带 Nugget 的空间球体中,重复效应协方差参数估计值表示各种半变差图功能:
极差 ρ 的估计值。
变差图估计值
经验半方差
要计算经验半方差,需要计算为变差图协方差选择的变量的所有点对之间的距离。距离的极差分为 10 个相等的区间。若数据不允许 10 个区间,则构造尽可能多的区间。
构造包含点对的距离分类。第 h 个距离分类包含其距离位于第 h 个区间的所有点对。
Ch
包含其距离位于第 h 个最大区间的点的距离分类
Z(x)
x 处的响应值,其中 x 是时间或空间坐标的向量
g(h)
距离分类 Ch 的半方差