General CMB and Primordial Trispectrum Estimation

We present trispectrum estimation methods which can be applied to general non-separable primordial and CMB trispectra. We present a general optimal estimator for the connected part of the trispectrum, for which we derive a quadratic term to incorporate the effects of inhomogeneous noise and masking. We describe a general algorithm for creating simulated maps with given arbitrary (and independent) power spectra, bispectra and trispectra. We propose a universal definition of the trispectrum parameter $T_{NL}$, so that the integrated bispectrum on the observational domain can be consistently compared between theoretical models. We define a shape function for the primordial trispectrum, together with a shape correlator and a useful parametrisation for visualizing the trispectrum. We derive separable analytic CMB solutions in the large-angle limit for constant and local models. We present separable mode decompositions which can be used to describe any primordial or CMB bispectra on their respective wavenumber or multipole domains. By extracting coefficients of these separable basis functions from an observational map, we are able to present an efficient estimator for any given theoretical model with a nonseparable trispectrum. The estimator has two manifestations, comparing the theoretical and observed coefficients at either primordial or late times. These mode decomposition methods are numerically tractable with order $l^5$ operations for the CMB estimator and approximately order $l^6$ for the general primordial estimator (reducing to order $l^3$ in both cases for a special class of models). We also demonstrate how the trispectrum can be reconstructed from observational maps using these methods.Comment: 38 pages, 9 figures. In v2 Figures 4-7 are altered slightly and some extra references are included in the bibliography. v3 matches version submitted to journal. Includes discussion of special case

Effect of folic acid supplementation in pregnancy on preeclampsia: The folic acid clinical trial study

Copyright © 2013 Shi Wu Wen et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.Preeclampsia (PE) is hypertension with proteinuria that develops during pregnancy and affects at least 5% of pregnancies. The Effect of Folic Acid Supplementation in Pregnancy on Preeclampsia: the Folic Acid Clinical Trial (FACT) aims to recruit 3,656 high risk women to evaluate a new prevention strategy for PE: supplementation of folic acid throughout pregnancy. Pregnant women with increased risk of developing PE presenting to a trial participating center between 80/7 and 166/7 weeks of gestation are randomized in a 1: 1 ratio to folic acid 4.0 mg or placebo after written consent is obtained. Intent-to-treat population will be analyzed. The FACT study was funded by the Canadian Institutes of Health Research in 2009, and regulatory approval from Health Canada was obtained in 2010. A web-based randomization system and electronic data collection system provide the platform for participating centers to randomize their eligible participants and enter data in real time. To date we have twenty participating Canadian centers, of which eighteen are actively recruiting, and seven participating Australian centers, of which two are actively recruiting. Recruitment in Argentina, UK, Netherlands, Brazil, West Indies, and United States is expected to begin by the second or third quarter of 2013. This trial is registered with NCT01355159. © 2013 Shi Wu Wen et al.The Canadian Institutes of Healt

Rapid Separable Analysis of Higher Order Correlators in Large Scale Structure

We present an efficient separable approach to the estimation and reconstruction of the bispectrum and the trispectrum from observational (or simulated) large scale structure data. This is developed from general CMB (poly-)spectra methods which exploit the fact that the bispectrum and trispectrum in the literature can be represented by a separable mode expansion which converges rapidly (with $n_\textrm{max}={\cal{O}}(30)$ terms). With an effective grid resolution $l_\textrm{max}$ (number of particles/grid points $N=l_\textrm{max}^3$), we present a bispectrum estimator which requires only ${\cal O}(n_\textrm{max} \times l_\textrm{max}^3)$ operations, along with a corresponding method for direct bispectrum reconstruction. This method is extended to the trispectrum revealing an estimator which requires only ${\cal O}(n_\textrm{max}^{4/3} \times l_\textrm{max}^3)$ operations. The complexity in calculating the trispectrum in this method is now involved in the original decomposition and orthogonalisation process which need only be performed once for each model. However, for non-diagonal trispectra these processes present little extra difficulty and may be performed in ${\cal O}(l_\textrm{max}^4)$ operations. A discussion of how the methodology may be applied to the quadspectrum is also given. An efficient algorithm for the generation of arbitrary nonGaussian initial conditions for use in N-body codes using this separable approach is described. This prescription allows for the production of nonGaussian initial conditions for arbitrary bispectra and trispectra. A brief outline of the key issues involved in parameter estimation, particularly in the non-linear regime, is also given

Universal Non-Gaussian Initial Conditions for N-body Simulations

In this paper we present the implementation of an efficient formalism for the generation of arbitrary non-Gaussian initial conditions for use in N-body simulations. The methodology involves the use of a separable modal approach for decomposing a primordial bispectrum or trispectrum. This approach allows for the far more efficient generation of the non-Gaussian initial conditions already described in the literature, as well as the generation for the first time of non-separable bispectra and the special class of diagonal-free trispectra. The modal approach also allows for the reconstruction of the spectra from given realisations, a fact which is exploited to provide an accurate consistency check of the simulations.Comment: 7 pages, 3 figure

The CMB Bispectrum

