Datos de precisión en la evaluación del ruido transmitido al interior de una habitación según el Real Decreto 1367/2007

El objetivo del presente trabajo es realizar una evaluación realista de la incertidumbre en los ensayos de ruido transmitido al interior de una sala según el Anexo IV del Real Decreto 1367/2007, contemplando todas las posibles causas que intervienen en la incertidumbre de forma que las imprecisiones observadas en los ejercicios de intercomparación se cubran con un único ensayo. Otra parte fundamental del trabajo es cuantificar la fuente de incertidumbre generada en la elección del punto de medición. En primer lugar se realiza la medición del nivel continuo equivalente ponderado A en habitaciones de distinto tamaño, a la que llegan los distintos tipos de ruido emitido desde el exterior de las mismas. El siguiente paso es realizar un análisis de los niveles medidos, tanto en su distribución espacial como en su evolución temporal, incidiendo en los valores máximos, LAeq,5s. A continuación se comprueba si la distribución de los niveles medidos se ajusta a una distribución normal mediante el software de análisis estadístico STATGRAPHICS. Determinando, en base al tamaño de las muestras escogidas, por medio de análisis estadísticos, la aportación a la incertidumbre generada en la elección del punto de medida, cuantificando su valor. Por último se realiza una medición del ruido de actividades según el Real Decreto 1367/2007, aportando una evaluación de la incertidumbre, teniendo en cuenta todas las fuentes que la generan, mediante el enfoque clásico de la GUM. ABSTRACT. The objective of this work is to carry out a realistic assessment of the uncertainty in the trials of noise transmitted into the interior of a room according to Annex IV of the Royal Decree 136772997, considering all of the causes involved in uncertainty in such a way that the inaccuracies observed in intercomparison exercises are covered with a single trial. Another fundamental part of the work is to quantify the source of uncertainty in the choice of the point of measurement. First is the weighted equivalent of the continuous level measurement in rooms of different sizes, which reach different types of noise emitted from outside of them. The next step is to perform an analysis of the measured levels, both in their spatial distribution and their temporal evolution, influencing the maximum values, LAeq,5s. Then it is checked to see whether the distribution of the measured levels conforms to a normal distribution using the STATGRAPHICS statistical analysis software. Determining therefore, based on the size of the selected samples, through statistical analysis, the contribution to the uncertainty generated in the choice of the measurement point, quantifying its value. At last is the measure of the noise of activities according to the Royal Decree 1367/2007, providing an assessment of the uncertainty, taking into account all sources that generate it, using the classical approach to the GUM.

Desarrollo de metodologías para mejorar la estimación de los hidrogramas de diseño para el cálculo de los órganos de desagüe de las presas

