A Poisson mixed model with nonnormal random effect distribution

In this paper, we propose a random intercept Poisson model in which the random effect is assumed to follow a generalized log-gamma (GLG) distribution. This random effect accommodates (or captures) the overdispersion in the counts and induces within-cluster correlation. We derive the first two moments for the marginal distribution as well as the intraclass correlation. Even though numerical integration methods are, in general, required for deriving the marginal models, we obtain the multivariate negative binomial model from a particular parameter setting of the hierarchical model. An iterative process is derived for obtaining the maximum likelihood estimates for the parameters in the multivariate negative binomial model. Residual analysis is proposed and two applications with real data are given for illustration. (C) 2011 Elsevier B.V. All rights reserved.

Segmental dataset and whole body expression data do not support the hypothesis that non-random movement is an intrinsic property of Drosophila retrogenes

Background: Several studies in Drosophila have shown excessive movement of retrogenes from the X chromosome to autosomes, and that these genes are frequently expressed in the testis. This phenomenon has led to several hypotheses invoking natural selection as the process driving male-biased genes to the autosomes. Metta and Schlotterer (BMC Evol Biol 2010, 10:114) analyzed a set of retrogenes where the parental gene has been subsequently lost. They assumed that this class of retrogenes replaced the ancestral functions of the parental gene, and reported that these retrogenes, although mostly originating from movement out of the X chromosome, showed female-biased or unbiased expression. These observations led the authors to suggest that selective forces (such as meiotic sex chromosome inactivation and sexual antagonism) were not responsible for the observed pattern of retrogene movement out of the X chromosome. Results: We reanalyzed the dataset published by Metta and Schlotterer and found several issues that led us to a different conclusion. In particular, Metta and Schlotterer used a dataset combined with expression data in which significant sex-biased expression is not detectable. First, the authors used a segmental dataset where the genes selected for analysis were less testis-biased in expression than those that were excluded from the study. Second, sex-biased expression was defined by comparing male and female whole-body data and not the expression of these genes in gonadal tissues. This approach significantly reduces the probability of detecting sex-biased expressed genes, which explains why the vast majority of the genes analyzed (parental and retrogenes) were equally expressed in both males and females. Third, the female-biased expression observed by Metta and Schltterer is mostly found for parental genes located on the X chromosome, which is known to be enriched with genes with female-biased expression. Fourth, using additional gonad expression data, we found that autosomal genes analyzed by Metta and Schlotterer are less up regulated in ovaries and have higher chance to be expressed in meiotic cells of spermatogenesis when compared to X-linked genes. Conclusions: The criteria used to select retrogenes and the sex-biased expression data based on whole adult flies generated a segmental dataset of female-biased and unbiased expressed genes that was unable to detect the higher propensity of autosomal retrogenes to be expressed in males. Thus, there is no support for the authors' view that the movement of new retrogenes, which originated from X-linked parental genes, was not driven by selection. Therefore, selection-based genetic models remain the most parsimonious explanations for the observed chromosomal distribution of retrogenes.

Bayesian analysis of random regression models using B-splines to model test-day milk yield of Holstein cattle in Brazil

The objective of this paper is to model variations in test-day milk yields of first lactations of Holstein cows by RR using B-spline functions and Bayesian inference in order to fit adequate and parsimonious models for the estimation of genetic parameters. They used 152,145 test day milk yield records from 7317 first lactations of Holstein cows. The model established in this study was additive, permanent environmental and residual random effects. In addition, contemporary group and linear and quadratic effects of the age of cow at calving were included as fixed effects. Authors modeled the average lactation curve of the population with a fourth-order orthogonal Legendre polynomial. They concluded that a cubic B-spline with seven random regression coefficients for both the additive genetic and permanent environment effects was to be the best according to residual mean square and residual variance estimates. Moreover they urged a lower order model (quadratic B-spline with seven random regression coefficients for both random effects) could be adopted because it yielded practically the same genetic parameter estimates with parsimony. (C) 2012 Elsevier B.V. All rights reserved.

Superdiffusion in a non-Markovian random walk model with a Gaussian memory profile

