Study on Bivariate Normal Distribution

Let (X, Y) be bivariate normal random vectors which represent the responses as a result of Treatment 1 and Treatment 2. The statistical inference about the bivariate normal distribution parameters involving missing data with both treatment samples is considered. Assuming the correlation coefficient ρ of the bivariate population is known, the MLE of population means and variance (ξ, η, and σ2) are obtained. Inferences about these parameters are presented. Procedures of constructing confidence interval for the difference of population means ξ – η and testing hypothesis about ξ – η are established. The performances of the new estimators and testing procedure are compared numerically with the method proposed in Looney and Jones (2003) on the basis of extensive Monte Carlo simulation. Simulation studies indicate that the testing power of the method proposed in this thesis study is higher.

Inferences about Parameters of Trivariate Normal Distribution with Missing Data

Multivariate normal distribution is commonly encountered in any field, a frequent issue is the missing values in practice. The purpose of this research was to estimate the parameters in three-dimensional covariance permutation-symmetric normal distribution with complete data and all possible patterns of incomplete data. In this study, MLE with missing data were derived, and the properties of the MLE as well as the sampling distributions were obtained. A Monte Carlo simulation study was used to evaluate the performance of the considered estimators for both cases when ρ was known and unknown. All results indicated that, compared to estimators in the case of omitting observations with missing data, the estimators derived in this article led to better performance. Furthermore, when ρ was unknown, using the estimate of ρ would lead to the same conclusion.

Consumer populations and nutritional transition in Spain in the 20th century: a methodology for their reconstruction

One feature of the modern nutrition transition is the growing consumption of animal proteins. The most common approach in the quantitative analysis of this change used to be the study of averages of food consumption. But this kind of analysis seems to be incomplete without the knowledge of the number of consumers. Data about consumers are not usually published in historical statistics. This article introduces a methodological approach for reconstructing consumer populations. This methodology is based on some assumptions about the diffusion process of foodstuffs and the modeling of consumption patterns with a log-normal distribution. This estimating process is illustrated with the specific case of milk consumption in Spain between 1925 and 1981. These results fit quite well with other data and indirect sources available showing that this dietary change was a slow and late process. The reconstruction of consumer population could shed a new light in the study of nutritional transitions.

Occupational exposure of truck drivers to dust and polynuclear aromatic hydrocarbons: a pilot study in Geneva, Switzerland.

The exposure to dust and polynuclear aromatic hydrocarbons (PAH) of 15 truck drivers from Geneva, Switzerland, was measured. The drivers were divided between "long-distance" drivers and "local" drivers and between smokers and nonsmokers and were compared with a control group of 6 office workers who were also divided into smokers and nonsmokers. Dust was measured on 1 workday both by a direct-reading instrument and by sampling. The local drivers showed higher exposure to dust (0.3 mg/m3) and PAH than the long-distance drivers (0.1 mg/m3), who showed no difference with the control group. This observation may be due to the fact that the local drivers spend more time in more polluted areas, such as streets with heavy traffic and construction sites, than do the long-distance drivers. Smoking does not influence exposure to dust and PAH of professional truck drivers, as measured in this study, probably because the ventilation rate of the truck cabins is relatively high even during cold days (11-15 r/h). The distribution of dust concentrations was shown in some cases to be quite different from the expected log-normal distribution. The contribution of diesel exhaust to these exposures could not be estimated since no specific tracer was used. However, the relatively low level of dust exposure dose not support the hypothesis that present day levels of diesel exhaust particulates play a significant role in the excess occurrence of lung cancer observed in professional truck drivers.

Calidad microbiológica de las aguas de baño del embalse Conde del Guadalhorce, de Málaga (España), 2000-2005.

Background. Determine the presence and evolution of indicators microorganisms of water pollution in “Conde del Guadalhorce” reservoir, Málaga city, Spain. A second objective was to analyze pollution degree and evaluate the sanitary quality of bathing water and compliance with European Directive 76/160/CE. Method. A total of 120 water samples were collected in two bathing freshwater sites during May to September sampling period between 2000 to 2005, and the numbers of total coliforms (CT), faecal coliforms (CF) and faecal streptococci (EF) were enumerated using the membrane filtration method. We used the log-normal distribution method and calculate the logarithmic means, percentile points, ratios CF:EF, ANOVA and Pearson correlations. Results. Only two samples overcome CF limit values at Camping sampling station during 2000 year. Ratios CF:EF values were higher (> 4) during 2000 to 2002, and lower (< 0,7) during 2003 to 2005. Significant differences (ANOVA F = 3,41, ∝ < 0,01) was only observed with EF during evaluated period. There was no significant difference between concentration means at bathing water sites (ANOVA, F = 3,395, ∝ < 0,01). The counts of CT and CF were significantly correlated in Kiosko water samples, while in Camping water, significant correlation (t = 0,632, p < 0,05) was only observed with EF at the Camping station during 2000, 2003 and 2005 years. Conclusions. “Conde del Guadalhorce” reservoir showed hygienic conditions for safety bathing. Globally, water bathing quality is good. CT, CF y EF indicators were agreed with UE Directive during 2000- 2005, with exception CF at Camping station in 2000 year. CT y CF concentrations at Camping were frecuently higher than Kiosko, it could be caused to swimmers abundance and recreational activities. There was a trend towards rising EF, it could be caused to faecal pollution source of animal origin, needed to research it.

