A Simple Method to Check the Reliability of Annual Sunspot Number in the Historical Period 1610-1847

A simple method to detect inconsistencies in low annual sunspot numbers based on the relationship between these values and the annual number of active days is described. The analysis allowed for the detection of problems in the annual sunspot number series clustered in a few specific periods and unambiguous, namely: i) before Maunder minimum, ii) the year 1652 during the Maunder minimum, iii) the year 1741 in Solar Cycle -1, and iv) the so-called "lost" solar cycle in 1790s and subsequent onset of the Dalton Minimum.Comment: 15 pages, 3 figures, to be published in Solar Physic

Multi-timescale Solar Cycles and the Possible Implications

Based on analysis of the annual averaged relative sunspot number (ASN) during 1700 -- 2009, 3 kinds of solar cycles are confirmed: the well-known 11-yr cycle (Schwabe cycle), 103-yr secular cycle (numbered as G1, G2, G3, and G4, respectively since 1700); and 51.5-yr Cycle. From similarities, an extrapolation of forthcoming solar cycles is made, and found that the solar cycle 24 will be a relative long and weak Schwabe cycle, which may reach to its apex around 2012-2014 in the vale between G3 and G4. Additionally, most Schwabe cycles are asymmetric with rapidly rising-phases and slowly decay-phases. The comparisons between ASN and the annual flare numbers with different GOES classes (C-class, M-class, X-class, and super-flare, here super-flare is defined as $\geq$ X10.0) and the annal averaged radio flux at frequency of 2.84 GHz indicate that solar flares have a tendency: the more powerful of the flare, the later it takes place after the onset of the Schwabe cycle, and most powerful flares take place in the decay phase of Schwabe cycle. Some discussions on the origin of solar cycles are presented.Comment: 8 pages, 4 figure

Joint Loop End Modeling Improves Covariance Model Based Non-coding RNA Gene Search

The effect of more detailed modeling of the interface between stem and loop in non-coding RNA hairpin structures on efficacy of covariance-model-based non-coding RNA gene search is examined. Currently, the prior probabilities of the two stem nucleotides and two loop-end nucleotides at the interface are treated the same as any other stem and loop nucleotides respectively. Laboratory thermodynamic studies show that hairpin stability is dependent on the identities of these four nucleotides, but this is not taken into account in current covariance models. It is shown that separate estimation of emission priors for these nucleotides and joint treatment of substitution probabilities for the two loop-end nucleotides leads to improved non-coding RNA gene search

Power-spectrum analysis of Super-Kamiokande solar neutrino data, taking into account asymmetry in the error estimates

The purpose of this article is to carry out a power-spectrum analysis (based on likelihood methods) of the Super-Kamiokande 5-day dataset that takes account of the asymmetry in the error estimates. Whereas the likelihood analysis involves a linear optimization procedure for symmetrical error estimates, it involves a nonlinear optimization procedure for asymmetrical error estimates. We find that for most frequencies there is little difference between the power spectra derived from analyses of symmetrized error estimates and from asymmetrical error estimates. However, this proves not to be the case for the principal peak in the power spectra, which is found at 9.43 yr-1. A likelihood analysis which allows for a "floating offset" and takes account of the start time and end time of each bin and of the flux estimate and the symmetrized error estimate leads to a power of 11.24 for this peak. A Monte Carlo analysis shows that there is a chance of only 1% of finding a peak this big or bigger in the frequency band 1 - 36 yr-1 (the widest band that avoids artificial peaks). On the other hand, an analysis that takes account of the error asymmetry leads to a peak with power 13.24 at that frequency. A Monte Carlo analysis shows that there is a chance of only 0.1% of finding a peak this big or bigger in that frequency band 1 - 36 yr-1. From this perspective, power spectrum analysis that takes account of asymmetry of the error estimates gives evidence for variability that is significant at the 99.9% level. We comment briefly on an apparent discrepancy between power spectrum analyses of the Super-Kamiokande and SNO solar neutrino experiments.Comment: 13 pages, 2 tables, 6 figure

Evaluation of regression models in metabolic physiology: predicting fluxes from isotopic data without knowledge of the pathway