Most superdiffusive Non-Markovian random walk models assume that correlations are maintained at all time scales, e. g., fractional Brownian motion, Levy walks, the Elephant walk and Alzheimer walk models. In the latter two models the random walker can always "remember" the initial times near t = 0. Assuming jump size distributions with finite variance, the question naturally arises: is superdiffusion possible if the walker is unable to recall the initial times? We give a conclusive answer to this general question, by studying a non-Markovian model in which the walker's memory of the past is weighted by a Gaussian centered at time t/2, at which time the walker had one half the present age, and with a standard deviation sigma t which grows linearly as the walker ages. For large widths we find that the model behaves similarly to the Elephant model, but for small widths this Gaussian memory profile model behaves like the Alzheimer walk model. We also report that the phenomenon of amnestically induced persistence, known to occur in the Alzheimer walk model, arises in the Gaussian memory profile model. We conclude that memory of the initial times is not a necessary condition for generating (log-periodic) superdiffusion. We show that the phenomenon of amnestically induced persistence extends to the case of a Gaussian memory profile.

Random regression analyses using B-spline functions to model growth of Nellore cattle

The objective of this study was to estimate (co)variance components using random regression on B-spline functions to weight records obtained from birth to adulthood. A total of 82 064 weight records of 8145 females obtained from the data bank of the Nellore Breeding Program (PMGRN/Nellore Brazil) which started in 1987, were used. The models included direct additive and maternal genetic effects and animal and maternal permanent environmental effects as random. Contemporary group and dam age at calving (linear and quadratic effect) were included as fixed effects, and orthogonal Legendre polynomials of age (cubic regression) were considered as random covariate. The random effects were modeled using B-spline functions considering linear, quadratic and cubic polynomials for each individual segment. Residual variances were grouped in five age classes. Direct additive genetic and animal permanent environmental effects were modeled using up to seven knots (six segments). A single segment with two knots at the end points of the curve was used for the estimation of maternal genetic and maternal permanent environmental effects. A total of 15 models were studied, with the number of parameters ranging from 17 to 81. The models that used B-splines were compared with multi-trait analyses with nine weight traits and to a random regression model that used orthogonal Legendre polynomials. A model fitting quadratic B-splines, with four knots or three segments for direct additive genetic effect and animal permanent environmental effect and two knots for maternal additive genetic effect and maternal permanent environmental effect, was the most appropriate and parsimonious model to describe the covariance structure of the data. Selection for higher weight, such as at young ages, should be performed taking into account an increase in mature cow weight. Particularly, this is important in most of Nellore beef cattle production systems, where the cow herd is maintained on range conditions. There is limited modification of the growth curve of Nellore cattle with respect to the aim of selecting them for rapid growth at young ages while maintaining constant adult weight.

Implant Stability Measurements of Two Immediate Loading Protocols for the Edentulous Mandible: Rigid and Semi-rigid Splinting of the Implants

Aim: Primary and secondary stabilities of immediately loaded mandibular implants restored with fixed prostheses (FP) using rigid or semirigid splinting systems were clinically and radiographically evaluated. Methods: Fifteen edentulous patients were rehabilitated using hybrid FP; each had 5 implants placed between the mental foramens. Two groups were randomly divided: group 1-FP with the conventional rigid bar splinting the implants and group 2-semi-rigid cantilever extension system with titanium bars placed in the 2 distal abutment cylinders. Primary stability was evaluated using resonance frequency analysis after installation of the implant abutments. The measurements were made at 3 times: T0, at baseline; T1, 4 months after implant placement; and T2, 8 months after implant placement. Presence of mobility and inflammation in the implant surrounding regions were checked. Stability data were submitted to statistical analysis for comparison between groups (P, 0.05). Results: Implant survival rate for the implants was of 100% in both groups. No significant differences in the mean implant stability quotient values were found for both groups from baseline and after the 8-month follow-up. Conclusion: The immediate loading of the implants was satisfactory, and both splinting conditions (rigid and semi-rigid) can be successfully used for the restoration of edentulous mandibles. (Implant Dent 2012;21:486-490)

The impact of healing time before loading on orthodontic mini-implant stability: a histomophometric study in minipigs

