Flood frequency analysis by a bivariate model based on copulas

La adecuada estimación de avenidas de diseño asociadas a altos periodos de retorno es necesaria para el diseño y gestión de estructuras hidráulicas como presas. En la práctica, la estimación de estos cuantiles se realiza normalmente a través de análisis de frecuencia univariados, basados en su mayoría en el estudio de caudales punta. Sin embargo, la naturaleza de las avenidas es multivariada, siendo esencial tener en cuenta características representativas de las avenidas, tales como caudal punta, volumen y duración del hidrograma, con el fin de llevar a cabo un análisis apropiado; especialmente cuando el caudal de entrada se transforma en un caudal de salida diferente durante el proceso de laminación en un embalse o llanura de inundación. Los análisis de frecuencia de avenidas multivariados han sido tradicionalmente llevados a cabo mediante el uso de distribuciones bivariadas estándar con el fin de modelar variables correlacionadas. Sin embargo, su uso conlleva limitaciones como la necesidad de usar el mismo tipo de distribuciones marginales para todas las variables y la existencia de una relación de dependencia lineal entre ellas. Recientemente, el uso de cópulas se ha extendido en hidrología debido a sus beneficios en relación al contexto multivariado, permitiendo superar los inconvenientes de las técnicas tradicionales. Una copula es una función que representa la estructura de dependencia de las variables de estudio, y permite obtener la distribución de frecuencia multivariada de dichas variables mediante sus distribuciones marginales, sin importar el tipo de distribución marginal utilizada. La estimación de periodos de retorno multivariados, y por lo tanto, de cuantiles multivariados, también se facilita debido a la manera en la que las cópulas están formuladas. La presente tesis doctoral busca proporcionar metodologías que mejoren las técnicas tradicionales usadas por profesionales para estimar cuantiles de avenida más adecuados para el diseño y la gestión de presas, así como para la evaluación del riesgo de avenida, mediante análisis de frecuencia de avenidas bivariados basados en cópulas. Las variables consideradas para ello son el caudal punta y el volumen del hidrograma. Con el objetivo de llevar a cabo un estudio completo, la presente investigación abarca: (i) el análisis de frecuencia de avenidas local bivariado centrado en examinar y comparar los periodos de retorno teóricos basados en la probabilidad natural de ocurrencia de una avenida, con el periodo de retorno asociado al riesgo de sobrevertido de la presa bajo análisis, con el fin de proporcionar cuantiles en una estación de aforo determinada; (ii) la extensión del enfoque local al regional, proporcionando un procedimiento completo para llevar a cabo un análisis de frecuencia de avenidas regional bivariado para proporcionar cuantiles en estaciones sin aforar o para mejorar la estimación de dichos cuantiles en estaciones aforadas; (iii) el uso de cópulas para investigar tendencias bivariadas en avenidas debido al aumento de los niveles de urbanización en una cuenca; y (iv) la extensión de series de avenida observadas mediante la combinación de los beneficios de un modelo basado en cópulas y de un modelo hidrometeorológico. Accurate design flood estimates associated with high return periods are necessary to design and manage hydraulic structures such as dams. In practice, the estimate of such quantiles is usually done via univariate flood frequency analyses, mostly based on the study of peak flows. Nevertheless, the nature of floods is multivariate, being essential to consider representative flood characteristics, such as flood peak, hydrograph volume and hydrograph duration to carry out an appropriate analysis; especially when the inflow peak is transformed into a different outflow peak during the routing process in a reservoir or floodplain. Multivariate flood frequency analyses have been traditionally performed by using standard bivariate distributions to model correlated variables, yet they entail some shortcomings such as the need of using the same kind of marginal distribution for all variables and the assumption of a linear dependence relation between them. Recently, the use of copulas has been extended in hydrology because of their benefits regarding dealing with the multivariate context, as they overcome the drawbacks of the traditional approach. A copula is a function that represents the dependence structure of the studied variables, and allows obtaining the multivariate frequency distribution of them by using their marginal distributions, regardless of the kind of marginal distributions considered. The estimate of multivariate return periods, and therefore multivariate quantiles, is also facilitated by the way in which copulas are formulated. The present doctoral thesis seeks to provide methodologies that improve traditional techniques used by practitioners, in order to estimate more appropriate flood quantiles for dam design, dam management and flood risk assessment, through bivariate flood frequency analyses based on the copula approach. The flood variables considered for that goal are peak flow and hydrograph volume. In order to accomplish a complete study, the present research addresses: (i) a bivariate local flood frequency analysis focused on examining and comparing theoretical return periods based on the natural probability of occurrence of a flood, with the return period associated with the risk of dam overtopping, to estimate quantiles at a given gauged site; (ii) the extension of the local to the regional approach, supplying a complete procedure for performing a bivariate regional flood frequency analysis to either estimate quantiles at ungauged sites or improve at-site estimates at gauged sites; (iii) the use of copulas to investigate bivariate flood trends due to increasing urbanisation levels in a catchment; and (iv) the extension of observed flood series by combining the benefits of a copula-based model and a hydro-meteorological model.

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.