En la presente Tesis se ha llevado a cabo el contraste y desarrollo de metodologías que permitan mejorar el cálculo de las avenidas de proyecto y extrema empleadas en el cálculo de la seguridad hidrológica de las presas. En primer lugar se ha abordado el tema del cálculo de las leyes de frecuencia de caudales máximos y su extrapolación a altos periodos de retorno. Esta cuestión es de gran relevancia, ya que la adopción de estándares de seguridad hidrológica para las presas cada vez más exigentes, implica la utilización de periodos de retorno de diseño muy elevados cuya estimación conlleva una gran incertidumbre. Es importante, en consecuencia incorporar al cálculo de los caudales de diseño todas la técnicas disponibles para reducir dicha incertidumbre. Asimismo, es importante hacer una buena selección del modelo estadístico (función de distribución y procedimiento de ajuste) de tal forma que se garantice tanto su capacidad para describir el comportamiento de la muestra, como para predecir de manera robusta los cuantiles de alto periodo de retorno. De esta forma, se han realizado estudios a escala nacional con el objetivo de determinar el esquema de regionalización que ofrece mejores resultados para las características hidrológicas de las cuencas españolas, respecto a los caudales máximos anuales, teniendo en cuenta el numero de datos disponibles. La metodología utilizada parte de la identificación de regiones homogéneas, cuyos límites se han determinado teniendo en cuenta las características fisiográficas y climáticas de las cuencas, y la variabilidad de sus estadísticos, comprobando posteriormente su homogeneidad. A continuación, se ha seleccionado el modelo estadístico de caudales máximos anuales con un mejor comportamiento en las distintas zonas de la España peninsular, tanto para describir los datos de la muestra como para extrapolar a los periodos de retorno más altos. El proceso de selección se ha basado, entre otras cosas, en la generación sintética de series de datos mediante simulaciones de Monte Carlo, y el análisis estadístico del conjunto de resultados obtenido a partir del ajuste de funciones de distribución a estas series bajo distintas hipótesis. Posteriormente, se ha abordado el tema de la relación caudal-volumen y la definición de los hidrogramas de diseño en base a la misma, cuestión que puede ser de gran importancia en el caso de presas con grandes volúmenes de embalse. Sin embargo, los procedimientos de cálculo hidrológico aplicados habitualmente no tienen en cuenta la dependencia estadística entre ambas variables. En esta Tesis se ha desarrollado un procedimiento para caracterizar dicha dependencia estadística de una manera sencilla y robusta, representando la función de distribución conjunta del caudal punta y el volumen en base a la función de distribución marginal del caudal punta y la función de distribución condicionada del volumen respecto al caudal. Esta última se determina mediante una función de distribución log-normal, aplicando un procedimiento de ajuste regional. Se propone su aplicación práctica a través de un procedimiento de cálculo probabilístico basado en la generación estocástica de un número elevado de hidrogramas. La aplicación a la seguridad hidrológica de las presas de este procedimiento requiere interpretar correctamente el concepto de periodo de retorno aplicado a variables hidrológicas bivariadas. Para ello, se realiza una propuesta de interpretación de dicho concepto. El periodo de retorno se entiende como el inverso de la probabilidad de superar un determinado nivel de embalse. Al relacionar este periodo de retorno con las variables hidrológicas, el hidrograma de diseño de la presa deja de ser un único hidrograma para convertirse en una familia de hidrogramas que generan un mismo nivel máximo en el embalse, representados mediante una curva en el plano caudal volumen. Esta familia de hidrogramas de diseño depende de la propia presa a diseñar, variando las curvas caudal-volumen en función, por ejemplo, del volumen de embalse o la longitud del aliviadero. El procedimiento propuesto se ilustra mediante su aplicación a dos casos de estudio. Finalmente, se ha abordado el tema del cálculo de las avenidas estacionales, cuestión fundamental a la hora de establecer la explotación de la presa, y que puede serlo también para estudiar la seguridad hidrológica de presas existentes. Sin embargo, el cálculo de estas avenidas es complejo y no está del todo claro hoy en día, y los procedimientos de cálculo habitualmente utilizados pueden presentar ciertos problemas. El cálculo en base al método estadístico de series parciales, o de máximos sobre un umbral, puede ser una alternativa válida que permite resolver esos problemas en aquellos casos en que la generación de las avenidas en las distintas estaciones se deba a un mismo tipo de evento. Se ha realizado un estudio con objeto de verificar si es adecuada en España la hipótesis de homogeneidad estadística de los datos de caudal de avenida correspondientes a distintas estaciones del año. Asimismo, se han analizado los periodos estacionales para los que es más apropiado realizar el estudio, cuestión de gran relevancia para garantizar que los resultados sean correctos, y se ha desarrollado un procedimiento sencillo para determinar el umbral de selección de los datos de tal manera que se garantice su independencia, una de las principales dificultades en la aplicación práctica de la técnica de las series parciales. Por otra parte, la aplicación practica de las leyes de frecuencia estacionales requiere interpretar correctamente el concepto de periodo de retorno para el caso estacional. Se propone un criterio para determinar los periodos de retorno estacionales de forma coherente con el periodo de retorno anual y con una distribución adecuada de la probabilidad entre las distintas estaciones. Por último, se expone un procedimiento para el cálculo de los caudales estacionales, ilustrándolo mediante su aplicación a un caso de estudio. The compare and develop of a methodology in order to improve the extreme flow estimation for dam hydrologic security has been developed. First, the work has been focused on the adjustment of maximum peak flows distribution functions from which to extrapolate values for high return periods. This has become a major issue as the adoption of stricter standards on dam hydrologic security involves estimation of high design return periods which entails great uncertainty. Accordingly, it is important to incorporate all available techniques for the estimation of design peak flows in order to reduce this uncertainty. Selection of the statistical model (distribution function and adjustment method) is also important since its ability to describe the sample and to make solid predictions for high return periods quantiles must be guaranteed. In order to provide practical application of previous methodologies, studies have been developed on a national scale with the aim of determining a regionalization scheme which features best results in terms of annual maximum peak flows for hydrologic characteristics of Spanish basins taking into account the length of available data. Applied methodology starts with the delimitation of regions taking into account basin’s physiographic and climatic characteristics and the variability of their statistical properties, and continues with their homogeneity testing. Then, a statistical model for maximum annual peak flows is selected with the best behaviour for the different regions in peninsular Spain in terms of describing sample data and making solid predictions for high return periods. This selection has been based, among others, on synthetic data series generation using Monte Carlo simulations and statistical analysis of results from distribution functions adjustment following different hypothesis. Secondly, the work has been focused on the analysis of the relationship between peak flow and volume and how to define design flood hydrographs based on this relationship which can be highly important for large volume reservoirs. However, commonly used hydrologic procedures do not take statistical dependence between these variables into account. A simple and sound method for statistical dependence characterization has been developed by the representation of a joint distribution function of maximum peak flow and volume which is based on marginal distribution function of peak flow and conditional distribution function of volume for a given peak flow. The last one is determined by a regional adjustment procedure of a log-normal distribution function. Practical application is proposed by a probabilistic estimation procedure based on stochastic generation of a large number of hydrographs. The use of this procedure for dam hydrologic security requires a proper interpretation of the return period concept applied to bivariate hydrologic data. A standard is proposed in which it is understood as the inverse of the probability of exceeding a determined reservoir level. When relating return period and hydrological variables the only design flood hydrograph changes into a family of hydrographs which generate the same maximum reservoir level and that are represented by a curve in the peak flow-volume two-dimensional space. This family of design flood hydrographs depends on the dam characteristics as for example reservoir volume or spillway length. Two study cases illustrate the application of the developed methodology. Finally, the work has been focused on the calculation of seasonal floods which are essential when determining the reservoir operation and which can be also fundamental in terms of analysing the hydrologic security of existing reservoirs. However, seasonal flood calculation is complex and nowadays it is not totally clear. Calculation procedures commonly used may present certain problems. Statistical partial duration series, or peaks over threshold method, can be an alternative approach for their calculation that allow to solve problems encountered when the same type of event is responsible of floods in different seasons. A study has been developed to verify the hypothesis of statistical homogeneity of peak flows for different seasons in Spain. Appropriate seasonal periods have been analyzed which is highly relevant to guarantee correct results. In addition, a simple procedure has been defined to determine data selection threshold on a way that ensures its independency which is one of the main difficulties in practical application of partial series. Moreover, practical application of seasonal frequency laws requires a correct interpretation of the concept of seasonal return period. A standard is proposed in order to determine seasonal return periods coherently with the annual return period and with an adequate seasonal probability distribution. Finally a methodology is proposed to calculate seasonal peak flows. A study case illustrates the application of the proposed methodology.

