R surf.gls 通过广义最小二乘法拟合趋势面

说明

通过广义最小二乘法拟合趋势面。

用法

surf.gls(np, covmod, x, y, z, nx = 1000, ...)

参数

np	多项式曲面的次数
covmod	评估协方差或相关函数的函数
x	x 坐标或包含 x 、 y 、 z 列的 DataFrame
y	y 坐标
z	z 坐标。将取代x$z
nx	协方差表的箱数。增加会增加准确性，并增加对象的大小。
...	covmod 的参数

值

包含组件的列表

beta	系数
x
y
z	以及其他仅供内部使用的内容。

例子

library(MASS)  # for eqscplot
data(topo, package="MASS")
topo.kr <- surf.gls(2, expcov, topo, d=0.7)
trsurf <- trmat(topo.kr, 0, 6.5, 0, 6.5, 50)
eqscplot(trsurf, type = "n")
contour(trsurf, add = TRUE)

prsurf <- prmat(topo.kr, 0, 6.5, 0, 6.5, 50)
contour(prsurf, levels=seq(700, 925, 25))
sesurf <- semat(topo.kr, 0, 6.5, 0, 6.5, 30)
eqscplot(sesurf, type = "n")
contour(sesurf, levels = c(22, 25), add = TRUE)

参考

Ripley, B. D. (1981) Spatial Statistics. Wiley.

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

也可以看看

trmat , surf.ls , prmat , semat , expcov , gaucov , sphercov