Deje que el modelo de la pregunta se escriba como
donde es un GP no observado con índice
y es un término de ruido normal con varianza . Por lo general, se supone que el GP es centrado, estacionario y no determinista. Tenga en cuenta que el término
puede considerarse como un GP (determinista) con kernel
donde
Yi=ϕ(xi)⊤β+h(xi)+εi(1)
h(x)x∈Rdεiσ2ϕ(x)⊤βϕ(x)⊤Bϕ(x)Bes una matriz de covarianza de `` valor infinito ''. De hecho, al tomar
con obtenemos las ecuaciones de kriging de la pregunta. Esto a menudo se denomina el
difuso anterior para
. Un posterior apropiado para
resulta solo cuando la matriz tiene rango completo. Entonces el modelo escribe tan bien como
donde es un GP . La misma interpretación de Bayes se puede usar con restricciones cuando ya no es un GP sino un
B:=ρIρ→∞ββΦYi=ζ(xi)+εi(2)
ζ(x)ζ(x)Función aleatoria intrínseca (IRF). La derivación se puede encontrar en el libro de G. Wahba. Las presentaciones legibles del concepto de IRF se encuentran, por ejemplo, en el libro de N. Cressie y el artículo de Mardia et al. Los IRF son similares a los procesos integrados bien conocidos en el contexto de tiempo discreto (como ARIMA): un IRF se transforma en un GP clásico mediante un tipo de operación de diferenciación.
Aquí hay dos ejemplos de IRF para . En primer lugar, considere un proceso de Wiener con su condición inicial reemplazada por una condición inicial difusa : es normal con una varianza infinita. Una vez que se conoce un valor , se puede predecir el IRF como es el GP de Wiener. En segundo lugar, considere un proceso Wiener integrado dado por la ecuación donde es un proceso de Wiener. Para obtener un GP ahora necesitamos dos parámetros escalares: dos valores
y parad=1ζ(x)ζ(0)=0ζ(0)ζ(x)
d2ζ(x)/dx2=dW(x)/dx
W(x)ζ(x)ζ(x′)x≠x′, o los valores
y en alguna elegida . Podemos considerar que los dos parámetros adicionales son conjuntamente gaussianos con una matriz de covarianza infinita . En ambos ejemplos, tan pronto como esté disponible un conjunto finito adecuado de observaciones, el IRF casi se enfrenta como un GP. Además, utilizamos un operador diferencial: y respectivamente. El espacio nulo es un espacio lineal de funciones
tal que . Contiene la función constante
ζ(x)dζ(x)/dxx2×2L:=d/dxL:=d2/dx2Fϕ(x)Lϕ=0ϕ1(x)=1en el primer caso y las funciones y
en el segundo caso. Tenga en cuenta que en el primer ejemplo
es GP para cualquier fijo en el primer ejemplo y de manera similar es un GP en el segundo caso.
ϕ1(x)=1ϕ2(x)=xζ(x)−ζ(x+δ)δζ(x−δ)−2ζ(x)+ζ(x+δ)
Para una dimensión general , considere un espacio lineal de funciones definidas en . Llamamos a un incremento
relativo a una colección finita de ubicaciones
y pesos reales tales que
Considere como el espacio nulo de nuestros ejemplos. Para el primer ejemplo podemos tomar, por ejemplo, con y
arbitraria ydFRdFsxi∈Rdsνi
∑i=1sνiϕ(xi)=0 for all ϕ∈F.
Fs=2x1x2[1,−1] . Para el segundo ejemplo, podemos tomar sy
espaciado igual . La definición de un IRF implica un espacio de funciones y una función que es
condicionalmente positivo wrt , lo que significa que
mantiene tan pronto como es un incremento wrt . De y
s=3xiν=[1,−2,1]Fg(x,x′)F∑i=1s∑j=1sνiνjg(xi,x′j)≥0
[νi,xi]si=1FFg(x,x′)
podemos hacer un núcleo de covarianza, por lo tanto, un GP como en Mardia et al. Podemos comenzar desde un operador diferencial lineal y usar el espacio nulo como ; la IRF tendrá conexión con la ecuación un ruido gaussiano.
LFLζ=
El cálculo de la predicción del IRF es casi el mismo que en la pregunta, con reemplazado por
, pero con ahora formando una base de . La restricción adicional
debe agregarse en el problema de optimización, lo que garantizará que
. Todavía podemos agregar más funciones básicas que no están en
si es necesario; esto tendrá el efecto de agregar un GP determinista, digamos
al IRF
k(x,x′)g(x,x′)ϕi(x)FΦ⊤α=0α⊤Kα≥0Fψ(x)⊤γζ(x) en (2).
La spline de placa delgada depende de un número entero tal que , el espacio contiene polinomios de bajo grado, con una dimensión depende de y . Se puede demostrar que si
es la siguiente función para luego
define un wrt condicionalmente positivo . La construcción se refiere a un operador diferencialmm>2dFp(m)mdE(r)r≥0
E(r):={(−1)m+1+d/2r2m−dlogrr2m−dd even,d odd,
g(x,x′):=E(∥x−x′∥)FL. Resulta que para y de la delgada spline placa hay nada que el spline cúbico natural, usual, que se refiere al ejemplo Wiener integrado anteriormente, con . Entonces (2) no es más que el modelo de spline de suavizado habitual. Cuando y el espacio nulo tiene dimensión
y se genera por las funciones , y .
d=1m=2g(x,x′)=|x−x′|3d=2m=2p(m)=31x1x2
Estadísticas de Cressie N para datos espaciales . Wiley 1993.
Mardia KV, Kent JT, Goodall CR y Little JA. Kriging y splines con información derivada. Biometrika (1996), 83,1, pp. 207-221.
Modelos Wahba G Spline para datos de observación . SIAM 1990.
Wang, Y Alisar Splines, Métodos y Aplicaciones . Chapman y Hall, 2011.