R fs.test FELSPLINE 测试函数

说明

实现基于 Tim Ramsay (2002) 提出的有限区域测试函数。

用法

fs.test(x,y,r0=.1,r=.5,l=3,b=1,exclude=TRUE)
fs.boundary(r0=.1,r=.5,l=3,n.theta=20)

参数

x , y	评估测试函数的点。
r0	测试域是一种弯曲的香肠。这是内弯的半径
r	香肠中心的曲线半径。
l	香肠臂的长度。
b	沿着香肠中心线的距离每增加一个单位，函数增加的速率。
exclude	外部点是否应该设置为NA？
n.theta	生成边界曲线时，在四分之一圆的分段线性表示中使用多少个点。

细节

源文章中没有给出函数的详细信息：但这非常接近。该函数是根据 Ramsay (2002) 进行修改的，它是凸出的，而不是平坦的：这可以更好地测试平滑器。

值

fs.test 返回函数求值，或 NA 返回边界外的点。 fs.boundary 返回要连接起来的 x,y 点列表，以便定义/绘制边界。

例子

require(mgcv)
## plot the function, and its boundary...
fsb <- fs.boundary()
m<-300;n<-150 
xm <- seq(-1,4,length=m);yn<-seq(-1,1,length=n)
xx <- rep(xm,n);yy<-rep(yn,rep(m,n))
tru <- matrix(fs.test(xx,yy),m,n) ## truth
image(xm,yn,tru,col=heat.colors(100),xlab="x",ylab="y")
lines(fsb$x,fsb$y,lwd=3)
contour(xm,yn,tru,levels=seq(-5,5,by=.25),add=TRUE)

作者

Simon N. Wood simon.wood@r-project.org

参考

Tim Ramsay (2002) "Spline smoothing over difficult regions" J.R.Statist. Soc. B 64(2):307-319