Análisis, desarrollo y optimización de un sistema para el diseño de redes topográficas valorando el modelo digital del terreno

El presente trabajo tiene como objetivo general el análisis de las técnicas de diseño y optimización de redes topográficas, observadas mediante topografía convencional (no satelital) el desarrollo e implementación de un sistema informático capaz de ayudar a la definición de la geometría más fiable y precisa, en función de la orografía del terreno donde se tenga que ubicar. En primer lugar se realizará un estudio de la metodología del ajuste mediante mínimos cuadrados y la propagación de varianzas, para posteriormente analizar su dependencia de la geometría que adopte la red. Será imprescindible determinar la independencia de la matriz de redundancia (R) de las observaciones y su total dependencia de la geometría, así como la influencia de su diagonal principal (rii), números de redundancia, para garantizar la máxima fiabilidad interna de la misma. También se analizará el comportamiento de los números de redundancia (rii) en el diseño de una red topográfica, la variación de dichos valores en función de la geometría, analizando su independencia respecto de las observaciones así como los diferentes niveles de diseño en función de los parámetros y datos conocidos. Ha de señalarse que la optimización de la red, con arreglo a los criterios expuestos, está sujeta a los condicionantes que impone la necesidad de que los vértices sean accesibles, y además sean visibles entre sí, aquellos relacionados por observaciones, situaciones que dependen esencialmente del relieve del terreno y de los obstáculos naturales o artificiales que puedan existir. Esto implica la necesidad de incluir en el análisis y en el diseño, cuando menos de un modelo digital del terreno (MDT), aunque lo más útil sería la inclusión en el estudio del modelo digital de superficie (MDS), pero esta opción no siempre será posible. Aunque el tratamiento del diseño esté basado en un sistema bidimensional se estudiará la posibilidad de incorporar un modelo digital de superficie (MDS); esto permitirá a la hora de diseñar el emplazamiento de los vértices de la red la viabilidad de las observaciones en función de la orografía y los elementos, tanto naturales como artificiales, que sobre ella estén ubicados. Este sistema proporcionaría, en un principio, un diseño óptimo de una red constreñida, atendiendo a la fiabilidad interna y a la precisión final de sus vértices, teniendo en cuenta la orografía, lo que equivaldría a resolver un planteamiento de diseño en dos dimensiones y media1; siempre y cuando se dispusiera de un modelo digital de superficie o del terreno. Dado que la disponibilidad de obtener de manera libre el MDS de las zonas de interés del proyecto, hoy en día es costoso2, se planteará la posibilidad de conjuntar, para el estudio del diseño de la red, de un modelo digital del terreno. Las actividades a desarrollar en el trabajo de esta tesis se describen en esta memoria y se enmarcan dentro de la investigación para la que se plantean los siguientes objetivos globales: 1. Establecer un modelo matemático del proceso de observación de una red topográfica, atendiendo a todos los factores que intervienen en el mismo y a su influencia sobre las estimaciones de las incógnitas que se obtienen como resultado del ajuste de las observaciones. 2. Desarrollar un sistema que permita optimizar una red topográfica en sus resultados, aplicando técnicas de diseño y simulación sobre el modelo anterior. 3. Presentar una formulación explícita y rigurosa de los parámetros que valoran la fiabilidad de una red topográfica y de sus relaciones con el diseño de la misma. El logro de este objetivo se basa, además de en la búsqueda y revisión de las fuentes, en una intensa labor de unificación de notaciones y de construcción de pasos intermedios en los desarrollos matemáticos. 4. Elaborar una visión conjunta de la influencia del diseño de una red, en los seis siguientes factores (precisiones a posteriori, fiabilidad de las observaciones, naturaleza y viabilidad de las mismas, instrumental y metodología de estacionamiento) como criterios de optimización, con la finalidad de enmarcar el tema concreto que aquí se aborda. 5. Elaborar y programar los algoritmos necesarios para poder desarrollar una aplicación que sea capaz de contemplar las variables planteadas en el apartado anterior en el problema del diseño y simulación de redes topográficas, contemplando el modelo digital de superficie. Podrían considerarse como objetivos secundarios, los siguientes apartados: Desarrollar los algoritmos necesarios para interrelacionar el modelo digital del terreno con los propios del diseño. Implementar en la aplicación informática la posibilidad de variación, por parte del usuario, de los criterios de cobertura de los parámetros (distribución normal o t de Student), así como los grados de fiabilidad de los mismos ABSTRACT The overall purpose of this work is the analysis of the techniques of design and optimization for geodetic networks, measured with conventional survey methods (not satellite), the development and implementation of a computational system capable to help on the definition of the most liable and accurate geometry, depending on the land orography where the network has to be located. First of all, a study of the methodology by least squares adjustment and propagation of variances will be held; then, subsequently, analyze its dependency of the geometry that the network will take. It will be essential to determine the independency of redundancy matrix (R) from the observations and its absolute dependency from the network geometry, as well as the influence of the diagonal terms of the R matrix (rii), redundancy numbers, in order to ensure maximum re liability of the network. It will also be analyzed first the behavior of redundancy numbers (rii) in surveying network design, then the variation of these values depending on the geometry with the analysis of its independency from the observations, and finally the different design levels depending on parameters and known data. It should be stated that network optimization, according to exposed criteria, is subject to the accessibility of the network points. In addition, common visibility among network points, which of them are connected with observations, has to be considered. All these situations depends essentially on the terrain relief and the natural or artificial obstacles that should exist. Therefore, it is necessary to include, at least, a digital terrain model (DTM), and better a digital surface model (DSM), not always available. Although design treatment is based on a bidimensional system, the possibility of incorporating a digital surface model (DSM) will be studied; this will allow evaluating the observations feasibility based on the terrain and the elements, both natural and artificial, which are located on it, when selecting network point locations. This system would provide, at first, an optimal design of a constrained network, considering both the internal reliability and the accuracy of its points (including the relief). This approach would amount to solving a “two and a half dimensional”3 design, if a digital surface model is available. As the availability of free DSM4 of the areas of interest of the project today is expensive, the possibility of combining a digital terrain model will arise. The activities to be developed on this PhD thesis are described in this document and are part of the research for which the following overall objectives are posed: 1. To establish a mathematical model for the process of observation of a survey network, considering all the factors involved and its influence on the estimates of the unknowns that are obtained as a result of the observations adjustment. 2. To develop a system to optimize a survey network results, applying design and simulation techniques on the previous model. 3. To present an explicit and rigorous formulation of parameters which assess the reliability of a survey network and its relations with the design. The achievement of this objective is based, besides on the search and review of sources, in an intense work of unification of notation and construction of intermediate steps in the mathematical developments. 4. To develop an overview of the influence on the network design of six major factors (posterior accuracy, observations reliability, viability of observations, instruments and station methodology) as optimization criteria, in order to define the subject approached on this document. 5. To elaborate and program the algorithms needed to develop an application software capable of considering the variables proposed in the previous section, on the problem of design and simulation of surveying networks, considering the digital surface model. It could be considered as secondary objectives, the following paragraphs: To develop the necessary algorithms to interrelate the digital terrain model with the design ones. To implement in the software application the possibility of variation of the coverage criteria parameters (normal distribution or Student t test) and therefore its degree of reliability.

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

