R trls.influence 趨勢麵的回歸診斷

說明

此函數提供了用於生成各種診斷的基本量，用於檢查 surf.ls 計算的趨勢麵的回歸擬合的質量。

用法

trls.influence(object)
## S3 method for class 'trls'
plot(x, border = "red", col = NA, pch = 4, cex = 0.6,
     add = FALSE, div = 8, ...)

參數

object , x	來自 surf.ls 的擬合趨勢麵模型
div	plot.trls 中影響圓半徑的比例因子
add	將影響圖添加到現有圖形 if TRUE
border , col , pch , cex , ...	附加圖形參數

值

trls.influence 返回包含組件的列表：

r	原始殘差由 residuals.trls 給出
hii	Hat 矩陣的對角元素
stresid	標準化殘差
Di	庫克的統計

例子

library(MASS)  # for eqscplot
data(topo, package = "MASS")
topo2 <- surf.ls(2, topo)
infl.topo2 <- trls.influence(topo2)
(cand <- as.data.frame(infl.topo2)[abs(infl.topo2$stresid) > 1.5, ])
cand.xy <- topo[as.integer(rownames(cand)), c("x", "y")]
trsurf <- trmat(topo2, 0, 6.5, 0, 6.5, 50)
eqscplot(trsurf, type = "n")
contour(trsurf, add = TRUE, col = "grey")
plot(topo2, add = TRUE, div = 3)
points(cand.xy, pch = 16, col = "orange")
text(cand.xy, labels = rownames(cand.xy), pos = 4, offset = 0.5)

參考

Unwin, D. J., Wrigley, N. (1987) Towards a general-theory of control point distribution effects in trend surface models. Computers and Geosciences, 13, 351-355.

Venables, W. N. and Ripley, B. D. (2002) Modern Applied Statistics with S. Fourth edition. Springer.

也可以看看

surf.ls , influence.measures , plot.lm