Symbolic Formulae for Linear Mixed Models

A statistical model is a mathematical representation of an often simplified or idealised data-generating process. In this paper, we focus on a particular type of statistical model, called linear mixed models (LMMs), that is widely used in many disciplines e.g.~agriculture, ecology, econometrics, psychology. Mixed models, also commonly known as multi-level, nested, hierarchical or panel data models, incorporate a combination of fixed and random effects, with LMMs being a special case. The inclusion of random effects in particular gives LMMs considerable flexibility in accounting for many types of complex correlated structures often found in data. This flexibility, however, has given rise to a number of ways by which an end-user can specify the precise form of the LMM that they wish to fit in statistical software. In this paper, we review the software design for specification of the LMM (and its special case, the linear model), focusing in particular on the use of high-level symbolic model formulae and two popular but contrasting R-packages in lme4 and asreml

Viral Load, Clinical Disease Severity and Cellular Immune Responses in Primary Varicella Zoster Virus Infection in Sri Lanka

BACKGROUND: In Sri Lanka, varicella zoster virus (VZV) is typically acquired during adulthood with significant associated disease morbidity and mortality. T cells are believed to be important in the control of VZV replication and in the prevention of reactivation. The relationship between viral load, disease severity and cellular immune responses in primary VZV infection has not been well studied. METHODOLOGY: We used IFNgamma ELISpot assays and MHC class II tetramers based on VZV gE and IE63 epitopes, together with quantitative real time PCR assays to compare the frequency and phenotype of specific T cells with virological and clinical outcomes in 34 adult Sri Lankan individuals with primary VZV infection. PRINCIPAL FINDINGS: Viral loads were found to be significantly higher in patients with moderate to severe infection compared to those with mild infection (p<0.001) and were significantly higher in those over 25 years of age (P<0.01). A significant inverse correlation was seen between the viral loads and the ex vivo IFNgamma ELISpot responses of patients (P<0.001, r = -0.85). VZV-specific CD4+ T cells expressed markers of intermediate differentiation and activation. CONCLUSIONS: Overall, these data show that increased clinical severity in Sri Lankan adults with primary VZV infection associates with higher viral load and reduced viral specific T cell responses

Identification and Characterization of Paramyosin from Cyst Wall of Metacercariae Implicated Protective Efficacy against Clonorchis sinensis Infection

Human clonorchiasis has been increasingly prevalent in recent years and results in a threat to the public health in epidemic regions, motivating current strategies of vaccines to combat Clonorchis sinensis (C. sinensis). In this study, we identified C. sinensis paramyosin (CsPmy) from the cyst wall proteins of metacercariae by proteomic approaches and characterized the expressed recombinant pET-26b-CsPmy protein (101 kDa). Bioinformatics analysis indicated that full-length sequences of paramyosin are conserved in helminthes and numerous B-cell/T-cell epitopes were predicted in amino acid sequence of CsPmy. Western blot analysis showed that CsPmy was expressed at four life stages of C. sinensis, both cyst wall proteins and soluble tegumental components could be probed by anti-CsPmy serum. Moreover, immunolocalization results revealed that CsPmy was specifically localized at cyst wall and excretory bladder of metacercaria, as well as the tegument, oral sucker and vitellarium of adult worm. Both immunoblot and immunolocalization results demonstrated that CsPmy was highly expressed at the stage of adult worm, metacercariae and cercaria, which could be supported by real-time PCR analysis. Both recombinant protein and nucleic acid of CsPmy showed strong immunogenicity in rats and induced combined Th1/Th2 immune responses, which were reflected by continuous high level of antibody titers and increased level of IgG1/IgG2a subtypes in serum. In vaccine trials, comparing with control groups, both CsPmy protein and DNA vaccine exhibited protective effect with significant worm reduction rate of 54.3% (p<0.05) and 36.1% (p<0.05), respectively. In consistence with immune responses in sera, elevated level of cytokines IFN-γ and IL-4 in splenocytes suggested that CsPmy could induce combined cellular immunity and humoral immunity in host. Taken together, CsPmy could be a promising vaccine candidate in the prevention of C. sinensis regarding its high immunogenicity and surface localization