Pulmonary malformation in transgenic mice expressing human keratinocyte growth factor in the lung.

Expression of human keratinocyte growth factor (KGF/FGF-7) was directed to epithelial cells of the developing embryonic lung of transgenic mice disrupting normal pulmonary morphogenesis during the pseudoglandular stage of development. By embryonic day 15.5(E15.5), lungs of transgenic surfactant protein C (SP-C)-KGF mice resembled those of humans with pulmonary cystadenoma. Lungs were cystic, filling the thoracic cavity, and were composed of numerous dilated saccules lined with glycogen-containing columnar epithelial cells. The normal distribution of SP-C proprotein in the distal regions of respiratory tubules was disrupted. Columnar epithelial cells lining the papillary structures stained variably and weakly for this distal respiratory cell marker. Mesenchymal components were preserved in the transgenic mouse lungs, yet the architectural relationship of the epithelium to the mesenchyme was altered. SP-C-KGF transgenic mice failed to survive gestation to term, dying before E17.5. Culturing mouse fetal lung explants in the presence of recombinant human KGF also disrupted branching morphogenesis and resulted in similar cystic malformation of the lung. Thus, it appears that precise temporal and spatial expression of KGF is likely to play a crucial role in the control of branching morphogenesis during fetal lung development.

Caracterização da anemia em cadelas com piometra

