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