Assessing dimensionality in dichotomous items when many subjects have all-zero responses: an example from psychiatry and a solution using mixture models

Common methods for determining the number of latent dimensions underlying an item set include eigenvalue analysis and examination of fit statistics for factor analysis models with varying number of factors. Given a set of dichotomous items, the authors demonstrate that these empirical assessments of dimensionality often incorrectly estimate the number of dimensions when there is a preponderance of individuals in the sample with all-zeros as their responses, for example, not endorsing any symptoms on a health battery. Simulated data experiments are conducted to demonstrate when each of several common diagnostics of dimensionality can be expected to under- or over-estimate the true dimensionality of the underlying latent variable. An example is shown from psychiatry assessing the dimensionality of a social anxiety disorder battery where 1, 2, 3, or more factors are identified, depending on the method of dimensionality assessment. An all-zero inflated exploratory factor analysis model (AZ-EFA) is introduced for assessing the dimensionality of the underlying subgroup corresponding to those possessing the measurable trait. The AZ-EFA approach is demonstrated using simulation experiments and an example measuring social anxiety disorder from a large nationally representative survey. Implications of the findings are discussed, in particular, regarding the potential for different findings in community versus patient populations

Temporal variability and statistics of the Strehl ratio in adaptive-optics images

We have investigated the temporal variability and statistics of the "instantaneous" Strehl ratio. The observations were carried out with the 3.63-m AEOS telescope equipped with a high-order adaptive optics system. In this paper Strehl ratio is defined as the peak intensity of a single short exposure. We have also studied the behaviour of the phase variance computed on the reconstructed wavefronts. We tested the Marechal approximation and used it to explain the observed negative skewness of the Strehl ratio distribution. The estimate of the phase variance is shown to fit a three-parameter Gamma distribution model. We show that simple scaling of the reconstructed wavefronts has a large impact on the shape of the Strehl ratio distribution.Comment: submitted to PAS

Infrared Excess and Molecular Clouds: A comparison of new suerveys of far-infrared and H I 21-cm emission at high galactic latitudes

We have created a map of the large-scale infrared surface brightness in excess of that associated with the atomic interstellar medium, using region-by-region correlations between the far-infrared and 21-cm line surface brightness. Our study updates and extends a previous attempt with IRAS and Berkeley/Parkes H I surveys. The far-infrared observations used here are from DIRBE, which extends far-infrared wavelength coverage to 240 um, so that we are reliably sampling the emission of large, thermal-equilibrium grains that dominate the dust mass. The H I data are from the combined Leiden-Dwingeloo and Parkes 21-cm line surveys. Using the maps of excess infrared emission at 100, 140, and 240 um, we created an atlas and identified the coherent structures. These infrared excess clouds can be caused both by dust that is warmer than average, or by dust associated with gas other than the atomic interstellar medium. We find very few warm clouds, such as the H II region around Spica. The majority of the infrared excess clouds are colder than the average atomic interstellar medium. These clouds are peaks of column density, and their excess infrared emission is due to dust associated with molecular gas. We identify essentially all known high-latitude molecular clouds in the infrared excess maps, and further identify a sample of new clouds with similar infrared properties. The infrared excess was correlated with CO line brightness, allowing us to measure the ratio of N(H2)/W(CO) for high-latitude clouds. The atlas of infrared excess may be a useful guide to regions of relatively high column density, which might cause high extinction toward extragalactic objects at optical and ultraviolet wavelengths and confusion for cosmic infrared and microwave backgrounds.Comment: * manuscript in 3 TeX files, 12 figures in 15 gif files submitted to the Astrophysical Journal on Feb. 12, 1998 * this version contains a new version of Figs. 7-8, with the correct cloud name

Alkaline-Silicate REE-HFSE Systems

