Stellar population gradients in Seyfert 2 galaxies. Northern sample

We use high signal-to-noise ratio long-slit spectra in the 3600-4700A range of the twenty brightest northern Seyfert 2 galaxies to study the variation of the stellar population properties as a function of distance from the nucleus. In order to characterize the stellar population and other continuum sources (e.g. featureless continuum FC) we have measured equivalent widths Ws of six absorption features, four continuum colours and their radial variations, and performed spectral population synthesis as a function of distance from the nucleus. About half the sample has CaIIK and G-band W values smaller at the nucleus than at 1 kpc from it, due to a younger population and/or FC. The stellar population synthesis shows that, while at the nucleus, 75% of the galaxies present contribution > 20% of ages younger or equal than 100Myr and/or of a FC, this proportion decreases to 45% at 3 kpc. In particular, 55% of the galaxies have contribution > 10% of the 3 Myr/FC component (a degenerate component in which one cannot separate what is due to a FC or to a 3 Myr stellar population) at the nucleus, but only 25% of them have this contribution at 3 kpc. As reference, the stellar population of 10 non-Seyfert galaxies, spanning the Hubble types of the Seyfert (from S0 to Sc) was also studied. A comparison between the stellar population of the Seyferts and that of the non-Seyferts shows systematic differences: the contribution of ages younger than 1 Gyr is in most cases larger in the Seyfert galaxies than in non-Seyferts, not only at the nucleus but up to 1 kpc from it.Comment: 23 pages, 18 figures, MNRAS in pres

Sample Attrition in the Presence of Population Attrition

This paper develops a method that accounts for non-ignorable sample attrition in the presence of population attrition for use with a non-representative panel sample. The method is applied to obtain attrition-correcting weights for the native and immigrant samples in the matched Current Population Survey (CPS). Of the two samples, the immigrant sample suffers from sample attrition due to changes in residence as well as population attrition caused by selective return migration. When there is population attrition, the second period cross-section is not representative of the first period population. Therefore, the existing sample attrition-correcting method developed by Hirano, Imbens, Ridder, and Rubin (2001) and Bhattacharya (2008) cannot be applied. We resolve this problem by generating a counterfactual, but representative cross-section prior to applying their method. The counterfactual sample can be obtained by weighting the second period cross-section by one minus the probability of population attrition. We show that the sample attrition and the population attrition processes are separately identified. This is useful because samples usually do not indicate which missing observations are due to sample attrition and which are due to population attrition. The attrition-correcting weights, once obtained, can be used in various studies of immigration using the CPS.

The prevalence of cam hip morphology in a general population sample

Objective Cam hip morphology is associated with femoroacetabular impingement (FAI) syndrome and causes hip osteoarthritis (OA). We aimed to assess the prevalence of cam hip morphology in a sample representative of the general population, using a measure with a predefined diagnostic accuracy. Design Patients aged 16–65, who were admitted to a major trauma centre and received a computed tomography (CT) pelvis were retrospectively screened for eligibility. Subjects with proximal femoral, acetabular or pelvic fractures and those who were deceased were excluded. Eligible subjects were divided into 10 groups based on gender and age. 20 subjects from each group were included. Subjects' index of multiple deprivation (IMD) and ethnicity were recorded. CT imaging was assessed and alpha angles (a measure of cam morphology) measured in the anterosuperior aspect of the femoral head neck junction. An alpha angle greater than 60° was considered to represent cam morphology. This measure and technique has a predefined sensitivity of 80% and specificity of 73% to detect cam morphology associated with FAI syndrome. The prevalence of cam morphology was reported as a proportion of subjects affected with 95% confidence intervals. Results 200 subjects were included. The sample was broadly representative of the UK general population in terms of IMD. 155 subjects (86%) identified as white. Cam morphology was present in 47% (95% CI 42,51) of subjects. Conclusions In this sample, broadly representative of the UK general population 47% of subjects had cam hip morphology; a hip shape associated with FAI syndrome and OA

Universality for the largest eigenvalue of sample covariance matrices with general population

