Event Weighted Tests for Detecting Periodicity in Photon Arrival Times

This paper treats the problem of detecting periodicity in a sequence of photon arrival times, which occurs, for example, in attempting to detect gamma-ray pulsars. A particular focus is on how auxiliary information, typically source intensity, background intensity, and incidence angles and energies associated with each photon arrival should be used to maximize the detection power. We construct a class of likelihood-based tests, score tests, which give rise to event weighting in a principled and natural way, and derive expressions quantifying the power of the tests. These results can be used to compare the efficacies of different weight functions, including cuts in energy and incidence angle. The test is targeted toward a template for the periodic lightcurve, and we quantify how deviation from that template affects the power of detection

Data-driven efficient score tests for deconvolution problems

We consider testing statistical hypotheses about densities of signals in deconvolution models. A new approach to this problem is proposed. We constructed score tests for the deconvolution with the known noise density and efficient score tests for the case of unknown density. The tests are incorporated with model selection rules to choose reasonable model dimensions automatically by the data. Consistency of the tests is proved

Effects of different energy intakes before and after calving on food intake, performance and blood hormones and metabolites in dairy cows

An experiment with 25 dairy cows was performed to investigate the effects of different energy intakes on food intake, performance and blood hormone and metabolite levels during the last 70 days of pregnancy and the first 125 days after parturition. Compared with animals fed ad libitum before parturition, cows fed only according to requirements during the same time showed no decrease of food intake at calving. Cows fed at a restricted level also showed a faster increase in food intake, a smaller energy deficiency at the onset of lactation and a smaller weight loss after parturition. Peak milk yield was also lower, but the peak was maintained for a longer period, and during the first 2 months of lactation these cows had higher concentrations of glucose and lactic acid, and lower concentrations of non-esterified fatty acids and ketone bodies in the blood. Weight losses, and decrease of milk production during the 1st weeks of lactation were higher in cows fed only 0·75 of requirements during the first 60 days of lactation than in animals whose energy intake was planned to cover energy requirements, and during the first 2 months of lactation levels of unesterified fatty acids, ketone bodies and urea were higher, whereas levels of insulin, thyroxine and triiodothyronine were lower. Endocrine and metabolic changes seen during early lactation suggested improved glucose homeostasis, diminished fat mobilization and ketogenesis in cows given reduced amounts of energy during the dry period compared to animals fed ad libitum. To prevent the development of acetonaemia an adequate energy supply during the first part of lactation is, however, of greater importanc

Finite-difference distributions for the Ginibre ensemble

The Ginibre ensemble of complex random matrices is studied. The complex valued random variable of second difference of complex energy levels is defined. For the N=3 dimensional ensemble are calculated distributions of second difference, of real and imaginary parts of second difference, as well as of its radius and of its argument (angle). For the generic N-dimensional Ginibre ensemble an exact analytical formula for second difference's distribution is derived. The comparison with real valued random variable of second difference of adjacent real valued energy levels for Gaussian orthogonal, unitary, and symplectic, ensemble of random matrices as well as for Poisson ensemble is provided.Comment: 8 pages, a number of small changes in the tex

Variable-frequency-train stimulation of skeletal muscle after spinal cord injury

Skeletal muscle, after spinal cord injury (SCI), becomes highly susceptible to fatigue. Variable-frequency trains (VFTs) enhance force in fatigued human skeletal muscle of able-bodied (AB) individuals. VFTs do this by taking advantage of the catch-like property of skeletal muscle. However, mechanisms responsible for fatigue in AB and SCI subjects may not be the same, and the efficacy of VFT stimulation after SCI is unknown. Accordingly, we tested the hypothesis that VFT stimulation would augment torque-time integral in SCI subjects. The quadriceps femoris muscle was stimulated with constant frequency trains (CFTs) (six 200 s square wave pulses separated by 70 ms) or VFTs (a train identical to the CFT, except that the first two pulses were separated by 5 ms) in SCI and AB subjects. After 180 contractions (50% duty cycle), isometric peak torque decreased 44, 56, and 67 percent, in the AB (n = 10), acute SCI (n = 10), and chronic SCI (n = 12) groups, respectively. In fatigued muscle, VFTs enhanced the torque-time integral by 18 percent in AB subjects and 6 percent in chronic SCI patients, and had no effect in acute SCI patients when compared to the corresponding CFT. The much faster rise times in SCI subjects (~80 ms vs. 120 ms in AB subjects) probably contributed to the inability of VFTs to enhance torque-time integrals in SCI patients. The results suggest that the use of VFT stimulation in patients with SCI may not be as efficacious as it is in AB persons

Patterson Function from Low-Energy Electron Diffraction Measured Intensities and Structural Discrimination

Surface Patterson Functions have been derived by direct inversion of experimental Low-Energy Electron Diffraction I-V spectra measured at multiple incident angles. The direct inversion is computationally simple and can be used to discriminate between different structural models. 1x1 YSi_2 epitaxial layers grown on Si(111) have been used to illustrate the analysis. We introduce a suitable R-factor for the Patterson Function to make the structural discrimination as objective as possible. From six competing models needed to complete the geometrical search, four could easily be discarded, achieving a very significant and useful reduction in the parameter space to be explored by standard dynamical LEED methods. The amount and quality of data needed for this analysis is discussed.Comment: 5 pages, 4 figure

