當前位置: 首頁>>代碼示例 >>用法及示例精選 >>正文


R Variogram.gls 計算 gls 對象殘差的半變異函數


R語言 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

繼承自類 "gls" 的對象,表示廣義最小二乘擬合模型。

distance

一個可選的數值向量,其中包含殘差對之間的距離。如果存在分組變量,則僅應給出同一組內殘差對之間的距離。如果缺失,則根據參數 formdatametric 的值計算距離,除非 object 包含 corSpatial 元素,在這種情況下,關聯的協變量(通過使用getCovariate 方法)。

form

一個可選的單邊公式,指定用於計算殘差對之間的距離的協變量,以及可選的用於對殘差進行分區的分組因子(必須出現在 form 中的 | 運算符的右側) )。默認值為 ~1 ,表示使用組內的觀察順序來獲取距離。

resType

一個可選字符串,指定要使用的殘差類型。如果 "response" ,則使用 "raw" 殘差(觀察到的擬合值);否則,如果 "pearson" ,則使用標準化殘差(原始殘差除以相應的標準誤差);否則,如果 "normalized" ,則使用歸一化殘差(標準化殘差預乘以估計誤差相關矩陣的逆平方根因子)。使用參數的部分匹配,因此隻需提供第一個字符。默認為 "pearson"

data

一個可選數據幀,用於解釋 form 中的變量。默認情況下,使用與 object 相同的數據。

na.action

一個函數,指示當數據包含 NA 時應該發生什麽。如果存在任何不完整的觀察結果,默認操作 ( na.fail ) 會導致打印錯誤消息並終止函數。

maxDist

用於計算兩個殘差之間的半變異函數的最大距離的可選數值。默認情況下,包括所有殘差對。

length.out

一個可選的整數值。當object 包含corSpatial 元素時,將計算其半變異函數值,並將該參數用作相應Variogram 方法的length.out 參數。默認為 50

collapse

一個可選字符串,指定要應用於各個半變異函數值的折疊類型。如果等於 "quantiles" ,則半變異函數值根據距離分布的分位數進行分割,每組的觀測值數量相同,距離間隔長度可能不同。否則,如果 "fixed" ,則根據相等長度的距離間隔劃分半變異函數值,每個間隔可能有不同數量的觀測值。否則,如果 "none" ,則不使用折疊並返回各個半變異函數值。默認為 "quantiles"

nint

一個可選整數,表示折疊半變異函數值時要使用的間隔數。默認為 20

robust

一個可選的邏輯值,指定在折疊各個值時是否應使用穩健的半變異函數估計量。如果TRUE使用魯棒估計器。默認為 FALSE

breaks

一個可選的數值向量,其中包含用於折疊半變異函數值的距離間隔的斷點。如果不丟失,則忽略 collapse 中指定的選項。

metric

指定要使用的距離度量的可選字符串。當前可用的選項是"euclidean",用於距離的平方和; "maximum"為最大差值; "manhattan" 為絕對差之和。使用參數的部分匹配,因此隻需提供前三個字符。默認為 "euclidean"

...

該泛型的某些方法需要額外的參數。此方法中沒有使用任何內容。

DataFrame ,其中 variogdist 列分別表示半變異函數值和相應的距離。如果半變異函數值折疊,則返回的 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-devel大神的英文原創作品 Calculate Semi-variogram for Residuals from a gls Object。非經特殊聲明,原始代碼版權歸原作者所有,本譯文未經允許或授權,請勿轉載或複製。