Critical ruptures in a bundle of slowly relaxing fibers

We study the damage enhanced creep rupture of disordered materials by means of a fiber bundle model. Broken fibers undergo a slow stress relaxation modeled by a Maxwell element whose stress exponent m can vary in a broad range. Under global load sharing we show that due to the strength disorder of fibers, the lifetime ʧ of the bundle has sample-to-sample fluctuations characterized by a log-normal distribution independent of the type of disorder. We determine the Monkman-Grant relation of the model and establish a relation between the rupture life tʄ and the characteristic time tm of the intermediate creep regime of the bundle where the minimum strain rate is reached, making possible reliable estimates of ʧ from short term measurements. Approaching macroscopic failure, the deformation rate has a finite time power law singularity whose exponent is a decreasing function of m. On the microlevel the distribution of waiting times is found to have a power law behavior with m-dependent exponents different below and above the critical load of the bundle. Approaching the critical load from above, the cutoff value of the distributions has a power law divergence whose exponent coincides with the stress exponent of Maxwell elements

Determining the best population-level alcohol consumption model and its impact on estimates of alcohol-attributable harms.

BACKGROUND: The goals of our study are to determine the most appropriate model for alcohol consumption as an exposure for burden of disease, to analyze the effect of the chosen alcohol consumption distribution on the estimation of the alcohol Population- Attributable Fractions (PAFs), and to characterize the chosen alcohol consumption distribution by exploring if there is a global relationship within the distribution. METHODS: To identify the best model, the Log-Normal, Gamma, and Weibull prevalence distributions were examined using data from 41 surveys from Gender, Alcohol and Culture: An International Study (GENACIS) and from the European Comparative Alcohol Study. To assess the effect of these distributions on the estimated alcohol PAFs, we calculated the alcohol PAF for diabetes, breast cancer, and pancreatitis using the three above-named distributions and using the more traditional approach based on categories. The relationship between the mean and the standard deviation from the Gamma distribution was estimated using data from 851 datasets for 66 countries from GENACIS and from the STEPwise approach to Surveillance from the World Health Organization. RESULTS: The Log-Normal distribution provided a poor fit for the survey data, with Gamma and Weibull distributions providing better fits. Additionally, our analyses showed that there were no marked differences for the alcohol PAF estimates based on the Gamma or Weibull distributions compared to PAFs based on categorical alcohol consumption estimates. The standard deviation of the alcohol distribution was highly dependent on the mean, with a unit increase in alcohol consumption associated with a unit increase in the mean of 1.258 (95% CI: 1.223 to 1.293) (R2 = 0.9207) for women and 1.171 (95% CI: 1.144 to 1.197) (R2 = 0. 9474) for men. CONCLUSIONS: Although the Gamma distribution and the Weibull distribution provided similar results, the Gamma distribution is recommended to model alcohol consumption from population surveys due to its fit, flexibility, and the ease with which it can be modified. The results showed that a large degree of variance of the standard deviation of the alcohol consumption Gamma distribution was explained by the mean alcohol consumption, allowing for alcohol consumption to be modeled through a Gamma distribution using only average consumption.

Size Effect and Notch Size Effect in Metal Fatigue

