Estimation and diagnostics for heteroscedastic nonlinear regression models based on scale mixtures of skew-normal distributions

An extension of some standard likelihood based procedures to heteroscedastic nonlinear regression models under scale mixtures of skew-normal (SMSN) distributions is developed. This novel class of models provides a useful generalization of the heteroscedastic symmetrical nonlinear regression models (Cysneiros et al., 2010), since the random term distributions cover both symmetric as well as asymmetric and heavy-tailed distributions such as skew-t, skew-slash, skew-contaminated normal, among others. A simple EM-type algorithm for iteratively computing maximum likelihood estimates of the parameters is presented and the observed information matrix is derived analytically. In order to examine the performance of the proposed methods, some simulation studies are presented to show the robust aspect of this flexible class against outlying and influential observations and that the maximum likelihood estimates based on the EM-type algorithm do provide good asymptotic properties. Furthermore, local influence measures and the one-step approximations of the estimates in the case-deletion model are obtained. Finally, an illustration of the methodology is given considering a data set previously analyzed under the homoscedastic skew-t nonlinear regression model. (C) 2012 Elsevier B.V. All rights reserved.

Subcritical crack growth and in vitro lifetime prediction of resin composites with different filler distributions

Objectives. Verify the influence of different filler distributions on the subcritical crack growth (SCG) susceptibility, Weibull parameters (m and sigma(0)) and longevity estimated by the strength-probability-time (SPT) diagram of experimental resin composites. Methods. Four composites were prepared, each one containing 59 vol% of glass powder with different filler sizes (d(50) = 0.5; 0.9; 1.2 and 1.9 mu m) and distributions. Granulometric analyses of glass powders were done by a laser diffraction particle size analyzer (Sald-7001, Shimadzu, USA). SCG parameters (n and sigma(f0)) were determined by dynamic fatigue (10(-2) to 10(2) MPa/s) using a biaxial flexural device (12 x 1.2 mm; n = 10). Twenty extra specimens of each composite were tested at 10(0) MPa/s to determine m and sigma(0). Specimens were stored in water at 37 degrees C for 24 h. Fracture surfaces were analyzed under SEM. Results. In general, the composites with broader filler distribution (C0.5 and C1.9) presented better results in terms of SCG susceptibility and longevity. C0.5 and C1.9 presented higher n values (respectively, 31.2 +/- 6.2(a) and 34.7 +/- 7.4(a)). C1.2 (166.42 +/- 0.01(a)) showed the highest and C0.5 (158.40 +/- 0.02(d)) the lowest sigma(f0) value (in MPa). Weibull parameters did not vary significantly (m: 6.6 to 10.6 and sigma(0): 170.6 to 176.4 MPa). Predicted reductions in failure stress (P-f = 5%) for a lifetime of 10 years were approximately 45% for C0.5 and C1.9 and 65% for C0.9 and C1.2. Crack propagation occurred through the polymeric matrix around the fillers and all the fracture surfaces showed brittle fracture features. Significance. Composites with broader granulometric distribution showed higher resistance to SCG and, consequently, higher longevity in vitro. (C) 2012 Academy of Dental Materials. Published by Elsevier Ltd. All rights reserved.

On scale-mixture Birnbaum-Saunders distributions

We present for the first time a justification on the basis of central limit theorems for the family of life distributions generated from scale-mixture of normals. This family was proposed by Balakrishnan et al. (2009) and can be used to accommodate unexpected observations for the usual Birnbaum-Saunders distribution generated from the normal one. The class of scale-mixture of normals includes normal, slash, Student-t, logistic, double-exponential, exponential power and many other distributions. We present a model for the crack extensions where the limiting distribution of total crack extensions is in the class of scale-mixture of normals. (C) 2012 Elsevier B.V. All rights reserved.

Atmospheric conditions, lunar phases, and childbirth: a multivariate analysis

Our objective was to assess extrinsic influences upon childbirth. In a cohort of 1,826 days containing 17,417 childbirths among them 13,252 spontaneous labor admissions, we studied the influence of environment upon the high incidence of labor (defined by 75th percentile or higher), analyzed by logistic regression. The predictors of high labor admission included increases in outdoor temperature (odds ratio: 1.742, P = 0.045, 95%CI: 1.011 to 3.001), and decreases in atmospheric pressure (odds ratio: 1.269, P = 0.029, 95%CI: 1.055 to 1.483). In contrast, increases in tidal range were associated with a lower probability of high admission (odds ratio: 0.762, P = 0.030, 95%CI: 0.515 to 0.999). Lunar phase was not a predictor of high labor admission (P = 0.339). Using multivariate analysis, increases in temperature and decreases in atmospheric pressure predicted high labor admission, and increases of tidal range, as a measurement of the lunar gravitational force, predicted a lower probability of high admission.