A piometra é uma condição mórbida caracterizada pela inflamação do útero com acúmulo de exsudatos, resultante de ações hormonais e geralmente associada à presença de bactérias no lúmen uterino. A anemia é a alteração hematológica mais frequentemente observada em cadelas com piometra e está associada à cronicidade da doença, diminuição da eritropoiese, devido ao efeito toxêmico na medula óssea, diminuição da disponibilidade de ferro ou perda de sangue para o útero. Adicionalmente, o efeito das toxinas bacterianas e os radicais livres gerados pelo metabolismo oxidativo dos neutrófilos podem resultar na modificação da estrutura antigênica da membrana do eritrócito, permitindo a ligação de imunoglobulinas em sua superfície e acelerando a destruição eritrocitária. Essa hipótese pode ser comprovada pela detecção de imunocomplexos na superfície eritrocitária de cadelas com piometra. O diagnóstico de piometra foi estabelecido em 33 cadelas atendidas no Serviço de Obstetrícia/Ginecologia do Hospital Veterinário da Universidade de São Paulo com base na anamnese, exame físico e exames subsidiários (ultrassonografia, hemograma e concentrações séricas de ureia e creatinina). As amostras sanguíneas foram coletadas em dois momentos. A primeira anterior a ovariosalpingohisterectomia (OSH) e a segunda, sete a dez dias após a OSH. A quantificação de hemácias com deposição de imunocomplexos IgG e IgM foi realizada utilizando-se anticorpos anti-IgG e anti-IgM (Bethyl®Laboratories) conjugadas a fluoresceína de isotiocianato (FITC), e a leitura realizada com citômetro de fluxo (FACS Calibur; Becton, Dickinson and Company© 2007 BD), sendo os resultados expressos em percentual de hemácias marcadas. Foram utilizados o Teste de Shapiro-Wilk para a avaliação da distribuição de dados e a comparação entre os grupos controle, pré e pós-OSH foi realizada valendo-se do Teste t ou Teste t pareado e Correlação de Pearson, e do Teste U de Mann-Whitney e Correlação de Spearman, para as variáveis com distribuição normal e não-normal, respectivamente. O valor de alfa estipulado foi de 0,05. Analisando os valores hematológicos de cada um dos cães incluídos no estudo, observa-se que 19 (57,6%) apresentavam anemia normocítica normocrômica não regenerativa no momento pré-OSH e cinco (15,2%) no momento pós-OSH. Em cães do grupo controle foram observadas 0,14 - 0,77% (0,43±0,18%) de hemácias marcadas com anticorpos anti-IgG FITC e 0,29 - 9,58% (0,68±0,29%) para anticorpos anti-IgM FITC. Já nos cães com piometra, foram encontradas 0,14 - 4,19% (0,96±0,86%) de hemácias marcadas com anticorpos anti-IgG FITC e 0,29 - 9,58% (1,37±1,71%) com anticorpos anti-IgM FITC, antecedendo a OSH. No momento pós-OSH observou-se 0,18 - 16,2% (2,77±3,67%) de hemácias marcadas para anticorpos anti-IgG FITC e 0,15 - 19,8% (4,01±4,46%) para anticorpos anti-IgM FITC. O percentual de hemácias marcadas com anticorpos anti-IgG FITC diferiu entre os grupos controle e piometra, pré-OSH (p<0,001) e pós-OSH (p<0,001). Em relação a anticorpos anti-IgM FITC, não foram observadas diferenças entre os grupos controle e pré-OSH (p=0,09), porém, após a OSH houve aumento na marcação de hemácias, quando comparado ao grupo controle (p<0,001). Apenas alguns animais apresentaram mais de 5% de hemácias marcadas, e isto ocorreu, principalmente, no momento pós-OSH. Entretanto, não resultou no agravamento da anemia, indicando que a piometra em cadelas está associada à deposição de imunoglobulinas G ou M na superfície das hemácias, sem, no entanto, promover hemólise ou agravamento da anemia

Avaliação em 3D da cabeça da mandíbula em crianças portadoras de mordida cruzada posterior unilateral por meio de tomografia computadorizada de feixe cônico

