Bias correction in a multivariate normal regression model with general parameterization

This paper derives the second-order biases Of maximum likelihood estimates from a multivariate normal model where the mean vector and the covariance matrix have parameters in common. We show that the second order bias can always be obtained by means of ordinary weighted least-squares regressions. We conduct simulation studies which indicate that the bias correction scheme yields nearly unbiased estimators. (C) 2009 Elsevier B.V. All rights reserved.

A robust Bayesian approach to null intercept measurement error model with application to dental data

Measurement error models often arise in epidemiological and clinical research. Usually, in this set up it is assumed that the latent variable has a normal distribution. However, the normality assumption may not be always correct. Skew-normal/independent distribution is a class of asymmetric thick-tailed distributions which includes the Skew-normal distribution as a special case. In this paper, we explore the use of skew-normal/independent distribution as a robust alternative to null intercept measurement error model under a Bayesian paradigm. We assume that the random errors and the unobserved value of the covariate (latent variable) follows jointly a skew-normal/independent distribution, providing an appealing robust alternative to the routine use of symmetric normal distribution in this type of model. Specific distributions examined include univariate and multivariate versions of the skew-normal distribution, the skew-t distributions, the skew-slash distributions and the skew contaminated normal distributions. The methods developed is illustrated using a real data set from a dental clinical trial. (C) 2008 Elsevier B.V. All rights reserved.

Albendazole-praziquantel interaction in healthy volunteers: kinetic disposition, metabolism and enantioselectivity

center dot Pharmacokinetic interactions between albendazole and praziquantel are based on plasma concentrations of the enantiomeric mixture of both drugs with contradictory data, although the antiparasitic activity arises from (-)-(R)-praziquantel and (+)-albendazole sulfoxide. WHAT THIS STUDY ADDS center dot The pharmacokinetic interaction between albendazole and praziquantel is enantioselective. Praziquantel increased the plasma concentrations of (+)-albendazole sulfoxide more than those of (-)-albendazole sulfoxide and the administration of albendazole did not change the kinetic disposition of (+)-(S)-praziquantel, but increased the plasma concentration of (-)-(R)-praziquantel. AIM This study investigated the kinetic disposition, metabolism and enantioselectivity of albendazole (ABZ) and praziquantel (PZQ) administered alone and in combination to healthy volunteers. METHODS A randomized crossover study was carried out in three phases (n = 9), in which some volunteers started in phase 1 (400 mg ABZ), others in phase 2 (1500 mg PZQ), and the remaining volunteers in phase 3 (400 mg ABZ + 1500 mg PZQ). Serial blood samples were collected from 0-48 h after drug administration. Pharmacokinetic parameters were calculated using a monocompartmental model with lag time and were analyzed using the Wilcoxon test; P < 0.05. RESULTS The administration of PZQ increased the plasma concentrations of (+)-ASOX (albendazole sulphoxide) by 264% (AUC 0.99 vs. 2.59 mu g ml-1 h), (-)-ASOX by 358% (0.14 vs. 0.50 mu g ml-1 h) and albendazole sulfone (ASON) by 187% (0.17 vs. 0.32 mu g ml-1 h). The administration of ABZ did not change the kinetic disposition of (+)-(S)-PZQ (-)-(R)-4-OHPZQ or (+)-(S)-4-OHPZQ, but increased the plasma concentration of (-)-(R)-PZQ by 64.77% (AUC 0.52 vs. 0.86 mu g ml-1 h). CONCLUSIONS The pharmacokinetic interaction between ABZ and PZQ in healthy volunteers was demonstrated by the observation of increased plasma concentrations of ASON, both ASOX enantiomers and (-)-(R)-PZQ. Clinically, the combination of ABZ and PZQ may improve the therapeutic efficacy as a consequence of higher concentration of both active drugs. On the other hand, the magnitude of this elevation may represent an increased risk of side effects, requiring, certainly, reduction of the dosage. However, further studies are necessary to evaluate the efficacy and safety of this combination.

Stable Model Predictive Control for Integrating Systems with Optimizing Targets