The new class of Kummer beta generalized distributions

Ng and Kotz (1995) introduced a distribution that provides greater flexibility to extremes. We define and study a new class of distributions called the Kummer beta generalized family to extend the normal, Weibull, gamma and Gumbel distributions, among several other well-known distributions. Some special models are discussed. The ordinary moments of any distribution in the new family can be expressed as linear functions of probability weighted moments of the baseline distribution. We examine the asymptotic distributions of the extreme values. We derive the density function of the order statistics, mean absolute deviations and entropies. We use maximum likelihood estimation to fit the distributions in the new class and illustrate its potentiality with an application to a real data set.

Quantitative Chemical Profile and Multivariate Statistical Analysis of Alembic Distilled Sugarcane Spirit Fractions

Concentrations of 39 organic compounds were determined in three fractions (head, heart and tail) obtained from the pot still distillation of fermented sugarcane juice. The results were evaluated using analysis of variance (ANOVA), Tukey's test, principal component analysis (PCA), hierarchical cluster analysis (HCA) and linear discriminant analysis (LDA). According to PCA and HCA, the experimental data lead to the formation of three clusters. The head fractions give rise to a more defined group. The heart and tail fractions showed some overlap consistent with its acid composition. The predictive ability of calibration and validation of the model generated by LDA for the three fractions classification were 90.5 and 100%, respectively. This model recognized as the heart twelve of the thirteen commercial cachacas (92.3%) with good sensory characteristics, thus showing potential for guiding the process of cuts.

IDENTIFICATION OF ETHANOLIC WOOD EXTRACTS USING ELECTRONIC ABSORPTION SPECTRUM AND MULTIVARIATE ANALYSIS

IDENTIFICATION OF ETHANOLIC WOOD EXTRACTS USING ELECTRONIC ABSORPTION SPECTRUM AND MULTIVARIATE ANALYSIS. The application of multivariate analysis to spectrophotometric (UV) data was explored for distinguishing extracts of cachaca woods commonly used in the manufacture of casks for aging cachacas (oak, cabretiva-parda, jatoba, amendoim and canela-sassafras). Absorbances close to 280 nm were more strongly correlated with oak and jatoba woods, whereas absorbances near 230 nm were more correlated with canela-sassafras and cabretiva-parda. A comparison between the spectrophotometric model and the model based on chromatographic (HPLC-DAD) data was carried out. The spectrophotometric model better explained the variance data (PC1 + PC2 = 91%) exhibiting potential as a routine method for checking aged spirits.

Size distributions and temporal variations of biological aerosol particles in the Amazon rainforest characterized by microscopy and real-time UV-APS fluorescence techniques during AMAZE-08

As a part of the AMAZE-08 campaign during the wet season in the rainforest of central Amazonia, an ultraviolet aerodynamic particle sizer (UV-APS) was operated for continuous measurements of fluorescent biological aerosol particles (FBAP). In the coarse particle size range (> 1 mu m) the campaign median and quartiles of FBAP number and mass concentration were 7.3x10(4) m(-3) (4.0-13.2x10(4) m(-3)) and 0.72 mu g m(-3) (0.42-1.19 mu g m(-3)), respectively, accounting for 24% (11-41%) of total particle number and 47% (25-65%) of total particle mass. During the five-week campaign in February-March 2008 the concentration of coarse-mode Saharan dust particles was highly variable. In contrast, FBAP concentrations remained fairly constant over the course of weeks and had a consistent daily pattern, peaking several hours before sunrise, suggesting observed FBAP was dominated by nocturnal spore emission. This conclusion was supported by the consistent FBAP number size distribution peaking at 2.3 mu m, also attributed to fungal spores and mixed biological particles by scanning electron microscopy (SEM), light microscopy and biochemical staining. A second primary biological aerosol particle (PBAP) mode between 0.5 and 1.0 mu m was also observed by SEM, but exhibited little fluorescence and no true fungal staining. This mode may have consisted of single bacterial cells, brochosomes, various fragments of biological material, and small Chromalveolata (Chromista) spores. Particles liquid-coated with mixed organic-inorganic material constituted a large fraction of observations, and these coatings contained salts likely from primary biological origin. We provide key support for the suggestion that real-time laser-induce fluorescence (LIF) techniques using 355 nm excitation provide size-resolved concentrations of FBAP as a lower limit for the atmospheric abundance of biological particles in a pristine environment. We also show some limitations of using the instrument for ambient monitoring of weakly fluorescent particles < 2 mu m. Our measurements confirm that primary biological particles, fungal spores in particular, are an important fraction of supermicron aerosol in the Amazon and that may contribute significantly to hydrological cycling, especially when coated by mixed inorganic material.