Durante o desenvolvimento da oclusão, a instalação de maloclusões podem resultar em desarmonias dento faciais de natureza e severidade diversas, podendo provocar alterações no desenvolvimento crânio facial, dentre as estruturas envolvidas as Articulações Temporo Mandibulares (ATM), podem sofrer alguma influência, dessa forma a avaliação desta região, no aspecto morfológico e funcional, constituí tema de interesse, sempre que levados em conta os aspectos funcionais da oclusão. A relação entre a forma e a função, tanto das cabeças da mandíbula, bem como o contorno da fossa mandibular com as maloclusões ainda é controversa e não está compreendida por completo, porém a literatura sobre o assunto, demonstra correlação entre a instalação de maloclusões e modificações neste sistema, mesmo que algumas alterações não sejam de ordem estatística e em amostras de indivíduos em tenra idade, as mesmas podem comprometer o desenvolvimento adequado em indivíduos adultos ou mesmo adultos jovens. Tendo como propósito nesse estudo a avaliação das cabeças da mandíbula quanto ao volume e superficíe dos lados direito e esquerdo, cruzado e não cruzado, a amostra selecionada foi de 20 indivíduos com mordida cruzada posterior unilateral, com idades entre 06 e 09 anos de idade, utilizando imagens de tomografia computadorizada por feixe cônico, imagens obtidas por um equipamento modelo i- Cat, sendo utilizado na reformatação e manipulação das imagens o programa computacional - NemoCeph 3D® versão 11.5. Nas medições propostas para esse estudo, utilizou-se o teste t pareado de Student para amostras com distribuição normal. Na observação das tabelas e seus respectivos gráficos, podemos verificar que na comparação entre os lados direito e esquerdo, e cruzado e não cruzado das cabeças da mandíbula, com relação ao volume e superfície, existem diferenças numéricas entre elas, porém não pode ser observado diferenças estatísticas significantes, nessa amostra especifica com a metodologia empregada para esse estudo. Assim foi possível concluir que nas Mordidas Cruzadas Posteriores Unilaterais as cabeças da mandíbula tanto em seu volume como em sua superfície não apresentaram diferenças estatisticamente significantes na amostra estudada.

Análise biomecânica de resistência à fratura em dentes caninos artificiais de cães, com restauração metálica fundida e retentores intrarradiculares retos e curvos com núcleo - Estudo comparativo com análise por elementos finitos

Os cães, por fatores diversos, acabam por apresentar dentes fraturados com ou sem exposição de polpa. Estas fraturas basicamente são identificadas como fraturas recuperáveis não complicadas, recuperáveis complicadas ou irrecuperáveis. As fraturas recuperáveis (localizadas apenas no esmalte e dentina) são tratadas com dentística restauradora. As recuperáveis complicadas (com lesões em esmalte, dentina e exposição do canal radicular) passam por tratamento endodôntico, podendo ser seguidas de restaurações metálicas. Os dentes mais comumente acometidos são os dentes caninos, superiores ou inferiores. Este trabalho em dentes artificiais simulando considerável destruição de sua porção coronal objetivou testar, após a adaptação da restauração metálica fundida, a resistência às fraturas no dente canino. Os dentes artificiais foram padronizados com uma técnica de replicação de raízes artificiais em molde de resina acrílica quimicamente ativada. Oitenta réplicas iguais de resina composta fotopolimerizável, padronizadas em tamanho e forma, foram construídas a partir desta técnica. Antes da reconstrução protética, aplicou-se o tratamento endodôntico, desobturação, preparo do canal radicular e moldagem. Proteticamente, um pino intrarradicular reto e outro curvo, ambos com núcleo para sustentar a coroa metálica fundida foram cimentados na porção coronal de cada raiz-réplica. Os núcleos e coroa metálica foram ambos ferulados ou estojados. Avaliou-se os dois tipos de restauração com pino intrarradicular curvos ou retos cimentados com cimento de fosfato de zinco ou resinoso para identificar o melhor conjunto restaurador. Os testes de resistência biomecânica de 80 raízes-réplicas foram divididos em 4 grupos com 20 corpos de prova para cada um dos grupos. Grupo 1: das raízes-réplicas com pino intrarradicular curvo cimentados com cimento resinoso. Grupo 2: das raízes-réplicas com pino intrarradicular curvo cimentados com cimento de fosfato de zinco. Grupo 3: das raízes-réplicas com pino intrarradicular reto cimentados com cimento resinoso. Grupo 4: das raízes-réplicas com pino intrarradicular reto cimentados com cimento de fosfato de zinco. Estes grupos foram submetidos a teste de força com pré-carga de 1,5 N, com velocidade de avanço constante de 0,05 mm por minuto em ponto pré- determinado (mésio-lateral vestibularizada) até ocorrência de fratura do conjunto ou parte dele em uma Máquina Universal Kratos. Com a avaliação biomecânica e estudo estatístico de Kruskall-Wallis, identificou-se que os dados obtidos não seguiram distribuição normal. Esta diferença mostrou-se com o p<0,05 na interpretação do teste. No caso de dados não paramétricos o post-hoc do Kruskal-Wallis foi o teste de U de Mann-Withney. Paralelamente, um estudo com análise de elementos finitos comparou os resultados obtidos. Não houve diferença significativa sobre o tipo de cimento utilizado ou que favorecesse o uso do pino reto ou do pino curvo, recaindo a escolha para o operador decidir de acordo com a melhor indicação para cada caso clínico