Hashimoto’s thyroiditis follow up along 10 years

Trabalho Final do Curso de Mestrado Integrado em Medicina, Faculdade de Medicina, Universidade de Lisboa, 2014

(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.

On the dimensionality of the Attachment Style Questionnaire in Italian clinical and nonclinical participants

This study examined the utility of the Attachment Style Questionnaire (ASQ) in an Italian sample of 487 consecutively admitted psychiatric participants and an independent sample of 605 nonclinical participants. Minimum average partial analysis of data from the psychiatric sample supported the hypothesized five-factor structure of the items; furthermore, multiple-group component analysis showed that this five-factor structure was not an artifact of differences in item distributions. The five-factor structure of the ASQ was largely replicated in the nonclinical sample. Furthermore, in both psychiatric and nonclinical samples, a two-factor higher order structure of the ASQ scales was observed. The higher order factors of Avoidance and Anxious Attachment showed meaningful relations with scales assessing parental bonding, but were not redundant with these scales. Multivariate normal mixture analysis supported the hypothesis that adult attachment patterns, as measured by the ASQ, are best considered as dimensional constructs.

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.

Issues of robustness and high dimensionality in cluster analysis

Finite mixture models are being increasingly used to model the distributions of a wide variety of random phenomena. While normal mixture models are often used to cluster data sets of continuous multivariate data, a more robust clustering can be obtained by considering the t mixture model-based approach. Mixtures of factor analyzers enable model-based density estimation to be undertaken for high-dimensional data where the number of observations n is very large relative to their dimension p. As the approach using the multivariate normal family of distributions is sensitive to outliers, it is more robust to adopt the multivariate t family for the component error and factor distributions. The computational aspects associated with robustness and high dimensionality in these approaches to cluster analysis are discussed and illustrated.

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.