Estimación Espacio-Temporal de Procesos Hilbert-Valuados
Estimación Espacio-Temporal de Procesos Hilbert-Valuados
T ECNOLOGÍAS E I NGENIERÍAS
TESIS DOCTORAL
Directoras de tesis:
María del Pilar Frías Bustamante
María Dolores Ruiz Medina
Tutora:
María Dolores Ruiz Medina
Coordinador:
Francisco De Asís Torres Ruiz
Granada, 2021
Editor: Universidad de Granada. Tesis Doctorales
Autor: Antoni Torres Signes
ISBN: 978-84-1306-891-6
URI: [Link]
Índice general
Portada I
Índice de figuras IX
Abstract XX
Resumen XXV
Introducción 2
Objetivos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
Objetivos generales . . . . . . . . . . . . . . . . . . . . . . . . . . 17
Estructura de la tesis . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
1. P RELIMINARES 27
1.1. Preliminares para el Capítulo 2 . . . . . . . . . . . . . . . . . . . 27
1.1.1. Estimadores de mínimo contraste y datos tapered . . . . . 27
1.1.2. Definición general de estimadores de mínimo contraste . . 30
1.1.3. Datos y periodograma tapered . . . . . . . . . . . . . . . . 33
Esperanza y covarianza del periodograma . . . . . . . . . 35
1.1.4. Ejemplo de estimación a partir de datos tapered . . . . . . 37
1.1.5. Sobre los estimadores Whittle con datos tapered . . . . . . 39
Caso de campos aleatorios débilmente dependientes . . . 39
Caso de campos aleatorios fuertemente dependientes . . . 41
1.2. Preliminares para los capítulos 3, 6 y 7 . . . . . . . . . . . . . . . 43
1.2.1. Procesos de Poisson doblemente estocásticos . . . . . . . . 43
V
Índice general
M ETODOLOGÍA Y R ESULTADOS 59
VI
Índice general
VII
Índice general
B IBLIOGRAFÍA 272
VIII
Índice de figuras
IX
Índice de figuras
X
Índice de figuras
5.1. Escala 10, N = 900. Datos curva sobre algunos nodos de una
rejilla regular espacial 30 × 30 . . . . . . . . . . . . . . . . . . . . 181
5.2. Escala 9, N = 900. Datos curva sobre algunos nodos de una
rejilla regular espacial 30 × 30 . . . . . . . . . . . . . . . . . . . . 181
5.3. Escala 8, N = 900. Datos curva sobre algunos nodos de una
rejilla regular espacial 30 × 30 . . . . . . . . . . . . . . . . . . . . 182
5.4. Escala 7, N = 900. Datos curva sobre algunos nodos de una
rejilla regular espacial 30 × 30 . . . . . . . . . . . . . . . . . . . . 182
5.5. nErrores cuadráticos medioso empíricos (EMSE), asociados con
λ bN,p,2 , p = 1, . . . , kN , N = 100, 900, 2500, 4900, 8100,
bN,p,1 , λ
12100, 16900, 22500 . . . . . . . . . . . . . . . . . . . . . . . . 184
XI
Índice de figuras
XII
Índice de figuras
XIII
Índice de figuras
XIV
Índice de figuras
XV
Índice de tablas
XVI
Índice de tablas
XVII
Índice de tablas
XVIII
Abstract
XX
Abstract
XXI
Abstract
Epidemiology and the study in general of the evolution, both spatial and
temporal, of several diseases has been the fundamental framework considered
for the contributions indicated. Specifically, real data have been used for the
estimation and functional prediction in time and space of prostate, breast and
brain cancer, as well as respiratory diseases, in Spanish provinces, from annual
or monthly observations, over periods ranging around thirty years. Furthermo-
re, by implementing the techniques presented throughout the thesis, an appli-
cation to real data has been carried out to analyze the incidence of a disease
in a foreign territory. In particular, the evolution of dengue fever in American
countries in recent years has been modeled. On the other hand, given the so-
cial emergency situation caused by the COVID-19 pandemic in the last stage of
development of the thesis, it has been considered pertinent to include a statis-
tical study on the estimation of the spatio-temporal evolution of the mortality
risk, as well as of the daily mortality cases caused by this disease in the Autono-
mous Communities, which allows reflecting, among other aspects, the effect of
the first state of alarm on the behavior of this evolution. In this way, the daily
mortality due to COVID-19 in the Spanish Communities during the first wave,
specifically from 8 March 2020 to 13 May 2020, has been modeled.
The latter practical applications have been developed, based on the infinite-
dimensional statistical techniques proposed in the development of the thesis,
under a classical and Bayesian approach, with modifications in the estimation
methodology. Subsequently, in both cases, an empirical comparison has been
made with other approaches. In the case of the risk of annual incidence of
dengue fever in American countries, a comparison has been made with tradi-
tional spatio-temporal models, including a Leroux model, an Intrinsic Condi-
XXII
Abstract
tional Autoregressive model and a Besag, York and Mollie model. In the case
of daily mortality risk by COVID–19 in the Spanish Autonomous Communities,
the proposed approaches have been compared with another methodology ba-
sed on the estimation by confidence intervals and probability densities using
bootstrap techniques, as well as with a battery of models in the context of Ma-
chine Learning, including Generalized Regression Neural Networks, Multilayer
Perceptron, Support Vector Regression, Bayesian Neural Networks, Neural Net-
works from Radial Function Bases, and Gaussian Processes. In addition, model
selection in the context of parametric non-linear regression is addressed.
XXIII
Resumen
XXV
Resumen
rrelados en el tiempo. Por otro lado, se desarrollan los procesos de Cox dirigidos
por log-intensidades espaciales infinito-dimensionales lineales, en este caso, es-
tos patrones espaciales se analizan desde una perspectiva de datos funcionales
correlados en el espacio.
En cuanto a los enfoques metodológicos adoptados para la estimación, en
el ámbito de procesos de Cox dirigidos por una log-intensidad O-U Hilbert-
valuada, aproximada mediante un proceso ARH(1), se ha utilizado el método
de los momentos empíricos. En el caso de procesos de Cox espaciales dirigidos
mediante una log-intensidad aleatoria lineal espacial infinito-dimensional, para
la estimación de los parámetros que modelizan la estructura paramétrica del
operador de densidad espectral, bajo la condición de estacionariedad espacial,
se han aplicado técnicas espectrales funcionales basadas en el operador perio-
dograma que extienden el funcional de Whittle. Como análisis preliminar, se
contribuye, en el caso de procesos espaciales real-valuados estacionarios en el
espacio, con la obtención de condiciones suficientes que garantizan la consis-
tencia y normalidad asintótica de estimadores de mínimo contraste, basados en
el periodograma tapered.
De forma concreta, en esta tesis, a partir de una perspectiva de proce-
sos de Cox infinito-dimensionales, o bien, procesos de Cox dirigidos por log-
intensidades lineales infinito-dimensionales, no necesariamente gaussianas, den-
tro del ámbito del análisis estadístico funcional de patrones puntuales en el
tiempo y/o espacio, se han establecido las siguientes contribuciones:
XXVI
Resumen
XXVII
Resumen
permite reflejar, entre otros aspectos, el efecto del primer estado de alarma so-
bre el comportamiento de dicha evolución. De esta manera, se ha modelizado
la mortalidad diaria por COVID–19, en las comunidades españolas, durante la
primera ola, en concreto, desde el 8 de marzo de 2020 hasta el 13 de mayo de
2020.
Estas últimas aplicaciones prácticas se han desarrollado, a partir de las téc-
nicas estadísticas infinito-dimensionales propuestas en el desarrollo de la tesis,
bajo un enfoque clásico y bayesiano, con modificaciones en la metodología de
estimación. Posteriormente, en ambos casos se ha realizado una comparativa
empírica con otros enfoques. En el caso del riesgo de incidencia anual por fie-
bre del dengue en países americanos, se ha realizado una comparativa con mo-
delos espacio-temporales tradicionales, incluyendo un modelo de Leroux, un
modelo Autorregresivo Condicionado Intrínseco y otro modelo de Besag, York
y Mollie. En el caso del riesgo de mortalidad diario por COVID–19 en las comu-
nidades autónomas españolas, los enfoques propuestos se han comparado con
otra metodología basada en la estimación por intervalos de confianza y densi-
dades de probabilidad mediante técnicas bootstrap, así como con una batería de
modelos en el contexto de Machine Learning, incluyendo Redes Neuronales de
Regresión Generalizada, Perceptrón Multicapa, Regresión de Soporte Vectorial,
Redes Neuronales Bayesianas, Redes Neuronales a partir de Bases de Funciones
Radiales, y Procesos Gaussianos. Además se aborda la selección de modelos en
el contexto de la regresión no lineal paramétrica.
XXVIII
Introducción
2
Introducción
narios, fractales y con dependencia de largo rango; véase, por ejemplo, Chung
(1996). En particular, WeiLin et al. (2012) aplicaron la estimación paramétrica
de mínimo contraste para aproximar el parámetro de desviación del proceso de
O-U cuando la correspondiente ecuación diferencial estocástica viene dada por
el movimiento browniano fraccionario.
La consistencia y normalidad asintótica del estimador de máxima verosimi-
litud de Whittle para los procesos autorregresivos de medias móviles integrados
fraccionalmente estacionales y estacionarios (SARFIMA) se ha demostrado en
Guo et al. (2009). La maximización de la verosimilitud de Whittle también se
ha considerado en diversos autores (Boissy et al., 2005; Chan y Tsai, 2012).
En Robinson y Sanz (2006), se consideran modificaciones a la estimación de
Whittle de modelos multilaterales sobre rejilla. En Leonenko y Taufer (2013),
se ha estudiado la estimación de la densidad espectral para campos aleatorios
gaussianos agregados, con posible dependencia fuerte, a partir de una expan-
sión en términos de polinomios ortogonales de Gegenbauer. Estos campos se
construyen a partir de la acumulación de campos aleatorios de memoria cor-
ta independientes y con idéntica distribución, a través de una densidad mixta
desconocida (a estimar).
Una alternativa a la familia de funcionales lineales de Whittle se propuso en
Ibragimov (1967), véase también Anh et al. (1999, 2004a,b, 2007a,c,d); Espejo
et al. (2015), entre otros. Un desarrollo de estos métodos de estimación en el
espectro se presenta en el Capítulo 2, a partir del periodograma tapared. Los
resultados de dicho capítulo sirven de base para la extensión al caso funcio-
nal que se estudia en el Capítulo 4. En el contexto de los modelos lineales en
espacios funcionales, motivado por el análisis estadístico de datos de grandes
dimensiones, se puede consultar Bosq (2000); Bosq y Blanke (2007); Ferraty y
Vieu (2006); Ramsay y Silverman (2005), entre otros. Un tratamiento completo
de la estimación basada en el método de los momentos y la proyección en bases
ortogonales de funciones se puede ver en Antoniadis y Spatinas (2003); Bosq
(2000); Bosq y Blanke (2007); Damon y Guillas (2005); Soltani y Hashemi
(2011). En el contexto de series autorregresivas espaciales hilbertianas desta-
can los trabajos de Ruiz-Medina y Espejo (2012, 2013).
En las dos últimas décadas, ha habido un desarrollo creciente de técnicas
que permiten el análisis estadístico de datos funcionales (curvas, superficies,
etc.) que están correlados en el tiempo o espacio, es decir, con una perspectiva
3
Introducción
4
Introducción
5
Introducción
6
Introducción
7
Introducción
8
Introducción
literatura financiera, después de la obra pionera de Cox et al. (1985). Una de las
propiedades del proceso de Feller es que se da en R+ , cumpliendo la condición
de no negatividad.
La simplicidad de un proceso de Cox radica en que las propiedades de los
momentos de un proceso de Cox son equivalentes a las propiedades de los
momentos de su correspondiente campo aleatorio de intensidad. De hecho, la
función de intensidad de un proceso de Cox es igual a la esperanza de la inten-
sidad en un campo aleatorio, y la densidad de covarianza de un proceso de Cox
es igual a la función de covarianza de su intensidad aleatoria.
Se pueden destacar los siguientes trabajos, que inicialmente promovieron
el desarrollo de la literatura en el ámbito de procesos puntuales y procesos de
Cox. En particular, en el contexto de la bioestadistica y fiabilidad, mencionamos
los trabajos iniciales de Gail et al. (1980), donde se aplica la teoría de procesos
puntuales para estudiar el tamaño de los tumores en ratas durante un período
de tiempo; Lando (1998), donde la descripción de ocurrencias de obtención de
créditos se logra a partir de los procesos puntuales; y Dalal y McIntosh (1994),
donde la prueba y la validación del software se realiza en la teoría de la fia-
bilidad. En Rathbun y Cressie (1994), se considera el contexto de procesos de
Cox log-gaussianos, dirigidos por una intensidad aleatoria con valores constan-
tes sobre los cuadrantes de una rejilla regular; véase también Besag (1974). El
muestreo de Gibbs se puede aplicar para explorar estos modelos discretizados,
véase Rathbun y Cressie (1994).
La literatura ofrece varias subclases de procesos de Cox de particular interés.
En la presente tesis se presta atención sobre los procesos de Cox log-gaussianos
(LGCP), véase, por ejemplo, Møller et al. (1998), definidos como procesos de
Cox con una intensidad aleatoria que es la exponencial de un campo aleatorio
gaussiano, o Rostami et al. (2012), donde se considera un proceso puntual de
Cox log-gaussiano para modelizar la distribución espacio-temporal de la mor-
talidad por abuso de sustancias en Irán. Este modelo de log-intensidad normal
proporciona un marco flexible en el análisis de patrones puntuales espaciales y
espacio-temporales, véase Diggle et al. (2013); González et al. (2016).
La caracterización completa de esta clase de procesos por la intensidad y la
densidad producto de segundo orden hace posible su aplicación en diferentes
campos, ver, por ejemplo, Rathbun y Cressie (1994) en bosques de pinos; Serra
et al. (2014) en incendios forestales. También se puede señalar algún ejem-
9
Introducción
plo sobre áreas de aplicación adicionales como ecología, con trabajos como
Serra et al. (2014); Waagepetersen et al. (2016); sismología (Møller y Tofta-
ker, 2014); o reconocimiento de texturas de escenas visuales (Simpson et al.,
2016). En la presente tesis se va más allá del caso real-valuado, introduciendo
procesos de Cox dirigidos por un proceso (o campo) aleatorio de log-intensidad
Hilbert-valuada.
En un contexto matemático y probabilístico, modelos lineales y no linea-
les de ecuaciones de estados en el espacio y tiempo para procesos funcionales
se han considerado en los trabajos de Kelbert et al. (2005); Leonenko y Ruiz-
Medina (2006); Leonenko et al. (2011), entre otros. La estimación paramétrica
de modelos de este tipo mediante mínimo contraste se puede ver, por ejem-
plo, en Anh et al. (2004a), en el caso espacial. Este tipo de modelos permi-
ten representar procesos espacio-temporales con difusión anómala espacial. En
particular, permiten representar la distribución espacial anómala de riesgos de
mortalidad (respectivamente incidencia) de enfermedades cuando se observan
patrones espaciales cuyas log-intensidades aleatorias son localmente singulares,
o bien, presentan discontinuidades, así como un comportamiento local errático
que corresponde a un modelo diferencial fraccionario, en los correspondientes
mapas de evolución de la enfermedad.
En el ámbito de la salud, la epidemiología investiga la evolución de los pa-
trones espaciales de distribución de enfermedades en la población humana, en-
tre otros aspectos. Es interesante ver las posibilidades de aplicar la Estadística
a la modelización de tales sucesos. Se consideran los registros de ocurrencias
para determinar los riesgos de mortalidad e incidencia asociados a una enfer-
medad cuando el objetivo es un análisis de la evolución. Así, la determinación
de patrones espaciales es de particular interés para identificar los territorios de
riesgo y, en base a los mapas de riesgo, se puede lograr un mejor desempeño de
los planes nacionales de emergencia de salud pública (Sánchez-Gómez et al.,
2017). La inversión en la vigilancia y respuesta a epidemias y pandemias, junto
con el mapeo de riesgos y su comunicación, son esenciales para la detección
temprana de los brotes. Lawson (2008) sugiere que los usuarios y lectores se
sienten más cómodos si los datos sanitarios agregados se analizan en mapas.
Tradicionalmente, los riesgos de mortalidad (o incidencia de una enferme-
dad) suelen estar representados en mapas que resumen la distribución espacial
o espacio-temporal de los patrones de mortalidad en diferentes regiones. La me-
10
Introducción
11
Introducción
12
Introducción
las referencias en ellos). Las predicciones a corto plazo pueden obtenerse adop-
tando el marco de los modelos compartimentados SIR (susceptible-infeccioso-
recuperado), basados en ecuaciones diferenciales ordinarias, véanse, por ejem-
plo, Angulo et al. (2013); Elhia et al. (2014); Ji et al. (2012); Kermack y Mc-
Kendrick (1927); Kuznetsov y Piccardi (1994); Milner y Zhao (2008); Pathak
et al. (2010); Tornatore et al. (2005); Yu et al. (2009); Zhang et al. (2008).
Existe una extensa literatura, que incluye diferentes versiones de modelos
compartimentados, como el SIR-susceptible, SIRS (Dushoff et al., 2004), y for-
mulaciones basadas en ecuaciones diferenciales de retardo, ver Beretta et al.
(2001); McCluskey (2010); Sekiguchi y Ishiwata (2010). Las extensiones espa-
ciales, basadas en modelos de reacción-difusión, que reflejan la propagación de
la enfermedad infecciosa en una región espacial pueden encontrarse, por ejem-
plo, en Guin y Mandal (2014); Webb (1981). También se analizan los modelos
SEIRD (susceptible, expuesto, infectado, recuperado, fallecido), que incorporan
la propagación espacial de la enfermedad con términos de difusión no homogé-
neos (Roques y Bonnefon, 2016; Roques et al., 2011).
La versión estocástica de los modelos de tipo SIR pretende cubrir varias li-
mitaciones detectadas en cuanto a la incertidumbre, en las observaciones, y
el proceso epidémico dinámico oculto. La modelización basada en cadenas de
Markov SIR (Anderson et al., 2017; Xu et al., 2007), y algunas formulaciones
estocásticas que implican redes complejas (Volz, 2008; Zhou et al., 2006) o la
gripe resistente a los medicamentos (Chao et al., 1996) constituyen algunas al-
ternativas. En Aalen et al. (2008); Abboud et al. (2019); Anderson et al. (2017);
Fleming y Harrington (1991), se adopta un marco de modelo estadístico baye-
siano SIRS teniendo en cuenta el error de observación en los recuentos, y la
incertidumbre en el espacio de parámetros. Más allá de la modelización SIR,
mencionamos aquí los enfoques basados en el análisis multivariante y de super-
vivencia para modelizar, por ejemplo, la infección, la incubación y los períodos
aleatorios de recuperación, que afectan a la contención del COVID-19, véase,
por ejemplo, Bolker y Grenfell (1996); Keeling et al. (1997); Pak et al. (2020);
Wasiur et al. (2019).
En una primera etapa, la mayoría de los modelos referidos han sido adapta-
dos y aplicados para aproximar la evolución espacio-temporal de la incidencia
y la mortalidad del COVID–19. Ese es el caso, por ejemplo, de los tres modelos
presentados en Roosa et al. (2020), que fueron validados con brotes de otras
13
Introducción
14
Introducción
Objetivos
La motivación fundamental para el desarrollo de esta tesis reside en el aná-
lisis estadístico de la evolución de patrones aleatorios espaciales de enferme-
dades, esencialmente en el marco epidemiológico. Se intenta proporcionar he-
rramientas estadísticas apropiadas para establecer la evolución en el tiempo así
como en el espacio, de mapas de riesgo o, en su defecto, mapas de mortalidad
o incidencia de una enfermedad.
La temática teórica se desarrolla en el contexto de modelos que utilizan los
procesos de Cox log-gaussianos en espacios infinito-dimensionales, en particu-
lar, se consideran espacios de Hilbert separables. Se analizan las propiedades
asintóticas de las estimaciones paramétricas para dichos modelos y también, de
una forma más general, para campos aleatorios lineales estacionarios a partir de
funcionales espectrales, basados en el periodograma y periodograma tapered,
así como su versión infinito-dimensional, es decir, el operador periodograma.
La teoría espectral de operadores compactos y autoadjuntos juega un pa-
pel fundamental en todos estos desarrollos, cuando se considera una diagona-
lización de la estructura de dependencia del proceso latente de interés, que
usualmente se trata de un proceso temporal funcional o espacial funcional.
En el contexto log-gaussiano, centramos la atención en las medidas gaussianas
Hilbert-valuadas, definidas mediante el producto infinito de medidas gaussia-
nas univariantes.
En el ámbito computacional y aplicado, se desarrollan e implementan nue-
vos algoritmos para el análisis estadístico y la construcción de mapas de riego
de enfermedades y mapas de mortalidad e incidencia causada por enfermedad,
15
Introducción
16
Introducción
Objetivos generales
(ii) Obtener los resultados teóricos que definen las condiciones que garanti-
zan un buen comportamiento asintótico de los estimadores y predictores
formulados (es decir, consistencia débil y fuerte y normalidad asintótica).
17
Introducción
(vi) Introducir una nueva clase de procesos espaciales puntuales definidos me-
diante la familia de procesos de Cox estudiados, dirigidos por una intensi-
dad lineal infinito-dimensional espacial. Se considerará, como ilustración
de los resultados de estimación derivados, la sub-clase de procesos de Cox
dirigidos mediante una intensidad espacial SARH(1).
(ix) Llevar a cabo aplicaciones prácticas reales con datos reales espacio-temporales
donde ilustrar el estudio, metodología y resultados indicados en los obje-
tivos anteriores.
Estructura de la tesis
Para el desarrollo de la tesis, se comienza con una introducción en la que
se realiza una motivación de la misma que incluye la literatura relacionada
destacada, los principales objetivos señalados y su estructura. Posteriormente,
en las secciones 1.1, 1.2 y 1.3, se establecen los preliminares necesarios para el
desarrollo en la obtención de los resultados de los capítulos posteriores.
De forma más concreta, en la Sección 1.1 se establecen los preliminares
sobre los estimadores de mínimo contraste a partir de datos tapered1 . Los apar-
tados de esta sección evidencian los problemas que surgen, en la estimación,
cuando se trabaja con campos aleatorios respecto a las metodologías de estima-
ción, basadas en el espectro unidimensional, en el caso de procesos estocásticos
temporales.
En la Sección 1.2 se presentan los preliminares y el fundamento teórico
de los procesos de Cox log-gaussianos. Primeramente se definen los procesos
de Poisson doblemente estocásticos, o procesos de Cox, y posteriormente se
1
Tapered se podría traducir como una disminución gradual del diámetro, anchura y grosor
de un objeto. En nuestro caso se trata de la disminución gradual en el borde de los campos
aleatorios a partir de una función taper, véase el Apartado 1.1.3.
18
Introducción
19
Introducción
20
Introducción
21
Introducción
22
Introducción
23
Introducción
24
Introducción
25
Capítulo 1
P RELIMINARES
C(h)
h → ρ(h) = ,
C(0)
27
Capítulo 1. P RELIMINARES
1 X 0
f (u) = d
C(h)e−ih u , u∈T
(2π) d
h∈Z
28
Capítulo 1. P RELIMINARES
1 X
Ĉ(h) = Z(h)Z(s + h),
N
s∈D(h)
1 X 0
I(λ) = 2
Ĉ(h)eih λ ,
(2π) h∈H
1 X 0
J(λ) = √ Z(s)eis λ ,
2π n1 n2 s∈D
29
Capítulo 1. P RELIMINARES
30
Capítulo 1. P RELIMINARES
en P0 –probabilidad.
Esta Condición (1.1) de subergodicidad transfiere el hecho de que el valor
UT (θ) − UT (θ0 ) que se obtiene al estimar el contraste de θ con θ0 en base
a Y (t), asintóticamente separa los parámetros. Esta Condición (1.1) se
puede fortalecer mediante la condición ergódica1 , obteniendo:
31
Capítulo 1. P RELIMINARES
Xi = m(Zi , θ0 ) + εi , i = 1, . . . , n,
A continuación se define K(θ, θ0 ) = lı́m inf n→∞ ni=1 (m(xi , θ) − m(xi , θ0 ))2 /n.
P
Para ver más características de este ejemplo puede consultarse Gaetan y Guyon
(2010).
Ejemplo 1.2. Otro ejemplo que puede resultar interesante es el contraste gaus-
siano en procesos de segundo orden.
Consideremos que X = (X(t), t ∈ Z) es una serie temporal centrada esta-
cionaria de segundo orden con densidad espectral fθ . El periodograma asocia-
do a las autocovarianzas empíricas, r̂n (k), a partir de las observaciones Xn =
(X(1), X(2), . . . , X(n)), es la estimación de la densidad espectral
n−1
1 X
In (λ) = r̂n (k)eiλk ,
2π k=−n+1
32
Capítulo 1. P RELIMINARES
Pn−|k|
donde r̂n (k) = r̂n (−k) = n−1 i=1 Xi Xi+k .
La estimación de f (λ) a partir del periodograma In (λ) resulta bastante pobre.
Sin embargo, se puede utilizar el contraste de Whittle, definido a partir de la
siguiente regularización de In (λ):
Z 2π
1 In (λ)
Un (θ) = log fθ (λ) + dλ. (1.3)
2π 0 fθ (λ)
1 X
ĈT (k) = Y (t)Y (t + k), k ∈ Zd ,
Td t, t+k ∈DT
33
Capítulo 1. P RELIMINARES
Se denota
T
X
H̃k,T (λ) = [h̃T (t)]k e−iλt ,
t=−T
X d
Y
k −i(λ,t)
Hk,T (λ) = [hT (t)] e = H̃k,T (λ(i) ), (1.4)
t∈LT i=1
34
Capítulo 1. P RELIMINARES
1
ITh (λ) = dh (λ)dhT (−λ). (1.5)
(2π)d H2,T (0) T
35
Capítulo 1. P RELIMINARES
véase, por ejemplo, Dahlhaus y Polonik (2002). Así pues, el periodograma tape-
red puede funcionar como estimador de la densidad espectral cuando T es sufi-
cientemente grande, ya que asintóticamente es un estimador insesgado, véase,
por ejemplo, Brillinger (1970), entre otros.
36
Capítulo 1. P RELIMINARES
d
Y ik − 1/2
hT (t) Y (t) = aT (t)Y (t) , donde aT (t) = h , t ∈ LT .
k=1
T
37
Capítulo 1. P RELIMINARES
1 X h iλ0 k
ITh (λ) = Ĉ T (k)e .
(2π)d k∈L
T
ITh (λ)
Z
1
UT (θ) = {log fθ (λ) + }dλ, T = [0, 2π)d ,
(2π)d T fθ (λ)
38
Capítulo 1. P RELIMINARES
39
Capítulo 1. P RELIMINARES
forma:
d
t − 1/2 Y ti − 1/2
hT (t) = h = h̃ ,
T i=1
T
y
g 2uρ
si 0 ≤ u ≤ 12 ρ
h̃(u) = 1 si 21 ρ ≤ u ≤ 12
h̃(1 − u) si 1 ≤ u ≤ 1
2
donde la función g : [0, 1] → [0, 1] es tal que g(0) = 0 y g(1) = 1. Esta función
es creciente y pertenece a C 2 , con 0 ≤ ρ ≤ 1. De esta forma, se obtiene el
denominado g–taper, el cual disminuye un 100(1 − ρ) % los valores del borde.
A continuación se formulan los resultados sobre el control del sesgo para los
estimadores tapered indicados por Guyon; véase, Guyon (1995).
R
Consideremos J(φ) = T f (λ)φ(λ) dλ, y el correspondiente estimador basa-
R
do en los datos tapered JT (φ) = T ITh (λ)φ(λ) dλ.
Proposición 1.1. Sea φ ∈ C(T), supongamos que h y la densidad espectral, f ,
pertenecen a C 2 . Entonces, cuando T → ∞,
donde el factor taper e(h) viene dado por (2.9), e(h) ≥ 1, y e(h) = 1, si no se
produce tapering. Se puede elegir un ρT tal que el factor taper tienda a 1.
40
Capítulo 1. P RELIMINARES
Nota 1.1. En vista de la anterior Proposición 1.1, con el fin de controlar el sesgo
para el estimador de Whittle, se pueden usar las siguientes condiciones:
la función taper h̃(t) está en C 2 [−1, 1]; la densidad espectral f (λ, θ) está en C 2 (T)
y la función ∂θ∂ i f (λ, θ) es continua con respecto a (λ, θ).
O(|λi |−αi (θ)−δ ), |λi | → 0 para cada i = 1, ..., d, y que f0 (λ, θ) es una función
de (λ, θ) positiva con segunda derivada continua.
A2 En este caso, suponemos que existe una función α : Θ → (0, d), tal que la
densidad espectral f (λ, θ) satisface, para todo δ > 0 : f (λ, θ) =
kλk−α(θ)+δ (1 + o(1)), donde el o(1) es uniforme con respecto a λ. Para
más detalles, véase Ludeña y Lavielle (1999).
wi (λi ) = |λi |ν , o bajo A2, w(λ) = 1T (λ), esto es, la función indicadora en el
d–toro.
La función taper hT,ε (t) = di=1 h̃ε tTi , donde h̃ε : [0, 1] → [0, 1] es de la
Q
siguiente forma:
u/ε si u ≤ ε
hε (u) = 1 si ε ≤ u ≤ 1 − ε
h (1 − u) si u > 1 − ε
ε
41
Capítulo 1. P RELIMINARES
con ε = T −γ .
En el curso de la deducción de la normalidad asintótica para el estimador
de Whittle, se obtuvo el siguiente resultado sobre el sesgo.
Proposición 1.2. Supongamos que se cumplen A1 o A2, y que w(λ) = di=1 wi (λi )
Q
bajo A1, o w(λ) = 1T (λ) bajo A2. Además, se considera que γ < (1 + v − d/2)/v
bajo A1, o γ < 2 − d/2 bajo A2. Entonces, cuando T → ∞
donde κ = ν bajo A1, o κ = 1 bajo A2, y por lo tanto, el sesgo es de orden o(T −d/2 ),
si d ≤ 3.
42
1.2. Preliminares para los capítulos 3, 6 y 7
Dentro de los procesos de Poisson hay una gran variedad de tipos. El caso
homogéneo es el más simple de ellos, sin embargo, resulta insuficiente para
modelizar la mayoría de casos reales. Para ampliar más las posibilidades del
proceso de Poisson se utilizan los procesos de Poisson doblemente estocásticos,
donde la intensidad no solo es no homogénea sino que a su vez es un proce-
so estocástico. A estos procesos también se les conoce como procesos de Cox,
quien los introdujo en su artículo, Cox (1955). Un proceso de Poisson doble-
mente estocástico puede ser visto como un procedimiento de aleatorización en
dos pasos. Un proceso λt , t > 0, se utiliza para generar otro proceso, Nt , ac-
tuando como su intensidad. En particular, se supone que Nt es un proceso de
Poisson condicionado a λt . Se pueden dar muchas posibles definiciones de un
proceso de Poisson doblemente estocástico. De todas las posibles, se ha seguido
la utilizada en la Sección 2.1, pp. 2–3, en Møller (2005), por ser muy general,
bien fundamentada y adaptada al marco espacial.
Primeramente, se define un proceso puntual espacial en Rd como un con-
junto aleatorio localmente finito X ⊂ Rd , lo que significa que para cualquier
conjunto de Borel acotado B ⊂ R2 , el número de puntos en XB = X ∩ B es una
variable aleatoria finita denotada por N (B).
Se dice que X es estacionario, o respectivamente isotrópico, si su distribu-
ción es invariante bajo traslaciones en Rd , o respectivamente rotaciones sobre
el origen en Rd .
Para cualquier entero n ≥ 1, la medida del momento factorial de orden n,
Capítulo 1. P RELIMINARES
6=
X
(n)
µ (A) = E 1[(u1 , . . . , un ) ∈ A]
u1 ,...,un ∈X
para conjuntos de Borel A ∈ Rdn , donde 6= sobre la suma significa que la suma
recorre todos los puntos distintos, dos a dos, u1 , . . . , un en X, y 1[·] es la función
indicadora. Se asume que µ(n) es localmente finita y tiene función de densidad
ρ(n) con respecto a la medida de Lebesgue en Rdn , donde llamamos a ρ(n) la
densidad producto de orden n. La medida de momento es µ(A) = µ(1) (A) =
E[N (A)], y ρ = ρ(1) es la función de intensidad. Intuitivamente, si u1 , . . . , un ∈ Rd
son distintos dos a dos, ρ(n) (u1 , . . . , un ) du1 · · · dun es la probabilidad de observar
la ocurrencia conjunta de n puntos de X en cada n regiones infinitesimales con
áreas du1 , . . . , dun , conteniendo u1 , . . . , un .
Una extensión natural del proceso de Poisson es considerar un proceso esto-
cástico no negativo Λ = (Λ(u))u∈Rd tal que X condicionado a Λ es un proceso de
Poisson con intensidad Λ, Cox (1955). Entonces, decimos que X es un proceso
de Cox (o doblemente estocástico) con intensidad Λ.
Como extensión inmediata de las propiedades de un proceso de Poisson X |
Λ, se tiene que X es estacionario o isotrópico si Λ es estacionaria o isotrópica
respectivamente. Además, la densidad producto
de manera que
E[Λ(u)Λ(v)]
ρ(u) = E[Λ(u)], g(u, v) = .
ρ(u)ρ(v)
44
Capítulo 1. P RELIMINARES
• Las propiedades teóricas son fácilmente extraídas. Por ejemplo, los mo-
mentos de orden superior se expresan simplemente por la intensidad y la
función de correlación par del proceso de Cox log-gaussiano.
45
Capítulo 1. P RELIMINARES
P 2
se sigue que j Nt,j < ∞ casi seguramente (c. s.). Entonces, Mt = {Nt,j , j ≥ 0}
define una variable aleatoria con valores en `2 , y satisface EkMt k2`2 < ∞. Así,
{Mt , t ≥ 0} es un proceso en tiempo continuo `2 -valuado.
Predicción clásica
46
Capítulo 1. P RELIMINARES
j0 (ω,λ,T )
Y (λj T )NT,j (ω)
L(MT (ω), λ) = exp (−λj T ) .
j=1
NT,j (ω)!
Claramente, ∞
P b
j=1 λT,j < ∞ (c. s.) y (λ)T es insesgado (considerando (λ)
b
como un parámetro con valores en `2 ). Se sigue que
T +h
f (MT ) = MT (1.7)
T
47
Capítulo 1. P RELIMINARES
Una función µ : C(H) → [0, ∞] se dice medible cilíndrica en C(H) si, para
cada subconjunto finito Γ ⊆ H ∗ , la restricción de µ sobre la σ–álgebra C(H, Γ)
es una medida. Una medida cilíndrica se dice finita si µ(H) < ∞.
Un proceso cilíndrico X en H es una familia {X(t), t > 0} de variables
aleatorias cilíndricas en H. Nótese, que la función característica ϕX : H ∗ → C,
con ϕX (f ∗ ) = E[exp(iX(f ∗ ))], f ∗ ∈ H ∗ , de una variable aleatoria cilíndrica es
definida positiva y continua en subespacios finitos. Entonces, existe una medi-
da cilíndrica µX con la misma función característica. Llamamos µX a la distri-
bución cilíndrica de X. El recíproco también es cierto. Una medida cilíndrica
fuertemente gaussiana µ, tiene como función característica ϕ dada por
∗ 1 ∗ ∗
ϕ(f ) = exp − hQf , f i) , ∀f ∗ ∈ H ∗ ,
2
48
Capítulo 1. P RELIMINARES
49
Capítulo 1. P RELIMINARES
Si existe un valor positivo j0 ≥ 1 tal que kρj0 kL(H) < 1, entonces, el proceso
ARH(1) X = {Xn , n ∈ Z} en (1.10) es estándar, y existe una única solución
estacionaria para la ecuación (1.10); véase la Proposición 10.5, p. 248, en Bosq
(2000).
Los operadores de autocovarianza y covarianza cruzada vienen dados res-
pectivamente por
CX = E[Xn ⊗ Xn ] = E[X0 ⊗ X0 ], n ∈ Z,
DX = E[Xn ⊗ Xn+1 ] = E[X0 ⊗ X1 ], n ∈ Z,
(1.11)
50
Capítulo 1. P RELIMINARES
para f, g ∈ H,
f ⊗ g(h) = f hg, hiH , ∀h ∈ H,
51
1.3. Preliminares para el Capítulo 4 y 5
Definición 1.3. (Ruiz-Medina, 2011). Un proceso funcional espacial
X = {Xi,j , (i, j) ∈ Z2 }, con valores en un espacio separable de Hilbert, H, se dice
que es un proceso autorregresivo espacial de orden uno, evaluado en un espacio de
Hilbert H, SARH(1), si es estacionario y cumple la siguiente ecuación:
lo que equivale a la existencia de un entero j0 ≥ 1 tal que kLjl 0 kL( S) } < 1, para
l = 1, 2, 3. De hecho, podemos considerar un vector, j0 = (j01 , j02 , j03 ), tal que
53
Capítulo 1. P RELIMINARES
p p p p
Xi,j = λp (L1 )Xi−1,j + λp (L2 )Xi,j−1 + λp (L3 )Xi−1,j−1 + pi,j , p ≥ 1, (1.15)
p
donde Xi,j = hXi,j , φp iH , y p (i, j) = hi,j , φp iH .
Nota 1.5. Bajo la suposición más restrictiva de que Ll , con l = 1, 2, 3, son operado-
res de Hilbert-Schmidt (positivos y autoadjuntos, como antes), se pueden obtener
resultados similares al desarrollo posterior para una base ortonormal arbitraria,
{ϕk }k≥1 de H, más allá de la restricción de un sistema común de autovectores pa-
ra Ll , l = 1, 2, 3, en la Suposición 1.1, ya que podemos considerar las siguientes
expansiones en serie de Ll , l = 1, 2, 3, en H,
∞ X
X ∞
Ll = Ll (ϕq )(ϕk )ϕk ⊗ ϕq , l = 1, 2, 3.
k=1 q=1
54
Capítulo 1. P RELIMINARES
I J I−1 J−1
1 XX 1 XX
R̂0,0 = Zi,j ⊗ Zi,j ; R̂1,1 = Zi+1,j+1 ⊗ Zi,j ;
IJ i=1 j=1 (I − 1)(J − 1) i=1 j=1
I J−1 I−1 XJ
1 XX 1 X
R̂0,1 = Zi,j+1 ⊗ Zi,j ; R̂1,0 = Zi+1,j ⊗ Zi,j ;
I(J − 1) i=1 j=1 (I − 1)J i=1 j=1
55
Capítulo 1. P RELIMINARES
ricos de autocovarianza:
−1
Λ̂(L1 ) C(R̂0,0 ) C(R̂1,1 ) C(R̂0,1 ) C(R̂1,0 )
Λ̂(L2 ) = C(R̂1,1 ) C(R̂0,0 ) C(R̂1,0 ) C(R̂0,1 )
donde,
para p, q ∈ {0, 1}, C(R̂p,q ) es una matriz M × M de entradas
R̂p,q (φm )(ψn ) y, para l = 1, 2, 3, Λ̂(Ll ) es una matriz diagonal, donde
m,n=1,...M
56
M ETODOLOGÍA Y R ESULTADOS
58
Metodología y Resultados
59
Capítulo 2
60
Capítulo 2. Análisis de las propiedades asintóticas de estimadores espaciales de
mínimo contraste basados en el periodograma tapered
61
Capítulo 2. Análisis de las propiedades asintóticas de estimadores espaciales de
mínimo contraste basados en el periodograma tapered
Restringimos nuestro estudio a los campos con d ≤ 3. Nótese que los re-
sultados sobre la consistencia de los estimadores se mantendrán para el caso
general, d ≥ 1. Sin embargo, para la normalidad asintótica se le impondrá la
restricción d ≤ 3, ya que sólo para estas dimensiones podemos controlar el
sesgo con la ayuda del tapering.
Nota 2.1. En lo que sigue, para la diferenciabilidad con respecto a θ, nos referimos
a la diferenciación en el interior de Θ.
adicionalmente suponemos:
ble fijo D ⊂ T , tal que, con probabilidad 1,
sup Y (t) = sup Y (t), ı́nf Y (t) = ı́nf Y (t), para todo conjunto abierto V,
t∈V ∩D t∈V t∈V ∩D t∈V
62
Capítulo 2. Análisis de las propiedades asintóticas de estimadores espaciales de
mínimo contraste basados en el periodograma tapered
B4. Las derivadas ∇θ ψ(λ, θ) existen y se verifican las condiciones que permiten
intercambiar la diferenciación bajo el signo de la integral en la ecuación
(2.3), es decir que
Z Z
∇θ ψ(λ, θ)w(λ) dλ = ∇θ ψ(λ, θ)w(λ) dλ = 0.
T T
Supongamos que existen las funciones αi : Θ → (0, 1), i = 1, ..., d, tales que
cumplen las condiciones siguientes (con la elección apropiada de una función
peso w(λ)). Para simplificar la notación en la formulación de las condiciones
siguientes, omitiremos el argumento θ en las funciones αi (θ) y las escribiremos
simplemente como αi .
B5. Para todo θ ∈ Θ, la densidad espectral f (λ, θ) = O( di=1 |λi |−αi ) cuando
Q
2 Qd
1. w(λ) ∂θ∂i ∂θj log ψ(λ, θ) = O( i=1 |λi |αi ) cuando λi → 0, para todo i, j,
θ ∈ Θ.
2. w(λ) ∂θ∂ i log ψ(λ, θ) = O( di=1 |λi |αi ) cuando λi → 0, para todo i, θ ∈
Q
Θ.
∂2
3. Las derivadas de segundo orden, ∂θi ∂θj
log ψ(λ, θ), con i, j = 1, . . . , q,
son continuas en θ.
B8. Las matrices S(θ) = (sij (θ))i,j=1,...,q y A(θ) = (aij (θ))i,j=1,...,q son definidas
63
Capítulo 2. Análisis de las propiedades asintóticas de estimadores espaciales de
mínimo contraste basados en el periodograma tapered
∂2
Z
(sij (θ)) = f (λ, θ)w(λ) log ψ(λ, θ) dλ
T ∂θi ∂θj
Z h ∂2
2
= σ (θ) w(λ) ψ(λ, θ)
T ∂θi ∂θj
1 ∂ ∂ i
− ψ(λ, θ) ψ(λ, θ) dλ,
ψ(λ, θ) ∂θi ∂θj
Z
∂ ∂
(aij (θ)) = 2(2π)d f 2 (λ, θ)w2 (λ) log ψ(λ, θ) log ψ(λ, θ) dλ
T ∂θi ∂θj
Z
d 2 2 ∂ ∂
= 2(2π) (σ (θ)) w2 (λ) ψ(λ, θ) ψ(λ, θ) dλ.
T ∂θi ∂θj
H2. La función taper h̃(t) es Lipschitz continua en [−1, 1], y h̃(−1) = h̃(1) = 0.
Nota 2.3. Las condiciones B7 y B8 se utilizan para lograr una convergencia apro-
piada del sesgo a 0.
Nota 2.4. Un ejemplo de un campo aleatorio con dependencia de largo rango para
el cual se satisfacen todas las condiciones anteriores (con la elección apropiada de
64
Capítulo 2. Análisis de las propiedades asintóticas de estimadores espaciales de
mínimo contraste basados en el periodograma tapered
y Z
U (θ) = − f (λ, θ0 )w(λ) log ψ(λ, θ) dλ. (2.6)
T
65
Capítulo 2. Análisis de las propiedades asintóticas de estimadores espaciales de
mínimo contraste basados en el periodograma tapered
es un estimador consistente del vector paramétrico θ. Esto es, existe una conver-
gencia en P0 probabilidad:
0P
θ̂T −→ θ0 , cuando T −→ ∞.
P
El estimador σ̂T2 −→ 0
σ 2 (θ0 ), cuando T −→ ∞, donde σ̂T2 es un estimador del
parámetro σ 2 (θ0 ) que viene dado por
Z
σ̂T2 = ITh (λ)w(λ) dλ.
T
D
T d/2 (θ̂T − θ0 ) −→ Nq (0, e(h)S−1 (θ0 )A(θ0 )S−1 (θ0 )), (2.8)
donde las entradas de las matrices S(θ) = (sij (θ)), y A(θ) = (aij (θ)) están defi-
nidas en la Condición B8, e(h) viene dado por la siguiente ecuación,
Z Z −2 d
4 2
e(h) = (h̃(t)) dt (h̃(t)) dt , (2.9)
66
Capítulo 2. Análisis de las propiedades asintóticas de estimadores espaciales de
mínimo contraste basados en el periodograma tapered
siendo válidos para un dominio de este tipo. Así, solamente tenemos que ajustar los
supuestos para una función taper h̃(t). Es decir, la Suposición H1 debe modificarse
de la siguiente manera: h̃(t) es una función medible positiva de variación acotada
con soporte en [0, 1] y h(0) = 0, h(1 − v) = h(v), para 0 ≤ v ≤ 21 .
Nota 2.6. Un ejemplo de una función taper h̃(t) que satisface las condiciones
asumidas es
1
h̃(t) = (1 + cos(4πt)), t ∈ [−1, 1].
2
Esto es simplemente una modificación de la conocida campana de coseno (o la
función taper de Tukey-Hanning):
1
h̃(t) = (1 − cos(2πt)), t ∈ [0, 1],
2
67
Capítulo 2. Análisis de las propiedades asintóticas de estimadores espaciales de
mínimo contraste basados en el periodograma tapered
En el anterior Teorema 2.2, las condiciones que ayudan a controlar el sesgo son
H1–H2 y B9.
Nota 2.8. Obsérvese que los resultados anteriores sobre estimadores de Ibragi-
mov2 se expresaron bajo las condiciones de integrabilidad de la densidad espectral
f (λ, θ) y la función w(λ) log ψ(λ, θ), así como de funcionales que involucran a di-
chas funciones (2.10). Cabe destacar también que la investigación del sesgo para
los estimadores de funcionales espectrales, en el ajuste no paramétrico, fue presen-
tada en Anh et al. (2007a); Sakhno (2007, 2014).
Las condiciones asumidas para los resultados de consistencia y normalidad
asintótica, se formulan de forma explícita, para facilitar su verificación y apli-
cación en la práctica.
Nota 2.9. Como se desprende de los resultados de Guyon, véase, por ejemplo,
Guyon (1995), resumidos en la siguiente sección: para controlar el sesgo de los
estimadores de mínimo contraste de Ibragimov en el caso de campos aleatorios
con dependencia de corto rango se pueden usar las siguientes condiciones: la fun-
ción taper h̃(t) pertenece a C 2 [−1, 1]; la densidad espectral f (λ, θ) y la función
ϕ(λ, θ) = w(λ) ∂θ∂ i log ψ(λ, θ) pertenecen a C 2 (T).
2.2. Demostraciones
Las demostraciones se basan en los resultados sobre propiedades asintóticas
de los funcionales JT (ϕ) obtenidos mediante la evaluación de sus cumulantes.
2
Véase, por ejemplo, Anh et al. (2004a,b, 2007b,c); Avram et al. (2010a).
68
Capítulo 2. Análisis de las propiedades asintóticas de estimadores espaciales de
mínimo contraste basados en el periodograma tapered
P
lı́m UT (θ) − UT (θ0 ) → U (θ) − U (θ0 ) = K(θ0 ; θ)
T →∞
69
Capítulo 2. Análisis de las propiedades asintóticas de estimadores espaciales de
mínimo contraste basados en el periodograma tapered
Por lo tanto, en vista del Teorema 3.4.1 de Guyon (1995), concluimos que
P0
el estimador θ̂T es consistente: θ̂T −→ θ0 , cuando T → ∞.
por lo tanto, K(θ0 ; θ) ≥ 0 y, además, K(θ0 ; θ) > 0 si ψ(λ, θ0 ) 6≡ ψ(λ, θ), para
θ0 6= θ, casi seguramente, respecto a la medida de Lebesgue.
La siguiente convergencia,
Z Z
P0
σ̂T2 = ITh (λ)w(λ) dλ 2
−→ σ (θ0 ) = f (λ, θ0 )w(λ) dλ (2.12)
T T
Demostración del Teorema 2.2: Aplicando el teorema del valor medio, se tiene
que
∇θ UT (θ̂T ) = ∇θ UT (θ0 ) + ∇θ ∇0θ UT (θT∗ )(θ̂T − θ0 ),
70
Capítulo 2. Análisis de las propiedades asintóticas de estimadores espaciales de
mínimo contraste basados en el periodograma tapered
donde A(θ) y S(θ) vienen dados en la Condición B8. Entonces, por el lema de
Slutsky, la relación (2.8) es una consecuencia de (2.13) y (i)–(ii).
Obsérvese que al cumplirse la Condición B4,
Z
f (λ, θ0 )w(λ)∇θ log ψ(λ, θ0 )dλ = 0,
T
por lo tanto, ∇θ UT (θ0 ) = JT (ϕi ) − J(ϕi ) i=1,...,q , y la convergencia (ii) se cum-
plirá si se produce la siguiente convergencia
D
T d/2 JT (ϕi ) − EJT (ϕi ) i=1,...,q → Nq (0, e(h)A(θ0 ))
(2.14)
y
T d/2 EJT (ϕi ) − J(ϕi ) → 0, i = 1, . . . , q,
(2.15)
donde ϕi (λ) = w(λ) ∂θ∂ i log ψ(λ, θ0 ). La convergencia (2.14) se cumple bajo las
condiciones B5 y B7(2), en vista del Teorema 2.4; la convergencia (i) se cum-
ple bajo las condiciones B5, B6, B7(1), B7(3), teniendo en cuenta la Nota 2.11
(y teniendo en cuenta la consistencia de θ̂T ). La convergencia (2.15) se cum-
71
Capítulo 2. Análisis de las propiedades asintóticas de estimadores espaciales de
mínimo contraste basados en el periodograma tapered
plirá bajo las condiciones H1–H2 y B9, en vista del Lema 2.1, que se indica a
continuación. Obsérvese que en los trabajos anteriores, la convergencia (2.15)
se impuso como hipótesis, véase, por ejemplo Anh et al. (2007b); Avram et al.
(2010b).
Lema 2.1. Sea la función taper h̃(t) tal que satisface los supuestos H1 y H2.
Supongamos además que una de las siguientes condiciones se cumple:
(i) f tiene segunda derivada acotada y ϕ ∈ L1 (T);
Demostración:
El análogo de este lema para los campos paramétricos continuos se obtuvo en
Sakhno (2007); los resultados sobre la evaluación del sesgo para los funcionales
espectrales de órdenes superiores en contextos discretos y continuos se estable-
cieron en Anh et al. (2007b) y Avram et al. (2010b). La presente demostración
utiliza las ideas de estos artículos. Tras realizar algunos cálculos sencillos, el
sesgo de JT (ϕ) se puede expresar como sigue:
Z Z Z
E [JT (ϕ) − J (ϕ)] = ϕ(λ)f (u) Φh2,T (u − λ) dudλ − ϕ(λ)f (λ)dλ
Z Z
= ϕ(λ) (f (λ + u) −f (λ)) Φh2,T (u) dudλ (2.17)
Z Z
= f (λ) (ϕ (λ + u) −ϕ(λ)) Φh2,T (λ) dudλ (2.18)
Z
= (g(u) − g(0)) Φh2,T (u) du, (2.19)
donde Φh2,T (u) = [(2π)d H2,T (0)]−1 |H1,T (u)|2 , es una función par que se in-
tegra a 1 (más específicamente, este es el kernel de tipo Fejér; véase la Sec-
ción 2.4), con H1,T (u) que viene dado en (1.4) para k = 1, y denotando
R
g(u) = f (λ)ϕ (λ + u) dλ.
72
Capítulo 2. Análisis de las propiedades asintóticas de estimadores espaciales de
mínimo contraste basados en el periodograma tapered
d d
X ∂f (λ) X
f (λ + u) − f (λ) − ui ≤ const · |ui |2 .
i=1
∂λi i=1
Consideremos la expresión
d Z
|H1,T (u)|2
Z
X ∂f (λ)
ϕ (λ) ui du dλ
i=1 T ∂λi T (2π)d H2,T (0)
Qd 2
d Z
X ∂f (λ)
Z
i=1 H̃1,T (ui )
= ϕ (λ) ui du dλ . (2.20)
i=1 T ∂λi T (2π)d H2,T (0)
2
Como H̃1,T (ui ) ·(2π H̃2,T (0))−1 es un kernel de tipo Fejér e integra 1, siguiendo
lo planteado en Dahlhaus (1983), las integrales dadas en (2.20) se reducen a
la expresión
Z π 2
H̃1,T (ui )
ui dui ,
−π 2π H̃2,T (0)
2
la cual es igual a cero, ya que H̃1,T (ui ) es una función par. Por lo tanto, (2.20)
es igual a cero y
2
d Z
X Z d
Y H̃1,T (ui )
|EJT (ϕ) − J (ϕ)| ≤ const · |ϕ (λ)| |ui |2 du dλ. (2.21)
i=1 T T i=1
2π H̃2,T (0)
π 2
H̃1,T (ui )
Z
2
|ui | dui , para i = 1, . . . , d.
−π 2π H̃2,T (0)
73
Capítulo 2. Análisis de las propiedades asintóticas de estimadores espaciales de
mínimo contraste basados en el periodograma tapered
Observamos que Z
H̃2,T (0) ∼ T h̃2 (t)dt (2.22)
cuando T → ∞.
A continuación, se procede como en Brillinger (1981), obteniéndose:
X X
H̃1,T (u) = h̃ (t/T ) e−iut = − ∆t (u) h̃ ((t + 1)/T ) − h̃ (t/T ) , (2.23)
t t
74
Capítulo 2. Análisis de las propiedades asintóticas de estimadores espaciales de
mínimo contraste basados en el periodograma tapered
donde di ∈ (−1/2, 1/2), i = 1, 2, φ(B1 , B2 , α, β)Yt1 ,t2 = Yt1 ,t2 − αYt1 −1,t2 −
βYt1 ,t2 −1 + αβYt1 −1,t2 −1 , y ∇d11 ∇d22 = (1 − B1 )d1 (1 − B2 )d2 , con Bi denotando el
operador de retardos para la coordenada ti , i = 1, 2, es decir, B1 Yt1 ,t2 = Yt1 −1,t2 ,
y B2 Yt1 ,t2 = Yt1 ,t2 −1 .
La densidad espectral del proceso definido por (2.26) viene dada por
75
Capítulo 2. Análisis de las propiedades asintóticas de estimadores espaciales de
mínimo contraste basados en el periodograma tapered
Mediante cálculos directos es fácil comprobar que se cumplen todas las con-
diciones de consistencia y normalidad asintótica del estimador θ̂T . En particular,
para la Condición B6, podemos usar la función v(λ) = v(λ1 , λ2 ) = |λ1 |2β |λ2 |2β ,
con β ∈ (0, 1/2). Para verificar la definición positiva de las matrices S(θ) y A(θ),
como se requiere en la Condición B8, podemos seguir análogamente lo consi-
derado de forma muy detallada en Espejo et al. (2015), ver verificación de la
Condición A8.
D
T d/2 (θ̂T − θ0 ) −→ Nq (0, e(h)S−1 (θ0 )A(θ0 )S−1 (θ0 )), (2.27)
donde, para la función taper escogida, tenemos e(h) ≈ 3.78. Los elementos de
las matrices S(θ) y A(θ) se calculan de forma análoga a pp. 671 y 672, en
Espejo et al. (2015). En particular, nótese que
∂ λi
f (λ, θ) = −2 log 2 sen f (λ, θ),
∂θi 2
Z
∂ 2 λi
σ (θ) = −2 log 2 sen w(λ)f (λ, θ)dλ.
∂θi [−π,π]2 2
Así se obtiene,
Z
w(λ) ∂ 2 ∂
sij = 3 2
σ (θ) f (λ, θ) dλ
[−π,π]2 σ (θ) ∂θj ∂θi
Z
w(λ) ∂ 2 ∂
− σ (θ) f (λ, θ) dλ
[−π,π]2 f (λ, θ) ∂θi ∂θi
3 ∂ 2 ∂ 2
= σ (θ) σ (θ)
σ 2 (θ) ∂θj ∂θi
Z
λi λj
−4 log 2 sen log 2 sen w(λ)f (λ, θ)dλ,
[−π,π]2 2 2
76
Capítulo 2. Análisis de las propiedades asintóticas de estimadores espaciales de
mínimo contraste basados en el periodograma tapered
y
Z
2 4 ∂ ∂
aij = 8π σ (θ) w2 (λ)ψ(λ, θ) ψ(λ, θ)dλ
[−π,π]2 ∂θi ∂θj
Z
2 2 ∂ 2 ∂ 2
= 8π w (λ) f (λ, θ) σ (θ) − σ (θ) f (λ, θ)
[−π,π]2 ∂θi ∂θi
∂ 2 ∂ 2
× f (λ, θ) σ (θ) − σ (θ) f (λ, θ) dλ.
∂θj ∂θj
77
Capítulo 2. Análisis de las propiedades asintóticas de estimadores espaciales de
mínimo contraste basados en el periodograma tapered
Para cada T , se usaron 100 valores simulados para estimar P0 (|θbT −θ0 | < ε).
El gráfico de las probabilidades muestrales P0 (|θbT − θ0 | < ε) en la Figura 2.3
confirma la convergencia de θbT = (dˆ1 , dˆ2 ) a θ0 = (0.2, 0.3) en probabilidad.
Finalmente, se investigan las propiedades asintóticas de σ̂T . De manera si-
milar a los casos anteriores, para cada T se construyen los diagramas de cajas
y el gráfico de las probabilidades muestrales para 100 simulaciones. Las figuras
78
Capítulo 2. Análisis de las propiedades asintóticas de estimadores espaciales de
mínimo contraste basados en el periodograma tapered
∆du11 ∆du22 Yt1 ,t2 = (I − 2u1 B1 + B12 )d1 (I − 2u2 B2 + B22 )d2 Yt1 ,t2 = t1 ,t2 ,
donde t1 ,t2 , con (t1 , t2 ) ∈ Z2 , es un campo de ruido blanco de media cero y
varianza común E[2t1 ,t2 ] = σ2 , y, para j = 1, 2,
79
Capítulo 2. Análisis de las propiedades asintóticas de estimadores espaciales de
mínimo contraste basados en el periodograma tapered
80
Capítulo 2. Análisis de las propiedades asintóticas de estimadores espaciales de
mínimo contraste basados en el periodograma tapered
donde di 6= 0, y |ui | ≤ 1, i = 1, 2.
(d)
El polinomio de Gegenbauer Cn (u) viene dado por
[n/2]
X (2u)n−2k Γ(d − k + n)
Cn(d) (u) = (−1)k .
k=0
k!(n − 2k)!Γ(d)
σ2
f (λ, θ) = |1 − 2u1 e−iλ1 + e−2iλ1 |−2d1 |1 − 2u2 e−iλ2 + e−2iλ2 |−2d2
(2π)2
σ2
= 2
|2 cos(λ1 ) − 2u1 |−2d1 |2 cos(λ2 ) − 2u2 |−2d2 .
(2π)
81
Capítulo 2. Análisis de las propiedades asintóticas de estimadores espaciales de
mínimo contraste basados en el periodograma tapered
Para la secuencia creciente de tamaños de muestra, T =10, 30, 50, 70, 90,
110, se generaron 200 repeticiones de campo aleatorio en (2.28), utilizando los
valores de los parámetros d1 =0.2, d2 =0.3, u1 =0.4, u2 =0.3, y σ2 =1.
Análogamente al ejemplo de la Sección 2.3.1, se aproximaron realizaciones
de Yt1 ,t2 por sumas truncadas con 100 términos en (2.28). El periodograma ITh
se calculó con hT (t1 , t2 )=h̃(t1 /T )h̃(t2 /T ), h̃(t) = 12 (1 − cos(2πt)), t ∈ [0, 1]. Se
utilizó la función peso ω(λ1 , λ2 ) = |2 cos(λ1 ) − 2u1 |2 |2 cos(λ2 ) − 2u2 |2 .
El análisis, las gráficas y las conclusiones sobre la ilustración de los resul-
tados obtenidos son análogas al ejemplo de la Sección 2.3.1. Para cada simu-
lación, se encontraron los argumentos de minimización θbT del funcional UT (θ)
en (2.4), ver Figura 2.6 y Tabla 2.2, en relación con los valores muestrales y
errores cuadráticos medios (ECM), respectivamente.
82
Capítulo 2. Análisis de las propiedades asintóticas de estimadores espaciales de
mínimo contraste basados en el periodograma tapered
Los gráficos Q-Q de normalidad de dˆ1 y dˆ2 en la Figura 2.7 coinciden con
las distribuciones normales teóricas. La gráfica de probabilidades muestrales
P0 (|θbT − θ0 | < ε) en la Figura 2.8 también confirma la convergencia en pro-
babilidad. Las figuras 2.9 y 2.10 respaldan la convergencia en probabilidad
P0 (|σ̂T2 − σ 2 (θ)| < ε) → 1 cuando T aumenta.
83
Capítulo 2. Análisis de las propiedades asintóticas de estimadores espaciales de
mínimo contraste basados en el periodograma tapered
84
Capítulo 2. Análisis de las propiedades asintóticas de estimadores espaciales de
mínimo contraste basados en el periodograma tapered
(|λ1 | − arc cos(u1 ))2 (|λ2 | − arc cos(u2 ))2 = |λ1 |2 |λ2 |2 . Resultados similares se ob-
tuvieron también para ω(λ1 , λ2 ) = |2 cos(λ1 ) − 2u1 |2 |2 cos(λ2 ) − 2u2 |2 .
El análisis, las gráficas y las son análogas al ejemplo de la Sección 2.3.1. Los
argumentos de minimización θbT del funcional UT (θ) en (2.4) se encontraron
numéricamente para cada simulación (ver Figura 2.11 y Tabla 2.3).
85
Capítulo 2. Análisis de las propiedades asintóticas de estimadores espaciales de
mínimo contraste basados en el periodograma tapered
Figura 2.11: Diagrama de cajas de los valores muestrales de dˆ1 y dˆ2 con ui = 1,
i = {1, 2}
86
Capítulo 2. Análisis de las propiedades asintóticas de estimadores espaciales de
mínimo contraste basados en el periodograma tapered
87
Capítulo 2. Análisis de las propiedades asintóticas de estimadores espaciales de
mínimo contraste basados en el periodograma tapered
Uno de los enfoques clásicos para llegar al teorema central del límite para
R
los funcionales espectrales, JT (ϕ) = T IT (λ) ϕ (λ) dλ, consiste en calcular y
evaluar sus cumulantes. Se presentan aquí algunos detalles de este enfoque, y
los resultados correspondientes para el caso de campos aleatorios gaussianos y
datos tapered, según se ha considerado en el desarrollo del presente capítulo.
Para establecer el teorema central del límite (TCL) para el funcional normaliza-
do J˜T (ϕ) = T d/2 (JT (ϕ) − EJT (ϕ)) bastará con proporcionar condiciones para
la convergencia al límite finito del cumulante de segundo orden de J˜T (ϕ) , y
proporcionar condiciones para la convergencia a cero del resto de cumulantes
de órdenes superiores.
88
Capítulo 2. Análisis de las propiedades asintóticas de estimadores espaciales de
mínimo contraste basados en el periodograma tapered
×f (λ3 )ϕ(λ4 )...f (λ2k−1 )ϕ(λ2k )H1,T (λ2 − λ1 )H1,T (λ3 − λ2 )...
Z Z
kd/2 k−1
= T 2 (k − 1)! f (λ)ϕ(λ + u1 )
u1 ,...,u2k−1 ∈T2k−1 λ∈T
2k−1
!
X
×f (λ + u1 + u2 )ϕ(λ + u1 + u2 + u3 ) . . . ϕ λ + ui
i=1
2k−1 2k−1
!
Y X
× H1,T (ui )H1,T − ui dλdu1 . . . du2k−1 ,
i=1 i=1
k−1 k−1
!
1 Y X
Φhk,T (λ1 , ..., λk−1 ) := H1,T (λj ) H1,T − λj (2.29)
(2π)d(k−1) Hk,T (0) j=1 j=1
89
Capítulo 2. Análisis de las propiedades asintóticas de estimadores espaciales de
mínimo contraste basados en el periodograma tapered
siempre que la función G sea acotada y continua en el punto (v1 , ..., vk−1 ) .
En particular, se tiene,
Z Z
3d
V ar(J˜T (ϕ)) = 2T (2π) H4,T (0)[H2,T (0)]
d −2
f (λ)ϕ(λ + u1 )
u1 ,u2 ,u3 ∈T3 λ∈T
Z
d 3d −2
= 2T (2π) H4,T (0)[H2,T (0)] f (λ1 )ϕ(λ2 )f (λ3 )ϕ(λ4 )
(λ1 ,λ2 ,λ3 ,λ4 )∈T4
Teorema 2.3. Sea X(t), t ∈ Zd , un campo aleatorio gaussiano de media cero con
densidad espectral f (λ) ∈ Lp y ϕ(λ) ∈ Lq , donde p1 + 1q ≤ 12 . Entonces,
D
T d/2 (JT (ϕ) − EJT (ϕ)) −→ N (0, σ 2 ) cuando T → ∞, (2.32)
donde Z
2 3d
σ = 2(2π) e(h) f 2 (λ)ϕ2 (λ)dλ, (2.33)
T
90
Capítulo 2. Análisis de las propiedades asintóticas de estimadores espaciales de
mínimo contraste basados en el periodograma tapered
Corolario 2.1. Bajo las condiciones del Teorema 2.3 y del Lema 2.1 se tiene,
D
T d/2 (JT (ϕ) − J (ϕ)) −→ N (0, σ 2 ) cuando T → ∞, (2.34)
91
Capítulo 2. Análisis de las propiedades asintóticas de estimadores espaciales de
mínimo contraste basados en el periodograma tapered
Proposición 2.1. Sea X(t), t ∈ Z, un proceso gaussiano con media cero y densi-
dad espectral f (λ).
Entonces,
D
T 1/2 (JT (ϕ) − EJT (ϕ)) −→ N (0, σ 2 ) cuando T → ∞, (2.35)
Rπ
donde σ 2 = 16π 3 −π
ϕ2 (λ)f 2 (λ)dλ.
II. Si existe α < 1 y β < 1 tales que α + β < 1/2, y f (λ) = |λ|−α L1 (λ) y
ϕ(λ) = |λ|−β L2 (λ) cuando λ → 0, donde L1 y L2 son funciones que varían
lentamente en cero, entonces, la Condición (ii) se cumple
D
T d/2 (JT (ϕ) − EJT (ϕ)) −→ N (0, σ 2 ) cuando T → ∞, (2.36)
A partir del Teorema 2.4 y el Lema 2.1 (véase la Sección 2.2) se obtiene el
siguiente corolario.
92
Capítulo 2. Análisis de las propiedades asintóticas de estimadores espaciales de
mínimo contraste basados en el periodograma tapered
Corolario 2.2. Supongamos que se cumplen las condiciones del Teorema 2.4 y del
Lema 2.1. Entonces,
D
T d/2 (JT (ϕ) − J (ϕ)) → N (0, σ 2 ) cuando T → ∞, (2.37)
1
donde ITY (λ) = 2πT | t∈LT hT (t)eiλt Y (t)|2 es el periodograma tapered que co-
P
rresponde a {Y (t), t ∈ LT }.
Dado que la densidad espectral fY (λ) del proceso Y (t) y la función ψ(λ)
satisfacen las condiciones del Teorema 2.3, para el funcional J˜TY (ψ) se tiene la
siguiente convergencia cuando T → ∞,
D
T 1/2 J˜TY (ψ) → N (0, σ 2 ), (2.38)
93
Capítulo 2. Análisis de las propiedades asintóticas de estimadores espaciales de
mínimo contraste basados en el periodograma tapered
donde
Z π Z π
2 3 2
σ = 16π e(h) ψ (λ)fY2 (λ)dλ 3
= 16π e(h) ϕ2 (λ)fX2 (λ)dλ. (2.39)
−π −π
Por lo tanto, para probar la afirmación del teorema, basta con demostrar que
donde Z π Z π
J˜TX (ϕ) = ϕ(λ)ITX (λ)dλ −E ϕ(λ)ITX (λ)dλ.
−π −π
Consideremos
T E|J˜TX (ϕ) − J˜TY (ψ)|2 = T E|J˜TX (ϕ)|2 + T E|J˜TY (ψ)|2 − 2T E J˜TX (ϕ)J˜TY (ψ). (2.41)
Z
T E|J˜TY (ψ)|2 = 2T d (2π)3d H4,T (0)[H2,T (0)]−2 fY (λ1 )ψ(λ2 )fY (λ3 )ψ(λ4 )
[−π,π]4
donde por fXY hemos denotado la densidad espectral cruzada de los procesos
94
Capítulo 2. Análisis de las propiedades asintóticas de estimadores espaciales de
mínimo contraste basados en el periodograma tapered
X(t) e Y (t).
Dado que Y (t) se obtiene a partir de X(t) con el uso del filtrado de la fun-
ción de transferencia D(iλ), entonces, fXY (λ) = D(−iλ)fX (λ), o fXY (λ) =
fY (λ) 2 fY (λ)
D(−iλ) |D(iλ)| 2 , recordemos también que ϕ(λ) = ψ(λ)|D(iλ)| , fX (λ) = |D(iλ)|2 .
95
Capítulo 2. Análisis de las propiedades asintóticas de estimadores espaciales de
mínimo contraste basados en el periodograma tapered
α /2 α /2
Y (t) = Y (t1 , ..., td ) = ∇1 1 ...∇d d X(t)
X∞ ∞ Y
X d
α /2
= ... Ckii X(t1 − k1 , ..., td − kd ),
k1 =0 kd =0 i=1
d
!
αi
Y λ
fY (λ1 , ..., λd ) = 2 sen fX (λ1 , ..., λd ) .
i=1
2
96
Capítulo 3
98
Capítulo 3. Procesos de Cox log-gaussianos temporales dirigidos mediante un
proceso de Ornstein-Uhlenbeck Hilbert-valuado
99
Capítulo 3. Procesos de Cox log-gaussianos temporales dirigidos mediante un
proceso de Ornstein-Uhlenbeck Hilbert-valuado
En particular,
∞
X
exp(tA)(f ) = exp(−λj (A)t) hf, φj iH φj , ∀f ∈ H. (3.3)
j=1
1/2
{ψj = λj (Q)φj }j≥1 (3.4)
define una base ortonormal del espacio de Hilbert con núcleo reproductor de
Q, W admite la siguiente representación:
∞
∗
X k·kQ ∗
W (t)(g ) = hiQ (ψj ), giH Bj (t), ∀g ∗ ∈ U ∗ = [Q1/2 (H) ], (3.5)
j=1
k·kQ
donde U ∗ denota el espacio dual de U = Q1/2 (H) , el cierre de Q1/2 (H),
con respecto a la norma k · kQ , inducida por Q, es decir, kf kQ = Q(f )(f ) =
hQ(f ∗ ), f iH , para todo f ∈ H, con, como anteriormente, f ∗ ∈ H ∗ denotando
el elemento dual de f, dado por el teorema de representación de Riesz. Aquí,
{Bj (t), t ≥ 0}j≥1 son procesos de Wiener estándar real-valuados independien-
tes y, como antes, iQ es la aplicación de inclusión continua que introduce el
rango de Q en H.
Demostración:
Bajo la Condición 3.1, a partir de las ecuaciones (3.3) y (3.5), la proyección de
la ecuación (1.8) en el sistema de autovectores ortonormales {φj }j≥1 nos lleva
100
Capítulo 3. Procesos de Cox log-gaussianos temporales dirigidos mediante un
proceso de Ornstein-Uhlenbeck Hilbert-valuado
hX(t, φj ), φp iH
Z t
= hexp(tA)(φj ), φp iH + exp((t − s)A)dW (s), φp
0 H
∞
X
= exp(−λk (A)t) hφk , φj iH hφk , φp iH
k=1
X∞ ∞ Z
X t
+ hφk , φp iH exp(−λk (A)(t − s)) hφk , iQ (ψl )iH Bl (ds)
k=1 l=1 0
Z t
= exp (−λp (A)t) + λp1/2 (Q) exp(−λp (A)(t − s))Bp (ds),
0
(3.6)
101
Capítulo 3. Procesos de Cox log-gaussianos temporales dirigidos mediante un
proceso de Ornstein-Uhlenbeck Hilbert-valuado
(3.9)
√
Z [t/τ ]
Yp,τ ([t/τ ]) = γτ[t/τ ] Yp (0) + λp1/2 (Q)γ −[s] dBp (s)
τγ [t/τ ]−1
0
Z t
1/2
∼ exp(−λp (A)t) + λp (Q) exp(−λp (A)t) exp(λp (A)s)dBp (s),
0
(3.11)
102
Capítulo 3. Procesos de Cox log-gaussianos temporales dirigidos mediante un
proceso de Ornstein-Uhlenbeck Hilbert-valuado
donde la serie (3.12) es finita para cualquier f ∈ H, y, como antes, {φj }j≥1 denota
el sistema de autovectores del operador de autocovarianza CX dado en (1.11).
A partir de las ecuaciones (3.9)–(3.11), la Suposición 3.1 implica la Supo-
sición 3.2, cuando el proceso considerado ARH(1) se construye como la apro-
ximación en tiempo discreto de un proceso O-U H–valuado, es decir, cuando
ambos procesos están caracterizados por una secuencia de procesos ARH(1)
independientes tal como se da en las ecuaciones (3.9)–(3.10), y la ecuación
(3.13) que aparece en el siguiente resultado:
Lema 3.1. Bajo la Condición 3.2, la proyección de la ecuación ARH(1), (1.10),
en {φj }j≥1 nos lleva a la siguiente expresión: para cada p ≥ 1
103
Capítulo 3. Procesos de Cox log-gaussianos temporales dirigidos mediante un
proceso de Ornstein-Uhlenbeck Hilbert-valuado
T −2
X
1
kT T −1
Xk (i)Xk (i + 1)
X i=0
ρbkT = bT,k (ρ)φk ⊗ φk ,
λ λ
bT,k (ρ) = ,
T −1
X
k=1 1
T
Xk2 (i)
i=0
k = 1, . . . kT ,
T −1
kT
!
bT,k (Rε ) = 1
X X
R
bε = bT,k (Rε )φk ⊗ φk ,
λ λ Xk2 (i) bT,k (ρ)]2 ,
1 − [λ
k=1
T i=0
k = 1, . . . kT ,
(3.14)
kT
X
A
bT = (1 − λ
bT,k (ρ))φk ⊗ φk (3.15)
k=1
XkT
Q
bT = bT,k (Rε )φk ⊗ φk ,
λ (3.16)
k=1
104
Capítulo 3. Procesos de Cox log-gaussianos temporales dirigidos mediante un
proceso de Ornstein-Uhlenbeck Hilbert-valuado
u2 dk − u2
Hk (u) = (−1)k e 2 e 2, k = 0, 1, . . . .
duk
105
Capítulo 3. Procesos de Cox log-gaussianos temporales dirigidos mediante un
proceso de Ornstein-Uhlenbeck Hilbert-valuado
E [Hk (Y (x))] = 0,
E [Hk (Y (x)) Hm (Y (y))] = δm,k m! (E[Y (x)Y (y)])m ; (3.18)
(3.20)
Demostración:
A la ecuación (3.20) se llega directamente a partir de (3.18), y de la defini-
ción de Yp en la ecuación (3.7), para cada p ≥ 1. Además, G(u) = exp(u) ∈
L2 (R, ϕ(u)du) implica
∞ 2
X Ck
kGk2L2 (R,ϕ(u)du) = < ∞. (3.21)
k=1
k!
106
Capítulo 3. Procesos de Cox log-gaussianos temporales dirigidos mediante un
proceso de Ornstein-Uhlenbeck Hilbert-valuado
h i2
λp (Q)
Condición (3.19), nos lleva a que 2λp (A)
< 1, con lo que se obtiene:
∞ k
Ck2 t t λp (Q)
X Z Z
exp (−λp (A)|u − v|) dudv
k=1
k! 0 0 2λp (A)
v v #2
u ∞ 2 2 u ∞
"Z Z
t t k
uX C k u X λ p (Q)
≤t t exp (−λp (A)|u − v|) dudv
k=1
k! k=1 0 0 2λ p (A)
v v #2
u ∞ 2 2 u ∞
"Z Z
t t k
uX C k u X λ p (Q)
≤t t exp (−λp (A)|u − v|) dudv
k=1
[k!]2 k=1 0 0 2λp (A)
v v
u ∞ 2 u ∞
uX Ck 2 uX λp (Q) 2k
≤Ket tt
k=1
k! k=1
2λp (A)
v
u ∞ 2
uX Ck t2
=t
k!
r h i2 < ∞, (3.22)
k=1 λ (Q)
1 − 2λpp (A)
donde
n K2es
e
o una constante positiva, sabiendo que (3.21) implica que la secuen-
Ck
cia k!
converge a cero y, en particular, para k suficientemente grande,
k≥1
C 4 C 2
se tiene que k!k < k!k < 1.
En el desarrollo subsiguiente se asume la siguiente condición:
2λp (A)
= O p−α ,
p p → ∞, α > 1.
[2λp (A)]2 − [λp (Q)]2
Demostración:
A partir de la ecuación (3.22), bajo la Condición (3.19), y las condiciones 3.1 y
107
Capítulo 3. Procesos de Cox log-gaussianos temporales dirigidos mediante un
proceso de Ornstein-Uhlenbeck Hilbert-valuado
3.3,
" ∞
# ∞
X X
Cp2 (t) E Cp2 (t)
E =
p=1 p=1
v
∞ uX
u ∞ Ck 2 t2
X
≤ K
e t r
p=1 k=1
k! h
λp (Q)
i2
1 − 2λp (A)
v
u ∞ 2 ∞ "
X Ck X 2 #−1/2
e 2t
u λ p (Q)
= Kt 1−
k=1
k! p=1
2λ p (A)
v
u ∞ 2 ∞
2
uX C k X 2λp (A)
= Kt
e t p < ∞.
k=1
k! p=1
[2λ p (A)] 2 − [λ (Q)]2
p
108
Capítulo 3. Procesos de Cox log-gaussianos temporales dirigidos mediante un
proceso de Ornstein-Uhlenbeck Hilbert-valuado
p0 (ω,Y
b T (ω),T,T +h) Z
T +h
Y
EYb T (ω) (C(T + h)|C(ω, T )) = exp Ybp (ω, s) ds + C(ω, T ),
p=1 T
donde Yb T (ω) = {Ybp (ω, s), p ≥ 1, T ≤ s ≤ T + h}, siendo Ybp (·), tal y como se
ha definido en (3.23), para cada p ≥ 1, el predictor de la intensidad aleatoria,
calculado a partir de las observaciones en el intervalo [0, T ], es decir,
109
Capítulo 3. Procesos de Cox log-gaussianos temporales dirigidos mediante un
proceso de Ornstein-Uhlenbeck Hilbert-valuado
110
Capítulo 3. Procesos de Cox log-gaussianos temporales dirigidos mediante un
proceso de Ornstein-Uhlenbeck Hilbert-valuado
\
Φ(C(T + h))(xi )
, i = 1, . . . , N,
ET +h (xi )
p0 (ω,Y
b T (ω),T,T +h) Z
T +h
Y
\
C(T + h) = bp (ω, s)ds + C(ω, T ),
λ (3.26)
p=1 T
Se generan tres muestras funcionales, bajo los supuestos 3.1, 3.2 y 3.3, con-
dicionadas a los registros disponibles sobre casos esperados de cáncer de mama,
próstata y encéfalo, en las provincias de España. Además, siguiendo las expre-
siones (3.24)–(3.25), nuestro modelo de observación se define a partir de la
siguiente ecuación, extraída de considerar el desarrollo de Taylor en el ln de
i)
la función ECtt(x
(xi )
1
en el valor observado λt (xi ), para cada t = 1, . . . , T, y para
1
En el caso de estudio de mortalidad, esta función se conoce como Razón de Mortalidad
Estandarizada (RME) y se trata de un estimador del Riesgo Relativo (RR), véase, por ejemplo,
Bivand et al. (2013), entre otros. Consideremos Ot,i y Et,i como el número de muertes ob-
servadas y esperadas (respectivamente) causadas por una enfermedad en el año t = 1, ..., T ,
111
Capítulo 3. Procesos de Cox log-gaussianos temporales dirigidos mediante un
proceso de Ornstein-Uhlenbeck Hilbert-valuado
cualquier i = 1, . . . , N,
Ct (xi )
Zet (xi ) = ln = ln(Ct (xi )) − ln Et (xi )
Et (xi )
Ct (xi ) − λt (xi )
= ln(λt (xi )) + − ln Et (xi )
i
λt (x )
i i i i
Ct (xi ) − λt (xi )
= ln Et (x ) + βt (x ) + Xt (x ) − ln Et (x ) +
λt (xi )
Ct (xi ) − λt (xi )
= βt (xi ) + Xt (xi ) + i
= βt (xi ) + Xt (xi ) + ξt (xi ),
λt (x )
(3.27)
Ct (xi ) − λt (xi )
i
ξt (x ) = , i = 1, . . . , N, t = 1, . . . , T (3.28)
λt (xi )
con t = 1, ..., T , i = 1, . . . , N .
112
Capítulo 3. Procesos de Cox log-gaussianos temporales dirigidos mediante un
proceso de Ornstein-Uhlenbeck Hilbert-valuado
113
Capítulo 3. Procesos de Cox log-gaussianos temporales dirigidos mediante un
proceso de Ornstein-Uhlenbeck Hilbert-valuado
−1
Kt = Pt|t−1 Φ∗M Cξt + ΦM Pt|t−1 Φ∗M
siendo Rεt la matriz de covarianzas del vector aleatorio εM ×1 (t), con las com-
ponentes εp (t) = hεt , φp iH , p = 1, . . . , M, y
Pt|t = E (X(t) − X(t b | t))∗ = Pt|t−1 − Kt ΦM Pt|t−1 .
b | t))(X(t) − X(t
X b | t − 1) = ΦM Λ(ρ)X(t
bt = ΦM X(t b − 1 | t − 1). (3.31)
b | 0) = 0
X(0
P0|0 = E XM ×1 (0)[XM ×1 (0)]T ,
rbt = exp(βbt + X
bt ), (3.32)
114
Capítulo 3. Procesos de Cox log-gaussianos temporales dirigidos mediante un
proceso de Ornstein-Uhlenbeck Hilbert-valuado
Paso 1. Partimos de los registros disponibles del logaritmo del Riesgo Relati-
vo (RR) de mortalidad asociado al cáncer de mama, próstata y encéfalo,
según se define en la ecuación (3.27), primera línea. A continuación se
procede a la eliminación de la tendencia lineal.
115
Capítulo 3. Procesos de Cox log-gaussianos temporales dirigidos mediante un
proceso de Ornstein-Uhlenbeck Hilbert-valuado
Tabla 3.1: Promedio espacial de los errores cuadráticos medios empíricos (ECMEs) por
año, obtenido a partir de 100 generaciones de la muestra funcional de tamaño T del
proceso de intensidad ARH(1), para los datos generados de los tres tipos de cáncer
116
Capítulo 3. Procesos de Cox log-gaussianos temporales dirigidos mediante un
proceso de Ornstein-Uhlenbeck Hilbert-valuado
para i = 1, . . . , N, y t ∈ R+ . Aquí,
117
Capítulo 3. Procesos de Cox log-gaussianos temporales dirigidos mediante un
proceso de Ornstein-Uhlenbeck Hilbert-valuado
Figura 3.1: Mapas de riesgo relativo de mortalidad por cáncer de próstata: para una
generación (figura superior) y su estimación (figura inferior), basada en la aproxima-
ción en tiempo discreto ARH(1) del modelo de intensidad log-gaussiana
118
Capítulo 3. Procesos de Cox log-gaussianos temporales dirigidos mediante un
proceso de Ornstein-Uhlenbeck Hilbert-valuado
Figura 3.2: Mapas de riesgo relativo de mortalidad por cáncer de mama: para una ge-
neración (figura superior) y su estimación (figura inferior), basada en la aproximación
en tiempo discreto ARH(1) del modelo de intensidad log-gaussiana
119
Capítulo 3. Procesos de Cox log-gaussianos temporales dirigidos mediante un
proceso de Ornstein-Uhlenbeck Hilbert-valuado
Figura 3.3: Mapas de riesgo relativo de mortalidad por cáncer de encéfalo: para una
generación (figura superior) y su estimación (figura inferior), basada en la aproxima-
ción en tiempo discreto ARH(1) del modelo de intensidad log-gaussiana
120
Capítulo 3. Procesos de Cox log-gaussianos temporales dirigidos mediante un
proceso de Ornstein-Uhlenbeck Hilbert-valuado
bt = R e b R−1
X e −βb ) (Zt − βt ), ∀t ∈ R+ , (3.34)
e b
Xt (Zt −βt ) (Z
e −βb )(Z
t t t t
donde
ν̂ $1 $2
Cáncer de próstata 0.6997 0.4103 0.5483
Cáncer de mama 0.6667 0.4441 0.4356
Cáncer de encéfalo 0.6883 0.4555 0.55435
121
Capítulo 3. Procesos de Cox log-gaussianos temporales dirigidos mediante un
proceso de Ornstein-Uhlenbeck Hilbert-valuado
Paso 1. Partimos de los registros disponibles sobre el logaritmo del riesgo re-
lativo de mortalidad asociado al cáncer de próstata, mama o encéfalo,
según se define en la ecuación (3.27), primera línea.
122
Capítulo 3. Procesos de Cox log-gaussianos temporales dirigidos mediante un
proceso de Ornstein-Uhlenbeck Hilbert-valuado
Tabla 3.3: Promedio espacial de los errores cuadráticos medios empíricos (ECMEs)
por año, obtenidos a partir de 100 generaciones de la muestra funcional de tamaño
T del modelo de intensidad aleatoria con DLR, para los tres tipos de datos de cáncer
generados
123
Capítulo 3. Procesos de Cox log-gaussianos temporales dirigidos mediante un
proceso de Ornstein-Uhlenbeck Hilbert-valuado
Paso 7. Para cada tiempo, t = 1, ..., T , y provincia, i = 1, ..., 50, se genera una
muestra de valores aleatorios de una distribución de Poisson con paráme-
tro, el elemento (t, i) de la matriz obtenida en el Paso 6.
124
Capítulo 3. Procesos de Cox log-gaussianos temporales dirigidos mediante un
proceso de Ornstein-Uhlenbeck Hilbert-valuado
Figura 3.4: Mapas de riesgo relativo de mortalidad por cáncer de próstata: para una
generación (figura superior) y su estimación (figura inferior), basada en el modelo de
intensidad log-gaussiana con DLR
125
Capítulo 3. Procesos de Cox log-gaussianos temporales dirigidos mediante un
proceso de Ornstein-Uhlenbeck Hilbert-valuado
Figura 3.5: Mapas de riesgo relativo de mortalidad por cáncer de mama: para una
generación (figura superior) y su estimación (figura inferior), basada en el modelo de
intensidad log-gaussiana con DLR
126
Capítulo 3. Procesos de Cox log-gaussianos temporales dirigidos mediante un
proceso de Ornstein-Uhlenbeck Hilbert-valuado
Figura 3.6: Mapas de riesgo relativo de mortalidad por cáncer de encéfalo: para una
generación (figura superior) y su estimación (figura inferior), basada en el modelo de
intensidad log-gaussiana con DLR
127
Capítulo 3. Procesos de Cox log-gaussianos temporales dirigidos mediante un
proceso de Ornstein-Uhlenbeck Hilbert-valuado
Figura 3.7: Enfoque O-U log-gaussiano. Mapas de España con el log-riesgo de morta-
lidad por tipo de cáncer: de próstata (año 1975), de mama (año 1975), y de encéfalo
(año 1986); en base a 1a generación de un proceso O-U (de la parte superior izquierda
a la inferior izquierda) y sus respectivas estimaciones con el filtrado de Kalman (de la
parte superior derecha a la inferior derecha)
128
Capítulo 3. Procesos de Cox log-gaussianos temporales dirigidos mediante un
proceso de Ornstein-Uhlenbeck Hilbert-valuado
Figura 3.8: Enfoque log-gaussiano con DLR. Mapas de España con el log-riesgo de
mortalidad por tipo de cáncer: de próstata (año 1975), de mama (año 1975), y de
encéfalo (año 1986); en base a 1 generación del proceso de Cox con intensidad log-
gaussiana con DLR (de la parte superior izquierda a la inferior izquierda) y sus res-
pectivas estimaciones de mínimos cuadrados de DLR (de la parte superior derecha a la
inferior derecha)
129
Capítulo 4
131
Capítulo 4. Procesos de Cox espaciales en un marco infinito-dimensional
132
Capítulo 4. Procesos de Cox espaciales en un marco infinito-dimensional
medibles, desde (Ω, A, P ) a la recta real o compleja, así como desde (Ω, A, P )
a un espacio de funciones. Sólo reflejamos dicha dependencia la primera vez,
cuando se introducen distribuciones de probabilidad condicionadas a la obser-
vación de valores muestrales. Normalmente, el subíndice denota la localización
del elemento de una familia de operadores o variables aleatorias dada, y el
argumento se refiere al elemento ubicado en el conjunto que define su soporte.
Sea H un espacio de Hilbert real separable de funciones, denotamos por H+
iH su versión compleja. En lo que sigue, h·, ·i y k · k denotan respectivamente el
producto interno y la norma en el espacio H+iH. La misma notación se utilizará
para el producto interno y la norma de H, considerado como un subespacio de
H + iH. A efectos prácticos nos referimos a H = L2 (T ), como el espacio de las
funciones de cuadrado integrable en el intervalo de tiempo T . En esta sección
B representa un conjunto acotado de Borel B ∈ B d .
Se denota por κz , z ∈ Rd a un proceso espacial con valores en H. Al apli-
car el Teorema de representación de Riesz, podemos definir un funcional alea-
torio Xz , sobre H utilizando la identidad
Xz (ϕ) = hκz , ϕi , ∀ϕ ∈ H, z ∈ Rd .
133
Capítulo 4. Procesos de Cox espaciales en un marco infinito-dimensional
Rd } definida como
∞
X [κz (t)]k
Λz (t) = exp (κz (t)) = , ∀t ∈ T , z ∈ Rd , (4.2)
k=1
k!
en la norma del espacio L2 (Ω, A, P ), dado por kXk2L2 (Ω,A,P ) = E[X 2 ], para cual-
quier variable aleatoria de segundo orden con media cero X en (Ω, A, P ).
134
Capítulo 4. Procesos de Cox espaciales en un marco infinito-dimensional
Z
= exp exp (xz (ϕ)) dz (exp(t) − 1) , t ∈ R, ∀ϕ ∈ H. (4.4)
B
135
Capítulo 4. Procesos de Cox espaciales en un marco infinito-dimensional
136
Capítulo 4. Procesos de Cox espaciales en un marco infinito-dimensional
(2)
ρϕ (zi − zj )
gϕ (zj − zi ) = = exp Rz i −z j
(ϕ)(ϕ) , ∀ϕ ∈ H.
[ρϕ ]2
(4.13)
d
!
1 X X
Fξ (h)(g) = exp −i ξj zj Rz (h)(g), ξ ∈ [−π, π]d , (4.14)
(2π)d d j=1
z∈Z
137
Capítulo 4. Procesos de Cox espaciales en un marco infinito-dimensional
P
1
P d
para todo h, g ∈ H+iH. Asumimos que la serie (2π) d z∈Zd exp −i j=1 j j Rz
ξ z
converge en la norma de la traza, véase, por ejemplo, Panaretos y Tavakoli
(2013b). Como antes, H + iH denota la versión complejo-valuada del espacio
de Hilbert real separable H, y el operador de covarianza Rz se ha introducido
en la ecuación (4.10). Asumimos la siguiente condición
NX 2
d −zd
" "N −z # #
1 1
X 1 X
··· Xy (ψk )Xy+z (ψl ) − Rz (ψk )(ψl ) →c. s. 0,
k,l≥1
N y1 =1 yd =1
Qd
donde se ha utilizado la notación y = (y1 , . . . , yd ) ∈ i=1 [1, Ni − zi ] ∩ Zd .
Nota 4.4. Como es habitual, la consistencia fuerte de una variable aleatoria eva-
luada en H se define a partir de la convergencia c. s. en la norma k · kH . Este hecho
se aplica en la Suposición 4.2, considerando el espacio de Hilbert S(H) de los ope-
radores de Hilbert–Schmidt sobre H. Así, la Suposición 4.2 define la consistencia
fuerte de los momentos de segundo orden funcionales empíricos de κz con respecto
a la norma S(H).
138
Capítulo 4. Procesos de Cox espaciales en un marco infinito-dimensional
139
Capítulo 4. Procesos de Cox espaciales en un marco infinito-dimensional
−1
Condición 4.6. La suma de Cesaro de las series de Fourier de Fξ,θ (ϕk )(ϕk ) con-
d
verge uniformemente en ξ ∈ [−π, π] , k ≥ 1, y θ ∈ Θ.
Nota 4.6. Como se demuestra en el siguiente Teorema 4.1, las condiciones 4.3–
4.6 proporcionan un escenario adecuado para la estimación consistente en sentido
fuerte en el dominio espectral funcional espacial del modelo paramétrico adoptado
para los autoelementos del operador de correlación espacial. Véanse también los
ejemplos 4.1 y 4.2 de la Sección 4.3.
N1 Nd d
! d
1 X X X Y
e ξN = p
κ ··· exp −i ξj yj κy , ξ ∈ [−π, π]d , N = Ni ,
N (2π)d y1 =1 yd =1 j=1 i=1
(4.17)
140
Capítulo 4. Procesos de Cox espaciales en un marco infinito-dimensional
N1 Nd X
N1 Nd
1 X X X
IξN (g)(h) = ··· ··· XN (h)
N (2π d ) y =1 y =1 z =1 z =1 y1 ,...,yd
1 d 1 d
d
!
X
×XzN1 ,...,zd (g) exp −i ξj (yj − zj ) , (4.18)
j=1
141
Capítulo 4. Procesos de Cox espaciales en un marco infinito-dimensional
n n
!
(n) 1 XX
ρϕ,θb (z1 , . . . , zn ) = [ρϕ,θbN ]n exp R b (ϕ)(ϕ) ,
N 2 i=1 j=1 zi −zj ,θN
para z1 , . . . , zn ∈ Zd , y ϕ ∈ H, con
!
(1) R0,θbN (ϕ)(ϕ)
ρϕ,θbN = ρd
ϕ,θ = ρϕ,θ (z) = exp , ∀z ∈ Zd . (4.22)
d
2
Demostración:
Consideramos aquí una versión compacta de la demostración. Algunos detalles
técnicos adicionales se pueden encontrar en la Nota 4.8. En concreto, la de-
mostración de este resultado se desprende de la hipótesis 4.2 que implica, para
cualquier base ortonormal, {ϕk , k ≥ 1} de H + iH, y para z = (z1 , . . . , zd ) ∈ Zd
con zi ∈ [−Ni + 1, Ni − 1], i = 1, . . . , d,
X
|C(z, k, l) − Rz,θ0 (ϕk )(ϕl )|2 →c. s. 0, N → ∞, (4.24)
k,l≥1
142
Capítulo 4. Procesos de Cox espaciales en un marco infinito-dimensional
NXd −zd
"N −z #
1 1
1 X
= ··· Xy (ϕk )Xy+z (ϕl ) , |zi | < Ni , i ∈ {1, . . . , d},
N y =1 y =1
1 d
d
!
M 1 X X
qξ,θ (ϕk )(ϕk ) = exp −i ξj zj
(2π)d j=1
zj ∈[−Mj +1,Mj −1]; j=1,...,d
d
Y zj
× 1− g(z1 , . . . , zd , θ, k)
j=1
Mj
Z
−1
= FM (ξ − $)F$,θ (ϕk )(ϕk )d$, (4.26)
[−π,π]d
−1
a Fξ,θ (ϕk )(ϕk ), uniformemente en k ≥ 1, ξ ∈ [−π, π]d , y θ ∈ Θ. Aquí, para
cada z = (z1 , . . . , zd ) ∈ dj=1 [−Mj + 1, Mj − 1], g(z, θ, k) = g(z1 , . . . , zd , θ, k)
Q
−1
denota el z–ésimo coeficiente de Fourier de F$,θ (ϕk )(ϕk ), véase la ecuación
(4.34) posterior.
Para un M suficientemente grande, aplicando la desigualdad de Hölder, a
partir de la Suposición 4.2 y las condiciones (4.3)–(4.6) se obtiene la conver-
gencia c. s. a cero, cuando N → ∞, para cualquier θ ∈ Θ, de
Z h i
M (N )
lı́m sup qξ,θ Iξ (ϕk )(ϕk ) − Fξ,θ0 (ϕk )(ϕk ) dξ
N →∞ k≥1 [−π,π]d
Z h i2
M (N )
≤c.s. lı́m sup qξ,θ Iξ (ϕk )(ϕk ) − Fξ,θ0 (ϕk )(ϕk ) dξ
N →∞ k≥1 [−π,π]d
(2π)d X M (N ) 2
= lı́m sup qξz ,θ Iξz (ϕk )(ϕk ) − qξMz ,θ Fξz ,θ0 (ϕk )(ϕk ) = 0,
N →∞ k≥1 N z
c.s.
(4.27)
143
Capítulo 4. Procesos de Cox espaciales en un marco infinito-dimensional
= S1 (N , M ) + S2 (N , M ) + S3 (M ), θ ∈ Θ. (4.28)
−1
La sumabilidad de Cesaro de la serie de Fourier de Fξ,θ en norma L(H + iH),
d
uniformemente en (ξ, θ) ∈ [−π, π] ×Θ en la Condición 4.6 implica que, cuando
M → ∞, ε(M ) → 0. Finalmente, bajo la Suposición 4.2, supk≥1 C(0, k, k) →
144
Capítulo 4. Procesos de Cox espaciales en un marco infinito-dimensional
S2 (N , M ) →c.s. 0 N → ∞
−1
ya que, según se deduce de la Nota 4.7, para cada k ≥ 1, Fξ,θ0 Fξ,θ (ϕk )(ϕk )
define el espectro de un proceso lineal estacionario cuya varianza de predicción
en un paso es la unidad. Específicamente, la varianza de ese proceso viene
−1
R
dada por [−π,π]d Fξ,θ0 Fξ,θ (ϕk )(ϕk )dξ, que debe ser mayor que la varianza de
predicción en un paso, a menos que el espectro sea constante para cada k ≥ 1.
−1
Este hecho no se cumple porque significa que la composición Fξ,θ0 Fξ,θ de los
−1
operadores Fξ,θ0 y Fξ,θ coincide con el operador identidad en H + iH, que no
es posible en vista de la condición de identificabilidad, 4.3, ya que θ 6= θ0 en
(4.31). Por lo tanto,
Z
1 −1
ı́nf sup Fξ,θ0 Fξ,θ (ϕk )(ϕk )dξ = 1. (4.32)
θ∈Θ k≥1 (2π)d [−π,π]d
e N (ϕk )X
IξN (ϕk )(ϕl ) = X ξ
e N (ϕl )
−ξ
d
!
1 X X
= C(z, k, l) exp −i ξj zj , k, l ≥ 1,
(2π)d Qd j=1
z∈ i=1 [−Ni +1,Ni −1]
−1
de la inversa F$,θ del operador de densidad espectral. Este truncamiento se ex-
presa de forma equivalente en la ecuación (4.26) en términos del kernel de Féjer
multidimensional,
M1 Md X
M1 Md d
!
1 X X X X
FM (ξ) = Qd ··· ··· exp −i (yj − xj )ξj . (4.35)
j=1 Mj y1 =1 yd =1 x1 =1 xd =1 j=1
(2π)d X M (N )
2
qξzj ,θ Iξz (ϕk )(ϕk ) − qξMzj ,θ Fξzj ,θ0 (ϕk )(ϕk )
N z j
j
d 2
(2π)d
X Y |uj |
= 1− g(u, θ, k)
N Mj
uj ∈[−Mj +1,Mj −1], j=1,...,d j=1
146
Capítulo 4. Procesos de Cox espaciales en un marco infinito-dimensional
147
Capítulo 4. Procesos de Cox espaciales en un marco infinito-dimensional
148
Capítulo 4. Procesos de Cox espaciales en un marco infinito-dimensional
donde ΦMi (θ) denota al operador de proyección sobre el subespacio generado por
los autovectores φ1 , . . . , φMi (θ) de (−∆)T (θ) , siendo Mi (θ) ≥ 1 una función acota-
da en θ ∈ Θ, para i = 1, 2. Aquí, βi , Ki ∈ R+ son conocidos, y Ki es continua
en θ ∈ Θ, para i = 1, 2. Aplicando el Teorema espectral sobre el cálculo espec-
tral para operadores autoadjuntos en un espacio separable de Hilbert, véase, por
ejemplo, Dautray y Lions (1985), los autovalores {λk,i (θ), k ≥ 1} del operador de
autocorrelación Li,θ vienen dados, para i = 1, 2, 3, por
−β1
λk,1 (θ) = 1[1,M1 (θ)] (k)K1 (θ, π) λk (−∆)T (θ)
−β2
λk,2 (θ) = 1[1,M2 (θ)] (k)K2 (θ, π) λk (−∆)T (θ)
λk,3 (θ) = −λk,1 (θ)λk,2 (θ), k ≥ 1, θ ∈ Θ, (4.40)
donde 1[1,Mi (θ)] (k) denota la función indicadora del conjunto [1, Mi (θ)] ⊂ N, con,
como antes, Mi (θ) ≥ 1, para θ ∈ Θ, y i = 1, 2. Para k ≥ 1, λk (−∆)T (θ) denota el
k–ésimo autovalor de (−∆)T (θ) . Consideramos el caso particular M1 (θ) = M2 (θ) =
M (θ), y asumimos que λk,i (θ), k ≥ 1, i = 1, 2, 3, son tales que las raíces del
polinomio Pk,θ (z1 , z2 ) = 1i=0 1j=0 pk,i,j,θ z1i z2j , están fuera del polidisco unitario
P P
(es decir, el producto de dos discos unitarios abiertos en el plano complejo). Aquí
pk,0,0,θ = 1, pk,1,0,θ = λk,1 (θ), pk,0,1,θ = λk,2 (θ), y pk,1,1,θ = −λk,1 (θ)λk,2 (θ), k ≥
1. A partir de las ecuaciones (4.37)–(4.40), calculamos el operador de densidad
149
Capítulo 4. Procesos de Cox espaciales en un marco infinito-dimensional
espectral
Bajo las condiciones asumidas sobre los parámetros del modelo (4.40), la Su-
posición 4.2 se satisface en nuestro marco SARH(1), como se deduce del Corolario
2.3 en (Bosq, 2000), aplicando resultados similares a los obtenidos en el Corola-
rio 4.1 y el Teorema 4.8 de dicha monografía. La condición de identificabilidad
4.3 se deduce de la caracterización de los autovectores y autovalores del operador
laplaciano negativo de Dirichlet a partir del soporte T (θ), definido unívocamente
por el parámetro θ. Las condiciones (4.15), 4.4 y 4.6 se obtienen a partir de las
ecuaciones (4.40)–(4.41), bajo las restricciones consideradas en los coeficientes
paramétricos {λk,i (θ), k ∈ [1, M (θ)], θ ∈ Θ, i = 1, 2} definiendo la familia de
polinomios {Pk,θ , k ∈ [1, M (θ)], θ ∈ Θ} . La Condición 4.5 también se cumple ya
que σk2 = 1, para cada k ≥ 1.
La consistencia fuerte del predictor paramétrico plug-in SARH(1)
N
κ
b i,j (θ) = L1,θbN κi−1,j + L2,θbN κi,j−1 + L3,θbN κi−1,j−1 , ∀(i, j) ∈ Z2 , (4.42)
se obtiene entonces bajo condiciones similares a las asumidas en Bosq (2000) para
el caso ARH(1).
Los resultados de la estimación se muestran en la Tabla 4.1 y en la Figura 4.1,
2−11/10 2−12/10
donde se ha considerado θ0 = 1, K1 (θ, π) = πθ 2−11/10 , K2 (θ, π) = πθ 2−12/10 , β1 = 11
20
,
12
β2 = 20 , M1 (θ) = M2 (θ) = [10(θ)], y θ ∈ Θ = [0.7, 4], con [·] denotando la fun-
ción parte entera. Se han probado los tamaños de muestra funcional espacial N =
40000, 62500, 90000, 122500, 160000, 202500, 250000, 302500. Específica-
mente, la Figura 4.1, en la parte izquierda, muestra los valores muestrales de θbN ,
basados en 100 generaciones de cada muestra funcional. Los errores cuadráticos
medios empíricos también se han calculado en la Tabla 4.1 y la Figura 4.1, en la
parte derecha.
150
Capítulo 4. Procesos de Cox espaciales en un marco infinito-dimensional
1.4 0.14
M.S.E.
S.D.
1.3 0.12
1.2 0.1
1.1
0.08
1
0.06
0.9
0.04
0.8
0.02
0.7
0
200 250 300 350 400 450 500 550 200 250 300 350 400 450 500 550
Figura 4.1: Diagramas de cajas de los valores muestrales de θbN , para θ0 = 1 basados
en 100 generaciones (a la izquierda). A la derecha, el error cuadrático medio empírico
(M.S.E), línea azul, y desviación típica (S.D.), línea roja continua
151
Capítulo 4. Procesos de Cox espaciales en un marco infinito-dimensional
donde, como antes, 1[1,M ] (k) es la función indicadora del conjunto [1, M ] ⊂ N, ΦM ,
M ≥ 1, denota el operador de proyección sobre el subespacio generado por los au-
tovectores φ1 , . . . , φM del operador laplaciano negativo de Dirichlet (−∆)[0,1] en el
intervalo [0, 1]. Aquí, {λk,1 (θ1 ), λk,2 (θ2 ), k ≥ 1, (θ1 , θ2 ) ∈ Θ} es la familia para-
métrica de autovalores de los operadores de autocorrelación
{L1,θ1 , L2,θ2 , (θ1 , θ2 ) ∈ Θ}.
Bajo condiciones similares a las del Ejemplo 4.1 sobre estas familias paramé-
tricas de autovalores de {L1,θ1 , L2,θ2 , (θ1 , θ2 ) ∈ Θ}, se cumplen las suposiciones
4.1–4.2, y las condiciones (4.15) y 4.3–4.6. Por lo tanto, la consistencia fuerte del
estimador paramétrico θbN se deduce del Teorema 4.1. En particular, se han con-
siderado los valores de los parámetros de escala θ1,1 = 1, θ2,1 = 1.5, y los valores
de los parámetros de localización θ1,2 = 1.6, y θ2,2 = 1.2, que se encuentran en el
interior de Θ = ([0.7, 1.3] × [1.3, 1.9]) × ([1.2, 1.8] × [0.9, 1.5]).
Véanse los resultados numéricos mostrados en los gráficos ubicados en la par-
te derecha de las figuras 4.2–4.5, a partir de las 100 generaciones de cada una
de las muestras funcionales espaciales de tamaños N = 40000, 62500, 90000,
122500, 160000, 202500, 250000, 302500. El error cuadrático medio empírico
también se ha calculado para estos tamaños de muestra (ver los gráficos ubicados
en la parte derecha de las figuras 4.2–4.5 y la Tabla 4.2). Se observa una buena
aproximación a partir de la técnica de estimación propuesta, para los tamaños de
muestra funcional espacial considerados, y para el valor del parámetro de trunca-
miento elegido, M = 10. Obsérvese que este valor del parámetro de truncamiento
proporciona la dimensión umbral para asegurar un comportamiento estable del
152
Capítulo 4. Procesos de Cox espaciales en un marco infinito-dimensional
1.1 0.035
M.S.E.
S.D.
0.03
1.05
0.025
1 0.02
0.015
0.95
0.01
0.9
0.005
0
200 250 300 350 400 450 500 550 200 250 300 350 400 450 500 550
Figura 4.2: Diagramas de cajas de los valores muestrales de θb1,1,N , para θ1,1,0 = 1,
basados en 100 generaciones (a la izquierda). A la derecha, el error cuadrático medio
empírico (M.S.E), línea azul, y desviación típica (S.D.), línea roja continua
1.6 0.025
M.S.E.
1.58 S.D.
1.56 0.02
1.54
1.52 0.015
1.5
1.48 0.01
1.46
1.44 0.005
1.42
1.4 0
200 250 300 350 400 450 500 550 200 250 300 350 400 450 500 550
Figura 4.3: Diagramas de cajas de los valores muestrales de θb2,1,N , para θ2,1,0 = 1.5,
basados en 100 generaciones (a la izquierda). A la derecha, el error cuadrático medio
empírico (M.S.E), línea azul, y desviación típica (S.D.), línea roja discontinua
153
Capítulo 4. Procesos de Cox espaciales en un marco infinito-dimensional
0.025
1.66 M.S.E.
S.D.
1.64
0.02
1.62
1.6
0.015
1.58
1.56
1.54 0.01
1.52
1.5 0.005
1.48
1.46
0
200 250 300 350 400 450 500 550 200 250 300 350 400 450 500 550
Figura 4.4: Diagramas de cajas de los valores muestrales de θb1,2,N , para θ1,2,0 = 1.6,
basados en 100 generaciones (a la izquierda). A la derecha, el error cuadrático medio
empírico (M.S.E), línea azul, y desviación típica (S.D.), línea roja continua
1.5 0.035
M.S.E.
1.45 S.D.
0.03
1.4
0.025
1.35
0.02
1.3
1.25 0.015
1.2
0.01
1.15
0.005
1.1
0
200 250 300 350 400 450 500 550 200 250 300 350 400 450 500 550
Figura 4.5: Diagramas de cajas de los valores muestrales de θb2,2,N , para θ2,2,0 = 1.2,
basados en 100 generaciones (a la izquierda). A la derecha, el error cuadrático medio
empírico (M.S.E), línea azul, y desviación típica (S.D.), línea roja discontinua
154
Capítulo 4. Procesos de Cox espaciales en un marco infinito-dimensional
156
Capítulo 4. Procesos de Cox espaciales en un marco infinito-dimensional
157
Capítulo 4. Procesos de Cox espaciales en un marco infinito-dimensional
n o
Las correspondientes estimaciones paramétricas λp,i (θbi,N ), p = 1, . . . , 10 ,
i = 1, 2, 3, de los espectros puntuales de Li,θbi,N , i = 1, 2, 3, se muestran en la
Tabla 4.3 y en la Figura 4.8. Finalmente, las curvas originales de log-intensidad
de enfermedades respiratorias (línea azul discontinua), y sus estimaciones (lí-
nea roja), en las 48 provincias españolas, se muestran en las figuras 4.9–4.10.
El peor y mejor ajuste del modelo (provincias de Soria y Pontevedra, respec-
tivamente) se amplían en la Figura 4.11. Asimismo, los mapas de riesgo de
mortalidad por enfermedades respiratorias originales y estimados, promedia-
dos anualmente, se muestran en la Figura 4.12.
158
Capítulo 4. Procesos de Cox espaciales en un marco infinito-dimensional
10 10 10 10 10
15 15 15 15 15
5 5 5 5 5
10 0 10 0 10 0 10 0 10 0
-5 -5 -5 -5 -5
5 5 5 5 5
-10 -10 -10 -10 -10
10 10 10 10 10
15 15 15 15 15
5 5 5 5 5
10 0 10 0 10 0 10 0 10 0
-5 -5 -5 -5 -5
5 5 5 5 5
-10 -10 -10 -10 -10
10 10 10 10 10
15 15 15 15 15
5 5 5 5 5
10 0 10 0 10 0 10 0 10 0
-5 -5 -5 -5 -5
5 5 5 5 5
-10 -10 -10 -10 -10
10 10 10 10 10
15 15 15 15 15
5 5 5 5 5
10 0 10 0 10 0 10 0 10 0
-5 -5 -5 -5 -5
5 5 5 5 5
-10 -10 -10 -10 -10
159
Capítulo 4. Procesos de Cox espaciales en un marco infinito-dimensional
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
-0.1
1 2 3 4 5 6 7 8 9 10
160
Capítulo 4. Procesos de Cox espaciales en un marco infinito-dimensional
1.5 2.6
3 2
2 2 2.4
1
0 1000 2000 0 1000 2000 0 1000 2000 0 1000 2000 0 1000 2000 0 1000 2000
Baleares Barcelona Burgos Cáceres Cádiz Castellón
3 4.5 2.5 2.5
3
2
2.5 4
2
2 3.5 1.5 2.5 2
0 1000 2000 0 1000 2000 0 1000 2000 0 1000 2000 0 1000 2000 0 1000 2000
Ciudad Real Córdoba Coruña (A) Cuenca Girona Granada
3 3.5 3
3 2.5
2.5
2.5 3 2
2.5 2 2.5
0 1000 2000 0 1000 2000 0 1000 2000 0 1000 2000 0 1000 2000 0 1000 2000
Guadalajara Guipúzcoa Huelva Huesca Jaén León
2 2.7 2.8
2.6 1.5 2.8
2.5 2.2
1.5 2.7 2.6
2.4
2.3 1 2.6
1 2 2.4
0 1000 2000 0 1000 2000 0 1000 2000 0 1000 2000 0 1000 2000 0 1000 2000
3.5
2 1.5 2.5 1.5 2.6
0 1000 2000 0 1000 2000 0 1000 2000 0 1000 2000 0 1000 2000 0 1000 2000
161
Capítulo 4. Procesos de Cox espaciales en un marco infinito-dimensional
Soria Pontevedra
1.3 3.1
1.2
3
1.1
2.9
1
2.8
0.9
2.7
0.8
2.6
0.7
2.5
0.6
0.5 2.4
0 500 1000 1500 2000 0 500 1000 1500 2000
Figura 4.11: Zoom. Provincias españolas con mayor y menor norma L1 ([0, 1725]) del
error relativo absoluto medio funcional asociado (Soria a la izquierda, Pontevedra a la
derecha, respectivamente). Las curvas originales aparecen en línea azul discontinua, y
las curvas estimadas se representan en línea roja continua
Tabla 4.4: ACVFARE. Promedio puntual anual del CVFARE obtenido a partir de la
validación
162
Capítulo 4. Procesos de Cox espaciales en un marco infinito-dimensional
163
Capítulo 5
165
Capítulo 5. Un modelo de recuento funcional espacial para el análisis de la
heterogeneidad en el tiempo
RX
z−y (f )(g) = E (Xz ⊗ Xy ) (f )(g) = hE (Xz ⊗ Xy ) (f ), giH , f, g ∈ H,
(5.1)
166
Capítulo 5. Un modelo de recuento funcional espacial para el análisis de la
heterogeneidad en el tiempo
satisfacen
∞ 1/2
X ?
kRX RX RX
z−y k`1 (H) = z−y z−y (ϕj ), ϕj < ∞,
j=1 H
A partir de (5.1), RX X
0 , con kernel r0 , es un operador traza autoadjunto (si-
métrico), que satisface
RX X
0 (φj ) = λj (R0 )φj , j ≥ 1, (5.2)
Es decir,
M 2
X
E Xz − hXz , φj iH φj → 0, M → ∞,
j=1 H
con E hXz , φj iH hXz , φp iH = δj,p λj (RX d
0 ), j ≥ 1, para cada z ∈ R . Aquí, δj,p
denota la función delta de Kronecker.
Supongamos que X es tal que, para cada z ∈ Rd , los valores funcionales de
Xz casi seguramente (c. s.) tienen su soporte en el intervalo temporal acotado
167
Capítulo 5. Un modelo de recuento funcional espacial para el análisis de la
heterogeneidad en el tiempo
R
Dadas las observaciones Ψz,ω0 = T Λz,ω0 (t)dt, z ∈ A ⊂ Rd , para cierto
ω0 ∈ Ω, el número de sucesos C(A), que ocurren durante el periodo T , en la
región A, sigue una distribución de probabilidad de Poisson de media Λ(A).
Nótese que, el predictor de mínimos cuadrados de C(A) viene dado por Λ(A),
introducido en (5.5), para cualquier conjunto de Borel acotado A ∈ B d . De
(5.4), la ecuación (5.3) lleva a la siguiente expresión de la variación de segundo
orden de Ψz :
∞ Z
"∞ #p
X Cp (t)Cp (s) X
E[Ψ2z ] = λj (RX
0 )φj ⊗ φj (t, s) dtds. (5.6)
p=0 T ×T p! j=1
168
Capítulo 5. Un modelo de recuento funcional espacial para el análisis de la
heterogeneidad en el tiempo
Demostración:
De la ecuación (5.6), aplicando la Proposición 4.9 en la p. 92 de Marinucci
y Peccati (2011), después de considerar la desigualdad de Cauchy-Schwarz,
en términos del producto interno introducido en la Fórmula (4.7) en la p.89
de Marinucci y Peccati (2011), las propiedades de expansión de Hermite nos
llevan a
∞ Z
"∞ #p
X C p (t)C p (s) X
E[Ψ2z ] = λj (RX0 )φj ⊗ φj (t, s) dtds
p=0 T ×T p! j=1
∞ Z
1
X p q
≤ E[Λz (t)]2 E[Λz (s)]2 E[Hp (Xz (t))]2 E[Hp (Xz (s))]2
p=0
p! T ×T
"∞ ∞
#p/2
X X
× λj (RX 0 )φj ⊗ φj (t, t) λj (RX0 )φj ⊗ φj (s, s) dtds
j=1 j=1
∞ Z
X 1
= exp (r0 (t, t)/2 + r0 (s, s)/2) [r0 (t, t)r0 (s, s)]p/2
p=0
p! T ×T
"∞ ∞
#p/2
X X
× λj (RX
0 )φj ⊗ φj (t, t) λj (RX0 )φj ⊗ φj (s, s) dtds
j=1 j=1
∞
(Z " ∞
#p )2
X 1 X
= exp (r0 (t, t)/2) λj (RX
0 )φj ⊗ φj (t, t) dt
p=0
p! T j=1
(∞ "∞ #p )2
X M2p X
≤ |T |2 exp 2M2 kRX
0 k`1 (H) λj (RX
0)
p=0
p! j=1
2
4kRX 2
= |T | exp 0 k`1 (H) M < ∞,
169
Capítulo 5. Un modelo de recuento funcional espacial para el análisis de la
heterogeneidad en el tiempo
× ψqy−z
1
(ψj:k ) · · · ψqy−z
p3
(ψj:k )ϕy−z y−z
q1 (ψj:k ) · · · ϕp3 (ψj:k )
"∞ #p1 " ∞ #p2 " ∞ #p3
X (1/2)p2 +p3 X X X
λh RX λl RX λq RX
≤ 0 z−y y−z
p1 ,p2 ,p3
p 1 !p 2 !p 3 ! h=1 l=1 q=1
1 X
= exp kR0 k`1 (H) + kRz−y k`1 (H) + kRy−z k`1 (H) < ∞, ∀z, y ∈ Rd ,
X X
2
(5.7)
170
Capítulo 5. Un modelo de recuento funcional espacial para el análisis de la
heterogeneidad en el tiempo
171
Capítulo 5. Un modelo de recuento funcional espacial para el análisis de la
heterogeneidad en el tiempo
1 X
Fω1 ,ω2 := RX
p,q exp (−i(pω1 + qω2 )) , (ω1 , ω2 ) ∈ [0, 2π) × [0, 2π),
(2π)2 2
(p,q)∈Z
(5.11)
que es un operador traza autoadjunto no negativo.
172
Capítulo 5. Un modelo de recuento funcional espacial para el análisis de la
heterogeneidad en el tiempo
o, equivalentemente por
eN ⊗ X
IωN1 ,ω2 := X ω1 ,ω2
eN .
ω1 ,ω2 (5.14)
Para una base ortonormal dada de wavelets con soporte compacto {ψj:k , k ∈
173
Capítulo 5. Un modelo de recuento funcional espacial para el análisis de la
heterogeneidad en el tiempo
S1 X
S2
eωN ,ω (ψj:k ) = 1 X
X 1 2
√ Xp,q (ψj:k ) exp (−i(pω1 + qω2 )) , k ∈ Γj , j ∈ Z
2π N p=1 q=1
(5.15)
eωN ,ω (ψj1 :k1 )X
IωN1 ,ω2 (ψj1 :k1 )(ψj2 :k2 ) = X 1 2
eωN ,ω (ψj2 :k2 )
1 2
S1 X
S2 X
S1 X
S2
X Xp,q (ψj 1 :k1 )Xp0 ,q0 (ψj2 :k2 ) exp (−i[(p − p0 )ω1 + (q − q 0 )ω2 ])
= ,
p=1 q=1 p0 =1 q 0 =1
(2π)2 N
(5.16)
y la aproximación multiresolución
{Xp,q (ψj:k ), p = 1, . . . , S1 , q = 1, . . . , S2 , k ∈ Γj , j ∈ Z}
174
Capítulo 5. Un modelo de recuento funcional espacial para el análisis de la
heterogeneidad en el tiempo
donde θ0,j:k denota el verdadero valor del parámetro, asociado con el nodo k a
escala j ∈ Γj . La densidad espectral espacial normalizada multiescalar
para cada $ = ($1 , $2 ) ∈ [0, 2π) × [0, 2π), siendo ηj:k ($) una función espacial
simétrica no negativa, para cada k ∈ Γj , y j ∈ Z, tal que ηj:k ($)fj:k ($, θj:k ) ∈
L1 ([0, 2π) × [0, 2π)) , el espacio de funciones integrables absolutas en [0, 2π) ×
[0, 2π), para cada θj:k ∈ Θj:k ⊂ Θ.
175
Capítulo 5. Un modelo de recuento funcional espacial para el análisis de la
heterogeneidad en el tiempo
N
donde IN,j:k ($) = I$ 1 ,$2
(ψj:k )(ψj:k ) denota, como antes, el periodograma mul-
tiescalar presentado en (5.16), para k ∈ Γj , y j ∈ Z.
para cada θj:k ∈ Θj:k ⊂ Θ, donde P0,j:k denota la medida asociada con función
de densidad fj:k ($, θ0,j:k ), para cada k ∈ Γj , y j ∈ Z. Para minimizar la di-
vergencia en (5.18), en la práctica, podemos calcular el mínimo de ÛN,j:k (θj:k )
sobre θj:k ∈ Θj:k , a través de los diferentes nodos k en cada escala j ∈ Z. Es
decir, consideraremos los estimadores paramétricos multiescalares
176
Capítulo 5. Un modelo de recuento funcional espacial para el análisis de la
heterogeneidad en el tiempo
parámetros
θj0 :ek,el = θj0 :ek,el,1 , θj0 :ek,el,2 , θj0 :ek,el,3
= L1 (ϕj0 :el )(ϕj0 :ek ), L2 (ϕj0 :el )(ϕj0 :ek ), L3 (ϕj0 :el )(ϕj0 :ek ) ∈ Θj0 :ek × Θj0 :el ,
θj0 :ek;j:k = θj0 :ek;j:k,1 , θj0 :ek;j:k,2 , θj0 :ek;j:k,3
= L1 (ψj:k )(ϕj0 :ek ), L2 (ψj:k )(ϕj0 :ek ), L3 (ψj:k )(ϕj0 :ek ) ∈ Θj0 :ek × Θj:k ,
θj:k;j0 :ek = θj:k;j0 :ek,1 , θj:k;j0 :ek,2 , θj:k;j0 :ek,3
= L1 (ϕj0 :ek )(ψj:k ), L2 (ϕj0 :ek )(ψj:k ), L3 (ϕj0 :ek )(ψj:k ) ∈ Θj:k × Θj0 :ek ,
177
Capítulo 5. Un modelo de recuento funcional espacial para el análisis de la
heterogeneidad en el tiempo
X X
X
bN,p,q (·) = θbN,j0 :ek,el,1 Xp−1,q (ϕj0 :el )ϕj0 :ek (·)
k∈Υj0 e
e l∈Υj0
X XX
+ θbN,j0 :ek;j:k,1 Xp−1,q (ψj:k )ϕj0 :ek (·)
k∈Υj0 j≥j0 k∈Γj
e
XX X
+ θbN,j:k;j0 :ek,1 Xp−1,q (ϕj0 :ek )ψj:k (·)
j≥j0 k∈Γj e
k∈Υj0
XXXX
+ θbN,j:k;l:h,1 Xp−1,q (ψl:h )ψj:k (·)
j≥j0 k∈Γj l≥j0 h∈Γl
X X
+ θbN,j0 :ek,el,2 Xp,q−1 (ϕj0 :el )ϕj0 :ek (·)
k∈Υj0 e
e l∈Υj0
X XX
+ θbN,j0 :ek;j:k,2 Xp,q−1 (ψj:k )ϕj0 :ek (·)
k∈Υj0 j≥j0 k∈Γj
e
XX X
+ θbN,j:k;j0 :ek,2 Xp,q−1 (ϕj0 :ek )ψj:k (·)
j≥j0 k∈Γj e
k∈Υj0
XXXX
+ θbN,j:k;l:h,2 Xp,q−1 (ψl:h )ψj:k (·)
j≥j0 k∈Γj l≥j0 h∈Γl
X X
+ θbN,j0 :ek,el,3 Xp−1,q−1 (ϕj0 :el )ϕj0 :ek (·)
k∈Υj0 e
e l∈Υj0
X XX
+ θbN,j0 :ek;j:k,3 Xp−1,q−1 (ψj:k )ϕj0 :ek (·)
k∈Υj0 j≥j0 k∈Γj
e
XX X
+ θbN,j:k;j0 :ek,3 Xp−1,q−1 (ϕj0 :ek )ψj:k (·)
j≥j0 k∈Γj e
k∈Υj0
XXXX
+ θbN,j:k;l:h,3 Xp−1,q−1 (ψl:h )ψj:k (·).
j≥j0 k∈Γj l≥j0 h∈Γl
(5.25)
178
Capítulo 5. Un modelo de recuento funcional espacial para el análisis de la
heterogeneidad en el tiempo
179
Capítulo 5. Un modelo de recuento funcional espacial para el análisis de la
heterogeneidad en el tiempo
180
Capítulo 5. Un modelo de recuento funcional espacial para el análisis de la
heterogeneidad en el tiempo
Node 2,2 Node 2,9 Node 2,16 Node 2,23 Node 2,30
5 4 5
2 4
0 2 2
0 0 0
-2 0
-4 -2
-5 -2 -5 -4
0 500 1000 0 500 1000 0 500 1000 0 500 1000 0 500 1000
Node 9,2 Node 9,9 Node 9,16 Node 9,23 Node 9,30
4 5
4 2 4
2
0 2 0 2
0 0 -2 0
-2
-4 -2 -4 -2
-6 -5 -4 -6 -4
0 500 1000 0 500 1000 0 500 1000 0 500 1000 0 500 1000
Node 16,2 Node 16,9 Node 16,16 Node 16,23 Node 16,30
5
5 5 5 5
0 0 0 0 0
-5 -5 -5 -5
-5
0 500 1000 0 500 1000 0 500 1000 0 500 1000 0 500 1000
Node 23,2 Node 23,9 Node 23,16 Node 23,23 Node 23,30
4 6
5 5 4
2 4
2 2
0
0 0 0 0
-2
-2 -2
-5 -4 -4
-5 -4
0 500 1000 0 500 1000 0 500 1000 0 500 1000 0 500 1000
Node 30,2 Node 30,9 Node 30,16 Node 30,23 Node 30,30
5 5 5
4 4
2 2 0
0 0 0 0
-2 -2 -5
-4 -4 -5
-10 -5
0 500 1000 0 500 1000 0 500 1000 0 500 1000 0 500 1000
Figura 5.1: Escala 10, N = 900. Datos curva sobre algunos nodos de una rejilla regular
espacial 30 × 30
Node 2,2 Node 2,9 Node 2,16 Node 2,23 Node 2,30
10 6 10
0 4 5
0 2 0
-5 0
0
-10 -10 -2 -10 -5
0 100 200 0 100 200 0 100 200 0 100 200 0 100 200
Node 9,2 Node 9,9 Node 9,16 Node 9,23 Node 9,30
5 5 10 5
5
0 0
0 0 0
-5 -5
-5
-10 -10 -10 -5
0 100 200 0 100 200 0 100 200 0 100 200 0 100 200
Node 16,2 Node 16,9 Node 16,16 Node 16,23 Node 16,30
10 10 10 10
5
5
0 0 0 0 0
-5
-5
-10 -10 -10
0 100 200 0 100 200 0 100 200 0 100 200 0 100 200
Node 23,2 Node 23,9 Node 23,16 Node 23,23 Node 23,30
10 10 10
5
10
5 0 5
0 0
0 0
-5
-5 -5
-10 -10 -10
0 100 200 0 100 200 0 100 200 0 100 200 0 100 200
Node 30,2 Node 30,9 Node 30,16 Node 30,23 Node 30,30
10 10 10 4
5
2
0 0
0 0 0 -2
-5
-4
-10 -6
-10 -10 -10
0 100 200 0 100 200 0 100 200 0 100 200 0 100 200
Figura 5.2: Escala 9, N = 900. Datos curva sobre algunos nodos de una rejilla regular
espacial 30 × 30
181
Capítulo 5. Un modelo de recuento funcional espacial para el análisis de la
heterogeneidad en el tiempo
Node 2,2 Node 2,9 Node 2,16 Node 2,23 Node 2,30
5 10 10 10
10
5 5
0 0
0 0
-5 0
-10 -5
-10 -10
0 50 100 0 50 100 0 50 100 0 50 100 0 50 100
Node 9,2 Node 9,9 Node 9,16 Node 9,23 Node 9,30
5 5 10 5 10
0
0 0
0 -5 0
-5 -5
-10
-10 -10 -10 -10
0 50 100 0 50 100 0 50 100 0 50 100 0 50 100
Node 16,2 Node 16,9 Node 16,16 Node 16,23 Node 16,30
10 10 10
10 10
5
5
0 0 0 0 0
-5 -5
-10 -10
-10
0 50 100 0 50 100 0 50 100 0 50 100 0 50 100
Node 23,2 Node 23,9 Node 23,16 Node 23,23 Node 23,30
10
5 10
5 10
0 5
0 0
0 0
-5 -5
-5
-10 -10 -10 -10
0 50 100 0 50 100 0 50 100 0 50 100 0 50 100
Node 30,2 Node 30,9 Node 30,16 Node 30,23 Node 30,30
10 10 10 5
5 0 0
0 0 0
-10 -5
-5
-10 -20 -10 -10
0 50 100 0 50 100 0 50 100 0 50 100 0 50 100
Figura 5.3: Escala 8, N = 900. Datos curva sobre algunos nodos de una rejilla regular
espacial 30 × 30
Node 2,2 Node 2,9 Node 2,16 Node 2,23 Node 2,30
20 10 20 20
5
0 5 10
0 0
-5 0
0
-10
-20 -20
0 20 40 60 0 20 40 60 0 20 40 60 0 20 40 60 0 20 40 60
Node 9,2 Node 9,9 Node 9,16 Node 9,23 Node 9,30
5 5 10
0 10
0 0 5
-5 0 0
-5 -10
-10 -5
-10 -10 -20
0 20 40 60 0 20 40 60 0 20 40 60 0 20 40 60 0 20 40 60
Node 16,2 Node 16,9 Node 16,16 Node 16,23 Node 16,30
20 20 10 20
10
0 0 0 0 0
Figura 5.4: Escala 7, N = 900. Datos curva sobre algunos nodos de una rejilla regular
espacial 30 × 30
182
Capítulo 5. Un modelo de recuento funcional espacial para el análisis de la
heterogeneidad en el tiempo
Tabla 5.2: Promedio por escala de los errores cuadráticos medios empíricos asocia-
dos a los estimadores paramétricos implicados en la aproximación de L1 , a partir de
100 generaciones de las muestras funcionales de tamaño N = 100, 900, 2500, 4900,
8100, 12100, 16900, 22500. Se consideran las escalas j = 6, 7, 8, 9, y el parámetro de
truncamiento kN = 10
Tabla 5.3: Promedio por escala de los errores cuadráticos medios empíricos asocia-
dos a los estimadores paramétricos implicados en la aproximación de L2 , a partir de
100 generaciones de las muestras funcionales de tamaño N = 100, 900, 2500, 4900,
8100, 12100, 16900, 22500. Se consideran las escalas j = 6, 7, 8, 9, y el parámetro de
truncamiento kN = 10
183
Capítulo 5. Un modelo de recuento funcional espacial para el análisis de la
heterogeneidad en el tiempo
0.03
p= 1
p= 2
0.025 p= 3
p= 4
p= 5
p= 6
0.02
p= 7
p= 8
p= 9
0.015 p= 10
0.01
0.005
0
10 30 50 70 90 110 130 150
S
0.04
p= 1
0.035 p= 2
p= 3
p= 4
0.03 p= 5
p= 6
p= 7
0.025
p= 8
p= 9
0.02 p= 10
0.015
0.01
0.005
0
10 30 50 70 90 110 130 150
S
Figura
n 5.5: Errores cuadráticos
o medios empíricos (EMSE), asociados con
λ bN,p,2 , p = 1, . . . , kN , N = 100, 900, 2500, 4900, 8100, 12100, 16900,
bN,p,1 , λ
22500
184
Capítulo 5. Un modelo de recuento funcional espacial para el análisis de la
heterogeneidad en el tiempo
0.7
0.7 0.7
0.6
0.6
0.6 0.5
0.5
0.4
0.5
0.4
0.3
0.4 0.3 0.2
0.1 0
0.2
0 -0.1
0.2 0.2
0.2
0.1
0.1 0.1
0
0 0
-0.1
-0.1 -0.1
-0.2
-0.2 -0.2
-0.3
0.4
0.3
0.2
0.1
-0.1
-0.2
-0.3
10 30 50 70 90 110 130 150
S
185
Capítulo 5. Un modelo de recuento funcional espacial para el análisis de la
heterogeneidad en el tiempo
186
Capítulo 5. Un modelo de recuento funcional espacial para el análisis de la
heterogeneidad en el tiempo
Figura 5.8: Campo X de log-intensidad espacial original (fila superior) y estimado (fi-
la inferior), para el tiempo t = 1/2, a través de las escalas j = 7, 8, 9, 10 (de izquierda
a derecha), sobre una rejilla regular espacial 10 × 10, a partir de los datos curva suavi-
zados
Figura 5.9: Campo X de log-intensidad espacial original (fila superior) y estimado (fi-
la inferior), para el tiempo t = 1/2, a través de las escalas j = 7, 8, 9, 10 (de izquierda
a derecha), sobre una rejilla regular espacial 30 × 30, a partir de los datos curva suavi-
zados
187
Capítulo 5. Un modelo de recuento funcional espacial para el análisis de la
heterogeneidad en el tiempo
Node 2,2 Node 2,6 Node 2,10 Node 2,14 Node 2,18
0.5 0.5 1
0.4 0.6
0.4
0 0.2 0.5
0 0.2
0 0
-0.5 -0.2 0
-0.5 -0.2
0 500 1000 0 500 1000 0 500 1000 0 500 1000 0 500 1000
Node 6,2 Node 6,6 Node 6,10 Node 6,14 Node 6,18
0.2 0.4 0.6
0.4 0.4
0 0.4
0.2 0.2 0.2
0.2
-0.2 0 0
0 0
-0.4 -0.2 -0.2
-0.2 -0.2
0 500 1000 0 500 1000 0 500 1000 0 500 1000 0 500 1000
Node 10,2 Node 10,6 Node 10,10 Node 10,14 Node 10,18
0.5 0.5 0.2
0.4 0.4
0
0.2 0
0 0.2 -0.2
0
0 -0.4
-0.2 -0.5
-0.5 -0.6
0 500 1000 0 500 1000 0 500 1000 0 500 1000 0 500 1000
Node 14,2 Node 14,6 Node 14,10 Node 14,14 Node 14,18
0.2 0.6 0.6 0.5
0.2 0.4 0.4
0 0 0.2
0.2 0
-0.2 -0.2 0
0
-0.4 -0.4 -0.2
-0.2 -0.5
-0.4
0 500 1000 0 500 1000 0 500 1000 0 500 1000 0 500 1000
Node 18,2 Node 18,6 Node 18,10 Node 18,14 Node 18,18
0.5 0.5 0.2
0.2
0 0
0
0 0
-0.2 -0.2 -0.2
-0.4 -0.4 -0.4
-0.5 -0.5
0 500 1000 0 500 1000 0 500 1000 0 500 1000 0 500 1000
Figura 5.10: Datos temporales y espaciales interpolados en una rejilla regular espacial
20 × 20
188
Capítulo 5. Un modelo de recuento funcional espacial para el análisis de la
heterogeneidad en el tiempo
Node 2,2 Node 2,6 Node 2,10 Node 2,14 Node 2,18
2 2
0.5 1
0
0 0.5 1 1
-0.5
-0.5 0
-1 0 0
0 20 40 60 0 20 40 60 0 20 40 60 0 20 40 60 0 20 40 60
Node 6,2 Node 6,6 Node 6,10 Node 6,14 Node 6,18
0 2 1.5 2
0.6
1
-0.5 0.4
1 0.5 1
0.2
-1 0 0
0
0
0 20 40 60 0 20 40 60 0 20 40 60 0 20 40 60 0 20 40 60
Node 10,2 Node 10,6 Node 10,10 Node 10,14 Node 10,18
1
0.5 0
1 1
0.5
0 0.5 -0.5
0.5
0 0
-0.5
0 -1
0 20 40 60 0 20 40 60 0 20 40 60 0 20 40 60 0 20 40 60
Node 14,2 Node 14,6 Node 14,10 Node 14,14 Node 14,18
0 1 2 1
0
-0.5 0.5 0.5
1
-0.5
0 0
-1
-1 0
0 20 40 60 0 20 40 60 0 20 40 60 0 20 40 60 0 20 40 60
Node 18,2 Node 18,6 Node 18,10 Node 18,14 Node 18,18
0
0 0 0 0
Figura 5.11: Datos temporales y espaciales interpolados en una rejilla regular espacial
20 × 20 a escala (nivel de resolución) 7
Node 2,2 Node 2,6 Node 2,10 Node 2,14 Node 2,18
0 0.1 0.2 0.4 0.4
0 0.2
-0.2 -0.1 0.2
0 500 1000 0 500 1000 0 500 1000 0 500 1000 0 500 1000
Node 6,2 Node 6,6 Node 6,10 Node 6,14 Node 6,18
-0.1 0.1 0.3 0.3 0.4
0.2 0.2
-0.15
0.05 0.1 0.3
0.1
-0.2
0 0
0 0.2
0 500 1000 0 500 1000 0 500 1000 0 500 1000 0 500 1000
Node 10,2 Node 10,6 Node 10,10 Node 10,14 Node 10,18
0.15
0.1 0.2 -0.06
0.05 0.2 0.1 -0.08
0 0.1 0.05 -0.1
-0.05 0.1 -0.12
0 0
0 500 1000 0 500 1000 0 500 1000 0 500 1000 0 500 1000
Node 14,2 Node 14,6 Node 14,10 Node 14,14 Node 14,18
0 0.4 0.2
0.05
0.1
0 0.3 0.1
-0.1
-0.05 0.05
0.2
-0.1 -0.2 0 0
0 500 1000 0 500 1000 0 500 1000 0 500 1000 0 500 1000
Node 18,2 Node 18,6 Node 18,10 Node 18,14 Node 18,18
0 0
0 0
0
-0.05
-0.1 -0.1 -0.1 -0.05
-0.1
-0.2 -0.1
-0.2 -0.2
0 500 1000 0 500 1000 0 500 1000 0 500 1000 0 500 1000
Figura 5.12: Datos curva suavizados mediante B–splines en una rejilla regular espacial
20 × 20
189
Capítulo 5. Un modelo de recuento funcional espacial para el análisis de la
heterogeneidad en el tiempo
190
Capítulo 5. Un modelo de recuento funcional espacial para el análisis de la
heterogeneidad en el tiempo
Figura 5.14: Gráficos de contorno del campo de log-intensidad observado para los
valores mensuales t = 108 (fila superior) y t = 216 (tercera fila), y el campo de log-
intensidad estimada para t = 108 (segunda fila) y t = 216 (fila inferior). Tanto los
valores observados como los estimados para los tiempos t = 108 y t = 216 se muestran
a través de las escalas j = 7, 8, 9, 10 (de izquierda a derecha), en el sistema wavelet de
Haar, a partir de los datos interpolados temporales y espaciales sobre una rejilla regular
20 × 20
191
Capítulo 5. Un modelo de recuento funcional espacial para el análisis de la
heterogeneidad en el tiempo
Figura 5.15: Gráficos de contorno del campo de log-intensidad observado para los
valores mensuales t = 108 (fila superior) y t = 216 (tercera fila), y el campo de log-
intensidad estimada para t = 108 (segunda fila) y t = 216 (fila inferior). Tanto los
valores observados como los estimados para los tiempos t = 108 y t = 216 se muestran
a través de las escalas j = 7, 8, 9, 10 (de izquierda a derecha), en el sistema wavelet
de Haar, a partir de los datos temporales interpolados y suavizados sobre una rejilla
regular 20 × 20
192
Capítulo 5. Un modelo de recuento funcional espacial para el análisis de la
heterogeneidad en el tiempo
Tabla 5.4: ALOOCVE. Promedio puntual anual de los errores de validación cruzada
193
Capítulo 5. Un modelo de recuento funcional espacial para el análisis de la
heterogeneidad en el tiempo
Figura 5.16: Promedio anual del número de muertes por enfermedades respiratorias
observadas y estimadas (arriba y abajo, respectivamente) en cada una de las 48 pro-
vincias españolas desde enero de 1980 hasta diciembre de 2015
194
A PLICACIONES A DATOS REALES
196
Aplicaciones a datos reales
197
Capítulo 6
198
Capítulo 6. Una metodología funcional bayesiana para el mapeo del riesgo de
dengue en América Latina y el Caribe
199
Capítulo 6. Una metodología funcional bayesiana para el mapeo del riesgo de
dengue en América Latina y el Caribe
6.1. Métodos
6.1.1. El modelo
200
Capítulo 6. Una metodología funcional bayesiana para el mapeo del riesgo de
dengue en América Latina y el Caribe
i = 1, . . . , N,
Ci,t
Zi,t = log = log Ci,t − log ei,t
ei,t
Ci,t − θi,t
= log θi,t + − log ei,t
θi,t
Ci,t − θi,t
= log ri,t +
θi,t
= log ri,t + vi,t , (6.2)
R0 φk = λk (R0 )φk .
201
Capítulo 6. Una metodología funcional bayesiana para el mapeo del riesgo de
dengue en América Latina y el Caribe
202
Capítulo 6. Una metodología funcional bayesiana para el mapeo del riesgo de
dengue en América Latina y el Caribe
6.1.3. Datos
203
Capítulo 6. Una metodología funcional bayesiana para el mapeo del riesgo de
dengue en América Latina y el Caribe
los países de América, que presentan hasta cinco ceros, desde 1998 hasta 2018.
Estos datos, que incluyen 672 valores con 24 ceros, están representados en la
Figura 6.1.
25000
10000
2005 2006 2007 2008 2009 2010 2011 5000
2500
1000
25
Figura 6.1: Evolución del número de DF observado desde 1998 hasta 2018 en los
países bajo estudio
204
Capítulo 6. Una metodología funcional bayesiana para el mapeo del riesgo de
dengue en América Latina y el Caribe
6.2. Resultados
Para estimar los autovalores del operador de covarianza hemos considerado
tanto el enfoque clásico como el bayesiano. Así, λ bk (R0 ) y λ
ek (R0 ) se calculan
para k = 1, . . . , M . Para obtener el estimador bayesiano, hemos utilizado la
clase de distribuciones a-priori gamma inversas, introducida en la Sección 6.1,
dando el valor = 0.05, es decir, una distribución IG(0.05, 0.05).
En la estimación de parámetros bayesianos implementada en la Sección 6.1,
se han considerado los valores empíricos de αk y βk , para k = 1, . . . , M = 19,
para la estimación bayesiana y clásica del operador de autocovarianza.
La elección del orden de truncamiento, M = 19, se ha hecho teniendo en
cuenta que los autovalores de un orden superior a T , cuando T < N , teórica-
mente son nulos, véase Bosq y Blanke (2007). Por otra parte, la determinación
del orden de truncamiento es compleja porque un valor demasiado pequeño
produce malas estimaciones pero, con un orden demasiado grande, el error de
predicción podría explotar, ver Bosq (2000). Hemos comprobado que ningu-
na de estas dos cosas es el caso de nuestra elección. Además, como criterio de
truncamiento, este orden de truncamiento lleva a un porcentaje muy alto de va-
riabilidad empírica explicada en términos de la traza aproximada del operador
de autocovarianza del proceso ARH(1), ver Ruiz-Medina et al. (2014).
Una medida de la variabilidad entre los valores observados y los estima-
dos se muestra en la Tabla 6.1. Se trata del Error Absoluto Empírico del log-
riesgo estimado en ambos casos, estimación bayesiana y clásica, | log(ri,t ) −
ri,t )|, respectivamente, para cada año, Promediado para todos los países
log(b
(i = 1, . . . , 32) bajo estudio (EAEA).
Para captar la heterogeneidad de las medidas de error desde otra pers-
pectiva, la Tabla 6.2 muestra el Error Cuadrático Empírico para el log-riesgo
tanto en la estimación bayesiana como en la clásica, (log(ri,t ) − log(b ri,t )2 y
(log(ri,t ) − N − rei,t )2 , respectivamente, para los 32 países considerados Prome-
205
Capítulo 6. Una metodología funcional bayesiana para el mapeo del riesgo de
dengue en América Latina y el Caribe
Tabla 6.1: Errores absolutos empíricos Promediados para el riesgo relativo de DF según
el promedio de los países bajo estudio, para la estimación bayesiana y clásica, (BEAEA)
y (CEAEA), respectivamente
206
Capítulo 6. Una metodología funcional bayesiana para el mapeo del riesgo de
dengue en América Latina y el Caribe
21 32
1 1 XX
EMAE = | log(ri,t ) − log(b
ri,t )|,
21 32 t=1 i=1
21 32
1 1 XX
EMSE = ri,t ))2 ,
(log(ri,t ) − log(b
21 32 t=1 i=1
207
Capítulo 6. Una metodología funcional bayesiana para el mapeo del riesgo de
dengue en América Latina y el Caribe
4
1998 1999 2000 2001 2002 2003 2004
2
−4
−8
−10
0 0 0
-2 -2 -2
log-risk
log-risk
log-risk
-4 -4 -4
-6 -6 -6
C-est. C-est. C-est.
B-est. B-est. B-est.
-8 Obs. -8 Obs. -8 Obs.
0 5 10 15 20 0 5 10 15 20 0 5 10 15 20
Years Years Years
French Guiana Costa Rica Venezuela
2 2 2
log-risk
log-risk
log-risk
1 1 1
0 0 0
-1 -1 -1
0 5 10 15 20 0 5 10 15 20 0 5 10 15 20
Years Years Years
208
Capítulo 6. Una metodología funcional bayesiana para el mapeo del riesgo de
dengue en América Latina y el Caribe
Tabla 6.3: Comparación del Error Absoluto Medio Empírico (EMAE) y el Error Cua-
drático Medio Empírico (EMSE) desde 1998 hasta 2018 para todos los países, según
las cinco metodologías aplicadas
Barbados Paraguay
2 2
0 0
-2 -2
log-risk
log-risk
-4 -4
-6 -6
Leroux Leroux
ICAR ICAR
-8 BYM -8 BYM
C-est. C-est.
B-est. B-est.
Obs. Obs.
0 2 4 6 8 10 12 14 16 18 20 22 0 2 4 6 8 10 12 14 16 18 20 22
Years Years
209
Capítulo 6. Una metodología funcional bayesiana para el mapeo del riesgo de
dengue en América Latina y el Caribe
6.3. Discusión
210
Capítulo 6. Una metodología funcional bayesiana para el mapeo del riesgo de
dengue en América Latina y el Caribe
Obsérvese que la Figura 6.2 informa sobre el riesgo en un país con respecto a
los demás en un año concreto. Podemos observar la evolución del riesgo en cada
país (comparado con el resto de los países) por año. Sin embargo, en general,
no podemos establecer una conexión directa entre esta figura y la Figura 6.1,
aunque pueda aparecer en algunos países. Por ejemplo, podemos observar que
el número de casos de DF en Argentina, a pesar de ser generalmente bajo,
exhibe una tendencia creciente a nivel global, además, el riesgo a nivel global
aumenta en comparación con el resto de los países del estudio (ver también la
Figura 6.3). Cabe destacar que en el caso argentino, el riesgo se ve mitigado por
la inclusión de toda la población en el estudio, cuando en realidad la parte sur
no es una población expuesta a la enfermedad. Por otro lado, Brasil presenta
tanto un alto número de casos de DF como un alto nivel de riesgo en todos los
años del estudio. Sin embargo, esta estabilidad del riesgo a lo largo del tiempo
(véanse las figuras 6.2 y 6.3) no refleja la gran variabilidad del número de casos
(véase la Figura 6.1).
De manera similar, Bhatt et al. (2013) predice que la transmisión del den-
gue es omnipresente en todos los trópicos, siendo los de América y Asia, los de
mayor riesgo. También muestran un mapa con niveles de riesgo durante 2010,
comparables a los calculados en nuestro trabajo. Sin embargo, Cattarino et al.
(2020) explica que los actuales mapas de riesgo del dengue proporcionan es-
timaciones de la carga de la enfermedad o de los límites de endemicidad en
lugar de la intensidad de la transmisión. Además, predicen una alta intensidad
de transmisión en todos los continentes, a ambos lados de los trópicos, con
puntos calientes en Sudamérica (Colombia, Venezuela y Brasil). Es pertinente
aclarar que los resultados de los trabajos citados no pueden compararse direc-
tamente. Emplean diferentes metodologías estadísticas y se centran en zonas y
años concretos. Sin embargo, en todos ellos, los países presentan un riesgo de
dengue similar a los encontrados con nuestro método.
211
Capítulo 6. Una metodología funcional bayesiana para el mapeo del riesgo de
dengue en América Latina y el Caribe
En la comparación con los enfoques previos (ver Tabla 6.3), se observa que
los errores de la metodología de estimación basada en un proceso ARH(1) tie-
nen un comportamiento más estable en las dos medidas propuestas, EMAE y
EMSE. Al comparar con los modelos espacio-temporales previos, los propuestos
aquí alcanzan un mejor ajuste con menores niveles de error. En particular, la
estimación bayesiana que proponemos reduce los errores encontrados con los
otros modelos y capta mejor la heterogeneidad presente en los datos.
Es importante destacar la información que dan las mediciones de la Tabla
6.3. El hecho de que, en los modelos espacio-temporales tradicionales conside-
rados, el EMSE sea grande comparado con el EMAE es indicativo de que estas
técnicas pueden presentar falsos negativos y positivos. El EMAE refleja que,
en la mayoría de los casos, estos modelos se ajustan correctamente a los datos
espacio-temporales, pero la gran diferencia con la EMSE refleja que, en momen-
tos específicos, hay grandes diferencias entre el riesgo estimado y el observado,
ver Figura 6.4 y también 6.3. Al comparar EMAE y EMSE (ver Tabla 6.3), bus-
camos valores pequeños y resultados similares, para que haya robustez en las
estimaciones.
En nuestra propuesta, EMAE y EMSE son más pequeños que en los modelos
espacio-temporales tradicionales analizados, por lo que la metodología de esti-
mación presentada muestra un mejor comportamiento y mayor estabilidad, ya
que las medidas de error son pequeñas y también cercanas. El EMSE es más pe-
queño en nuestros enfoques que en los otros modelos y en particular, en nuestra
propuesta, el EMSE es más pequeño que el EMAE, lo que corrobora la bondad
del ajuste ya que, al ser una medida cuadrática, penaliza grandes errores y sólo
puede ser más pequeña si los errores están por debajo de uno.
En el proceso de elaboración del presente trabajo han surgido varios proble-
mas que limitan el estudio. El número de casos de dengue notificados contiene
incoherencias debido a los casos asintomáticos y mal clasificados en la base de
datos primaria. Los problemas de la base de datos suelen ser las discrepancias
en las tasas de confirmación de los laboratorios (sensibilidad y especificidad) y,
a veces, los países revisan las definiciones y clasificaciones de los casos confor-
me pasa el tiempo. Estas diferencias en la calidad de los datos de los distintos
países causan dificultades para el análisis de los datos. La modelización de zo-
nas con grandes diferencias en sus extensiones da lugar a estimaciones un poco
peores, especialmente en países pequeños como las islas del Caribe, donde la
212
Capítulo 6. Una metodología funcional bayesiana para el mapeo del riesgo de
dengue en América Latina y el Caribe
213
Capítulo 6. Una metodología funcional bayesiana para el mapeo del riesgo de
dengue en América Latina y el Caribe
214
Capítulo 7
216
Capítulo 7. Análisis de la mortalidad por COVID–19 mediante curvas
multivariantes de regresión a partir de datos blandos y mediante aprendizaje
automático
217
Capítulo 7. Análisis de la mortalidad por COVID–19 mediante curvas
multivariantes de regresión a partir de datos blandos y mediante aprendizaje
automático
218
Capítulo 7. Análisis de la mortalidad por COVID–19 mediante curvas
multivariantes de regresión a partir de datos blandos y mediante aprendizaje
automático
219
Capítulo 7. Análisis de la mortalidad por COVID–19 mediante curvas
multivariantes de regresión a partir de datos blandos y mediante aprendizaje
automático
donde C0∞ (D) denota el espacio de funciones infinitas diferenciables, con sopor-
te compacto contenido en D. Para cada conjunto acotado D ∈ B d , se define
donde
Z
gt (ψp,$p , θ(p)) = gt (z, θ(p))ψp,$p (z)dz
Dp
220
Capítulo 7. Análisis de la mortalidad por COVID–19 mediante curvas
multivariantes de regresión a partir de datos blandos y mediante aprendizaje
automático
Z
hf, giH = f (z)g(z)dz, (7.4)
Rd
N
X
gt (θ(p)) = (Ak (p) cos(ϕk (p)t) + Bk (p) sen(ϕk (p)t)) , p = 1, . . . , P, t ∈ R+ ,
k=1
(7.5)
cuyos parámetros dependientes del espacio están dados por las escalas tempora-
les (ϕ1 (·), . . . , ϕN (·)) , y los coeficientes de Fourier (A1 (·), B1 (·), . . . , AN (·), BN (·)) .
Para simplificar, consideraremos que los parámetros de escala ϕk , k = 1, . . . , N,
son conocidos, y fijos sobre las regiones espaciales P . Además, Ck2 (·) = A2k (·) +
Bk2 (·) > 0, para k = 1, . . . , N, donde N denota el parámetro de truncamien-
to, que se seleccionará según el factor penalizado propuesto en Chapelle et al.
(2002), como explicamos con más detalle en la Sección 7.2. Así,
221
Capítulo 7. Análisis de la mortalidad por COVID–19 mediante curvas
multivariantes de regresión a partir de datos blandos y mediante aprendizaje
automático
P
X
εt (ψp,$p ) = ρ(ψq,$q )(ψp,$p )εt−1 (ψq,$q ) + νt (ψp,$p ), (7.6)
q=1
Aquí, νt (ψp,$p ), p = 1, . . . , P , t ∈ R+ , se supone que son vectores gaussianos
independientes de media cero P –dimensional. Para p, q ∈ {1, . . . P }, la proyec-
ción ρ(ψp,$p )(ψq,$q ) mantiene la autocorrelación lineal temporal en cada región
espacial para p = q, y la correlación cruzada lineal temporal entre regiones para
p 6= q del error de regresión {εt (·), t ∈ R+ }, véase Bosq (2000).
T
1X 2
LT (θbT (p)) = ı́nf c LT (θ(p)) = ı́nf ln (Λt ) (ψp,$p ) − gt (θ(p)) . (7.7)
θ(p)∈Θ θ(p)∈Θc T t=1
222
Capítulo 7. Análisis de la mortalidad por COVID–19 mediante curvas
multivariantes de regresión a partir de datos blandos y mediante aprendizaje
automático
ción entre la función de pérdida cuadrática esperada en θbT (p), y el valor espe-
rado del riesgo mínimo empírico a partir de la identidad
h i2
E ln (Λt ) (ψp,$p ) − gt (ψp,$p , θT (p)
b
−1 PN !
h i N 1/λi
= E LT (θbT (p)) 1 − 1 + i=1 , (7.8)
T T
1 eT
T
−1 T
LT (θT (p)) = R (p) IT ×T − Φ Φ Φ
b Φ R(p),
e (7.9)
T
con, para p = 1, . . . , P ,
∞
!
X
R(p)
e = (Ak (p) cos(ϕk t) + Bk (p) sen(ϕk t)) + εt (ψp,$p ), t = 1, . . . , T .
k=N +1
\
ln (Λt )(ψp,$p ) = gt (θbT (p)), p = 1, . . . , P, (7.10)
véase el Teorema 1 en Ivanov et al. (2015) sobre las condiciones para la consis-
tencia débil de (7.10).
223
Capítulo 7. Análisis de la mortalidad por COVID–19 mediante curvas
multivariantes de regresión a partir de datos blandos y mediante aprendizaje
automático
T
bY (ψp,$p )(ψq,$q ) = 1X
R0,T Yt (ψp,$p )Yt (ψq,$q ),
T t=1
T −1
bY (ψp,$p )(ψq,$q ) = 1 X
R1,T Yt (ψq,$q )Yt+1 (ψp,$p ), p, q = 1, . . . , P,
T − 1 t=1
(7.11)
k(T )
X ψp,$p , φk,T ψq,$q , φl,T
ρbk(T ) (ψp,$p )(ψq,$q ) = H H bY (φk,T )(φl,T ),
R1,T
bY )
λk,T (R
k,l=1 0,T
(7.12)
T
X
bY =
R bY )[φk,T ⊗ φk,T ],
λk,T (R (7.13)
0,T 0,T
k=1
P
X
k(T )
Ybt (ψp,$p ) = ρbk(T ) (ψq,$q )(ψp,$p )Yt−1 (ψq,$q ), t ≥ 1. (7.14)
q=1
224
Capítulo 7. Análisis de la mortalidad por COVID–19 mediante curvas
multivariantes de regresión a partir de datos blandos y mediante aprendizaje
automático
L
ep (ε1p , . . . , εT p , ε0q , . . . , ε(T −1)q ρ(ψq,$q )(ψp,$p ), q = 1, . . . , P )
2
1
PT PP
exp − 2σ2 t=1 εt (ψp,$p ) − q=1 εt−1 (ψq,$q )ρ(ψq,$q )(ψp,$p )
p
= √ T
σp 2π
P
Y apq −1 bpq −1
× ρ(ψq,$q )(ψp,$p ) 1 − ρ(ψq,$q )(ψp,$p )
q=1
I{0<ρ(ψq,$q )(ψp,$p )<1}
×
B(apq , bpq )
T
!
1 1 X 2
= √ T exp − 2 νt (ψp,$p )
σp 2π 2σp t=1
P
Y apq −1 bpq −1
× ρ(ψq,$q )(ψp,$p ) 1 − ρ(ψq,$q )(ψp,$p )
q=1
I{0<ρ(ψq,$q )(ψp,$p )<1}
× ,
B(apq , bpq )
(7.15)
Γ(apq )Γ(bpq )
B(apq , bpq ) = .
Γ(apq + bpq )
P
X
εet (ψp,$p ) = ρe(ψq,$q )(ψp,$p )εt−1 (ψq,$q ), t ≥ 1, (7.16)
q=1
con ρe(ψ1,$1 )(ψp,$p ), . . . , ρe(ψP,$P )(ψp,$p ) que se calcula maximizando (7.15),
para encontrar la moda de la distribución a-posteriori (véase Bosq y Ruiz-
Medina (2014), donde se introduce la estimación bayesiana en un marco
225
Capítulo 7. Análisis de la mortalidad por COVID–19 mediante curvas
multivariantes de regresión a partir de datos blandos y mediante aprendizaje
automático
226
Capítulo 7. Análisis de la mortalidad por COVID–19 mediante curvas
multivariantes de regresión a partir de datos blandos y mediante aprendizaje
automático
SC/PE A
b1 (·) A
b2 (·) A
b3 (·) A
b4 (·) A
b5 (·) A
b6 (·)
C1 3.6343 -0.4814 -0.0075 -0.0258 0.0189 0.0193
C2 3.4345 -0.3923 0.0416 0.0265 -0.0709 -0.0572
C3 3.2031 -0.1364 -0.0088 0.0221 0.0430 0.0289
C4 3.1445 -0.1118 0.0041 0.0337 0.0062 0.0072
C5 3.1015 -0.0693 -0.0345 0.0352 -0.0112 0.0003
C6 3.1347 -0.1397 0.0020 0.0300 -0.0061 -0.0002
C7 4.0591 -0.5487 -0.0907 0.0951 0.0992 0.0842
C8 3.8032 -0.5500 -0.1007 0.0633 0.0139 0.0277
C9 4.5095 -0.7435 -0.1134 0.1809 0.2231 0.2026
C10 3.6321 -0.4685 -0.0540 0.0384 -0.0152 0.0011
C11 3.2967 -0.2274 -0.0083 0.0553 0.0250 0.0240
C12 3.3454 -0.2122 -0.0927 -0.0330 0.0724 0.0679
C13 4.8419 -0.6790 -0.2455 0.0311 0.0554 0.0667
C14 3.0941 -0.1037 0.0210 0.0141 -0.0016 0.0041
C15 3.2877 -0.2598 -0.0524 0.0842 -0.0423 -0.0348
C16 3.6870 -0.4302 -0.0086 0.0078 -0.0027 -0.0017
C17 3.2197 -0.2071 0.0162 0.0079 0.0206 0.0110
SC/PE B
b1 (·) B
b2 (·) B
b3 (·) B
b4 (·) B
b5 (·) B
b6 (·)
C1 0 -0.0052 -0.1330 -0.0123 0.0064 -0.0195
C2 0 -0.0367 -0.0998 -0.0462 -0.0343 -0.0107
C3 0 -0.0531 -0.0074 -0.0142 -0.0003 0.0020
C4 0 -0.0074 -0.0284 -0.0151 -0.0092 0.0012
C5 0 0.0433 -0.0438 -0.0116 -0.0118 0.0046
C6 0 0.0018 -0.0174 -0.0068 -0.0089 0.0000
C7 0 -0.0365 -0.2451 -0.1791 -0.0820 0.0026
C8 0 0.0953 -0.2389 -0.0431 -0.0313 -0.0045
C9 0 -0.1587 -0.4054 -0.2269 -0.1010 0.0047
C10 0 0.1118 -0.1579 -0.0458 -0.0418 -0.0220
C11 0 0.0754 -0.1138 -0.0166 -0.0048 0.0072
C12 0 -0.1104 -0.1338 0.1330 0.0761 -0.0017
C13 0 0.4654 -0.1302 -0.1602 -0.1061 -0.0038
C14 0 0.0355 -0.0560 0.0119 0.0025 -0.0044
C15 0 -0.0187 -0.0021 -0.0897 -0.0562 0.0134
C16 0 0.0025 -0.0707 -0.0638 -0.0439 -0.0267
C17 0 0.0389 -0.0270 -0.0174 -0.0006 0.0019
227
Capítulo 7. Análisis de la mortalidad por COVID–19 mediante curvas
multivariantes de regresión a partir de datos blandos y mediante aprendizaje
automático
3500 3000
3000
2500
2500
2000
2000
1500
1500
1000
1000
500 500
0 0
−500
0 50 100 150 200 250 300 −500
0 50 100 150 200 250 300
120 100
100
80
80
60
60
40
40
20
20
0 0
−20 −20
0 50 100 150 200 250 300 0 50 100 150 200 250 300
4 4
3 3
2 2
1 1
0 0
−1 −1
−2 −2
0 50 100 150 200 250 300 0 50 100 150 200 250 300
228
Capítulo 7. Análisis de la mortalidad por COVID–19 mediante curvas
multivariantes de regresión a partir de datos blandos y mediante aprendizaje
automático
P
1 X
L265 (ωi ) = L265 (ωi , θb265 (p)), ωi ∈ Ω, i = 1, . . . , 1000,
P p=1
229
Capítulo 7. Análisis de la mortalidad por COVID–19 mediante curvas
multivariantes de regresión a partir de datos blandos y mediante aprendizaje
automático
0.14
0.13
0.12
0.11
0.1
0.09
0.08
0.07
0.06
0.05
0.04
0 100 200 300 400 500 600 700 800 900 1000
BOOTSTRAP HISTOGRAM
BOOTSTRAP PROBABILITY DENSITY
250
25
200
20
150
15
100 10
50 5
0 0
0.04 0.05 0.06 0.07 0.08 0.09 0.1 0.11 0.12 0.13 0.14 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16
Figura 7.2: Valores muestrales de la media espacial del riesgo mínimo empírico en
la regresión trigonométrica, en la parte superior. Histograma bootstrap (parte inferior
izquierda) y densidad de probabilidad bootstrap (parte inferior derecha) de la media
espacial del riesgo mínimo empírico de regresión, basado en 1000 muestras bootstrap
230
Capítulo 7. Análisis de la mortalidad por COVID–19 mediante curvas
multivariantes de regresión a partir de datos blandos y mediante aprendizaje
automático
0.7 6
0.65 5
0.6 4
BASED ON 100 BOOTSRAT SAMPLES
0.55 3
0.5 2
0.45 1
0.4 0
0 50 100 150 200 250 300 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
231
Capítulo 7. Análisis de la mortalidad por COVID–19 mediante curvas
multivariantes de regresión a partir de datos blandos y mediante aprendizaje
automático
295 295
day 10 day 13 day 16 day 19 day 22 day 10 day 13 day 16 day 19 day 22
80 80
40 40
day 25 day 28 day 31 day 34 day 37 day 25 day 28 day 31 day 34 day 37
20 20
10 10
day 40 day 43 day 46 day 49 day 52 day 40 day 43 day 46 day 49 day 52
5 5
day 55 day 58 day 61 day 64 day 55 day 58 day 61 day 64
0 0
295 295
day 10 day 13 day 16 day 19 day 22 day 10 day 13 day 16 day 19 day 22
80 80
40 40
day 25 day 28 day 31 day 34 day 37 day 25 day 28 day 31 day 34 day 37
20 20
10 10
day 40 day 43 day 46 day 49 day 52 day 40 day 43 day 46 day 49 day 52
5 5
day 55 day 58 day 61 day 64 day 55 day 58 day 61 day 64
0 0
Figura 7.4: Mapas de riesgo de mortalidad por COVID–19, desde el 8 de marzo hasta
el 13 de mayo de 2020. Mapas observados (arriba) y estimados (abajo), calculados a
partir de la regresión trigonométrica, combinados con predictores residuales clásicos
(abajo izquierda) y bayesianos (abajo derecha)
232
Capítulo 7. Análisis de la mortalidad por COVID–19 mediante curvas
multivariantes de regresión a partir de datos blandos y mediante aprendizaje
automático
T \
ln(Λ t )(ψp,ρp ) − ln(Λt )(ψp,ρp )
1X
, p = 1, . . . , P. (7.18)
T t=1 ln(Λ )(ψ ) + ln(Λ \
t p,ρp t )(ψp,ρp ) /2
233
Capítulo 7. Análisis de la mortalidad por COVID–19 mediante curvas
multivariantes de regresión a partir de datos blandos y mediante aprendizaje
automático
Hemos calculado la media de los SMAPE obtenidos en cada una de las k ite-
raciones del procedimiento de validación cruzada aleatoria k–fold. Esta técnica
de validación consiste en dividir aleatoriamente la muestra funcional en una
muestra de entrenamiento y otra de validación en cada una de las iteraciones
k. El ajuste del modelo se realiza a partir de la muestra de entrenamiento, y
las salidas objetivo se definen a partir de la muestra de validación o de prue-
ba. Al ejecutar cada modelo diez veces y promediar los SMAPE, eliminamos las
fluctuaciones debidas a los pesos iniciales aleatorios (para los modelos MLP y
BNN), y las diferencias en la estimación de los parámetros en todos los mé-
todos, debido a la especificación aleatoria de la división de la muestra en el
procedimiento de validación cruzada aleatoria k–fold.
Los SMAPE de validación cruzada aleatoria 10–fold basados en diez ejecu-
ciones se muestran en la Tabla 7.6, para las seis técnicas ML probadas, GRNN,
MLP, SVR, BNN, RBF y GP, cuando se consideran datos duros (véase también
la Tabla 7.13 de la Sección 7.5 sobre los resultados de validación cruzada alea-
toria 5–fold). La Tabla 7.7 proporciona los resultados de validación cruzada
aleatoria 10–fold, de la categoría de datos blandos (véase también la Tabla 7.14
de la Sección 7.5 sobre los resultados de validación cruzada aleatoria 5–fold).
Los correspondientes resultados de validación cruzada del enfoque presentado
a partir de datos blandos se muestran en la Tabla 7.8.
La selección de hiperparámetros del modelo ML se ha realizado aplicando
una validación cruzada aleatoria k–fold (k = 5, 10). Nuestra selección se ha
realizado a partir de un conjunto adecuado de candidatos. En concreto, los nú-
meros óptimos de nodos ocultos (NH) en la implementación de MLP y BNN se
han seleccionado de los conjuntos de candidatos [0, 1, 3, 5, 7, 9] y [1, 3, 5, 7, 9],
respectivamente. Los resultados de la validación cruzada aleatoria en ambos
casos, k = 5, 10, llevan a la misma elección del valor óptimo de NH, a saber,
NH= 1 para MLP, y NH= 5 para BNN. Este último muestra ligeras diferencias
con respecto a los valores NH= 3, 7, en la implementación de validación cruza-
da aleatoria 10–fold. Del mismo modo, hemos seleccionado los respectivos pará-
metros de dispersión β y ancho de banda h en los procedimientos RBF y GRNN.
Así, después de aplicar una validación cruzada aleatoria k–fold, con k = 5, 10,
se obtienen los valores óptimos β = 2.5, y h = 0.05, a partir de los conjun-
tos de candidatos [2.5, 5, 7.5, 10, 12.5, 15, 17.5, 20] y [0.05, 0.1, 0.2, 0.3, 0.5, 0.6, 0.7],
respectivamente (véase la Sección 7.5).
234
Capítulo 7. Análisis de la mortalidad por COVID–19 mediante curvas
multivariantes de regresión a partir de datos blandos y mediante aprendizaje
automático
Tabla 7.6: Categoría de datos duros. Promedio de los SMAPE, basado en 10 ejecu-
ciones de validación cruzada aleatoria 10–fold. Como se indica, los valores mostrados
deben ser multiplicados por 10−2
Tabla 7.7: Categoría de datos blandos. Promedio de los SMAPE, basado en 10 ejecu-
ciones de validación cruzada aleatoria 10–fold. Como se indica, los valores mostrados
deben ser multiplicados por 10−2
235
Capítulo 7. Análisis de la mortalidad por COVID–19 mediante curvas
multivariantes de regresión a partir de datos blandos y mediante aprendizaje
automático
SC C. k10 B. k10
C1 0.0024 0.7106(10)−3
C2 0.0019 0.4003(10)−3
C3 0.0016 0.6797(10)−3
C4 0.0017 0.4367(10)−3
C5 0.0023 0.6530(10)−3
C6 0.0018 0.5854(10)−3
C7 0.0017 0.6341(10)−3
C8 0.0016 0.6593(10)−3
C9 0.0013 0.5979(10)−3
C10 0.0019 0.6954(10)−3
C11 0.0017 0.5444(10)−3
C12 0.0016 0.5016(10)−3
C13 0.0020 0.4832(10)−3
C14 0.0026 0.6544(10)−3
C15 0.0023 0.6616(10)−3
C16 0.0015 0.7134(10)−3
C17 0.0022 0.6781(10)−3
M. 0.0019 0.60524(10)−3
T. 0.0321 0.0103
236
Capítulo 7. Análisis de la mortalidad por COVID–19 mediante curvas
multivariantes de regresión a partir de datos blandos y mediante aprendizaje
automático
237
Capítulo 7. Análisis de la mortalidad por COVID–19 mediante curvas
multivariantes de regresión a partir de datos blandos y mediante aprendizaje
automático
238
Capítulo 7. Análisis de la mortalidad por COVID–19 mediante curvas
multivariantes de regresión a partir de datos blandos y mediante aprendizaje
automático
250 800
700
200
600
500
150
400
100 300
200
50
100
0 0
1.5 2 2.5 3 3.5 4 4.5 5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5
−3
x 10
−3 x 10
250 3500
3000
200
2500
150 2000
1500
100
1000
50
500
0 0
4 5 6 7 8 9 10 11 12 13 2 4 6 8 10 12 14
−4 −4
x 10 x 10
4500
200 4000
3500
150 3000
2500
100 2000
1500
50 1000
500
0 0
2 3 4 5 6 7 8 2 3 4 5 6 7 8 9 10
−4 −4
x 10 x 10
Figura 7.5: Categoría de datos duros. A partir de 1000 muestras bootstrap, se repre-
sentan las densidades de probabilidad y los histogramas de los SMAPE promediados
espacialmente, para GRNN (arriba), MLP (centro) y SVR lineal (abajo)
239
Capítulo 7. Análisis de la mortalidad por COVID–19 mediante curvas
multivariantes de regresión a partir de datos blandos y mediante aprendizaje
automático
5000
200
4000
150
3000
100
2000
50 1000
0
0 2 3 4 5 6 7 8
2 3 4 5 6 7 8
−4
−4 x 10
x 10
7000
200
6000
150 5000
4000
100
3000
2000
50
1000
0 0
3 4 5 6 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5
−4 −4
x 10 x 10
BOOTSTRAPPED HISTOGRAM GP
BOOTSTRAPPED PROBABILITY DENSITY GP
250
7000
200 6000
5000
150
4000
100 3000
2000
50
1000
0
3 4 5 6 0
−4 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5
x 10 −4
x 10
Figura 7.6: Categoría de datos duros. A partir de 1000 muestras bootstrap, se repre-
sentan las densidades de probabilidad y los histogramas de los SMAPE promediados
espacialmente, para BNN (arriba), RBF (centro) y GP (abajo)
240
Capítulo 7. Análisis de la mortalidad por COVID–19 mediante curvas
multivariantes de regresión a partir de datos blandos y mediante aprendizaje
automático
1000
200
800
150
600
100
400
50
200
0 0
0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 3.5
−3 −3
x 10 x 10
BOOTSTRAPPED HISTOGRAM MLP (SD) BOOTSTRAPPED PROBABILITY DENSITY MLP (SD)
250 2500
200 2000
150 1500
100 1000
50 500
0 0
0.4 0.6 0.8 1 1.2 1.4 1.6 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8
−3 −3
x 10 x 10
4000
250
3500
200
3000
2500
150
2000
100
1500
1000
50
500
0
3 4 5 6 7 8 9 10 11 0
−4
2 3 4 5 6 7 8 9 10 11 12
x 10 −4
x 10
Figura 7.7: Categoría de datos blandos. A partir de 1000 muestras bootstrap, se re-
presentan las densidades de probabilidad y los histogramas de los SMAPE promediados
espacialmente, para GRNN (arriba), MLP (centro) y SVR no lineal (abajo)
241
Capítulo 7. Análisis de la mortalidad por COVID–19 mediante curvas
multivariantes de regresión a partir de datos blandos y mediante aprendizaje
automático
5000
200
4000
150
3000
100
2000
50
1000
0 0
2 3 4 5 6 7 8 2 3 4 5 6 7 8
−4 −4
x 10 x 10
250
2
200
1.5
150
1
100
0.5
50
0 0
1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 3
−4 −4
x 10 x 10
200 12000
10000
150
8000
100 6000
4000
50
2000
0
1.6 1.8 2 2.2 2.4 2.6 2.8 3 3.2 3.4 3.6 0
−4 1 2 3
x 10 −4
x 10
Figura 7.8: Categoría de datos blandos. A partir de 1000 muestras bootstrap, se re-
presentan las densidades de probabilidad y los histogramas de los SMAPE promediados
espacialmente, para BNN (arriba), RBF (centro) y GP (abajo)
242
Capítulo 7. Análisis de la mortalidad por COVID–19 mediante curvas
multivariantes de regresión a partir de datos blandos y mediante aprendizaje
automático
250 1200
1000
200
800
150
600
100
400
50
200
0 0
1.5 2 2.5 3 3.5 4 4.5 5 1.5 2 2.5 3 3.5 4 4.5 5
−3 −3
x 10 x 10
BAYESIAN BOOTSTRAPPED HISTOGRAM (SD) BAYESIAN BOOTSTRAPPED PROBABILITY DENSITY (SD)
300 8000
7000
250
6000
200
5000
150 4000
3000
100
2000
50
1000
0 0
2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 1 2 3 4 5 6 7
−4
x 10 x 10
−4
Figura 7.9: Categoría de datos blandos. A partir de 1000 muestras bootstrap, se re-
presentan las densidades de probabilidad y los histogramas de los SMAPE promediados
espacialmente, para la regresión trigonométrica combinada con la predicción residual,
basada en los momento empíricos, clásica (arriba) y bayesiana (abajo)
243
Capítulo 7. Análisis de la mortalidad por COVID–19 mediante curvas
multivariantes de regresión a partir de datos blandos y mediante aprendizaje
automático
de la respuesta:
NH
X
yb = η0 + ηk g(βkT x), (7.19)
k=1
NH
X
ybm = yb(xm ) = ln (λ\
t+1 ) (φm ) = η0 + ηk g(βkT xm ), (7.20)
k=1
donde, para m = 1, . . . , M (T ),
244
Capítulo 7. Análisis de la mortalidad por COVID–19 mediante curvas
multivariantes de regresión a partir de datos blandos y mediante aprendizaje
automático
NH
kxm − cj k2
X
ybt+1 (xm ) = ln (λ\
t+1 ) (φm ) = ηj exp (7.22)
j=1
β2
M
1 X
L = kβk2 + C |ym − f (xm )| . (7.24)
2 m=1
245
Capítulo 7. Análisis de la mortalidad por COVID–19 mediante curvas
multivariantes de regresión a partir de datos blandos y mediante aprendizaje
automático
Esta distribución a-priori pone más peso en los valores pequeños de los pará-
metros de red cercanos a cero. La densidad de probabilidad a-posteriori, dados
los datos observados O = o, y el valor ν del parámetro de regularización, se
246
Capítulo 7. Análisis de la mortalidad por COVID–19 mediante curvas
multivariantes de regresión a partir de datos blandos y mediante aprendizaje
automático
define como
p(o/w, ν)p(w/ν)
p(w/o, ν) = . (7.28)
p(o/ν)
Considerando que los errores también se distribuyen de forma gaussiana,
la probabilidad condicionada de los datos O dados los parámetros ν y w, se
obtiene como ν M/2
p(o/w, ν) = exp (−νEO (o)) , (7.29)
π
donde M denota el número de puntos de datos de entrenamiento. A partir de
las ecuaciones (7.27)–(7.29),
247
Capítulo 7. Análisis de la mortalidad por COVID–19 mediante curvas
multivariantes de regresión a partir de datos blandos y mediante aprendizaje
automático
Y = Z + , (7.34)
248
Capítulo 7. Análisis de la mortalidad por COVID–19 mediante curvas
multivariantes de regresión a partir de datos blandos y mediante aprendizaje
automático
249
Capítulo 7. Análisis de la mortalidad por COVID–19 mediante curvas
multivariantes de regresión a partir de datos blandos y mediante aprendizaje
automático
250
Capítulo 7. Análisis de la mortalidad por COVID–19 mediante curvas
multivariantes de regresión a partir de datos blandos y mediante aprendizaje
automático
Tabla 7.12: Categoría de datos blandos. Promedio de los SMAPE, basado en 10 eje-
cuciones de validación cruzada aleatoria 5–fold. Como se indica, los valores mostrados
deben ser multiplicados por 10−2
251
Capítulo 7. Análisis de la mortalidad por COVID–19 mediante curvas
multivariantes de regresión a partir de datos blandos y mediante aprendizaje
automático
252
Capítulo 7. Análisis de la mortalidad por COVID–19 mediante curvas
multivariantes de regresión a partir de datos blandos y mediante aprendizaje
automático
0 0 0 0 0
0 20 40 60 0 20 40 60 0 20 40 60 0 20 40 60 0 20 40 60 0 20 40 60
Figura 7.12: Categoría de datos duros. Curvas de mortalidad acumulada por COVID–
19 observadas y estimadas, a partir de la implementación de Redes Neuronales de
Regresión Generalizada (GRNN), Perceptrón Multicapa (MLP), Regresión de Soporte
Vectorial (SVR), Redes Neuronales Bayesianas (BNN), Redes Neuronales a partir de
Bases de Funciones Radiales (RBF), y Procesos Gaussianos (GP)
Andalucía Aragón Asturias Islas Baleares Canarias Cantabria
800 300 150 200
200
1000 600 150
150 100
200
400 100 100
500
100 50
200 50 50
0 0 0 0 0
0 20 40 60 0 20 40 60 0 20 40 60 0 20 40 60 0 20 40 60 0 20 40 60
0 0 0 0 0
0 20 40 60 0 20 40 60 0 20 40 60 0 20 40 60 0 20 40 60 0 20 40 60
Figura 7.13: Categoría de datos blandos. Curvas de mortalidad acumulada por por
COVID–19 observadas y estimadas, a partir de la implementación de Redes Neuronales
de Regresión Generalizada (GRNN), Perceptrón Multicapa (MLP), Regresión de Sopor-
te Vectorial (SVR), Redes Neuronales Bayesianas (BNN), Redes Neuronales a partir de
Bases de Funciones Radiales (RBF), Procesos Gaussianos (GP), y regresión trigonomé-
trica combinada con predicción residual clásica y bayesiana
253
Capítulo 7. Análisis de la mortalidad por COVID–19 mediante curvas
multivariantes de regresión a partir de datos blandos y mediante aprendizaje
automático
Tabla 7.13: Categoría de datos duros. Promedio de los SMAPE, basado en 10 ejecu-
ciones de validación cruzada aleatoria 5–fold. (Como se indica, los valores mostrados
deben multiplicarse por 10−2 )
254
Capítulo 7. Análisis de la mortalidad por COVID–19 mediante curvas
multivariantes de regresión a partir de datos blandos y mediante aprendizaje
automático
Tabla 7.14: Categoría de datos blandos. Promedio de los SMAPE, basado en 10 ejecu-
ciones de validación cruzada aleatoria 5–fold. (Como se indica, los valores mostrados
deben multiplicarse por 10−2 )
255
Capítulo 7. Análisis de la mortalidad por COVID–19 mediante curvas
multivariantes de regresión a partir de datos blandos y mediante aprendizaje
automático
256
C ONCLUSIONES Y LÍNEAS ABIERTAS
258
Conclusiones y líneas abiertas
259
Conclusiones y líneas abiertas
260
Conclusiones y líneas abiertas
Uno de los principales problemas que habrá que afrontar en un futuro próximo
es la definición de una versión multivariante de la clase de procesos puntuales
introducidos, así como su estimación funcional espectral en un marco multiva-
riante infinito-dimensional.
En el Capítulo 5, la metodología de predicción funcional espacial multies-
calar presentada permite el análisis de la heterogeneidad en diferentes escalas
temporales del campo de log-intensidad. Es bien sabido que las técnicas de
preprocesamiento de FDA (por ejemplo, suavizado mediante B–splines) suelen
ocultar o eliminar la variación local a niveles de alta resolución (véase, por
ejemplo, Goia y Vieu (2016); Horváth y Kokoszka (2012); Müller y Stadtmü-
ller (2005), entre otros). El enfoque de estimación adoptado en este capítulo
supera esta limitación, proporcionando un marco más flexible. Así, se puede
realizar una elección adecuada de la escala en la que se debe analizar el campo
de log-intensidad, según los objetivos del estudio, y las incertidumbres en los
recuentos asociadas a la falta de información de la muestra.
El enfoque basado en la entropía relativa a la frecuencia funcional espacial,
que se ha presentado en el Capítulo 5, no requiere información adicional sobre
la distribución a-priori de probabilidad de los parámetros, como en el marco
bayesiano. En particular, en Vicente et al. (2020), los autores presentan un
novedoso enfoque sobre la modelización de recuento multivariante bayesiano
espacio-temporal a partir de datos de zonas (véanse también las referencias en
el mismo para una visión general del estado de la técnica).
La técnica de reducción de dimensiones basada en FDA que aquí se conside-
ra evita los graves problemas de cálculo que suelen surgir cuando las matrices
de covarianza de alta dimensión, asociadas a variables gaussianas latentes e hi-
perparámetros, están involucradas en la metodología de estimación bayesiana.
De hecho, los modelos autorregresivos condicionados multivariantes (MCARs)
se han desarrollado ampliamente en un marco plenamente bayesiano. Los en-
foques basados en Cadenas de Markov Monte Carlo (MCMC) y la Aproxima-
ción Anidada Integrada de Laplace (INLA), véase, por ejemplo, Lindgren y Rue
(2015), y Rue et al. (2009), suelen aplicarse para el ajuste de los modelos en
este contexto. Sin embargo, siguen abiertos varios problemas en relación con
la separabilidad y la modelización de covarianza espacio-temporal, cuando se
adopta un marco de FDA.
El enfoque de estimación presentado en dicho Capítulo 5 proporciona una
261
Conclusiones y líneas abiertas
262
Conclusiones y líneas abiertas
263
C ONCLUSIONS AND OPEN LINES OF
RESEARCH
265
Conclusions and open lines
266
Conclusions and open lines
the fit of dependence models with different rank, an excess of smoothing is ob-
served while, for the log-Gaussian ARH(1) estimation, the smoothing is weaker,
allowing for more heterogeneity in the spatial patterns.
In relation to the contributions of Chapter 4, the extension of the Cox family
of processes introduced in Chapter 3 to the spatial context stands out. In par-
ticular, this extension makes possible the application of estimation techniques,
based on second-order statistics, under homogeneity in space and heterogeneity
in time, contributing to the analysis of point patterns and counting processes,
by means of more flexible models.
On the other hand, recent advances in computer science have made possible
the application of feasible probabilistic prediction techniques for the statistical
analysis of complex data in Hilbert, Banach, or metric/semi-metric spaces, in-
cluding the case of data that are not in a vector space.
FDA techniques allow the analysis of local regularity properties of high-
dimensional data of a continuous nature, or else, coming from an underlying
phenomenon or process continuous in space and/or time, avoiding a bad ap-
proach to classical multivariate estimation problems. In particular, we refer the
reader to Dai y Müller (2018); Goia y Vieu (2016); Horváth y Kokoszka (2012);
Ruiz-Medina (2012), among others. This Chapter 4 puts more complex mathe-
matical structures to work, dealing with elements of Hilbert spaces in the fra-
mework of point pattern analysis.
In summary, the approach presented in this chapter is a relevant contribu-
tion in the field of spectral analysis of infinite-dimensional spatial linear proces-
ses under spatial homogeneity, as well as a remarkable contribution to the field
of statistical analysis of counting processes, in the context of Spatio-Temporal
Statistics.
Since the spatial functional estimation approach presented here is applied
in the spectral domain, its continuous spatial parameter version can be for-
mulated directly from the results developed in Chapter 4, with an appropriate
reformulation of the conditions given. One of the main problems to be faced
in the near future is the definition of a multivariate version of the introduced
class of point processes, as well as their spectral functional estimation in an
infinite-dimensional multivariate framework.
In Chapter 5, the presented multiscale spatial functional prediction metho-
dology allows the analysis of heterogeneity at different time scales of the log-
267
Conclusions and open lines
268
Conclusions and open lines
269
B IBLIOGRAFÍA
271
O. O. Aalen, O. Borgan, y H. K. Gjessing. Survival and event history analysis: a
process point of view. Springer Science & Business Media, New-York, 2008.
272
J. M. Angulo y M. D. Ruiz-Medina. On the orthogonal representation of gene-
ralized random fields. Statistics & Probability Letters, 31:145–153, 1997.
273
A. Aswi, S. M. Cramb, P. Moraga, y K. Mengersen. Bayesian spatial and spatio-
temporal approaches to modelling dengue fever: a systematic review. Epi-
demiology and Infection, 33:1–14, 2019. URL [Link]
S0950268818002807.
J-M. Azaïs y M. Wschebor. Level Sets and Extrema of Random Processes and
Fields. John Wiley & Sons, Inc., 2009.
274
S. Basu y G. C. Reinsel. Properties of the spatial unilateral first-order ARMA
model. Advances in Applied Probability, 25(3):631–648, 1993.
275
B. M. Bolker y B. Grenfell. Impact of vaccination on the spatial correlation
and persistence of measles dynamics. Proceedings of the National Academy of
Sciences, 93:12648–12653, 1996.
D. R. Brillinger. Time Series: Data Analysis and Theory. Holden Day, San Fran-
cisco, 1981.
276
A. Chakraborty y A. E. Gelfand. Analyzing spatial point patterns subject to
measurement error. Bayesian Analysis, 5:97–122, 2001.
277
C. F. Chung. Estimating a generalized long memory process. Journal of Econo-
metrics, 73:237–259, 1996.
R. Dahlhaus. Spectral analysis with tapered data. Journal of Time Series Analy-
sis, 4(3):163–175, 1983.
S. Dalal y A. McIntosh. When to stop testing for large software systems with
changing code. IEEE Transactions on Software Engineering, 20:318–323, 1994.
278
D. Daley y D. Vere-Jones. An Introduction to the Theory of Point Processes: Volume
II: General Theory and Structure. New York: Springer-Verlag, 2008.
279
P. J. Diggle, P. Moraga, B. Rowlingson, y B. M. Taylor. Spatial and spatio-
temporal log-Gaussian Cox processes: Extending the geostatistical paradigm.
Statistica Science, 28:542–563, 2013.
P. Doukhan, J. León, y P. Soulier. Central and non central limit theorems for
strongly dependent stationary Gaussian field. Rebrape, 10:205–223, 1996.
280
J. P. Fouque, G. Papanicolaou, y K. R. Sircar. Derivatives in Financial Markets
with Stochastic Volatility. Cambridge University Press, 2000.
R. Fox y M. S. Taqqu. Central limit theorems for quadratic forms in random va-
riables having long-range dependence. Probability Theory and Related Fields,
74:213–240, 1987.
281
A. Goia y P. Vieu. An introduction to recent advances in high/infinite dimensio-
nal statistics. Journal of Multivariate Analysis, 146:1–6, 2016.
282
T. Hastie, R. Tibshirani, y J. Friedman. The elements of statistical learning.
Springer-Verlag, New-York, 2001.
S. Hörmann y P. Kokoszka. Functional time series. En C.R. Rao, ed., Time Series
Analysis: Methods and Applications, tomo 30, págs. 157–186. Springer, 2012.
283
I. A. Ibragimov. On maximum likelihood estimation of parameters of the spec-
tral density of stationary time series. Theory of Probability and Its Applications,
12(1):115–119, 1967.
284
A. Jalilian, Y. Guan, y R. Waagepetersen. Orthogonal series estimation of the
pair correlation function of a spatial point process. Statistica Sinica, 29:769–
787, 2019.
C. Ji, D. Jiang, y N. Shi. The behavior of an sir epidemic model with stochastic
perturbation. Stochastic Analysis and Applications, 30:755–773, 2012.
285
P. Kokoszka, R. Gabrys, y L. Horváth. Tests for error correlation in the functional
linear model. Journal of the American Statistical Association, 105:1113–1125,
2010.
N. Leonenko. Limit Theorems for Random Fields with Singular Spectrum. Sprin-
ger Science + Business Media, B. V., 1999.
286
N. N. Leonenko y M. D. Ruiz-Medina. Scaling laws for the multidimensional
burgers equation with quadratic external potentials. Journal of Statistical
Physics, 124:191–205, 2006.
287
D. Marinucci y G. Peccati. Random fields on the sphere. Representation, limit
theorems and cosmological applications. London Mathematical Society Lecture
Note Series 389. Cambridge University Press, Cambridge, 2011.
C. C. McCluskey. Complete global stability for an sir epidemic model with delay
distributed or discrete. Nonlinear Analysis: Real World Applications, 11:55–59,
2010.
288
J. Møller y H. Toftaker. Geometric anisotropic spatial point pattern analysis and
cox processes. Scandinavian Journal of Statistics, 41(2):414–435, 2014.
289
V. M. Panaretos y S. Tavakoli. Fourier analysis of stationary time series in fun-
ction space. Annals of Statistics, 41:568–603, 2013b.
A. Remuzzi y G. Remuzzi. Covid-19 and italy: what next? The Lancet, 2020.
URL [Link]/10.1016/S0140-6736(20)30690-5.
290
A. C. Restrepo, P. Baker, y A. C. Clements. National spatial and temporal pat-
terns of notified dengue cases, colombia 2007–2010. Tropical Medicine &
International Health, 19:863–871, 2014. URL [Link]
tmi.12325.
291
M. D. Ruiz-Medina. Spatial autoregressive and moving average Hilbertian pro-
cesses. Journal of Multivariate Analysis, 102:292–305, 2011.
292
M. Sekiguchi y E. Ishiwata. Global dynamics of a discretized sirs epidemic
model with time delay. Journal of Mathematical Analysis and Applications,
371:195–202, 2010.
293
E. Tornatore, S. M. Buccellato, y P. Vetro. Stability of a stochastic sir system.
Physica A: Statistical Mechanics and its Applications, 354:111–126, 2005.
294
areas. Stochastic Environmental Research and Risk Assessment, 34:1421–1440,
2020.
295
H. Yu, J. M. Angulo, M. Cheng, J. Wu, y G. Christakos. An online spatio tempo-
ral prediction model for dengue fever epidemic in kaohsiung (taiwan). Bio-
metrical Journal, 56:428–440, 2014. URL [Link]
201200270.
J. Yu, D. Jiang, y N. Shi. Global stability of two-group sir model with random
perturbation. Journal of Mathematical Analysis and Applications, 360:235–
244, 2009.
W.-b. Zhang, Y. Ge, M. Liu, P. M. Atkinson, J. Wang, Zhang X., y Z. Tian. Risk
assessment of the step-by-step return-to-work policy in beijing following the
covid-19 epidemic peak. Stochastic Environmental Research and Risk Assess-
ment, 2020. URL [Link]
F. Zhou, T. Yu, R. Du, G. Fan, Y. Liu, Z. Liu, J. Xiang, Y. Wang, B. Song, X. Gu,
L. Guan, Y. Wei, H. Li, X. Wu, Ji. Xu, S. Tu, Y. Zhang, H. Chen, y B. Cao.
Clinical course and risk factors for mortality of adult inpatients with covid–
19 in wuhan, china: a retrospective cohort study. The Lancet, 2020. URL
[Link]/10.1016/S0140-6736(20)30566-3.
296