This paper concern the development of a stable model predictive controller (MPC) to be integrated with real time optimization (RTO) in the control structure of a process system with stable and integrating outputs. The real time process optimizer produces Optimal targets for the system inputs and for Outputs that Should be dynamically implemented by the MPC controller. This paper is based oil a previous work (Comput. Chem. Eng. 2005, 29, 1089) where a nominally stable MPC was proposed for systems with the conventional control approach where only the outputs have set points. This work is also based oil the work of Gonzalez et at. (J. Process Control 2009, 19, 110) where the zone control of stable systems is studied. The new control for is obtained by defining ail extended control objective that includes input targets and zone controller the outputs. Additional decision variables are also defined to increase the set of feasible solutions to the control problem. The hard constraints resulting from the cancellation of the integrating modes Lit the end of the control horizon are softened,, and the resulting control problem is made feasible to a large class of unknown disturbances and changes of the optimizing targets. The methods are illustrated with the simulated application of the proposed,approaches to a distillation column of the oil refining industry.

Time course of haemodynamic, respiratory and inflammatory disturbances induced by experimental acute pulmonary polystyrene microembolism

Background and objective The time course of cardiopulmonary alterations after pulmonary embolism has not been clearly demonstrated and nor has the role of systemic inflammation on the pathogenesis of the disease. This study aimed to evaluate over 12 h the effects of pulmonary embolism caused by polystyrene microspheres on the haemodynamics, lung mechanics and gas exchange and on interleukin-6 production. Methods Ten large white pigs (weight 35-42 kg) had arterial and pulmonary catheters inserted and pulmonary embolism was induced in five pigs by injection of polystyrene microspheres (diameter similar to 300 mu mol l(-1)) until a value of pulmonary mean arterial pressure of twice the baseline was obtained. Five other animals received only saline. Haemodynamic and respiratory data and pressure-volume curves of the respiratory system were collected. A bronchoscopy was performed before and 12 h after embolism, when the animals were euthanized. Results The embolism group developed hypoxaemia that was not corrected with high oxygen fractions, as well as higher values of dead space, airway resistance and lower respiratory compliance levels. Acute haemodynamic alterations included pulmonary arterial hypertension with preserved systemic arterial pressure and cardiac index. These derangements persisted until the end of the experiments. The plasma interleukin-6 concentrations were similar in both groups; however, an increase in core temperature and a nonsignificant higher concentration of bronchoalveolar lavage proteins were found in the embolism group. Conclusion Acute pulmonary embolism induced by polystyrene microspheres in pigs produces a 12-h lasting hypoxaemia and a high dead space associated with high airway resistance and low compliance. There were no plasma systemic markers of inflammation, but a higher central temperature and a trend towards higher bronchoalveolar lavage proteins were found. Eur J Anaesthesiol 27:67-76 (C) 2010 European Society of Anaesthesiology.

Higgs pair production at the LHC in models with universal extra dimensions

In this Letter we study the process of gluon fusion into a pair of Higgs bosons in a model with one universal extra dimension. We find that the contributions from the extra top quark Kaluza-Klem excitations lead to a Higgs pair production cross section at the LHC that can be significantly altered compared to the Standard Model value for small values of the compactification scale. (C) 2007 Elsevier B.V. All rights reserved.

Kinetics of autocatalytic acid hydrolysis of cellulose with crystalline and amorphous fractions

The acid hydrolysis of cellulose with crystalline and amorphous fractions is analyzed on the basis of autocatalytic model with a positive feedback of acid production from the degraded biopolymer. In the condition of low acid rate production compared with hydrolysis rate, both fraction of cellulose decrease exponentially with linear and cubic time dependence, and the normalized number of scissions per cellulose chain follows a sigmoid behavior with reaction time. The model predicts that self generated acidic compounds from cellulose accelerate the degradation of the biopolymer. However, if the acidic compounds produced are volatile species, then their release under low pressure will reduce the global rate of degradation of cellulose toward its intrinsic rate value determined by the residual acid catalyst present in the starting material.

Aging and fluctuation-dissipation ratio in a nonequilibrium q-state lattice model

A generalized version of the nonequilibrium linear Glauber model with q states in d dimensions is introduced and analyzed. The model is fully symmetric, its dynamics being invariant under all permutations of the q states. Exact expressions for the two-time autocorrelation and response functions on a d-dimensional lattice are obtained. In the stationary regime, the fluctuation-dissipation theorem holds, while in the transient the aging is observed with the fluctuation-dissipation ratio leading to the value predicted for the linear Glauber model.

