Linear quantile mixed models

Dependent data arise in many studies. For example, children with the same parents or living in neighbouring geographic areas tend to be more alike in many characteristics than individuals chosen at random from the population at large; observations taken repeatedly on the same individual are likely to be more similar than observations from different individuals. Frequently adopted sampling designs, such as cluster, multilevel, spatial, and repeated measures (or longitudinal or panel), may induce this dependence, which the analysis of the data needs to take into due account. In a previous publication (Geraci and Bottai, Biostatistics 2007), we proposed a conditional quantile regression model for continuous responses where a random intercept was included along with fixed-coefficient predictors to account for between-subjects dependence in the context of longitudinal data analysis. Conditional on the random intercept, the response was assumed to follow an asymmetric Laplace distribution. The approach hinged upon the link existing between the minimization of weighted least absolute deviations, typically used in quantile regression, and the maximization of Laplace likelihood. As a follow up to that study, here we consider an extension of those models to more complex dependence structures in the data, which are modelled by including multiple random effects in the linear conditional quantile functions. Differently from the Gibbs sampling expectation-maximization approach proposed previously, the estimation of the fixed regression coefficients and of the random effects covariance matrix is based on a combination of Gaussian quadrature approximations and optimization algorithms. The former include Gauss-Hermite and Gauss-Laguerre quadratures for, respectively, normal and double exponential (i.e., symmetric Laplace) random effects; the latter include a gradient search algorithm and general purpose optimizers. As a result, some of the computational burden associated with large Gibbs sample sizes is avoided. We also discuss briefly an estimation approach based on generalized Clarke derivatives. Finally, a simulation study is presented and some preliminary results are shown

Simulation of UHE muons propagation for GEANT3

A simulation package for the transport of high energy muons has been developed. It has been conceived to replace the muon propagation software modules implemented in the detector simulation program GEANT3. Here we discuss the results achieved with our package and we check the agreement with numerical calculations up to 10**8 GeV.Comment: 21 pages, 6 figures, 1 Table. AMSTeX document, acknowledgments adde

MPC performance monitoring of a rigorously simulated industrial process

We address in this paper the application of a recently proposed MPC performance monitoring method to a rigorously simulated industrial process. The methodology aims at detecting possible sources of suboptimal performance of linear offset-free MPC algorithms by analysis of the prediction error sequence, discriminating between the presence of plant/model mismatch and incorrect disturbance/state estimation, and proposing for each scenario an appropriate corrective action. We focus on the applicability of the method to large-scale industrial systems, which typically comprise a block structure, devising efficient and scalable diagnosis and correction procedures. We also discuss and support the application of this method when the controlled plant shows a mild nonlinear behavior mainly associated with operating point changes. A high-fidelity dynamic simulation model of a crude distillation unit was developed in UniSim® Design and used as representative test bench. Results show the efficacy of the method and indicate possible research directions for further improvements. © IFAC

Application of a method to diagnose the source of performance degradation in MPC systems

Model Predictive Control systems may suffer from performance degradation mainly for two reasons: (i) external unmeasured disturbances are not estimated correctly, (ii) the (linear) dynamic model used by the MPC does not match (any longer) the actual process response. In this work we present the application of a method to detect when performance is not optimal, to diagnose the source of performance degradation and to propose appropriate corrections. In the simplest situation (i), optimal performance can be restored by recomputing the estimator parameters; in the other case (ii), re-identification becomes necessary. The method is based on analysis of the prediction error, i.e. the difference between the actual measured output and the corresponding model prediction, and uses three main tools: a statistical (whiteness) test on the prediction error sequence, a subspace identification method to detect the order of the input-to-prediction error system, and a nonlinear optimization algorithm to recompute optimal estimator parameters. We illustrate the effectiveness of the method on a large-scale rigorously simulated industrial process. Copyright © 2013, AIDIC Servizi S.r.l

Differential Age-Related Declines in Cardiorespiratory Fitness Between People With and Without Type 2 Diabetes Mellitus

Objective To assess the extent to which the established age-related decline in cardiorespiratory fitness (CRF) is augmented in adult men with type 2 diabetes mellitus (T2DM). Participants and Methods This study used data from the Aerobics Center Longitudinal Study, conducted between September 18, 1974, and August 3, 2006, in primarily non-Hispanic white, middle-to-upper class adults. The analyses were restricted to adult men with complete data on age, CRF, and T2DM (35,307 participants). Quantile regression models were used to estimate age-related differences in CRF, estimated using a maximal treadmill test, between persons with and without T2DM. Smoking status and birth cohort served as covariates. Results Age-related declines in CRF were observed in men with and without T2DM. For men younger than 60 years, at low-mid percentiles of the CRF distribution the magnitude of the age-related decline in CRF was significantly higher (P-values=.00, .02) in men with T2DM than in those without T2DM. At upper percentiles, the decline with age between the 2 groups was virtually identical. Significant declines in CRF in men 45 years or younger were observed only at high levels of CRF for those without T2DM and at low levels of CRF for those with T2DM (P-values .00, .04). Conclusion This study reported that men younger than 60 years with T2DM at the low-mid CRF percentiles experience an accelerated age-related decline in CRF. Men younger than 60 years with T2DM exhibiting high levels of CRF experienced a decline in CRF comparable to men without T2DM. This study highlights the importance of incorporating sufficient levels of exercise or activity to maintain high CRF in men with T2DM

Splitting neutrino masses and showering into Sky

Neutrino masses might be as light as a few time the atmospheric neutrino mass splitting. High Energy ZeV cosmic neutrinos (in Z-Showering model) might hit relic ones at each mass in different resonance energies in our nearby Universe. This non-degenerated density and energy must split UHE Z-boson secondaries (in Z-Burst model) leading to multi injection of UHECR nucleons within future extreme AUGER energy. Secondaries of Z-Burst as neutral gamma, below a few tens EeV are better surviving local GZK cut-off and they might explain recent Hires BL-Lac UHECR correlations at small angles. A different high energy resonance must lead to Glashow's anti-neutrino showers while hitting electrons in matter. In air, Glashow's anti-neutrino showers lead to collimated and directional air-showers offering a new Neutrino Astronomy. At greater energy around PeV, Tau escaping mountains and Earth and decaying in flight are effectively showering in air sky. These Horizontal showering is splitting by geomagnetic field in forked shapes. Such air-showers secondaries release amplified and beamed gamma bursts (like observed TGF), made also by muon and electron pair bundles, with their accompanying rich Cherenkov flashes. Also planet' s largest (Saturn, Jupiter) atmosphere limbs offer an ideal screen for UHE GZK and Z-burst tau neutrino, because their largest sizes. Titan thick atmosphere and small radius are optimal for discovering up-going resonant Glashow resonant showers. Earth detection of Neutrino showering by twin Magic Telescopes on top mountains, or by balloons and satellites arrays facing the limbs are the simplest and cheapest way toward UHE Neutrino Astronomy .Comment: 4 pages, 7 figures; an author's name correction and Journal Referenc