Objective: To evaluate healing time before loading, areas compression and tension and location of insertion on mini-implant stability. Design: Six minipigs were used. Each animal received 3 mini-implants in each quadrant: 1 mini-implant was used as an unloaded control (G1, n = 24); the other 2 were loaded with 150 g-force at three time intervals (G2: immediate loading, G3: after 15 days and G4: after 30 days), with 16 mini-implant in each experimental group. After 120 days, tissue blocks of the areas of interest were harvested. Clinical analysis (exact Fisher test) determined the survival rate. Histological analysis (Kontron KS 300TM, Zeiss) quantified the fractional bone-toimplant contact (%BIC) and bone area (%BA) at each healing time point, areas of interest, and insertion site (ANOVA and t tests for dependent and independent samples). Results: The mini-implant survival rates were G1: 71%, G2: 50%, G3: 75% and G4: 63%, with no statistical differences between them. The groups presented similar %BIC and %BA. There were no differences between the compression and tension sides or maxillary and mandibular insertion sites. Conclusions: These results suggest that low-intensity immediate or early orthodontic loading does not affect mini-implant stability, because similar histomorphometric results were observed for all the groups, with partial osseointegration of the mini-implants present.

The effect of CCL3 and CCR1 in bone remodeling induced by mechanical loading during orthodontic tooth movement in mice

Bone remodeling is affected by mechanical loading and inflammatory mediators, including chemokines. The chemokine (C–C motif) ligand 3 (CCL3) is involved in bone remodeling by binding to C–C chemokine receptors 1 and 5 (CCR1 and CCR5) expressed on osteoclasts and osteoblasts. Our group has previously demonstrated that CCR5 down-regulates mechanical loading-induced bone resorption. Thus, the present study aimed to investigate the role of CCR1 and CCL3 in bone remodeling induced by mechanical loading during orthodontic tooth movement in mice. Our results showed that bone remodeling was significantly decreased in CCL3−/− and CCR1−/− mice and in animals treated with Met-RANTES (an antagonist of CCR5 and CCR1). mRNA levels of receptor activator of nuclear factor kappa-B (RANK), its ligand RANKL, tumor necrosis factor alpha (TNF-α) and RANKL/osteoprotegerin (OPG) ratio were diminished in the periodontium of CCL3−/− mice and in the group treated with Met-RANTES. Met-RANTES treatment also reduced the levels of cathepsin K and metalloproteinase 13 (MMP13). The expression of the osteoblast markers runt-related transcription factor 2 (RUNX2) and periostin was decreased, while osteocalcin (OCN) was augmented in CCL3−/− and Met-RANTES-treated mice. Altogether, these findings show that CCR1 is pivotal for bone remodeling induced by mechanical loading during orthodontic tooth movement and these actions depend, at least in part, on CCL3.

Study of an anisotropic polymeric cellular material under compression loading

This paper emphasizes the influence of micro mechanisms of failure of a cellular material on its phenomenological response. Most of the applications of cellular materials comprise a compression loading. Thus, the study focuses on the influence of the anisotropy in the mechanical behavior of cellular material under cyclic compression loadings. For this study, a Digital Image Correlation (DIC) technique (named Correli) was applied, as well as SEM (Scanning Electron Microscopy) images were analyzed. The experimental results are discussed in detail for a closed-cell rigid poly (vinyl chloride) (PVC) foam, showing stress-strain curves in different directions and why the material can be assumed as transversely isotropic. Besides, the present paper shows elastic and plastic Poisson's ratios measured in different planes, explaining why the plastic Poisson's ratios approach to zero. Yield fronts created by the compression loadings in different directions and the influence of spring-back phenomenon on hardening curves are commented, also.

Effects of specimen geometry and loading mode on crack growth resistance curves of a high-strength pipeline girth weld

This work presents an investigation of the ductile tearing properties for a girth weld made of an API 5L X80 pipeline steel using experimentally measured crack growth resistance curves. Use of these materials is motivated by the increasing demand in the number of applications for manufacturing high strength pipes for the oil and gas industry including marine applications and steel catenary risers. Testing of the pipeline girth welds employed side-grooved, clamped SE(T) specimens and shallow crack bend SE(B) specimens with a weld centerline notch to determine the crack growth resistance curves based upon the unloading compliance (UC) method using the single specimen technique. Recently developed compliance functions and η-factors applicable for SE(T) and SE(B) fracture specimens with homogeneous material and overmatched welds are introduced to determine crack growth resistance data from laboratory measurements of load-displacement records.

Elastic Propagation in random media: applications to the imaging of volcano structures