Sur les tests lisses d'ajustement dans le context des series chronologiques

La plupart des modèles en statistique classique repose sur une hypothèse sur la distribution des données ou sur une distribution sous-jacente aux données. La validité de cette hypothèse permet de faire de l’inférence, de construire des intervalles de confiance ou encore de tester la fiabilité du modèle. La problématique des tests d’ajustement vise à s’assurer de la conformité ou de la cohérence de l’hypothèse avec les données disponibles. Dans la présente thèse, nous proposons des tests d’ajustement à la loi normale dans le cadre des séries chronologiques univariées et vectorielles. Nous nous sommes limités à une classe de séries chronologiques linéaires, à savoir les modèles autorégressifs à moyenne mobile (ARMA ou VARMA dans le cas vectoriel). Dans un premier temps, au cas univarié, nous proposons une généralisation du travail de Ducharme et Lafaye de Micheaux (2004) dans le cas où la moyenne est inconnue et estimée. Nous avons estimé les paramètres par une méthode rarement utilisée dans la littérature et pourtant asymptotiquement efficace. En effet, nous avons rigoureusement montré que l’estimateur proposé par Brockwell et Davis (1991, section 10.8) converge presque sûrement vers la vraie valeur inconnue du paramètre. De plus, nous fournissons une preuve rigoureuse de l’inversibilité de la matrice des variances et des covariances de la statistique de test à partir de certaines propriétés d’algèbre linéaire. Le résultat s’applique aussi au cas où la moyenne est supposée connue et égale à zéro. Enfin, nous proposons une méthode de sélection de la dimension de la famille d’alternatives de type AIC, et nous étudions les propriétés asymptotiques de cette méthode. L’outil proposé ici est basé sur une famille spécifique de polynômes orthogonaux, à savoir les polynômes de Legendre. Dans un second temps, dans le cas vectoriel, nous proposons un test d’ajustement pour les modèles autorégressifs à moyenne mobile avec une paramétrisation structurée. La paramétrisation structurée permet de réduire le nombre élevé de paramètres dans ces modèles ou encore de tenir compte de certaines contraintes particulières. Ce projet inclut le cas standard d’absence de paramétrisation. Le test que nous proposons s’applique à une famille quelconque de fonctions orthogonales. Nous illustrons cela dans le cas particulier des polynômes de Legendre et d’Hermite. Dans le cas particulier des polynômes d’Hermite, nous montrons que le test obtenu est invariant aux transformations affines et qu’il est en fait une généralisation de nombreux tests existants dans la littérature. Ce projet peut être vu comme une généralisation du premier dans trois directions, notamment le passage de l’univarié au multivarié ; le choix d’une famille quelconque de fonctions orthogonales ; et enfin la possibilité de spécifier des relations ou des contraintes dans la formulation VARMA. Nous avons procédé dans chacun des projets à une étude de simulation afin d’évaluer le niveau et la puissance des tests proposés ainsi que de les comparer aux tests existants. De plus des applications aux données réelles sont fournies. Nous avons appliqué les tests à la prévision de la température moyenne annuelle du globe terrestre (univarié), ainsi qu’aux données relatives au marché du travail canadien (bivarié). Ces travaux ont été exposés à plusieurs congrès (voir par exemple Tagne, Duchesne et Lafaye de Micheaux (2013a, 2013b, 2014) pour plus de détails). Un article basé sur le premier projet est également soumis dans une revue avec comité de lecture (Voir Duchesne, Lafaye de Micheaux et Tagne (2016)).

Sur les tests lisses d'ajustement dans le context des series chronologiques