Impact of the solvent capacity constraint on E. coli metabolism

<p>Abstract</p> <p>Background</p> <p>Obtaining quantitative predictions for cellular metabolic activities requires the identification and modeling of the physicochemical constraints that are relevant at physiological growth conditions. Molecular crowding in a cell's cytoplasm is one such potential constraint, as it limits the solvent capacity available to metabolic enzymes.</p> <p>Results</p> <p>Using a recently introduced flux balance modeling framework (FBAwMC) here we demonstrate that this constraint determines a metabolic switch in <it>E. coli </it>cells when they are shifted from low to high growth rates. The switch is characterized by a change in effective optimization strategy, the excretion of acetate at high growth rates, and a global reorganization of <it>E. coli </it>metabolic fluxes, the latter being partially confirmed by flux measurements of central metabolic reactions.</p> <p>Conclusion</p> <p>These results implicate the solvent capacity as an important physiological constraint acting on <it>E. coli </it>cells operating at high metabolic rates and for the activation of a metabolic switch when they are shifted from low to high growth rates. The relevance of this constraint in the context of both the aerobic ethanol excretion seen in fast growing yeast cells (Crabtree effect) and the aerobic glycolysis observed in rapidly dividing cancer cells (Warburg effect) should be addressed in the future.</p

Search for CP violation in D+→ϕπ+ and D+s→K0Sπ+ decays

A search for CP violation in D + → ϕπ + decays is performed using data collected in 2011 by the LHCb experiment corresponding to an integrated luminosity of 1.0 fb−1 at a centre of mass energy of 7 TeV. The CP -violating asymmetry is measured to be (−0.04 ± 0.14 ± 0.14)% for candidates with K − K + mass within 20 MeV/c 2 of the ϕ meson mass. A search for a CP -violating asymmetry that varies across the ϕ mass region of the D + → K − K + π + Dalitz plot is also performed, and no evidence for CP violation is found. In addition, the CP asymmetry in the D+s→K0Sπ+ decay is measured to be (0.61 ± 0.83 ± 0.14)%

Performance of CMS muon reconstruction in pp collision events at sqrt(s) = 7 TeV

The performance of muon reconstruction, identification, and triggering in CMS has been studied using 40 inverse picobarns of data collected in pp collisions at sqrt(s) = 7 TeV at the LHC in 2010. A few benchmark sets of selection criteria covering a wide range of physics analysis needs have been examined. For all considered selections, the efficiency to reconstruct and identify a muon with a transverse momentum pT larger than a few GeV is above 95% over the whole region of pseudorapidity covered by the CMS muon system, abs(eta) < 2.4, while the probability to misidentify a hadron as a muon is well below 1%. The efficiency to trigger on single muons with pT above a few GeV is higher than 90% over the full eta range, and typically substantially better. The overall momentum scale is measured to a precision of 0.2% with muons from Z decays. The transverse momentum resolution varies from 1% to 6% depending on pseudorapidity for muons with pT below 100 GeV and, using cosmic rays, it is shown to be better than 10% in the central region up to pT = 1 TeV. Observed distributions of all quantities are well reproduced by the Monte Carlo simulation.Comment: Replaced with published version. Added journal reference and DO

Expression of miR-21 and its targets (PTEN, PDCD4, TM1) in flat epithelial atypia of the breast in relation to ductal carcinoma in situ and invasive carcinoma

