Variogram.gls
位于 nlme
包(package)。 说明
此方法函数根据 gls
拟合计算残差的半变异函数。如果存在分组因子,则计算同一组水平内的残差对的半变异函数值。如果 collapse
与 "none"
不同,则使用 Cressie (1993) 中定义的稳健估计器 (robust = TRUE
) 或相同距离间隔内的值的平均值来折叠各个半变异函数值。半变异函数对于误差项相关结构建模非常有用。
用法
## S3 method for class 'gls'
Variogram(object, distance, form, resType, data,
na.action, maxDist, length.out, collapse, nint, breaks,
robust, metric, ...)
参数
object |
继承自类 |
distance |
一个可选的数值向量,其中包含残差对之间的距离。如果存在分组变量,则仅应给出同一组内残差对之间的距离。如果缺失,则根据参数 |
form |
一个可选的单边公式,指定用于计算残差对之间的距离的协变量,以及可选的用于对残差进行分区的分组因子(必须出现在 |
resType |
一个可选字符串,指定要使用的残差类型。如果 |
data |
一个可选数据帧,用于解释 |
na.action |
一个函数,指示当数据包含 |
maxDist |
用于计算两个残差之间的半变异函数的最大距离的可选数值。默认情况下,包括所有残差对。 |
length.out |
一个可选的整数值。当 |
collapse |
一个可选字符串,指定要应用于各个半变异函数值的折叠类型。如果等于 |
nint |
一个可选整数,表示折叠半变异函数值时要使用的间隔数。默认为 |
robust |
一个可选的逻辑值,指定在折叠各个值时是否应使用稳健的半变异函数估计量。如果 |
breaks |
一个可选的数值向量,其中包含用于折叠半变异函数值的距离间隔的断点。如果不丢失,则忽略 |
metric |
指定要使用的距离度量的可选字符串。当前可用的选项是 |
... |
该泛型的某些方法需要额外的参数。此方法中没有使用任何内容。 |
值
DataFrame ,其中 variog
和 dist
列分别表示半变异函数值和相应的距离。如果半变异函数值折叠,则返回的 DataFrame 中将包含一个额外的列 n.pairs
,其中包含每个半变异函数计算中使用的残差对的数量。如果 object
包含 corSpatial
元素,则返回值中将包含一个数据帧及其相应的半变异函数,作为属性 "modelVariog"
。返回值继承自类 Variogram
。
例子
fm1 <- gls(weight ~ Time * Diet, BodyWeight)
Vm1 <- Variogram(fm1, form = ~ Time | Rat)
print(head(Vm1), digits = 3)
作者
José Pinheiro and Douglas Bates bates@stat.wisc.edu
参考
Cressie, N.A.C. (1993), "Statistics for Spatial Data", J. Wiley & Sons.
也可以看看
gls
, Variogram
, Variogram.default
, Variogram.lme
, plot.Variogram
相关用法
- R Variogram.corSpher 计算 corSpher 对象的半变异函数
- R Variogram.corGaus 计算 corGaus 对象的半变异函数
- R Variogram.lme 计算 lme 对象残差的半变异函数
- R Variogram.corSpatial 计算 corSpatial 对象的半变异函数
- R Variogram.corLin 计算 corLin 对象的半变异函数
- R Variogram.corExp 计算 corExp 对象的半变异函数
- R Variogram.default 计算半变异函数
- R Variogram.corRatio 计算 corRatio 对象的半变异函数
- R Variogram 计算半变异函数
- R VarCorr 提取方差和相关分量
- R Pixel X 射线像素强度随时间的变化
- R corARMA ARMA(p,q) 相关结构
- R getGroupsFormula 提取分组公式
- R corRatio 有理二次相关结构
- R logLik.glsStruct glsStruct 对象的对数似然
- R intervals.lmList lmList 系数的置信区间
- R corLin 线性相关结构
- R plot.augPred 绘制 augPred 对象
- R print.varFunc 打印 varFunc 对象
- R recalc 重新计算压缩线性模型对象
- R getGroups.lme 提取 lme 对象组
- R nlmeStruct 非线性混合效应结构
- R predict.nlme 来自 nlme 对象的预测
- R corSymm 一般相关结构
- R qqnorm.gls gls 对象残差的正态图
注:本文由纯净天空筛选整理自R-devel大神的英文原创作品 Calculate Semi-variogram for Residuals from a gls Object。非经特殊声明,原始代码版权归原作者所有,本译文未经允许或授权,请勿转载或复制。