The Overlooked Potential of Generalized Linear Models in Astronomy - I: Binomial Regression

Revealing hidden patterns in astronomical data is often the path to fundamental scientific breakthroughs; meanwhile the complexity of scientific inquiry increases as more subtle relationships are sought. Contemporary data analysis problems often elude the capabilities of classical statistical techniques, suggesting the use of cutting edge statistical methods. In this light, astronomers have overlooked a whole family of statistical techniques for exploratory data analysis and robust regression, the so-called Generalized Linear Models (GLMs). In this paper -- the first in a series aimed at illustrating the power of these methods in astronomical applications -- we elucidate the potential of a particular class of GLMs for handling binary/binomial data, the so-called logit and probit regression techniques, from both a maximum likelihood and a Bayesian perspective. As a case in point, we present the use of these GLMs to explore the conditions of star formation activity and metal enrichment in primordial minihaloes from cosmological hydro-simulations including detailed chemistry, gas physics, and stellar feedback. We predict that for a dark mini-halo with metallicity $\approx 1.3 \times 10^{-4} Z_{\bigodot}$, an increase of $1.2 \times 10^{-2}$ in the gas molecular fraction, increases the probability of star formation occurrence by a factor of 75%. Finally, we highlight the use of receiver operating characteristic curves as a diagnostic for binary classifiers, and ultimately we use these to demonstrate the competitive predictive performance of GLMs against the popular technique of artificial neural networks.Comment: 20 pages, 10 figures, 3 tables, accepted for publication in Astronomy and Computin

Evaluación de la maloclusión, alteraciones funcionales y hábitos orales en una población escolar: Tarragona y Barcelona

Fundamento: La etiología de la maloclusión puede ser multifactorial y difícil de clasificar durante el desarrollo del individuo. Los objetivos del estudio son conocer las características de la maloclusión y su posible relación con las alteraciones funcionales y los hábitos orales. Métodos: Se estudió a 1.270 individuos representativos de la población escolar de Cataluña; 596 niños y 674 niñas de entre 6 y 14 años. Se obtuvo información mediante registros clínicos y cuestionarios. Se trata de un estudio observacional descriptivo mixto transversal, llevado a cabo entre 2006 y 2007 a través de un muestreo no probabilística de conveniencia. Se ha realizado el calibrado de los examinadores. Resultados: Fueron incluidos 1051 escolares con una edad media de 9,32 años. El 72,8% presenta según Angle Clase I, 19,0% y 5,2% Clase II/1 y II/2 respectivamente y 2,9% Clase III. El hábito onicofagia fue el más frecuente con un 46,4%. El mayor nivel de limitación funcional fue la hipertrofia amigdalar con un 21,2%. Existe relación estadísticamente significativa entre maloclusión sagital y tipo de respiración y movilidad lingual y entre maloclusión horizontal y tiempo de succión digital y movilidad lingual (p<0,05). Conclusiones: El diagnóstico de alteraciones funcionales y de hábitos orales puede advertir de la presencia de maloclusión

A probabilistic approach to emission-line galaxy classification

