R: Двумерная непараметрическая регрессия

Question

Какие пакеты и функции в R могут выполнять двумерную неаддитивную локальную регрессию / сглаживание. Например, рассмотрим

b<-seq(-6*pi,6*pi,length=100)
xy<-expand.grid(b,b) 
x=xy[[1]]
y=xy[[2]]
z= sin(x)+cos(y) + 2*sin(x)*cos(y)
contour(b,b,matrix(z,100,100))

Какие функции могут это оценить?

Answer 1

вы можете сделать это с лессом:

fit <- loess( z ~ x+ y, span=0.01 )
dev.new()

contour( b, b, matrix( predict(fit), 100, 100 ) )

Answer 2

mgcv имеет множество вариантов 2-D сплайна.