<p>Abstract</p> <p>Background</p> <p>Flat epithelial atypia (FEA) of the breast is characterised by a few layers of mildly atypical luminal epithelial cells. Genetic changes found in ductal carcinoma in situ (DCIS) and invasive ductal breast cancer (IDC) are also found in FEA, albeit at a lower concentration. So far, miRNA expression changes associated with invasive breast cancer, like miR-21, have not been studied in FEA.</p> <p>Methods</p> <p>We performed miRNA in-situ hybridization (ISH) on 15 cases with simultaneous presence of normal breast tissue, FEA and/or DCIS and 17 additional cases with IDC. Expression of the miR-21 targets PDCD4, TM1 and PTEN was investigated by immunohistochemistry.</p> <p>Results</p> <p>Two out of fifteen cases showed positive staining for miR-21 in normal breast ductal epithelium, seven out of fifteen cases were positive in the FEA component and nine out of twelve cases were positive in the DCIS component. A positive staining of miR-21 was observed in 15 of 17 IDC cases. In 12 cases all three components were present in one tissue block and an increase of miR-21 from normal breast to FEA and to DCIS was observed in five cases. In three cases the FEA component was negative, whereas the DCIS component was positive for miR-21. In three other cases, normal, FEA and DCIS components were negative for miR-21 and in the last case all three components were positive. Overall we observed a gradual increase in percentage of miR-21 positive cases from normal, to FEA, DCIS and IDC. Immunohistochemical staining for PTEN revealed no obvious changes in staining intensities in normal, FEA, DCIS and IDC. Cytoplasmic staining of PDCD4 increased from normal to IDC, whereas, the nuclear staining decreased. TM1 staining decreased from positive in normal breast to negative in most DCIS and IDC cases. In FEA, the staining pattern for TM1 was similar to normal breast tissue.</p> <p>Conclusion</p> <p>Upregulation of miR-21 from normal ductal epithelial cells of the breast to FEA, DCIS and IDC parallels morphologically defined carcinogenesis. No clear relation was observed between the staining pattern of miR-21 and its previously reported target genes.</p

Search for new physics with same-sign isolated dilepton events with jets and missing transverse energy

A search for new physics is performed in events with two same-sign isolated leptons, hadronic jets, and missing transverse energy in the final state. The analysis is based on a data sample corresponding to an integrated luminosity of 4.98 inverse femtobarns produced in pp collisions at a center-of-mass energy of 7 TeV collected by the CMS experiment at the LHC. This constitutes a factor of 140 increase in integrated luminosity over previously published results. The observed yields agree with the standard model predictions and thus no evidence for new physics is found. The observations are used to set upper limits on possible new physics contributions and to constrain supersymmetric models. To facilitate the interpretation of the data in a broader range of new physics scenarios, information on the event selection, detector response, and efficiencies is provided.Comment: Published in Physical Review Letter

Performance of CMS muon reconstruction in pp collision events at sqrt(s) = 7 TeV

The performance of muon reconstruction, identification, and triggering in CMS has been studied using 40 inverse picobarns of data collected in pp collisions at sqrt(s) = 7 TeV at the LHC in 2010. A few benchmark sets of selection criteria covering a wide range of physics analysis needs have been examined. For all considered selections, the efficiency to reconstruct and identify a muon with a transverse momentum pT larger than a few GeV is above 95% over the whole region of pseudorapidity covered by the CMS muon system, abs(eta) < 2.4, while the probability to misidentify a hadron as a muon is well below 1%. The efficiency to trigger on single muons with pT above a few GeV is higher than 90% over the full eta range, and typically substantially better. The overall momentum scale is measured to a precision of 0.2% with muons from Z decays. The transverse momentum resolution varies from 1% to 6% depending on pseudorapidity for muons with pT below 100 GeV and, using cosmic rays, it is shown to be better than 10% in the central region up to pT = 1 TeV. Observed distributions of all quantities are well reproduced by the Monte Carlo simulation.Comment: Replaced with published version. Added journal reference and DO

Azimuthal anisotropy of charged particles at high transverse momenta in PbPb collisions at sqrt(s[NN]) = 2.76 TeV

The azimuthal anisotropy of charged particles in PbPb collisions at nucleon-nucleon center-of-mass energy of 2.76 TeV is measured with the CMS detector at the LHC over an extended transverse momentum (pt) range up to approximately 60 GeV. The data cover both the low-pt region associated with hydrodynamic flow phenomena and the high-pt region where the anisotropies may reflect the path-length dependence of parton energy loss in the created medium. The anisotropy parameter (v2) of the particles is extracted by correlating charged tracks with respect to the event-plane reconstructed by using the energy deposited in forward-angle calorimeters. For the six bins of collision centrality studied, spanning the range of 0-60% most-central events, the observed v2 values are found to first increase with pt, reaching a maximum around pt = 3 GeV, and then to gradually decrease to almost zero, with the decline persisting up to at least pt = 40 GeV over the full centrality range measured.Comment: Replaced with published version. Added journal reference and DO