La plupart des modèles en statistique classique repose sur une hypothèse sur la distribution des données ou sur une distribution sous-jacente aux données. La validité de cette hypothèse permet de faire de l’inférence, de construire des intervalles de confiance ou encore de tester la fiabilité du modèle. La problématique des tests d’ajustement vise à s’assurer de la conformité ou de la cohérence de l’hypothèse avec les données disponibles. Dans la présente thèse, nous proposons des tests d’ajustement à la loi normale dans le cadre des séries chronologiques univariées et vectorielles. Nous nous sommes limités à une classe de séries chronologiques linéaires, à savoir les modèles autorégressifs à moyenne mobile (ARMA ou VARMA dans le cas vectoriel). Dans un premier temps, au cas univarié, nous proposons une généralisation du travail de Ducharme et Lafaye de Micheaux (2004) dans le cas où la moyenne est inconnue et estimée. Nous avons estimé les paramètres par une méthode rarement utilisée dans la littérature et pourtant asymptotiquement efficace. En effet, nous avons rigoureusement montré que l’estimateur proposé par Brockwell et Davis (1991, section 10.8) converge presque sûrement vers la vraie valeur inconnue du paramètre. De plus, nous fournissons une preuve rigoureuse de l’inversibilité de la matrice des variances et des covariances de la statistique de test à partir de certaines propriétés d’algèbre linéaire. Le résultat s’applique aussi au cas où la moyenne est supposée connue et égale à zéro. Enfin, nous proposons une méthode de sélection de la dimension de la famille d’alternatives de type AIC, et nous étudions les propriétés asymptotiques de cette méthode. L’outil proposé ici est basé sur une famille spécifique de polynômes orthogonaux, à savoir les polynômes de Legendre. Dans un second temps, dans le cas vectoriel, nous proposons un test d’ajustement pour les modèles autorégressifs à moyenne mobile avec une paramétrisation structurée. La paramétrisation structurée permet de réduire le nombre élevé de paramètres dans ces modèles ou encore de tenir compte de certaines contraintes particulières. Ce projet inclut le cas standard d’absence de paramétrisation. Le test que nous proposons s’applique à une famille quelconque de fonctions orthogonales. Nous illustrons cela dans le cas particulier des polynômes de Legendre et d’Hermite. Dans le cas particulier des polynômes d’Hermite, nous montrons que le test obtenu est invariant aux transformations affines et qu’il est en fait une généralisation de nombreux tests existants dans la littérature. Ce projet peut être vu comme une généralisation du premier dans trois directions, notamment le passage de l’univarié au multivarié ; le choix d’une famille quelconque de fonctions orthogonales ; et enfin la possibilité de spécifier des relations ou des contraintes dans la formulation VARMA. Nous avons procédé dans chacun des projets à une étude de simulation afin d’évaluer le niveau et la puissance des tests proposés ainsi que de les comparer aux tests existants. De plus des applications aux données réelles sont fournies. Nous avons appliqué les tests à la prévision de la température moyenne annuelle du globe terrestre (univarié), ainsi qu’aux données relatives au marché du travail canadien (bivarié). Ces travaux ont été exposés à plusieurs congrès (voir par exemple Tagne, Duchesne et Lafaye de Micheaux (2013a, 2013b, 2014) pour plus de détails). Un article basé sur le premier projet est également soumis dans une revue avec comité de lecture (Voir Duchesne, Lafaye de Micheaux et Tagne (2016)).

(Table 1) Mean grain sizes and standard deviations of gas hydrate samples

Methane hydrates are present in marine seep systems and occur within the gas hydrate stability zone. Very little is known about their crystallite sizes and size distributions because they are notoriously difficult to measure. Crystal size distributions are usually considered as one of the key petrophysical parameters because they influence mechanical properties and possible compositional changes, which may occur with changing environmental conditions. Variations in grain size are relevant for gas substitution in natural hydrates by replacing CH4 with CO2 for the purpose of carbon dioxide sequestration. Here we show that crystallite sizes of gas hydrates from some locations in the Indian Ocean, Gulf of Mexico and Black Sea are in the range of 200-400 µm; larger values were obtained for deeper-buried samples from ODP Leg 204. The crystallite sizes show generally a log-normal distribution and appear to vary sometimes rapidly with location.

Variability in Modified Estimators of VaR and ES