The Pea Aphid Uses a Version of the Terminal System during Oviparous, but not Viviparous, Development

Background: In most species of aphid, female nymphs develop into either sexual or asexual adults depending on the length of the photoperiod to which their mothers were exposed. The progeny of these sexual and asexual females, in turn, develop in dramatically different ways. The fertilized oocytes of sexual females begin embryogenesis after being deposited on leaves (oviparous development) while the oocytes of asexual females complete embryogenesis within the mother (viviparous development). Compared with oviparous development, viviparous development involves a smaller transient oocyte surrounded by fewer somatic epithelial cells and a smaller early embryo that comprises fewer cells. To investigate whether patterning mechanisms differ between the earliest stages of the oviparous and viviparous modes of pea aphid development, we examined the expression of pea aphid orthologs of genes known to specify embryonic termini in other insects. Results: Here we show that pea aphid oviparous ovaries express torso-like in somatic posterior follicle cells and activate ERK MAP kinase at the posterior of the oocyte. In addition to suggesting that some posterior features of the terminal system are evolutionarily conserved, our detection of activated ERK in the oocyte, rather than in the embryo, suggests that pea aphids may transduce the terminal signal using a mechanism distinct from the one used in Drosophila. In contrast with oviparous development, the pea aphid version of the terminal system does not appear to be used during viviparous development, since we did not detect expression of torso-like in the somatic epithelial cells that surround either the oocyte or the blastoderm embryo and we did not observe restricted activated ERK in the oocyte. Conclusions: We suggest that while oviparous oocytes and embryos may specify posterior fate through an aphid terminal system, viviparous oocytes and embryos employ a different mechanism, perhaps one that does not rely on an interaction between the oocyte and surrounding somatic cells. Together, these observations provide a striking example of a difference in the fundamental events of early development that is both environmentally induced and encoded by the same genome

Validation of differential gene expression algorithms: Application comparing fold-change estimation to hypothesis testing

<p>Abstract</p> <p>Background</p> <p>Sustained research on the problem of determining which genes are differentially expressed on the basis of microarray data has yielded a plethora of statistical algorithms, each justified by theory, simulation, or ad hoc validation and yet differing in practical results from equally justified algorithms. Recently, a concordance method that measures agreement among gene lists have been introduced to assess various aspects of differential gene expression detection. This method has the advantage of basing its assessment solely on the results of real data analyses, but as it requires examining gene lists of given sizes, it may be unstable.</p> <p>Results</p> <p>Two methodologies for assessing predictive error are described: a cross-validation method and a posterior predictive method. As a nonparametric method of estimating prediction error from observed expression levels, cross validation provides an empirical approach to assessing algorithms for detecting differential gene expression that is fully justified for large numbers of biological replicates. Because it leverages the knowledge that only a small portion of genes are differentially expressed, the posterior predictive method is expected to provide more reliable estimates of algorithm performance, allaying concerns about limited biological replication. In practice, the posterior predictive method can assess when its approximations are valid and when they are inaccurate. Under conditions in which its approximations are valid, it corroborates the results of cross validation. Both comparison methodologies are applicable to both single-channel and dual-channel microarrays. For the data sets considered, estimating prediction error by cross validation demonstrates that empirical Bayes methods based on hierarchical models tend to outperform algorithms based on selecting genes by their fold changes or by non-hierarchical model-selection criteria. (The latter two approaches have comparable performance.) The posterior predictive assessment corroborates these findings.</p> <p>Conclusions</p> <p>Algorithms for detecting differential gene expression may be compared by estimating each algorithm's error in predicting expression ratios, whether such ratios are defined across microarray channels or between two independent groups.</p> <p>According to two distinct estimators of prediction error, algorithms using hierarchical models outperform the other algorithms of the study. The fact that fold-change shrinkage performed as well as conventional model selection criteria calls for investigating algorithms that combine the strengths of significance testing and fold-change estimation.</p

Coherent frequentism

By representing the range of fair betting odds according to a pair of confidence set estimators, dual probability measures on parameter space called frequentist posteriors secure the coherence of subjective inference without any prior distribution. The closure of the set of expected losses corresponding to the dual frequentist posteriors constrains decisions without arbitrarily forcing optimization under all circumstances. This decision theory reduces to those that maximize expected utility when the pair of frequentist posteriors is induced by an exact or approximate confidence set estimator or when an automatic reduction rule is applied to the pair. In such cases, the resulting frequentist posterior is coherent in the sense that, as a probability distribution of the parameter of interest, it satisfies the axioms of the decision-theoretic and logic-theoretic systems typically cited in support of the Bayesian posterior. Unlike the p-value, the confidence level of an interval hypothesis derived from such a measure is suitable as an estimator of the indicator of hypothesis truth since it converges in sample-space probability to 1 if the hypothesis is true or to 0 otherwise under general conditions.Comment: The confidence-measure theory of inference and decision is explicitly extended to vector parameters of interest. The derivation of upper and lower confidence levels from valid and nonconservative set estimators is formalize