A Homonuclear Rotational Echo Double-Resonance Method for Measuring Site-Resolved Distance Distributions in I=1/2 Spin Pairs, Clusters, and Multispin Systems

Deutsche Forschungsgemeinschaft [SFB 858]

Levy stable distributions via associated integral transform

We present a method of generation of exact and explicit forms of one-sided, heavy-tailed Levy stable probability distributions g(alpha)(x), 0 <= x < infinity, 0 < alpha < 1. We demonstrate that the knowledge of one such a distribution g a ( x) suffices to obtain exactly g(alpha)p ( x), p = 2, 3, .... Similarly, from known g(alpha)(x) and g(beta)(x), 0 < alpha, beta < 1, we obtain g(alpha beta)( x). The method is based on the construction of the integral operator, called Levy transform, which implements the above operations. For a rational, alpha = l/k with l < k, we reproduce in this manner many of the recently obtained exact results for g(l/k)(x). This approach can be also recast as an application of the Efros theorem for generalized Laplace convolutions. It relies solely on efficient definite integration. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4709443]

Effect of Genetic Variants Related to Lipid Metabolism as Risk Factors for Cholelithiasis After Bariatric Surgery in Brazilian Population

The manifestation of cholelithiasis after bariatric surgery may depend on genetic factors related to lipid metabolism, including apolipoprotein E (APOE) and cholesteryl ester transfer protein (CETP) gene polymorphisms. We investigated the association between APOE HhaI and CETP TaqIB polymorphisms [PCR-RFLP] and occurrence of cholelithiasis over up to 8 months of follow-up after gastroplasty to Roux-en-Y gastric bypass in 220 patients distributed in Group 1 (G1) 114 with cholelithiasis postoperatively and Group 2 (G2) 106 without cholelithiasis, including biochemical and anthropometric profiles analyses. In our series, the allelic and genotypic distributions of CETP TaqIB and APOE HhaI polymorphisms were similar in both groups (P > 0.05). The subgroup analysis evidenced that 54% of the patients from G1, APOE*4 allele carriers compared with APOE*3/3 carriers, presented altered low-density lipoprotein cholesterol (LDL cholesterol) serum levels (P = 0.022) before bariatric surgery. The B1 allele for CETP was associated to more quickly elevation of HDL cholesterol levels just in individuals without cholelitiasis (P < 0.0001). The multivariate logistic regression analysis demonstrates correlation between APOE*4 allele, higher total cholesterol (TC) serum levels and prediposition to cholelitiasis in preoperative period. However, the presence of postoperative cholelithiasis was not associated with altered lipid profile. The CETP TaqIB and APOE HhaI polymorphisms do not seem to have association with gallstones in the late postoperative bariatric surgery, considering that these genetic variants do not differ subgroups of patients who are eligible to routine prophylactic cholecystectomy, at least in Brazilian population.

Multivariate analysis in the maximum strength performance