Thesis (Master's)--University of Washington, 2016-06

Copulas in QTL mapping

The standard variance components method for mapping quantitative trait loci is derived on the assumption of normality. Unsurprisingly, statistical tests based on this method do not perform so well if this assumption is not satisfied. We use the statistical concept of copulas to relax the assumption of normality and derive a test that can perform well under any distribution of the continuous trait. In particular, we discuss bivariate normal copulas in the context of sib-pair studies. Our approach is illustrated by a linkage analysis of lipoprotein(a) levels, whose distribution is highly skewed. We demonstrate that the asymptotic critical levels of the test can still be calculated using the interval mapping approach. The new method can be extended to more general pedigrees and multivariate phenotypes in a similar way as the original variance components method.

Nutrient partitioning during treatment of tuberculosis: gain in body fat mass but not in protein mass

Background: Tuberculosis is an important cause of wasting. The functional consequences of wasting and recovery may depend on the distribution of lost and gained nutrient stores between protein and fat masses. Objective: The goal was to study nutrient partitioning, ie, the proportion of weight change attributable to changes in fat mass (FM) versus protein mass (PM), during anti mycobacterial treatment. Design: Body-composition measures were made of 21 men and 9 women with pulmonary tuberculosis at baseline and after 1 and 6 mo of treatment. All subjects underwent dual-energy X-ray absorptiometry and deuterium bromide dilution tests, and a four-compartment model of FM, total body water (TBW), bone minerals (BM), and PM was derived. The ratio of PM to FM at any time was expressed as the energy content (p-ratio). Changes in the p-ratio were related to disease severity as measured by radiologic criteria. Results: Patients gained 10% in body weight (P < 0.001) from baseline to month 6. This was mainly due to a 44% gain in FM (P < 0.001); PM, BM, and TBW did not change significantly. Results were similar in men and women. The p-ratio decreased from baseline to month 1 and then fell further by month 6. Radiologic disease severity was not correlated with changes in the p-ratio. Conclusions: Microbiological cure of tuberculosis does not restore PM within 6 mo, despite a strong anabolic response. Change in the p-ratio is a suitable parameter for use in studying the effect of disease on body composition because it allows transformation of such effects into a normal distribution across a wide range of baseline proportion between fat and protein mass.

Towards a framework for combining stochastic and deterministic descriptions of nonstationary financial time series

We present in this paper ideas to tackle the problem of analysing and forecasting nonstationary time series within the financial domain. Accepting the stochastic nature of the underlying data generator we assume that the evolution of the generator's parameters is restricted on a deterministic manifold. Therefore we propose methods for determining the characteristics of the time-localised distribution. Starting with the assumption of a static normal distribution we refine this hypothesis according to the empirical results obtained with the methods anc conclude with the indication of a dynamic non-Gaussian behaviour with varying dependency for the time series under consideration.

Non-normally distributed data and non-parametric statistics

Different types of numerical data can be collected in a scientific investigation and the choice of statistical analysis will often depend on the distribution of the data. A basic distinction between variables is whether they are ‘parametric’ or ‘non-parametric’. When a variable is parametric, the data come from a symmetrically shaped distribution known as the ‘Gaussian’ or ‘normal distribution’ whereas non-parametric variables may have a distribution which deviates markedly in shape from normal. This article describes several aspects of the problem of non-normality including: (1) how to test for two common types of deviation from a normal distribution, viz., ‘skew’ and ‘kurtosis’, (2) how to fit the normal distribution to a sample of data, (3) the transformation of non-normally distributed data and scores, and (4) commonly used ‘non-parametric’ statistics which can be used in a variety of circumstances.

Statnote 14: the correlation of two variables (Pearson's 'r')

Pearson's correlation coefficient (‘r’) is one of the most widely used of all statistics. Nevertheless, care needs to be used in interpreting the results because with large numbers of observations, quite small values of ‘r’ become significant and the X variable may only account for a small proportion of the variance in Y. Hence, ‘r squared’ should always be calculated and included in a discussion of the significance of ‘r’. The use of ‘r’ also assumes that the data follow a bivariate normal distribution (see Statnote 17) and this assumption should be examined prior to the study. If the data do not conform to such a distribution, the use of a non-parametric correlation coefficient should be considered. A significant correlation should not be interpreted as indicating ‘causation’ especially in observational studies, in which the two variables may be correlated because of their mutual correlations with other confounding variables.