Transition to complete synchronization in phase-coupled oscillators with nearest neighbor coupling

We investigate synchronization in a Kuramoto-like model with nearest neighbor coupling. Upon analyzing the behavior of individual oscillators at the onset of complete synchronization, we show that the time interval between bursts in the time dependence of the frequencies of the oscillators exhibits universal scaling and blows up at the critical coupling strength. We also bring out a key mechanism that leads to phase locking. Finally, we deduce forms for the phases and frequencies at the onset of complete synchronization.

Respiratory allergy to Blomia tropicalis: Immune response in four syngeneic mouse strains and assessment of a low allergen-dose, short-term experimental model

Background: The dust mite Blomia tropicalis is an important source of aeroallergens in tropical areas. Although a mouse model for B. tropicalis extract (BtE)-induced asthma has been described, no study comparing different mouse strains in this asthma model has been reported. The relevance and reproducibility of experimental animal models of allergy depends on the genetic background of the animal, the molecular composition of the allergen and the experimental protocol. Objectives: This work had two objectives. The first was to study the anti-B. tropicalis allergic responses in different mouse strains using a short-term model of respiratory allergy to BtE. This study included the comparison of the allergic responses elicited by BtE with those elicited by ovalbumin in mice of the strain that responded better to BtE sensitization. The second objective was to investigate whether the best responder mouse strain could be used in an experimental model of allergy employing relatively low BtE doses. Methods: Groups of mice of four different syngeneic strains were sensitized subcutaneously with 100 mu g of BtE on days 0 and 7 and challenged four times intranasally, at days 8, 10, 12, and 14, with 10 mu g of BtE. A/J mice, that were the best responders to BtE sensitization, were used to compare the B. tropicalis-specific asthma experimental model with the conventional experimental model of ovalbumin (OVA)-specific asthma. A/J mice were also sensitized with a lower dose of BtE. Results: Mice of all strains had lung inflammatory-cell infiltration and increased levels of anti-BtE IgE antibodies, but these responses were significantly more intense in A/J mice than in CBA/J, BALB/c or C57BL/6J mice. Immunization of A/J mice with BtE induced a more intense airway eosinophil influx, higher levels of total IgE, similar airway hyperreactivity to methacholine but less intense mucous production, and lower levels of specific IgE, IgG1 and IgG2 antibodies than sensitization with OVA. Finally, immunization with a relatively low BtE dose (10 mu g per subcutaneous injection per mouse) was able to sensitize A/J mice, which were the best responders to high-dose BtE immunization, for the development of allergy-associated immune and lung inflammatory responses. Conclusions: The described short-term model of BtE-induced allergic lung disease is reproducible in different syngeneic mouse strains, and mice of the A/J strain was the most responsive to it. In addition, it was shown that OVA and BtE induce quantitatively different immune responses in A/J mice and that the experimental model can be set up with low amounts of BtE.

Heavy quark production in deep inelastic electron-nucleus scattering

Heavy quark production has been very well studied over the last years both theoretically and experimentally. Theory has been used to study heavy quark production in ep collisions at HERA, in pp collisions at Tevatron and RHIC, in pA and dA collisions at RHIC, and in AA collisions at CERN-SPS and RHIC. However, to the best of our knowledge, heavy quark production in eA has received almost no attention. With the possible construction of a high energy electron-ion collider, updated estimates of heavy quark production are needed. We address the subject from the perspective of saturation physics and compute the heavy quark production cross section with the dipole model. We isolate shadowing and nonlinear effects, showing their impact on the charm structure function and on the transverse momentum spectrum.

POLLING SYSTEMS WITH PARAMETER REGENERATION, THE GENERAL CASE

We consider a polling model with multiple stations, each with Poisson arrivals and a queue of infinite capacity. The service regime is exhaustive and there is Jacksonian feedback of served customers. What is new here is that when the server comes to a station it chooses the service rate and the feedback parameters at random; these remain valid during the whole stay of the server at that station. We give criteria for recurrence, transience and existence of the sth moment of the return time to the empty state for this model. This paper generalizes the model, when only two stations accept arriving jobs, which was considered in [Ann. Appl. Probab. 17 (2007) 1447-1473]. Our results are stated in terms of Lyapunov exponents for random matrices. From the recurrence criteria it can be seen that the polling model with parameter regeneration can exhibit the unusual phenomenon of null recurrence over a thick region of parameter space.

