Assessment of variance components in nonlinear mixed-effects elliptical models

The issue of assessing variance components is essential in deciding on the inclusion of random effects in the context of mixed models. In this work we discuss this problem by supposing nonlinear elliptical models for correlated data by using the score-type test proposed in Silvapulle and Silvapulle (1995). Being asymptotically equivalent to the likelihood ratio test and only requiring the estimation under the null hypothesis, this test provides a fairly easy computable alternative for assessing one-sided hypotheses in the context of the marginal model. Taking into account the possible non-normal distribution, we assume that the joint distribution of the response variable and the random effects lies in the elliptical class, which includes light-tailed and heavy-tailed distributions such as Student-t, power exponential, logistic, generalized Student-t, generalized logistic, contaminated normal, and the normal itself, among others. We compare the sensitivity of the score-type test under normal, Student-t and power exponential models for the kinetics data set discussed in Vonesh and Carter (1992) and fitted using the model presented in Russo et al. (2009). Also, a simulation study is performed to analyze the consequences of the kurtosis misspecification.

Influence diagnostics in heteroscedastic and/or autoregressive nonlinear elliptical models for correlated data

In this paper, we propose nonlinear elliptical models for correlated data with heteroscedastic and/or autoregressive structures. Our aim is to extend the models proposed by Russo et al. [22] by considering a more sophisticated scale structure to deal with variations in data dispersion and/or a possible autocorrelation among measurements taken throughout the same experimental unit. Moreover, to avoid the possible influence of outlying observations or to take into account the non-normal symmetric tails of the data, we assume elliptical contours for the joint distribution of random effects and errors, which allows us to attribute different weights to the observations. We propose an iterative algorithm to obtain the maximum-likelihood estimates for the parameters and derive the local influence curvatures for some specific perturbation schemes. The motivation for this work comes from a pharmacokinetic indomethacin data set, which was analysed previously by Bocheng and Xuping [1] under normality.

A theoretical model for estimating the margination constant of leukocytes

Abstract Background Blood leukocytes constitute two interchangeable sub-populations, the marginated and circulating pools. These two sub-compartments are found in normal conditions and are potentially affected by non-normal situations, either pathological or physiological. The dynamics between the compartments is governed by rate constants of margination (M) and return to circulation (R). Therefore, estimates of M and R may prove of great importance to a deeper understanding of many conditions. However, there has been a lack of formalism in order to approach such estimates. The few attempts to furnish an estimation of M and R neither rely on clearly stated models that precisely say which rate constant is under estimation nor recognize which factors may influence the estimation. Results The returning of the blood pools to a steady-state value after a perturbation (e.g., epinephrine injection) was modeled by a second-order differential equation. This equation has two eigenvalues, related to a fast- and to a slow-component of the dynamics. The model makes it possible to identify that these components are partitioned into three constants: R, M and SB; where SB is a time-invariant exit to tissues rate constant. Three examples of the computations are worked and a tentative estimation of R for mouse monocytes is presented. Conclusions This study establishes a firm theoretical basis for the estimation of the rate constants of the dynamics between the blood sub-compartments of white cells. It shows, for the first time, that the estimation must also take into account the exit to tissues rate constant, SB.

Color-space distortions following long-term occupational exposure to mercury vapor