We invoke a Gaussian mixture model (GMM) to jointly analyse two traditional emission-line classification schemes of galaxy ionization sources: the Baldwin-Phillips-Terlevich (BPT) and $\rm W_{H\alpha}$ vs. [NII]/H$\alpha$ (WHAN) diagrams, using spectroscopic data from the Sloan Digital Sky Survey Data Release 7 and SEAGal/STARLIGHT datasets. We apply a GMM to empirically define classes of galaxies in a three-dimensional space spanned by the $\log$ [OIII]/H$\beta$, $\log$ [NII]/H$\alpha$, and $\log$ EW(H${\alpha}$), optical parameters. The best-fit GMM based on several statistical criteria suggests a solution around four Gaussian components (GCs), which are capable to explain up to 97 per cent of the data variance. Using elements of information theory, we compare each GC to their respective astronomical counterpart. GC1 and GC4 are associated with star-forming galaxies, suggesting the need to define a new starburst subgroup. GC2 is associated with BPT's Active Galaxy Nuclei (AGN) class and WHAN's weak AGN class. GC3 is associated with BPT's composite class and WHAN's strong AGN class. Conversely, there is no statistical evidence -- based on four GCs -- for the existence of a Seyfert/LINER dichotomy in our sample. Notwithstanding, the inclusion of an additional GC5 unravels it. The GC5 appears associated to the LINER and Passive galaxies on the BPT and WHAN diagrams respectively. Subtleties aside, we demonstrate the potential of our methodology to recover/unravel different objects inside the wilderness of astronomical datasets, without lacking the ability to convey physically interpretable results. The probabilistic classifications from the GMM analysis are publicly available within the COINtoolbox (https://cointoolbox.github.io/GMM\_Catalogue/).Comment: Accepted for publication in MNRA

Exploring the spectroscopic diversity of type Ia supernovae with DRACULA: a machine learning approach

The existence of multiple subclasses of type Ia supernovae (SNeIa) has been the subject of great debate in the last decade. One major challenge inevitably met when trying to infer the existence of one or more subclasses is the time consuming, and subjective, process of subclass definition. In this work, we show how machine learning tools facilitate identification of subtypes of SNeIa through the establishment of a hierarchical group structure in the continuous space of spectral diversity formed by these objects. Using Deep Learning, we were capable of performing such identification in a 4 dimensional feature space (+1 for time evolution), while the standard Principal Component Analysis barely achieves similar results using 15 principal components. This is evidence that the progenitor system and the explosion mechanism can be described by a small number of initial physical parameters. As a proof of concept, we show that our results are in close agreement with a previously suggested classification scheme and that our proposed method can grasp the main spectral features behind the definition of such subtypes. This allows the confirmation of the velocity of lines as a first order effect in the determination of SNIa subtypes, followed by 91bg-like events. Given the expected data deluge in the forthcoming years, our proposed approach is essential to allow a quick and statistically coherent identification of SNeIa subtypes (and outliers). All tools used in this work were made publicly available in the Python package Dimensionality Reduction And Clustering for Unsupervised Learning in Astronomy (DRACULA) and can be found within COINtoolbox (https://github.com/COINtoolbox/DRACULA).Comment: 16 pages, 12 figures, accepted for publication in MNRA

On the possibility of magneto-structural correlations: detailed studies of di-nickel carboxylate complexes

A series of water-bridged dinickel complexes of the general formula [Ni&lt;sub&gt;2&lt;/sub&gt;(μ&lt;sub&gt;2&lt;/sub&gt;-OH&lt;sub&gt;2&lt;/sub&gt;)(μ2- O&lt;sub&gt;2&lt;/sub&gt;C&lt;sup&gt;t&lt;/sup&gt;Bu)&lt;sub&gt;2&lt;/sub&gt;(O&lt;sub&gt;2&lt;/sub&gt;C&lt;sup&gt;t&lt;/sup&gt;Bu)2(L)(L0)] (L = HO&lt;sub&gt;2&lt;/sub&gt;C&lt;sup&gt;t&lt;/sup&gt;Bu, L0 = HO&lt;sub&gt;2&lt;/sub&gt;C&lt;sup&gt;t&lt;/sup&gt;Bu (1), pyridine (2), 3-methylpyridine (4); L = L0 = pyridine (3), 3-methylpyridine (5)) has been synthesized and structurally characterized by X-ray crystallography. The magnetic properties have been probed by magnetometry and EPR spectroscopy, and detailed measurements show that the axial zero-field splitting, D, of the nickel(ii) ions is on the same order as the isotropic exchange interaction, J, between the nickel sites. The isotropic exchange interaction can be related to the angle between the nickel centers and the bridging water molecule, while the magnitude of D can be related to the coordination sphere at the nickel sites

Porous GaN and high-k MgO-GaN MOS diode layers grown in a single step on silicon

Porous GaN polycrystalline layers with n-type conduction characteristics were catalytically grown from Mg films formed by decomposition of a Mg2N3 precursor typically employed for activating p-type conduction in GaN. After being exposed to oxygen, the Mg film oxidized to a polycrystalline high-κ oxide between the ohmic alloy interlayer contact and the porous GaN, while maintaining a clean interface. Electrical measurements on devices coupled to composition analysis and electron microscopy of the component layers confirm that a MOS-type porous GaN diode on silicon can be formed by chemical vapor deposition in a single growth regime

Knee osteoarthritis and periarticular structure quantified by ultrasound. A case-control study

Background Assessment of pain and physical function is complex in patients with knee osteoarthritis (OA), as standard criteria are lacking.A previous study examining correlations between functional capacity and pain (WOMAC) and anthropometric characteristics and periarticular knee structure (quantified by ultrasound imaging) in females with knee OA found increased quadriceps muscle density was associated with higher functional disability and pain scores, suggesting that not only joint wear and symptom severity are involved and more objective measures are necessary..

Metal catalyzed porous n-type GaN layers: low resistivity ohmic contacting and single-step MgO/GaN diode formation

Porous GaN crystals have been successfully grown and electrically contacted simultaneously on Pt- and Au-coated silicon substrates as porous crystals and as porous layers. By the direct reaction of metallic Ga and NH3 gas through chemical vapor deposition, intermetallic metal-Ga alloys form at the GaN-metal interface, allowing vapoursolid- solid seeding and subsequent growth of porous GaN. Currentvoltage and capacitance-voltage measurements confirm that the intermetallic seed layers prevent interface oxidation and give a highquality reduced workfunction contact that allows exceptionally low contact resistivity. Additionally, the simultaneous formation of a lower workfunction intermetallic permits ohmic electron transport to n-type GaN grown using high workfunction metals that best catalyze the formation of porous GaN layers and may be employed to seed and ohmically contact a range of III-N compounds and alloys for broadband absorption and emission. Additionally, we show how a porous GaN rectifying diode can be formed by oxidatively crystallizing Mg typically employed for p-doping GaN, as a layer formed under porous structure resulting in a high-k polycrystalline MgO dielectric

Glioblastoma Therapy with Cytotoxic Mesenchymal Stromal Cells Optimized by Bioluminescence Imaging of Tumor and Therapeutic Cell Response

Genetically modified adipose tissue derived mesenchymal stromal cells (hAMSCs) with tumor homing capacity have been proposed for localized therapy of chemo- and radiotherapy resistant glioblastomas. We demonstrate an effective procedure to optimize glioblastoma therapy based on the use of genetically modified hAMSCs and in vivo non invasive monitoring of tumor and therapeutic cells. Glioblastoma U87 cells expressing Photinus pyralis luciferase (Pluc) were implanted in combination with hAMSCs expressing a trifunctional Renilla reniformis luciferase-red fluorescent protein-thymidine kinase reporter in the brains of SCID mice that were subsequently treated with ganciclovir (GCV). The resulting optimized therapy was effective and monitoring of tumor cells by bioluminescence imaging (BLI) showed that after 49 days GCV treatment reduced significantly the hAMSC treated tumors; by a factor of 104 relative to controls. Using a Pluc reporter regulated by an endothelial specific promoter and in vivo BLI to image hAMSC differentiation we gained insight on the therapeutic mechanism. Implanted hAMSCs homed to tumor vessels, where they differentiated to endothelial cells. We propose that the tumor killing efficiency of genetically modified hAMSCs results from their association with the tumor vascular system and should be useful vehicles to deliver localized therapy to glioblastoma surgical borders following tumor resection