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。非經特殊聲明,原始代碼版權歸原作者所有,本譯文未經允許或授權,請勿轉載或複製。