Development of renewable energy infrastructure requires critical raw materials, such as the rare earth elements (REEs, including scandium) and niobium, and is driving expansion and diversification in their supply chains. Although alternative sources are being explored, the majority of the world’s resources of these elements are found in alkaline-silicate rocks and carbonatites. These magmatic systems also represent major sources of fluorine and phosphorus. Exploration models for critical raw materials are comparatively less well developed than those for major and precious metals, such as iron, copper, and gold, where most of the mineral exploration industry continues to focus. The diversity of lithologic relationships and a complex nomenclature for many alkaline rock types represent further barriers to the exploration and exploitation of REE-high field strength element (HFSE) resources that will facilitate the green revolution. We used a global review of maps, cross sections, and geophysical, geochemical, and petrological observations from alkaline systems to inform our description of the alkaline-silicate REE + HFSE mineral system from continental scale (1,000s km) down to deposit scale (~1 km lateral). Continental-scale targeting criteria include a geodynamic trigger for low-degree mantle melting at high pressure and a mantle source enriched in REEs, volatile elements, and alkalies. At the province and district scales, targeting criteria relate to magmatic-system longevity and the conditions required for extensive fractional crystallization and the residual enrichment of the REEs and HFSEs. A compilation of maps and geophysical data were used to construct an interactive 3-D geologic model (25-km cube) that places mineralization within a depth and horizontal reference frame. It shows typical lithologic relationships surrounding orthomagmatic REE-Nb-Ta-Zr-Hf mineralization in layered agpaitic syenites, roof zone REE-Nb-Ta mineralization, and mineralization of REE-Nb-Zr associated with peralkaline granites and pegmatites. The resulting geologic model is presented together with recommended geophysical and geochemical approaches for exploration targeting, as well as mineral processing and environmental factors pertinent for the development of mineral resources hosted by alkaline-silicate magmatic systems

The Dense Gas Mass Fraction and the Relationship to Star Formation in M51

Observations of 12CO J = 1 - 0 and HCN J = 1 - 0 emission from NGC 5194 (M51) made with the 50 m Large Millimeter Telescope and the SEQUOIA focal plane array are presented. Using the HCN-to-CO ratio, we examine the dense gas mass fraction over a range of environmental conditions within the galaxy. Within the disk, the dense gas mass fraction varies along the spiral arms but the average value over all spiral arms is comparable to the mean value of interarm regions. We suggest that the near-constant dense gas mass fraction throughout the disk arises from a population of density-stratified, self-gravitating molecular clouds and the required density threshold to detect each spectral line. The measured dense gas fraction significantly increases in the central bulge in response to the effective pressure, P e , from the weight of the stellar and gas components. This pressure modifies the dynamical state of the molecular cloud population and, possibly, the HCN-emitting regions in the central bulge from self-gravitating to diffuse configurations in which P e is greater than the gravitational energy density of individual clouds. Diffuse molecular clouds comprise a significant fraction of the molecular gas mass in the central bulge, which may account for the measured sublinear relationships between the surface densities of the star formation rate and molecular and dense gas

The 1.6 micron near infrared nuclei of 3C radio galaxies: Jets, thermal emission or scattered light?

Using HST NICMOS 2 observations we have measured 1.6-micron near infrared nuclear luminosities of 100 3CR radio galaxies with z<0.3, by modeling and subtracting the extended emission from the host galaxy. We performed a multi-wavelength statistical analysis (including optical and radio data) of the properties of the nuclei following classification of the objects into FRI and FRII, and LIG (low-ionization galaxies), HIG (high-ionization galaxies) and BLO (broad-lined objects) using the radio morphology and optical spectra, respectively. The correlations among near infrared, optical, and radio nuclear luminosity support the idea that the near infrared nuclear emission of FRIs has a non-thermal origin. Despite the difference in radio morphology, the multi-wavelength properties of FRII LIG nuclei are statistically indistinguishable from those of FRIs, an indication of a common structure of the central engine. All BLOs show an unresolved near infrared nucleus and a large near infrared excess with respect to FRII LIGs and FRIs of equal radio core luminosity. This requires the presence of an additional (and dominant) component other than the non-thermal light. Considering the shape of their spectral energy distribution, we ascribe the origin of their near infrared light to hot circumnuclear dust. A near infrared excess is also found in HIGs, but their nuclei are substantially fainter than those of BLO. This result indicates that substantial obscuration along the line-of-sight to the nuclei is still present at 1.6 micron. Nonetheless, HIGs nuclei cannot simply be explained in terms of dust obscuration: a significant contribution from light reflected in a circumnuclear scattering region is needed to account for their multiwavelength properties.Comment: 20 pages, 16 figures. Accepted for publication on Ap

Hierarchical Spherical Model from a Geometric Point of View