We use a separable mode expansion estimator with WMAP data to estimate the bispectrum for all the primary families of non-Gaussian models. We review the late-time mode expansion estimator methodology which can be applied to any non-separable primordial and CMB bispectrum model, and we demonstrate how the method can be used to reconstruct the CMB bispectrum from an observational map. We extend the previous validation of the general estimator using local map simulations. We apply the estimator to the coadded WMAP 5-year data, reconstructing the WMAP bispectrum using $l<500$ multipoles and $n=31$ orthonormal 3D eigenmodes. We constrain all popular nearly scale-invariant models, ensuring that the theoretical bispectrum is well-described by a convergent mode expansion. Constraints from the local model \fnl=54.4\pm 29.4 and the equilateral model \fnl=143.5\pm 151.2 (\Fnl = 25.1\pm 26.4) are consistent with previously published results. (Here, we use a nonlinearity parameter \Fnl normalised to the local case, to allow more direct comparison between different models.) Notable new constraints from our method include those for the constant model \Fnl = 35.1 \pm 27.4 , the flattened model \Fnl = 35.4\pm 29.2, and warm inflation \Fnl = 10.3\pm 27.2. We investigate feature models surveying a wide parameter range in both the scale and phase, and we find no significant evidence of non-Gaussianity in the models surveyed. We propose a measure \barFnl for the total integrated bispectrum and find that the measured value is consistent with the null hypothesis that CMB anisotropies obey Gaussian statistics. We argue that this general bispectrum survey with the WMAP data represents the best evidence for Gaussianity to date and we discuss future prospects, notably from the Planck satellite

Functional Imaging of the Outer Retinal Complex using High Fidelity Imaging Retinal Densitometry

We describe a new technique, high fidelity Imaging Retinal Densitometry (IRD), which probes the functional integrity of the outer retinal complex. We demonstrate the ability of the technique to map visual pigment optical density and synthesis rates in eyes with and without macular disease. A multispectral retinal imaging device obtained precise measurements of retinal reflectance over space and time. Data obtained from healthy controls and 5 patients with intermediate AMD, before and after photopigment bleaching, were used to quantify visual pigment metrics. Heat maps were plotted to summarise the topography of rod and cone pigment kinetics and descriptive statistics conducted to highlight differences between those with and without AMD. Rod and cone visual pigment synthesis rates in those with AMD (v = 0.043 SD 0.019 min-1 and v = 0.119 SD 0.046 min-1, respectively) were approximately half those observed in healthy controls (v = 0.079 SD 0.024 min-1 for rods and v = 0.206 SD 0.069 min-1 for cones). By mapping visual pigment kinetics across the central retina, high fidelity IRD provides a unique insight into outer retinal complex function. This new technique will improve the phenotypic characterisation, diagnosis and treatment monitoring of various ocular pathologies, including AMD

Optimal bispectrum constraints on single-field models of inflation

We use WMAP 9-year bispectrum data to constrain the free parameters of an 'effective field theory' describing fluctuations in single-field inflation. The Lagrangian of the theory contains a finite number of operators associated with unknown mass scales. Each operator produces a fixed bispectrum shape, which we decompose into partial waves in order to construct a likelihood function. Based on this likelihood we are able to constrain four linearly independent combinations of the mass scales. As an example of our framework we specialize our results to the case of 'Dirac-Born-Infeld' and 'ghost' inflation and obtain the posterior probability for each model, which in Bayesian schemes is a useful tool for model comparison. Our results suggest that DBI-like models with two or more free parameters are disfavoured by the data by comparison with single parameter models in the same class

Oscillations in the bispectrum

There exist several models of inflation that produce primordial bispectra that contain a large number of oscillations. In this paper we discuss these models, and aim at finding a method of detecting such bispectra in the data. We explain how the recently proposed method of mode expansion of bispectra might be able to reconstruct these spectra from separable basis functions. Extracting these basis functions from the data might then lead to observational constraints on these models.Comment: 6 pages, 2 figures, submitted to JOP: Conference Series, PASCOS 201

The characteristics of sexual abuse in sport: A multidimensional scaling analysis of events described in media reports

Most research on sexual abuse has been conducted within family settings (Fergusson & Mullen, 1999). In recent years, following several high profile convictions and scandals, research into sexual abuse has also encompassed institutional and community settings such as sport and the church (Gallagher, 2000; Wolfe et al., 2003). Research into sexual abuse in sport, for example, began with both prevalence studies (Kirby & Greaves, 1996; Leahy, Pretty & Tenenbaum, 2002) and qualitative analyses of the processes and experiences of athlete sexual abuse (Brackenridge, 1997; Cense & Brackenridge, 2001, Toftegaard Nielsen, 2001). From such work, descriptions of the modus operandi of abusers in sport, and the experiences and consequences for athlete victims, have been provided, informing both abuse prevention work and coach education. To date, however, no study has provided empirical support for multiple associations or identified patterns of sex offending in sport in ways that might allow comparisons with research-generated models of offending outside sport. This paper reports on an analysis of 159 cases of criminally defined sexual abuse, reported in the print media over a period of 15 years. The main aim of the study was to identify the nature of sex offending in sport focusing on the methods and locations of offences. The data were analysed using multidimensional scaling (MDS), as a data reduction method, in order to identify the underlying themes within the abuse and explore the inter-relationships of behaviour, victim and context variables. The findings indicate that there are specific themes that can be identified within the perpetrator strategies that include ‘intimate’, ‘aggressive’, and ‘’dominant’ modes of interaction. The same patterns that are described here within the specific context of sport are consistent with themes that emerge from similar behavioural analyses of rapists (Canter & Heritage, 1990; Bishopp, 2003) and child molester groups (Canter, Hughes & Kirby, 1998). These patterns show a correspondence to a broader behavioural model – the interpersonal circumplex (e.g., Leary 1957). Implications for accreditation and continuing professional education of sport psychologists are noted