用 si 表示 yi 的位置,其中通过反映空间或时间的坐标指定 si。通常通过假定协方差是欧几里得距离 dij(si 和 sj 之间)的函数,限制空间或时间结构。该协方差可以表示为 ,其中 表示观测值 yi 和 yj 之间的相关性。
下面显示 JMP 中可用的空间模型的相关性结构。这些结构由参数 ρ 表示,除非另行约束,否则该参数值为正数。
•
|
•
|
•
|
•
|
若过程是各向同性的,则半方差仅依赖于两个点之间的距离 h 且函数可以表示为:
定义为截距。它表示 h = 0 处的跳跃不连续性。
定义为在达到更大距离的平稳值处的半变差图的值。它对应于观测值的方差。在不带 nugget 效应的模型中,sill 为 。在带 nugget 效应的模型中,sill 为 ,其中 c1 表示 nugget。偏 sill 定义为 。
定义为半变差图达到 sill 处的距离。在小于极差的距离处,观测值在空间上是相关的。对于大于等于极差的距离,空间相关性实际上为零。在球形模型中,ρ 为极差。在指数模型中,3ρ 为实际极差。在高斯模型中, 为实际极差。实际极差定义为协方差下降到 95% 的 sill 处的距离。
在“拟合混合模型”报表 - 带 Nugget 的空间球体中,重复效应协方差参数估计值表示各种半变差图功能:
极差 ρ 的估计值。
对于给定的各向同性空间结构,使用非线性最小二乘法将观测数据拟合为空间相关性结构中的合适函数,以得到估计的变差图。
要计算经验半方差,需要计算为变差图协方差选择的变量的所有点对之间的距离。距离的极差分为 10 个相等的区间。若数据不允许 10 个区间,则构造尽可能多的区间。
Ch
包含其距离位于第 h 个最大区间的点的距离分类
Z(x)
x 处的响应值,其中 x 是时间或空间坐标的向量
g(h)
按以下方式定义半方差函数 γ: