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