This paper is aimed at deriving the universality of the largest eigenvalue of a class of high-dimensional real or complex sample covariance matrices of the form $\mathcal{W}_N=\Sigma^{1/2}XX^*\Sigma ^{1/2}$. Here, $X=(x_{ij})_{M,N}$ is an $M\times N$ random matrix with independent entries $x_{ij},1\leq i\leq M,1\leq j\leq N$ such that $\mathbb{E}x_{ij}=0$, $\mathbb{E}|x_{ij}|^2=1/N$. On dimensionality, we assume that $M=M(N)$ and $N/M\rightarrow d\in(0,\infty)$ as $N\rightarrow\infty$. For a class of general deterministic positive-definite $M\times M$ matrices $\Sigma$, under some additional assumptions on the distribution of $x_{ij}$'s, we show that the limiting behavior of the largest eigenvalue of $\mathcal{W}_N$ is universal, via pursuing a Green function comparison strategy raised in [Probab. Theory Related Fields 154 (2012) 341-407, Adv. Math. 229 (2012) 1435-1515] by Erd\H{o}s, Yau and Yin for Wigner matrices and extended by Pillai and Yin [Ann. Appl. Probab. 24 (2014) 935-1001] to sample covariance matrices in the null case ($\Sigma=I$). Consequently, in the standard complex case ($\mathbb{E}x_{ij}^2=0$), combing this universality property and the results known for Gaussian matrices obtained by El Karoui in [Ann. Probab. 35 (2007) 663-714] (nonsingular case) and Onatski in [Ann. Appl. Probab. 18 (2008) 470-490] (singular case), we show that after an appropriate normalization the largest eigenvalue of $\mathcal{W}_N$ converges weakly to the type 2 Tracy-Widom distribution $\mathrm{TW}_2$. Moreover, in the real case, we show that when $\Sigma$ is spiked with a fixed number of subcritical spikes, the type 1 Tracy-Widom limit $\mathrm{TW}_1$ holds for the normalized largest eigenvalue of $\mathcal {W}_N$, which extends a result of F\'{e}ral and P\'{e}ch\'{e} in [J. Math. Phys. 50 (2009) 073302] to the scenario of nondiagonal $\Sigma$ and more generally distributed $X$.Comment: Published in at http://dx.doi.org/10.1214/14-AOS1281 the Annals of Statistics (http://www.imstat.org/aos/) by the Institute of Mathematical Statistics (http://www.imstat.org

A Genome-Wide Association Study of Neuroticism in a Population-Based Sample

Neuroticism is a moderately heritable personality trait considered to be a risk factor for developing major depression, anxiety disorders and dementia. We performed a genome-wide association study in 2,235 participants drawn from a population-based study of neuroticism, making this the largest association study for neuroticism to date. Neuroticism was measured by the Eysenck Personality Questionnaire. After Quality Control, we analysed 430,000 autosomal SNPs together with an additional 1.2 million SNPs imputed with high quality from the Hap Map CEU samples. We found a very small effect of population stratification, corrected using one principal component, and some cryptic kinship that required no correction. NKAIN2 showed suggestive evidence of association with neuroticism as a main effect (p<10(-6)) and GPC6 showed suggestive evidence for interaction with age (p approximate to 10(-7)). We found support for one previously-reported association (PDE4D), but failed to replicate other recent reports. These results suggest common SNP variation does not strongly influence neuroticism. Our study was powered to detect almost all SNPs explaining at least 2% of heritability, and so our results effectively exclude the existence of loci having a major effect on neuroticism

Statistical guarantees for the EM algorithm: From population to sample-based analysis

We develop a general framework for proving rigorous guarantees on the performance of the EM algorithm and a variant known as gradient EM. Our analysis is divided into two parts: a treatment of these algorithms at the population level (in the limit of infinite data), followed by results that apply to updates based on a finite set of samples. First, we characterize the domain of attraction of any global maximizer of the population likelihood. This characterization is based on a novel view of the EM updates as a perturbed form of likelihood ascent, or in parallel, of the gradient EM updates as a perturbed form of standard gradient ascent. Leveraging this characterization, we then provide non-asymptotic guarantees on the EM and gradient EM algorithms when applied to a finite set of samples. We develop consequences of our general theory for three canonical examples of incomplete-data problems: mixture of Gaussians, mixture of regressions, and linear regression with covariates missing completely at random. In each case, our theory guarantees that with a suitable initialization, a relatively small number of EM (or gradient EM) steps will yield (with high probability) an estimate that is within statistical error of the MLE. We provide simulations to confirm this theoretically predicted behavior

Star formation histories and evolution of 35 brightest E+A galaxies from SDSS DR5

We pick out the 35 brightest galaxies from Goto's E+A galaxies catalogue which are selected from the Sloan Digital Sky Survey Data Release 5. The spectra of E+As are prominently characterized by the strong Balmer absorption lines but little [Oii] or H_alpha emission lines. In this work we study the stellar populations of the sample galaxies by fitting their spectra using ULySS, which is a robust full spectrum fitting method. We fit each of the sample with 1-population (a single stellar population-a SSP) and 3-population (3 SSPs) models, separately. By 1-population fits, we obtain SSP-equivalent ages and metallicities which correspond to the `luminosity-weighted' averages. By 3-population fits, we divide components into three groups in age (old stellar population-OSP, intermediate-age stellar population-ISP, and young stellar population-YSP), and then get the optimal age, metallicity and population fractions in both mass and light for OSP, ISP and YSP. During the fits, both Pegase.HR/Elodie3.1 and Vazdekis/Miles are used as two independent population models. The two models result in generally consistent conclusions as follows: for all the sample galaxies, YSPs (< 1Gyr) make important contributions to the light. However, the dominant contributors to mass are OSPs. We also reconstruct the smoothing star formation histories (SFHs) by giving star formation rate (SFR) versus evolutionary age. In addition, we fit the E+A sample and 34 randomly selected elliptical galaxies with 2-population (2 SSPs) model. We obtain the equivalent age of old components for each of the E+A sample and elliptical galaxies. By comparison, the old components of E+As are statistically much younger than those of ellipticals. From the standpoint of the stellar population age, this probably provides an evidence for the proposed evolutionary link from E+As to early-types (E/S0s).Comment: 16 pages, 9 figures, Accepted for publication on MNRA