Simultaneous saccharification and ethanol fermentation of oxalic acid pretreated corncob assessed with response surface methodology

Response surface methodology was used to evaluate optimal time, temperature and oxalic acid concentration for simultaneous saccharification and fermentation (SSF) of corncob particles by Pichia stipitis CBS 6054. Fifteen different conditions for pretreatment were examined in a 2(3) full factorial design with six axial points. Temperatures ranged from 132 to 180 degrees C, time from 10 to 90 min and oxalic acid loadings from 0.01 to 0.038 g/g solids. Separate maxima were found for enzymatic saccharification and hemicellulose fermentation, respectively, with the condition for maximum saccharification being significantly more severe. Ethanol production was affected by reaction temperature more than by oxalic acid and reaction time over the ranges examined. The effect of reaction temperature was significant at a 95% confidence level in its effect on ethanol production. Oxalic acid and reaction time were statistically significant at the 90% level. The highest ethanol concentration (20 g/l) was obtained after 48 h with an ethanol volumetric production rate of 0.42 g ethanol l(-1) h(-1). The ethanol yield after SSF with P. stipitis was significantly higher than predicted by sequential saccharification and fermentation of substrate pretreated under the same condition. This was attributed to the secretion of beta-glucosidase by P. stipitis. During SSF, free extracellular beta-glucosidase activity was 1.30 pNPG U/g with P. stipitis, while saccharification without the yeast was 0.66 pNPG U/g. Published by Elsevier Ltd.

Effects of Feed Time, Organic Loading and Shock Loads in Anaerobic Whey Treatment by an AnSBBR with Circulation

The aim of this work was to investigate the effect of different feeding times (2, 4, and 6 h) and organic loading rates (3, 6 and 12 gCOD l(-1) day(-1)) on the performance of an anaerobic sequencing batch reactor containing immobilized biomass, as well as to verify the minimum amount of alkalinity that can be added to the influent. The reactor, in which mixing was achieved by recirculation of the liquid phase, was maintained at 30 +/- 1A degrees C, possessed 2.5 l reactional volume and treated 1.5 l cheese whey in 8-h cycles. Results showed that the effect of feeding time on reactor performance was more pronounced at higher values of organic loading rates (OLR). During operation at an OLR of 3 gCOD l(-1) day(-1), change in feeding time did not affect efficiency of organic matter removal from the reactor. At an OLR of 6 gCOD l(-1) day(-1), reactor efficiency improved in relation to the lower loading rate and tended to drop at longer feeding times. At an OLR of 12 gCOD l(-1) day(-1) the reactor showed to depend more on feeding time; higher feeding times resulted in a decrease in reactor efficiency. Under all conditions shock loads of 24 gCOD l(-1) day(-1) caused an increase in acids concentration in the effluent. However, despite this increase, the reactor regained stability readily and alkalinity supplied to the influent showed to be sufficient to maintain pH close to neutral during operation. Regardless of applied OLR, operation with feeding time of 2 h was which provided improved stability and rendered the process less susceptible to shock loads.

Dynamic Imaging in Electrical Impedance Tomography of the Human Chest With Online Transition Matrix Identification

One of the electrical impedance tomography objectives is to estimate the electrical resistivity distribution in a domain based only on electrical potential measurements at its boundary generated by an imposed electrical current distribution into the boundary. One of the methods used in dynamic estimation is the Kalman filter. In biomedical applications, the random walk model is frequently used as evolution model and, under this conditions, poor tracking ability of the extended Kalman filter (EKF) is achieved. An analytically developed evolution model is not feasible at this moment. The paper investigates the identification of the evolution model in parallel to the EKF and updating the evolution model with certain periodicity. The evolution model transition matrix is identified using the history of the estimated resistivity distribution obtained by a sensitivity matrix based algorithm and a Newton-Raphson algorithm. To numerically identify the linear evolution model, the Ibrahim time-domain method is used. The investigation is performed by numerical simulations of a domain with time-varying resistivity and by experimental data collected from the boundary of a human chest during normal breathing. The obtained dynamic resistivity values lie within the expected values for the tissues of a human chest. The EKF results suggest that the tracking ability is significantly improved with this approach.