It is a well known phenomenon that the constant amplitude fatigue limit of a large component is lower than the fatigue limit of a small specimen made of the same material. In notched components the opposite occurs: the fatigue limit defined as the maximum stress at the notch is higher than that achieved with smooth specimens. These two effects have been taken into account in most design handbooks with the help of experimental formulas or design curves. The basic idea of this study is that the size effect can mainly be explained by the statistical size effect. A component subjected to an alternating load can be assumed to form a sample of initiated cracks at the end of the crack initiation phase. The size of the sample depends on the size of the specimen in question. The main objective of this study is to develop a statistical model for the estimation of this kind of size effect. It was shown that the size of a sample of initiated cracks shall be based on the stressed surface area of the specimen. In case of varying stress distribution, an effective stress area must be calculated. It is based on the decreasing probability of equally sized initiated cracks at lower stress level. If the distribution function of the parent population of cracks is known, the distribution of the maximum crack size in a sample can be defined. This makes it possible to calculate an estimate of the largest expected crack in any sample size. The estimate of the fatigue limit can now be calculated with the help of the linear elastic fracture mechanics. In notched components another source of size effect has to be taken into account. If we think about two specimens which have similar shape, but the size is different, it can be seen that the stress gradient in the smaller specimen is steeper. If there is an initiated crack in both of them, the stress intensity factor at the crack in the larger specimen is higher. The second goal of this thesis is to create a calculation method for this factor which is called the geometric size effect. The proposed method for the calculation of the geometric size effect is also based on the use of the linear elastic fracture mechanics. It is possible to calculate an accurate value of the stress intensity factor in a non linear stress field using weight functions. The calculated stress intensity factor values at the initiated crack can be compared to the corresponding stress intensity factor due to constant stress. The notch size effect is calculated as the ratio of these stress intensity factors. The presented methods were tested against experimental results taken from three German doctoral works. Two candidates for the parent population of initiated cracks were found: the Weibull distribution and the log normal distribution. Both of them can be used successfully for the prediction of the statistical size effect for smooth specimens. In case of notched components the geometric size effect due to the stress gradient shall be combined with the statistical size effect. The proposed method gives good results as long as the notch in question is blunt enough. For very sharp notches, stress concentration factor about 5 or higher, the method does not give sufficient results. It was shown that the plastic portion of the strain becomes quite high at the root of this kind of notches. The use of the linear elastic fracture mechanics becomes therefore questionable.

Start-up rheometry of highly polydisperse magnetorheological fluids: experiments and simulations

An extensive experimental and simulation study is carried out in conventional magnetorheological fluids formulated by dispersion of mixtures of carbonyl iron particles having different sizes in Newtonian carriers. Apparent yield stress data are reported for a wide range of polydispersity indexes (PDI) from PDI = 1.63 to PDI = 3.31, which for a log-normal distribution corresponds to the standard deviation ranging from to . These results demonstrate that the effect of polydispersity is negligible in this range in spite of exhibiting very different microstructures. Experimental data in the magnetic saturation regime are in quantitative good agreement with particle-level simulations under the assumption of dipolar magnetostatic forces. The insensitivity of the yield stresses to the polydispersity can be understood from the interplay between the particle cluster size distribution and the packing density of particles inside the clusters.

Análise estatística de impedância da pele

In this work we analyze the skin bioimpedance statistical distribution. We focus on the study of two distinct samples: the statistics of impedance of several points in the skin of a single individual and the statistics over a population (many individuals) but in a single skin point. The impedance data was obtained from the literature (Pearson, 2007). Using the Shapiro-Wilk test and the assymmetry test we conclude that the impedance of a population is better described by an assymetric and non-normal distribution. On the other side, the data concerning the individual impedance seems to follow a normal distribution. We have performed a goodnes of fitting test and the better distribution to fit the data of a population is the log-normal distribution. It is interesting to note that our result for skin impedance is in simtony with body impedance from the literature of electrical engeneering. Our results have an impact over the statistical planning and modelling of skin impedance experiments. Special attention we should drive to the treatment of outliers in this kind of dataset. The results of this work are important in the general discussion of low impedance of points of acupuncture and also in the problem of skin biopotentials used in equipments like the Electrodermal Screen Tests.

Factors determining the structure of fish communities in Pantanal lagoons (MS, Brazil)

The fish communities of lagoons in the Nhecolandia Pantanal were studied to determine the factors which are responsible for the composition and abundance of species. Fishes were collected in 19 lagoons during August 1997, after their isolation from the River Negro, using beach seines (15 x 1.5 m; 2 mm mesh). A total of 51 species were collected. In the lagoons, or in parts with dense macrophytes, a screened box trap was used. Fishing was also accomplished with hooks of several sizes. Species richness was estimated by the jack-knife procedure, after adjustment to the log-normal distribution and with von Bertalanffy's equation (asymptotic). The most important factors in the community organization were macrophyte cover, piscivore abundance and depth of the lagoons. The role of these habitats in the Pantanal ecosystem was discussed.

Improving methods in gap ecology: Revisiting size and shape distributions using a model selection approach