This study explores the ability of regression models, with no knowledge of the underlying physiology, to estimate physiological parameters relevant for metabolism and endocrinology. Four regression models were compared: multiple linear regression (MLR), principal component regression (PCR), partial least-squares regression (PLS) and regression using artificial neural networks (ANN). The pathway of mammalian gluconeogenesis was analyzed using [U−(13)C]glucose as tracer. A set of data was simulated by randomly selecting physiologically appropriate metabolic fluxes for the 9 steps of this pathway as independent variables. The isotope labeling patterns of key intermediates in the pathway were then calculated for each set of fluxes, yielding 29 dependent variables. Two thousand sets were created, allowing independent training and test data. Regression models were asked to predict the nine fluxes, given only the 29 isotopomers. For large training sets (>50) the artificial neural network model was superior, capturing 95% of the variability in the gluconeogenic flux, whereas the three linear models captured only 75%. This reflects the ability of neural networks to capture the inherent non-linearities of the metabolic system. The effect of error in the variables and the addition of random variables to the data set was considered. Model sensitivities were used to find the isotopomers that most influenced the predicted flux values. These studies provide the first test of multivariate regression models for the analysis of isotopomer flux data. They provide insight for metabolomics and the future of isotopic tracers in metabolic research where the underlying physiology is complex or unknown

A new calibrated sunspot group series since 1749: statistics of active day fractions

Although the sunspot-number series have existed since the mid-19th century, they are still the subject of intense debate, with the largest uncertainty being related to the "calibration" of the visual acuity of individual observers in the past. Daisy-chain regression methods are applied to inter-calibrate the observers which may lead to significant bias and error accumulation. Here we present a novel method to calibrate the visual acuity of the key observers to the reference data set of Royal Greenwich Observatory sunspot groups for the period 1900-1976, using the statistics of the active-day fraction. For each observer we independently evaluate their observational thresholds [S_S] defined such that the observer is assumed to miss all of the groups with an area smaller than S_S and report all the groups larger than S_S. Next, using a Monte-Carlo method we construct, from the reference data set, a correction matrix for each observer. The correction matrices are significantly non-linear and cannot be approximated by a linear regression or proportionality. We emphasize that corrections based on a linear proportionality between annually averaged data lead to serious biases and distortions of the data. The correction matrices are applied to the original sunspot group records for each day, and finally the composite corrected series is produced for the period since 1748. The corrected series displays secular minima around 1800 (Dalton minimum) and 1900 (Gleissberg minimum), as well as the Modern grand maximum of activity in the second half of the 20th century. The uniqueness of the grand maximum is confirmed for the last 250 years. It is shown that the adoption of a linear relationship between the data of Wolf and Wolfer results in grossly inflated group numbers in the 18th and 19th centuries in some reconstructions

Ecological Complex Systems

Main aim of this topical issue is to report recent advances in noisy nonequilibrium processes useful to describe the dynamics of ecological systems and to address the mechanisms of spatio-temporal pattern formation in ecology both from the experimental and theoretical points of view. This is in order to understand the dynamical behaviour of ecological complex systems through the interplay between nonlinearity, noise, random and periodic environmental interactions. Discovering the microscopic rules and the local interactions which lead to the emergence of specific global patterns or global dynamical behaviour and the noises role in the nonlinear dynamics is an important, key aspect to understand and then to model ecological complex systems.Comment: 13 pages, Editorial of a topical issue on Ecological Complex System to appear in EPJ B, Vol. 65 (2008

Variation in Structure and Process of Care in Traumatic Brain Injury: Provider Profiles of European Neurotrauma Centers Participating in the CENTER-TBI Study.

INTRODUCTION: The strength of evidence underpinning care and treatment recommendations in traumatic brain injury (TBI) is low. Comparative effectiveness research (CER) has been proposed as a framework to provide evidence for optimal care for TBI patients. The first step in CER is to map the existing variation. The aim of current study is to quantify variation in general structural and process characteristics among centers participating in the Collaborative European NeuroTrauma Effectiveness Research in Traumatic Brain Injury (CENTER-TBI) study. METHODS: We designed a set of 11 provider profiling questionnaires with 321 questions about various aspects of TBI care, chosen based on literature and expert opinion. After pilot testing, questionnaires were disseminated to 71 centers from 20 countries participating in the CENTER-TBI study. Reliability of questionnaires was estimated by calculating a concordance rate among 5% duplicate questions. RESULTS: All 71 centers completed the questionnaires. Median concordance rate among duplicate questions was 0.85. The majority of centers were academic hospitals (n = 65, 92%), designated as a level I trauma center (n = 48, 68%) and situated in an urban location (n = 70, 99%). The availability of facilities for neuro-trauma care varied across centers; e.g. 40 (57%) had a dedicated neuro-intensive care unit (ICU), 36 (51%) had an in-hospital rehabilitation unit and the organization of the ICU was closed in 64% (n = 45) of the centers. In addition, we found wide variation in processes of care, such as the ICU admission policy and intracranial pressure monitoring policy among centers. CONCLUSION: Even among high-volume, specialized neurotrauma centers there is substantial variation in structures and processes of TBI care. This variation provides an opportunity to study effectiveness of specific aspects of TBI care and to identify best practices with CER approaches