A continuous version of the hierarchical spherical model at dimension d=4 is investigated. Two limit distribution of the block spin variable X^{\gamma}, normalized with exponents \gamma =d+2 and \gamma =d at and above the critical temperature, are established. These results are proven by solving certain evolution equations corresponding to the renormalization group (RG) transformation of the O(N) hierarchical spin model of block size L^{d} in the limit L to 1 and N to \infty . Starting far away from the stationary Gaussian fixed point the trajectories of these dynamical system pass through two different regimes with distinguishable crossover behavior. An interpretation of this trajectories is given by the geometric theory of functions which describe precisely the motion of the Lee--Yang zeroes. The large--$N$ limit of RG transformation with L^{d} fixed equal to 2, at the criticality, has recently been investigated in both weak and strong (coupling) regimes by Watanabe \cite{W}. Although our analysis deals only with N=\infty case, it complements various aspects of that work.Comment: 27 pages, 6 figures, submitted to Journ. Stat. Phy

Population genetic structure, antibiotic resistance, capsule switching and evolution of invasive pneumococci before conjugate vaccination in Malawi

INTRODUCTION: Pneumococcal infections cause a high death toll in Sub Saharan Africa (SSA) but the recently rolled out pneumococcal conjugate vaccines (PCV) will reduce the disease burden. To better understand the population impact of these vaccines, comprehensive analysis of large collections of pneumococcal isolates sampled prior to vaccination is required. Here we present a population genomic study of the invasive pneumococcal isolates sampled before the implementation of PCV13 in Malawi. MATERIALS AND METHODS: We retrospectively sampled and whole genome sequenced 585 invasive isolates from 2004 to 2010. We determine the pneumococcal population genetic structure and assessed serotype prevalence, antibiotic resistance rates, and the occurrence of serotype switching. RESULTS: Population structure analysis revealed 22 genetically distinct sequence clusters (SCs), which consisted of closely related isolates. Serotype 1 (ST217), a vaccine-associated serotype in clade SC2, showed highest prevalence (19.3%), and was associated with the highest MDR rate (81.9%) followed by serotype 12F, a non-vaccine serotype in clade SC10 with an MDR rate of 57.9%. Prevalence of serotypes was stable prior to vaccination although there was an increase in the PMEN19 clone, serotype 5 ST289, in clade SC1 in 2010 suggesting a potential undetected local outbreak. Coalescent analysis revealed recent emergence of the SCs and there was evidence of natural capsule switching in the absence of vaccine induced selection pressure. Furthermore, majority of the highly prevalent capsule-switched isolates were associated with acquisition of vaccine-targeted capsules. CONCLUSIONS: This study provides descriptions of capsule-switched serotypes and serotypes with potential to cause serotype replacement post-vaccination such as 12F. Continued surveillance is critical to monitor these serotypes and antibiotic resistance in order to design better infection prevention and control measures such as inclusion of emerging replacement serotypes in future conjugate vaccines

Alkaline-Silicate REE-HFSE Systems

Development of renewable energy infrastructure requires critical raw materials, such as the rare earth elements (REEs, including scandium) and niobium, and is driving expansion and diversification in their supply chains. Although alternative sources are being explored, the majority of the world’s resources of these elements are found in alkaline-silicate rocks and carbonatites. These magmatic systems also represent major sources of fluorine and phosphorus. Exploration models for critical raw materials are comparatively less well developed than those for major and precious metals, such as iron, copper, and gold, where most of the mineral exploration industry continues to focus. The diversity of lithologic relationships and a complex nomenclature for many alkaline rock types represent further barriers to the exploration and exploitation of REE-high field strength element (HFSE) resources that will facilitate the green revolution. We used a global review of maps, cross sections, and geophysical, geochemical, and petrological observations from alkaline systems to inform our description of the alkaline-silicate REE + HFSE mineral system from continental scale (1,000s km) down to deposit scale (~1 km lateral). Continental-scale targeting criteria include a geodynamic trigger for low-degree mantle melting at high pressure and a mantle source enriched in REEs, volatile elements, and alkalies. At the province and district scales, targeting criteria relate to magmatic-system longevity and the conditions required for extensive fractional crystallization and the residual enrichment of the REEs and HFSEs. A compilation of maps and geophysical data were used to construct an interactive 3-D geologic model (25-km cube) that places mineralization within a depth and horizontal reference frame. It shows typical lithologic relationships surrounding orthomagmatic REE-Nb-Ta-Zr-Hf mineralization in layered agpaitic syenites, roof zone REE-Nb-Ta mineralization, and mineralization of REE-Nb-Zr associated with peralkaline granites and pegmatites. The resulting geologic model is presented together with recommended geophysical and geochemical approaches for exploration targeting, as well as mineral processing and environmental factors pertinent for the development of mineral resources hosted by alkaline-silicate magmatic systems