Questions: We assess gap size and shape distributions, two important descriptors of the forest disturbance regime, by asking: which statistical model best describes gap size distribution; can simple geometric forms adequately describe gap shape; does gap size or shape vary with forest type, gap age or the method used for gap delimitation; and how similar are the studied forests and other tropical and temperate forests? Location: Southeastern Atlantic Forest, Brazil. Methods: Analysing over 150 gaps in two distinct forest types (seasonal and rain forests), a model selection framework was used to select appropriate probability distributions and functions to describe gap size and gap shape. The first was described using univariate probability distributions, whereas the latter was assessed based on the gap area-perimeter relationship. Comparisons of gap size and shape between sites, as well as size and age classes were then made based on the likelihood of models having different assumptions for the values of their parameters. Results: The log-normal distribution was the best descriptor of gap size distribution, independently of the forest type or gap delimitation method. Because gaps became more irregular as they increased in size, all geometric forms (triangle, rectangle and ellipse) were poor descriptors of gap shape. Only when small and large gaps (> 100 or 400m2 depending on the delimitation method) were treated separately did the rectangle and isosceles triangle become accurate predictors of gap shape. Ellipsoidal shapes were poor descriptors. At both sites, gaps were at least 50% longer than they were wide, a finding with important implications for gap microclimate (e.g. light entrance regime) and, consequently, for gap regeneration. Conclusions: In addition to more appropriate descriptions of gap size and shape, the model selection framework used here efficiently provided a means by which to compare the patterns of two different types of forest. With this framework we were able to recommend the log-normal parameters μ and σ for future comparisons of gap size distribution, and to propose possible mechanisms related to random rates of gap expansion and closure. We also showed that gap shape varied highly and that no single geometric form was able to predict the shape of all gaps, the ellipse in particular should no longer be used as a standard gap shape. © 2012 International Association for Vegetation Science.

Processo de attachment, densidade de raios e a ocorrência de flashover em linhas de transmissão de eletricidade

Esta tese refere-se ao uso de modelos de attachment de lideres (Leader Progression Model – LPM) para estimativa da distancia de salto, que, junto com dados de densidade de raios como coletado pelo sistema de detecção de raios do SIPAM, são usados para estimar a taxa de flashover em linhas de transmissão de eletricidade na Amazônia. O modelo de progressão de líder desenvolvido nesta tese é denominado ModSalto, que, para estimar a distancia de salto: 1- integra a densidade linear de cargas elétricas devida ao líder, proporcional a prospectiva corrente de primeira descarga (1º stroke) Ip, para determinar o campo elétrico produzido pelo líder descendente na estrutura sob estudo (para-raios, arestas, condutores, etc.); 2 – integra a distribuição imagem do líder descendente, fazendo uso da característica de poder das pontas como fator de estimulo e intensificação do campo elétrico devido a streamers nas partes aguçadas das estruturas sob estudo, para temporizar o momento do processo de attachment. O gatilho para o líder descendente, por hipótese, se deve a ejeção aleatória de pacotes de cargas elétricas em domínios turbulentos no interior das nuvens que recebem energia por processos de cascata da turbulência geral, e o comportamento do líder descendente, deve obedecer à equação da força de Lorentz, no espaço de campos cruzados elétrico devido às nuvens a cima e o campo magnético da Terra, que obriga as cargas do líder a desenvolverem movimentos cicloidais que podem explicar a natureza tortuosa do trajeto do líder descendente. Com o objetivo de formalizar dados consistentes de densidade de raios é feita uma reanalise do conjunto de dados coletado pelo LLS SIPAM de outubro de 2006 a julho de 2008 na região amazônica, com cerca de 3 milhões de eventos, comparando-os com dados de torres instrumentadas para evidenciar-se sua qualidade e usabilidade. Dados de elevação de terreno do SRTM da NASA são usados para gerar formula do raio de atração (Ra) das estruturas passiveis de serem atingidas por raios e para gerar formulas de área de atração, usadas para quantificar o numero de raios que provavelmente atingirão determinada estrutura, baseado no valor de densidade de raios (raios/km2/ano) na área em estudo.

Using complex networks to quantify consistency in the use of words

In this paper we have quantified the consistency of word usage in written texts represented by complex networks, where words were taken as nodes, by measuring the degree of preservation of the node neighborhood. Words were considered highly consistent if the authors used them with the same neighborhood. When ranked according to the consistency of use, the words obeyed a log-normal distribution, in contrast to Zipf's law that applies to the frequency of use. Consistency correlated positively with the familiarity and frequency of use, and negatively with ambiguity and age of acquisition. An inspection of some highly consistent words confirmed that they are used in very limited semantic contexts. A comparison of consistency indices for eight authors indicated that these indices may be employed for author recognition. Indeed, as expected, authors of novels could be distinguished from those who wrote scientific texts. Our analysis demonstrated the suitability of the consistency indices, which can now be applied in other tasks, such as emotion recognition.