Color vision was examined in subjects with long-term occupational exposure to mercury (Hg) vapor. The color vision impairment was assessed by employing a quantitative measure of distortion of individual and group perceptual color spaces. Hg subjects (n = 18; 42.1 +/- 6.5 years old; exposure time = 10.4 +/- 5.0 years; time away from the exposure source = 6.8 +/- 4.6 years) and controls (n = 18; 46.1 +/- 8.4 years old) were examined using two arrangement tests, D-15 and D-15d, in the traditional way, and also in a triadic procedure. From each subject`s `odd-one-out` choices, matrices of inter-cap subjective dissimilarities were derived and processed by non-metric multidimensional scaling (MDS). D-15d results differed significantly between the Hg-group and the control group (p < 0.05), with the impairment predominantly along the tritan axis. 2D perceptual color spaces, individual and group, were reconstructed, with the dimensions interpreted as the red-green (RG) and the blue-yellow (BY) systems. When color configurations from the Hg-group were compared to those of the controls, they presented more fluctuations along both chromatic dimensions, indicating a statistically significant difference along the BY axis. In conclusion, the present findings confirm that color vision impairments persist in subjects that have received long-term occupational exposure to Hg-vapor although, at the time of testing, they were presenting mean urinary concentration within the normal range for non-exposed individuals. Considering the advantages of the triadic procedure in clinical evaluation of acquired color vision deficiencies, further studies should attempt to verify and/or improve its efficacy.

Permian non-marine bivalves of the Falkland Islands and their palaeoenvironmental significance

We describe the occurrence of non-marine bivalves in exposures of the Middle Permian (Capitanian) Brenton Loch Formation on the southern shore of Choiseul Sound, East Falklands. The bivalves are associated with ichnofossils and were collected from a bed in the upper part of the formation, within a 25 cm thick interval of dark siltstones and mudstones with planar lamination, overlain by massive sandstones. The shells are articulated, with the valves either splayed open or closed. At the top of the succession, mudstone beds nearly 1.5 m above the bivalve-bearing layers yielded well-preserved Glossopteris sp. cf. G. communis leaf fossils. The closed articulated condition of some shells indicates preservation under high sedimentation rates with low residence time of bioclasts at the sediment/water interface. However, the presence of specimens with splayed shells is usually correlated to the slow decay of the shell ligament in oxygen-deficient bottom waters. The presence of complete carbonized leaves of Glossopteris associated with the bivalve-bearing levels also suggests a possibly dysoxic-anoxic bottom environment. Overall, our data suggest that the bivalves were preserved by abrupt burial, possibly by distal sediment flows into a Brenton Loch lake, and may represent autochthonous to parautochthonous fossil accumulations. The shells resemble those of anthracosiids and are herein assigned to Palaeanodonta sp. aff. P. dubia, a species also found in the Permian succession of the Karoo Basin, South Africa. Our results confirm that (a) the true distributions in space and time of all Permian non-marine (freshwater) bivalves are not yet well known, and (b) there is no evidence for marine conditions in the upper part of the Brenton Loch Formation.

GROUND STATE AND NON-GROUND STATE SOLUTIONS OF SOME STRONGLY COUPLED ELLIPTIC SYSTEMS

We study an elliptic system of the form Lu = vertical bar v vertical bar(p-1) v and Lv = vertical bar u vertical bar(q-1) u in Omega with homogeneous Dirichlet boundary condition, where Lu := -Delta u in the case of a bounded domain and Lu := -Delta u + u in the cases of an exterior domain or the whole space R-N. We analyze the existence, uniqueness, sign and radial symmetry of ground state solutions and also look for sign changing solutions of the system. More general non-linearities are also considered.

Identifying regulational alterations in gene regulatory networks by state space representation of vector autoregressive models and variational annealing

Background: In the analysis of effects by cell treatment such as drug dosing, identifying changes on gene network structures between normal and treated cells is a key task. A possible way for identifying the changes is to compare structures of networks estimated from data on normal and treated cells separately. However, this approach usually fails to estimate accurate gene networks due to the limited length of time series data and measurement noise. Thus, approaches that identify changes on regulations by using time series data on both conditions in an efficient manner are demanded. Methods: We propose a new statistical approach that is based on the state space representation of the vector autoregressive model and estimates gene networks on two different conditions in order to identify changes on regulations between the conditions. In the mathematical model of our approach, hidden binary variables are newly introduced to indicate the presence of regulations on each condition. The use of the hidden binary variables enables an efficient data usage; data on both conditions are used for commonly existing regulations, while for condition specific regulations corresponding data are only applied. Also, the similarity of networks on two conditions is automatically considered from the design of the potential function for the hidden binary variables. For the estimation of the hidden binary variables, we derive a new variational annealing method that searches the configuration of the binary variables maximizing the marginal likelihood. Results: For the performance evaluation, we use time series data from two topologically similar synthetic networks, and confirm that our proposed approach estimates commonly existing regulations as well as changes on regulations with higher coverage and precision than other existing approaches in almost all the experimental settings. For a real data application, our proposed approach is applied to time series data from normal Human lung cells and Human lung cells treated by stimulating EGF-receptors and dosing an anticancer drug termed Gefitinib. In the treated lung cells, a cancer cell condition is simulated by the stimulation of EGF-receptors, but the effect would be counteracted due to the selective inhibition of EGF-receptors by Gefitinib. However, gene expression profiles are actually different between the conditions, and the genes related to the identified changes are considered as possible off-targets of Gefitinib. Conclusions: From the synthetically generated time series data, our proposed approach can identify changes on regulations more accurately than existing methods. By applying the proposed approach to the time series data on normal and treated Human lung cells, candidates of off-target genes of Gefitinib are found. According to the published clinical information, one of the genes can be related to a factor of interstitial pneumonia, which is known as a side effect of Gefitinib.

Asymptotic integral kernel for ensembles of random normal matrices with radial potentials

The method of steepest descent is used to study the integral kernel of a family of normal random matrix ensembles with eigenvalue distribution P-N (z(1), ... , z(N)) = Z(N)(-1)e(-N)Sigma(N)(i=1) V-alpha(z(i)) Pi(1 <= i<j <= N) vertical bar z(i) - z(j)vertical bar(2), where V-alpha(z) = vertical bar z vertical bar(alpha), z epsilon C and alpha epsilon inverted left perpendicular0, infinity inverted right perpendicular. Asymptotic formulas with error estimate on sectors are obtained. A corollary of these expansions is a scaling limit for the n-point function in terms of the integral kernel for the classical Segal-Bargmann space. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.3688293]

Transiting exoplanets from the CoRoT space mission XXII. CoRoT-16b: a hot Jupiter with a hint of eccentricity around a faint solar-like star

Aims. We report the discovery of CoRoT-16b, a low density hot jupiter that orbits a faint G5V star (mV = 15.63) in 5.3523 +/- 0.0002 days with slight eccentricity. A fit of the data with no a priori assumptions on the orbit leads to an eccentricity of 0.33 +/- 0.1. We discuss this value and also derive the mass and radius of the planet. Methods. We analyse the photometric transit curve of CoRoT-16 given by the CoRoT satellite, and radial velocity data from the HARPS and HIRES spectrometers. A combined analysis using a Markov chain Monte Carlo algorithm is used to get the system parameters. Results. CoRoT-16b is a 0.535 -0.083/+0.085 M-J, 1.17 -0.14/+0.16 R-J hot Jupiter with a density of 0.44 -0.14/+0.21 g cm(-3). Despite its short orbital distance (0.0618 +/- 0.0015 AU) and the age of the parent star (6.73 +/- 2.8 Gyr), the planet orbit exhibits significantly non-zero eccentricity. This is very uncommon for this type of objects as tidal effects tend to circularise the orbit. This value is discussed taking into account the characteristics of the star and the observation accuracy.

Generalized non-commutative inflation

Non-commutative geometry indicates a deformation of the energy-momentum dispersion relation f (E) = E/pc (not equal 1) for massless particles. This distorted energy-momentum relation can affect the radiation-dominated phase of the universe at sufficiently high temperature. This prompted the idea of non-commutative inflation by Alexander et al (2003 Phys. Rev. D 67 081301) and Koh and Brandenberger (2007 JCAP06(2007) 021 and JCAP11(2007) 013). These authors studied a one-parameter family of a non-relativistic dispersion relation that leads to inflation: the a family of curves f (E) = 1 + (lambda E)(alpha). We show here how the conceptually different structure of symmetries of non-commutative spaces can lead, in a mathematically consistent way, to the fundamental equations of non-commutative inflation driven by radiation. We describe how this structure can be considered independently of (but including) the idea of non-commutative spaces as a starting point of the general inflationary deformation of SL(2, C). We analyze the conditions on the dispersion relation that leads to inflation as a set of inequalities which plays the same role as the slow-roll conditions on the potential of a scalar field. We study conditions for a possible numerical approach to obtain a general one-parameter family of dispersion relations that lead to successful inflation.

Expression of Mitochondrial Non-coding RNAs (ncRNAs) Is Modulated by High Risk Human Papillomavirus (HPV) Oncogenes

The study of RNA and DNA oncogenic viruses has proved invaluable in the discovery of key cellular pathways that are rendered dysfunctional during cancer progression. An example is high risk human papillomavirus (HPV), the etiological agent of cervical cancer. The role of HPV oncogenes in cellular immortalization and transformation has been extensively investigated. We reported the differential expression of a family of human mitochondrial non-coding RNAs (ncRNAs) between normal and cancer cells. Normal cells express a sense mitochondrial ncRNA (SncmtRNA) that seems to be required for cell proliferation and two antisense transcripts (ASncmtRNAs). In contrast, the ASncmtRNAs are down-regulated in cancer cells. To shed some light on the mechanisms that trigger down-regulation of the ASncmtRNAs, we studied human keratinocytes (HFK) immortalized with HPV. Here we show that immortalization of HFK with HPV-16 or 18 causes down-regulation of the ASncmtRNAs and induces the expression of a new sense transcript named SncmtRNA-2. Transduction of HFK with both E6 and E7 is sufficient to induce expression of SncmtRNA-2. Moreover, E2 oncogene is involved in down-regulation of the ASncmtRNAs. Knockdown of E2 in immortalized cells reestablishes in a reversible manner the expression of the ASncmtRNAs, suggesting that endogenous cellular factors(s) could play functions analogous to E2 during non-HPV-induced oncogenesis.

Comparing non-stationary and irregularly spaced time series

In this paper, we present approximate distributions for the ratio of the cumulative wavelet periodograms considering stationary and non-stationary time series generated from independent Gaussian processes. We also adapt an existing procedure to use this statistic and its approximate distribution in order to test if two regularly or irregularly spaced time series are realizations of the same generating process. Simulation studies show good size and power properties for the test statistic. An application with financial microdata illustrates the test usefulness. We conclude advocating the use of these approximate distributions instead of the ones obtained through randomizations, mainly in the case of irregular time series. (C) 2012 Elsevier B.V. All rights reserved.

Induction of chronic non-inflammatory widespread pain increases cardiac sympathetic modulation in rats

Fibromyalgia (FM) is characterized by chronic non-inflammatory widespread pain (CWP) and changes in sympathetic function. In attempt to elucidate the pathophysiological mechanisms of FM we used a well-established CWP animal model. We aimed to evaluate changes in cardiac autonomic balance and baroreflex function in response to CWP induction in rats. CWP was induced by two injections of acidic saline (pH 4.0, n = 8) five days apart into the left gastrocnemius muscle. Control animals were injected twice with normal saline (pH 7.2, n = 6). One day after the second injection of acidic saline or normal saline, the animals had pulse interval (PI) and systolic arterial pressure (SAP) variability, and spontaneous baroreflex sensitivity (BRS) evaluated. After induction of CWP, there was an increase of power in the low frequency (LF) band of PI spectrum (12.75 +/- 1.04 nu), a decrease in the high frequency (HF) band (87.25 +/- 1.04 nu) and an increase of LF/HF ratio (0.16 +/- 0.01), when compared to control animals (7.83 +/- 1.13 nu LF; 92.16 +/- 1.13 nu HF; 0.08 +/- 0.01 LF/HF). In addition, there was an increase of power in the LF band of SAP spectrum (7.93 +/- 1.39 mmHg(2)) when compared to control animals (2.97 +/- 0.61 mmHg(2)). BRS was lower in acidic saline injected rats (0.59 +/- 0.06 ms/mmHg) when compared to control animals (0.71 +/- 0.03 ms/mmHg). Our results showed that induction of CWP in rats shifts cardiac sympathovagal balance towards sympathetic predominance and decreases BRS. These data corroborate findings in humans with FM. (C) 2011 Elsevier B.V. All rights reserved.

Differences in Evolution of Children with Non-severe Acute Lower Respiratory Tract Infection With and Without Radiographically Diagnosed Pneumonia

Objective: To identify differences in the evolution of children with non-severe acute lower respiratory tract infection between those with and without radiographically diagnosed pneumonia. Design: Prospective cohort study. Setting: A public university pediatric hospital in Salvador, Northeast Brazil. Patients: Children aged 2-59 months. Methods: By active surveillance, the pneumonia cases were prospectively identified in a 2-year period. Each case was followed-up for changes in various clinical symptoms and signs. Demographic, clinical and radiographic data were recorded in standardized forms. Exclusion was due to antibiotic use in the previous 48 hours, signs of severe disease, refusal to give informed consent, underlying chronic illness, hospitalization in the previous 7 days or amoxicillin allergy. Chest X-ray (CXR) was later read by at least 2 independent pediatric radiologists. Main Outcome Measures: Radiographic diagnosed pneumonia based on agreed detection of pulmonary infiltrate or pleural effusion in 2 assessments. Results: A total of 382 patients receiving amoxicillin were studied, of whom, 372 (97.4%) had concordant radiographic diagnosis which was pneumonia (52%), normal CXR (41%). and others (7%). By multivariate analysis, age (OR=1.03; 95% CI: 1.02-1.05), disease >= 5days (OR = 1.04; 95% CI: 1.001-1.08), reduced pulmonary expansion (OR = 3.3; 95% CI: 1.4-8.0), absence of wheezing (OR = 0.5; 95% CI: 0.3-0.9), crackles on admission (OR = 2.0; 95% CI: 1.2-3.5), inability to drink on day 1 (OR = 4.2; 95% CI: 1.05-17.3), consolidation percussion sign (OR = 7.0; 95% CI: 1.5-32.3), tachypnea (OR = 2.0; 95% CI: 1.09-3.6) and fever (OR = 3.6; 95% CI: 1.4-9.4) on day 2 were independently associated with pneumonia. The highest positive predictive value was at the 2nd day of evolution for tachypnea (71.0%) and fever (81.1%). Conclusion: Persistence of fever or tachypnea up to the second day of amoxicillin treatment is predictive of radiographically diagnosed pneumonia among children with non-severe lower respiratory tract diseases.

The concept of quasi-integrability for modified non-linear Schrodinger models

We consider modifications of the nonlinear Schrodinger model (NLS) to look at the recently introduced concept of quasi-integrability. We show that such models possess an in finite number of quasi-conserved charges which present intriguing properties in relation to very specific space-time parity transformations. For the case of two-soliton solutions where the fields are eigenstates of this parity, those charges are asymptotically conserved in the scattering process of the solitons. Even though the charges vary in time their values in the far past and the far future are the same. Such results are obtained through analytical and numerical methods, and employ adaptations of algebraic techniques used in integrable field theories. Our findings may have important consequences on the applications of these models in several areas of non-linear science. We make a detailed numerical study of the modified NLS potential of the form V similar to (vertical bar psi vertical bar(2))(2+epsilon), with epsilon being a perturbation parameter. We perform numerical simulations of the scattering of solitons for this model and find a good agreement with the results predicted by the analytical considerations. Our paper shows that the quasi-integrability concepts recently proposed in the context of modifications of the sine-Gordon model remain valid for perturbations of the NLS model.