This study performed an exploratory analysis of the anthropometrical and morphological muscle variables related to the one-repetition maximum (1RM) performance. In addition, the capacity of these variables to predict the force production was analyzed. 50 active males were submitted to the experimental procedures: vastus lateralis muscle biopsy, quadriceps magnetic resonance imaging, body mass assessment and 1RM test in the leg-press exercise. K-means cluster analysis was performed after obtaining the body mass, sum of the left and right quadriceps muscle cross-sectional area (Sigma CSA), percentage of the type II fibers and the 1RM performance. The number of clusters was defined a priori and then were labeled as high strength performance (HSP1RM) group and low strength performance (LSP1RM) group. Stepwise multiple regressions were performed by means of body mass, Sigma CSA, percentage of the type II fibers and clusters as predictors' variables and 1RM performance as response variable. The clusters mean +/- SD were: 292.8 +/- 52.1 kg, 84.7 +/- 17.9 kg, 19249.7 +/- 1645.5 mm(2) and 50.8 +/- 7.2% for the HSP1RM and 254.0 +/- 51.1 kg, 69.2 +/- 8.1 kg, 15483.1 +/- 1 104.8 mm(2) and 51.7 +/- 6.2 %, for the LSP1RM in the 1RM, body mass, Sigma CSA and muscle fiber type II percentage, respectively. The most important variable in the clusters division was the Sigma CSA. In addition, the Sigma CSA and muscle fiber type II percentage explained the variance in the 1RM performance (Adj R-2 = 0.35, p = 0.0001) for all participants and for the LSP1RM (Adj R-2 = 0.25, p = 0.002). For the HSP1RM, only the Sigma CSA was entered in the model and showed the highest capacity to explain the variance in the 1RM performance (Adj R-2 = 0.38, p = 0.01). As a conclusion, the muscle CSA was the most relevant variable to predict force production in individuals with no strength training background.

Understanding the correlates of adolescents' dietary intake patterns. A multivariate analysis

We investigated dietary intake patterns (DIP) in adolescents (14-18 year-olds) and the association with demographic and socioeconomic characteristics and lifestyle variables. This school-based survey was carried out among high school students from the city of Maringa in the state of Parana (PR), Brazil (2007). The sample included 991 students (54.5% girls) from high schools. DIPs were investigated by the frequency of weekly consumption of each food group: vegetables, fruit, rice, beans, fried food, sweet food, milk, soda, meat, eggs, alcoholic drinks. Independent variables were: demographic and socioeconomic characteristics and lifestyle variables. DIPS were identified using principal component analysis with orthogonal rotation (varimax). Three components were extracted. Component 1 (fried foods, sweets and soft drinks) was positively associated with not having breakfast for girls and dinner for boys. Moreover, component 2 (consumption of fruit and vegetables) was positively associated with having breakfast at home for boys and number of meals for girls. Component 3 (beans, eggs and meat) was positively associated with having lunch, employment and sedentary behavior level for girls. However, it was negatively associated with having lunch and dinner for boys. Adolescents who have healthier eating patterns also had other healthier behaviors regardless of gender. However, factors associated with dietary patterns differ between boys and girls. (C) 2012 Elsevier Ltd. All rights reserved.

Signal propagation in Bayesian networks and its relationship with intrinsically multivariate predictive variables

A set of predictor variables is said to be intrinsically multivariate predictive (IMP) for a target variable if all properly contained subsets of the predictor set are poor predictors of the. target but the full set predicts the target with great accuracy. In a previous article, the main properties of IMP Boolean variables have been analytically described, including the introduction of the IMP score, a metric based on the coefficient of determination (CoD) as a measure of predictiveness with respect to the target variable. It was shown that the IMP score depends on four main properties: logic of connection, predictive power, covariance between predictors and marginal predictor probabilities (biases). This paper extends that work to a broader context, in an attempt to characterize properties of discrete Bayesian networks that contribute to the presence of variables (network nodes) with high IMP scores. We have found that there is a relationship between the IMP score of a node and its territory size, i.e., its position along a pathway with one source: nodes far from the source display larger IMP scores than those closer to the source, and longer pathways display larger maximum IMP scores. This appears to be a consequence of the fact that nodes with small territory have larger probability of having highly covariate predictors, which leads to smaller IMP scores. In addition, a larger number of XOR and NXOR predictive logic relationships has positive influence over the maximum IMP score found in the pathway. This work presents analytical results based on a simple structure network and an analysis involving random networks constructed by computational simulations. Finally, results from a real Bayesian network application are provided. (C) 2012 Elsevier Inc. All rights reserved.

Sets of probability distributions, independence, and convexity

This paper analyzes concepts of independence and assumptions of convexity in the theory of sets of probability distributions. The starting point is Kyburg and Pittarelli's discussion of "convex Bayesianism" (in particular their proposals concerning E-admissibility, independence, and convexity). The paper offers an organized review of the literature on independence for sets of probability distributions; new results on graphoid properties and on the justification of "strong independence" (using exchangeability) are presented. Finally, the connection between Kyburg and Pittarelli's results and recent developments on the axiomatization of non-binary preferences, and its impact on "complete" independence, are described.