High-frequency seismograms contain features that reflect the random inhomogeneities of the earth. In this work I use an imaging method to locate the high contrast small- scale heterogeneity respect to the background earth medium. This method was first introduced by Nishigami (1991) and than applied to different volcanic and tectonically active areas (Nishigami, 1997, Nishigami, 2000, Nishigami, 2006). The scattering imaging method is applied to two volcanic areas: Campi Flegrei and Mt. Vesuvius. Volcanic and seismological active areas are often characterized by complex velocity structures, due to the presence of rocks with different elastic properties. I introduce some modifications to the original method in order to make it suitable for small and highly complex media. In particular, for very complex media the single scattering approximation assumed by Nishigami (1991) is not applicable as the mean free path becomes short. The multiple scattering or diffusive approximation become closer to the reality. In this thesis, differently from the ordinary Nishigami’s method (Nishigami, 1991), I use the mean of the recorded coda envelope as reference curve and calculate the variations from this average envelope. In this way I implicitly do not assume any particular scattering regime for the "average" scattered radiation, whereas I consider the variations as due to waves that are singularly scattered from the strongest heterogeneities. The imaging method is applied to a relatively small area (20 x 20 km), this choice being justified by the small length of the analyzed codas of the low magnitude earthquakes. I apply the unmodified Nishigami’s method to the volcanic area of Campi Flegrei and compare the results with the other tomographies done in the same area. The scattering images, obtained with frequency waves around 18 Hz, show the presence of high scatterers in correspondence with the submerged caldera rim in the southern part of the Pozzuoli bay. Strong scattering is also found below the Solfatara crater, characterized by the presence of densely fractured, fluid-filled rocks and by a strong thermal anomaly. The modified Nishigami’s technique is applied to the Mt. Vesuvius area. Results show a low scattering area just below the central cone and a high scattering area around it. The high scattering zone seems to be due to the contrast between the high rigidity body located beneath the crater and the low rigidity materials located around it. The central low scattering area overlaps the hydrothermal reservoirs located below the central cone. An interpretation of the results in terms of geological properties of the medium is also supplied, aiming to find a correspondence of the scattering properties and the geological nature of the material. A complementary result reported in this thesis is that the strong heterogeneity of the volcanic medium create a phenomenon called "coda localization". It has been verified that the shape of the seismograms recorded from the stations located at the top of the volcanic edifice of Mt. Vesuvius is different from the shape of the seismograms recorded at the bottom. This behavior is justified by the consideration that the coda energy is not uniformly distributed within a region surrounding the source for great lapse time.

Probabilistic approaches in civil engineering: generation of random fields and structural identification with genetic algorithms

The inherent stochastic character of most of the physical quantities involved in engineering models has led to an always increasing interest for probabilistic analysis. Many approaches to stochastic analysis have been proposed. However, it is widely acknowledged that the only universal method available to solve accurately any kind of stochastic mechanics problem is Monte Carlo Simulation. One of the key parts in the implementation of this technique is the accurate and efficient generation of samples of the random processes and fields involved in the problem at hand. In the present thesis an original method for the simulation of homogeneous, multi-dimensional, multi-variate, non-Gaussian random fields is proposed. The algorithm has proved to be very accurate in matching both the target spectrum and the marginal probability. The computational efficiency and robustness are very good too, even when dealing with strongly non-Gaussian distributions. What is more, the resulting samples posses all the relevant, welldefined and desired properties of “translation fields”, including crossing rates and distributions of extremes. The topic of the second part of the thesis lies in the field of non-destructive parametric structural identification. Its objective is to evaluate the mechanical characteristics of constituent bars in existing truss structures, using static loads and strain measurements. In the cases of missing data and of damages that interest only a small portion of the bar, Genetic Algorithm have proved to be an effective tool to solve the problem.

Generating multiple random exams using mathematica package and LATEX

[EN]In the framework of the European Higher Education Area, the assessment has been one of the most important aspects considered. In the Spanish Universities, one of the main differences with regard to the previous system is the incorporation of the continuous assessment to the evaluation process that is understood in several ways depending on the Universities, the courses and the lecturers. Focus on our context, a course of Mathematics of the first academic year in the Faculty of Business Administration at the University of Las Palmas de Gran Canaria (Spain), the continuous assessment has brought the preparation of a large amount of different tests to evaluate the students enrolled in it, therefore the incorporation of new tools and skills in order to make the teaching–learning process easier and more dynamic has become a need mainly in those degrees with a large number of students as the case we consider here. In this work we provide an efficient and effective way to elaborate random multiple-choice examina tion tests (although essay exams are also possible) by using Mathematica package and LATEXin order to make easier to the lectures the preparation of a large number of mid-term tests for a large number of students.