Australian National University

Summary Trace-element data for mid-ocean ridge basalts (MORBs) and ocean island basalts (OIB) are used to formulate chemical systematics for oceanic basalts. The data suggest that the order of trace-element incompatibility in oceanic basalts is Cs ≈ Rb ≈ (≈ Tl) ≈ Ba(≈ W) > Th > U ≈ Nb = Ta ≈ K > La > Ce ≈ Pb > Pr (≈ Mo) ≈ Sr > P ≈ Nd (> F) > Zr = Hf ≈ Sm > Eu ≈ Sn (≈ Sb) ≈ Ti > Dy ≈ (Li) > Ho = Y > Yb. This rule works in general and suggests that the overall fractionation processes operating during magma generation and evolution are relatively simple, involving no significant change in the environment of formation for MORBs and OIBs. In detail, minor differences in element ratios correlate with the isotopic characteristics of different types of OIB components (HIMU, EM, MORB). These systematics are interpreted in terms of partial-melting conditions, variations in residual mineralogy, involvement of subducted sediment, recycling of oceanic lithosphere and processes within the low velocity zone. Niobium data indicate that the mantle sources of MORB and OIB are not exact complementary reservoirs to the continental crust. Subduction of oceanic crust or separation of refractory eclogite material from the former oceanic crust into the lower mantle appears to be required. The negative europium anomalies observed in some EM-type OIBs and the systematics of their key element ratios suggest the addition of a small amount (⩽1% or less) of subducted sediment to their mantle sources. However, a general lack of a crustal signature in OIBs indicates that sediment recycling has not been an important process in the convecting mantle, at least not in more recent times (⩽2 Ga). Upward migration of silica-undersaturated melts from the low velocity zone can generate an enriched reservoir in the continental and oceanic lithospheric mantle. We propose that the HIMU type ( eg St Helena) OIB component can be generated in this way. This enriched mantle can be re-introduced into the convective mantle by thermal erosion of the continental lithosphere and by the recycling of the enriched oceanic lithosphere back into the mantle.

A basic problem in computer vision is to understand the structure of a real world scene given several images of it. Techniques for solving this problem are taken from projective geometry and photogrammetry. Here, the authors cover the geometric principles and their algebraic representation in terms of camera projection matrices, the fundamental matrix and the trifocal tensor. The theory and methods of computation of these entities are discussed with real examples, as is their use in the reconstruction of scenes from multiple images. The new edition features an extended introduction covering the key ideas in the book (which itself has been updated with additional examples and appendices) and significant new results which have appeared since the first edition. Comprehensive background material is provided, so readers familiar with linear algebra and basic numerical methods can understand the projective geometry and estimation algorithms presented, and implement the algorithms directly from the book.

Abstract genalex is a user‐friendly cross‐platform package that runs within Microsoft Excel, enabling population genetic analyses of codominant, haploid and binary data. Allele frequency‐based analyses include heterozygosity, F statistics, Nei's genetic distance, population assignment, probabilities of identity and pairwise relatedness. Distance‐based calculations include amova , principal coordinates analysis (PCA), Mantel tests, multivariate and 2D spatial autocorrelation and twogener . More than 20 different graphs summarize data and aid exploration. Sequence and genotype data can be imported from automated sequencers, and exported to other software. Initially designed as tool for teaching, genalex 6 now offers features for researchers as well. Documentation and the program are available at http://www.anu.edu.au/BoZo/GenAlEx/

IQ-TREE (http://www.iqtree.org, last accessed February 6, 2020) is a user-friendly and widely used software package for phylogenetic inference using maximum likelihood. Since the release of version 1 in 2014, we have continuously expanded IQ-TREE to integrate a plethora of new models of sequence evolution and efficient computational approaches of phylogenetic inference to deal with genomic data. Here, we describe notable features of IQ-TREE version 2 and highlight the key advantages over other software.

We present the first public version (v0.2) of the open-source and community-developed Python package, Astropy. This package provides core astronomy-related functionality to the community, including support for domain-specific file formats such as flexible image transport system (FITS) files, Virtual Observatory (VO) tables, and common ASCII table formats, unit and physical quantity conversions, physical constants specific to astronomy, celestial coordinate and time transformations, world coordinate system (WCS) support, generalized containers for representing gridded as well as tabular data, and a framework for cosmological transformations and conversions. Significant functionality is under activedevelopment, such as a model fitting framework, VO client and server tools, and aperture and point spread function (PSF) photometry tools. The core development team is actively making additions and enhancements to the current code base, and we encourage anyone interested to participate in the development of future Astropy versions.

On September 14, 2015 at 09:50:45 UTC the two detectors of the Laser Interferometer Gravitational-Wave Observatory simultaneously observed a transient gravitational-wave signal. The signal sweeps upwards in frequency from 35 to 250 Hz with a peak gravitational-wave strain of 1.0×10(-21). It matches the waveform predicted by general relativity for the inspiral and merger of a pair of black holes and the ringdown of the resulting single black hole. The signal was observed with a matched-filter signal-to-noise ratio of 24 and a false alarm rate estimated to be less than 1 event per 203,000 years, equivalent to a significance greater than 5.1σ. The source lies at a luminosity distance of 410(-180)(+160) Mpc corresponding to a redshift z=0.09(-0.04)(+0.03). In the source frame, the initial black hole masses are 36(-4)(+5)M⊙ and 29(-4)(+4)M⊙, and the final black hole mass is 62(-4)(+4)M⊙, with 3.0(-0.5)(+0.5)M⊙c(2) radiated in gravitational waves. All uncertainties define 90% credible intervals. These observations demonstrate the existence of binary stellar-mass black hole systems. This is the first direct detection of gravitational waves and the first observation of a binary black hole merger.

SUMMARY: GenAlEx: Genetic Analysis in Excel is a cross-platform package for population genetic analyses that runs within Microsoft Excel. GenAlEx offers analysis of diploid codominant, haploid and binary genetic loci and DNA sequences. Both frequency-based (F-statistics, heterozygosity, HWE, population assignment, relatedness) and distance-based (AMOVA, PCoA, Mantel tests, multivariate spatial autocorrelation) analyses are provided. New features include calculation of new estimators of population structure: G'(ST), G''(ST), Jost's D(est) and F'(ST) through AMOVA, Shannon Information analysis, linkage disequilibrium analysis for biallelic data and novel heterogeneity tests for spatial autocorrelation analysis. Export to more than 30 other data formats is provided. Teaching tutorials and expanded step-by-step output options are included. The comprehensive guide has been fully revised. AVAILABILITY AND IMPLEMENTATION: GenAlEx is written in VBA and provided as a Microsoft Excel Add-in (compatible with Excel 2003, 2007, 2010 on PC; Excel 2004, 2011 on Macintosh). GenAlEx, and supporting documentation and tutorials are freely available at: http://biology.anu.edu.au/GenAlEx. CONTACT: rod.peakall@anu.edu.au.

Abstract The initial stimulus for the theorizing presented here was provided by certain experimental investigations of intergroup behavior. The laboratory analog of real-world ethnocentrism is in-group bias—that is, the tendency to favor the in-group over the out-group in evaluations and behavior. Not only are incompatible group interests not always sufficient to generate conflicts but there is a good deal of experimental evidence that these conditions are not always necessary for the development of competition and discrimination between groups (for example, Ferguson & Kelley, 1964; Rabbie & Wilkens, 1971; Doise & Sinclair, 1973; Doise & Weinberger, 1973).

The planetary boundaries framework defines a safe operating space for humanity based on the intrinsic biophysical processes that regulate the stability of the Earth system. Here, we revise and update the planetary boundary framework, with a focus on the underpinning biophysical science, based on targeted input from expert research communities and on more general scientific advances over the past 5 years. Several of the boundaries now have a two-tier approach, reflecting the importance of cross-scale interactions and the regional-level heterogeneity of the processes that underpin the boundaries. Two core boundaries—climate change and biosphere integrity—have been identified, each of which has the potential on its own to drive the Earth system into a new state should they be substantially and persistently transgressed.

On August 17, 2017 at 12∶41:04 UTC the Advanced LIGO and Advanced Virgo gravitational-wave detectors made their first observation of a binary neutron star inspiral. The signal, GW170817, was detected with a combined signal-to-noise ratio of 32.4 and a false-alarm-rate estimate of less than one per <a:math xmlns:a="http://www.w3.org/1998/Math/MathML" display="inline"><a:mrow><a:mrow><a:mn>8.0</a:mn><a:mo>×</a:mo><a:msup><a:mrow><a:mn>10</a:mn></a:mrow><a:mrow><a:mn>4</a:mn></a:mrow></a:msup></a:mrow><a:mtext> </a:mtext><a:mtext> </a:mtext><a:mi>years</a:mi></a:mrow></a:math>. We infer the component masses of the binary to be between 0.86 and <c:math xmlns:c="http://www.w3.org/1998/Math/MathML" display="inline"><c:mrow><c:mn>2.26</c:mn><c:mtext> </c:mtext><c:mtext> </c:mtext><c:msub><c:mrow><c:mi>M</c:mi></c:mrow><c:mrow><c:mo stretchy="false">⊙</c:mo></c:mrow></c:msub></c:mrow></c:math>, in agreement with masses of known neutron stars. Restricting the component spins to the range inferred in binary neutron stars, we find the component masses to be in the range <f:math xmlns:f="http://www.w3.org/1998/Math/MathML" display="inline"><f:mrow><f:mn>1.17</f:mn><f:mi>–</f:mi><f:mn>1.60</f:mn><f:mtext> </f:mtext><f:mtext> </f:mtext><f:msub><f:mrow><f:mi>M</f:mi></f:mrow><f:mrow><f:mo stretchy="false">⊙</f:mo></f:mrow></f:msub></f:mrow></f:math>, with the total mass of the system <i:math xmlns:i="http://www.w3.org/1998/Math/MathML" display="inline"><i:mrow><i:mn>2.7</i:mn><i:msubsup><i:mrow><i:mn>4</i:mn></i:mrow><i:mrow><i:mo>−</i:mo><i:mn>0.01</i:mn></i:mrow><i:mrow><i:mo>+</i:mo><i:mn>0.04</i:mn></i:mrow></i:msubsup><i:msub><i:mrow><i:mi>M</i:mi></i:mrow><i:mrow><i:mo stretchy="false">⊙</i:mo></i:mrow></i:msub></i:mrow></i:math>. The source was localized within a sky region of <l:math xmlns:l="http://www.w3.org/1998/Math/MathML" display="inline"><l:mrow><l:mn>28</l:mn><l:mtext> </l:mtext><l:mtext> </l:mtext><l:mrow><l:msup><l:mrow><l:mi>deg</l:mi></l:mrow><l:mrow><l:mn>2</l:mn></l:mrow></l:msup></l:mrow></l:mrow></l:math> (90% probability) and had a luminosity distance of <n:math xmlns:n="http://www.w3.org/1998/Math/MathML" display="inline"><n:mrow><n:mrow><n:mn>4</n:mn><n:msubsup><n:mrow><n:mn>0</n:mn></n:mrow><n:mrow><n:mo>−</n:mo><n:mn>14</n:mn></n:mrow><n:mrow><n:mo>+</n:mo><n:mn>8</n:mn></n:mrow></n:msubsup><n:mtext> </n:mtext><n:mtext> </n:mtext></n:mrow><n:mrow><n:mi>Mpc</n:mi></n:mrow></n:mrow></n:math>, the closest and most precisely localized gravitational-wave signal yet. The association with the <p:math xmlns:p="http://www.w3.org/1998/Math/MathML" display="inline"><p:mi>γ</p:mi></p:math>-ray burst GRB 170817A, detected by Fermi-GBM 1.7 s after the coalescence, corroborates the hypothesis of a neutron star merger and provides the first direct evidence of a link between these mergers and short <r:math xmlns:r="http://www.w3.org/1998/Math/MathML" display="inline"><r:mi>γ</r:mi></r:math>-ray bursts. Subsequent identification of transient counterparts across the electromagnetic spectrum in the same location further supports the interpretation of this event as a neutron star merger. This unprecedented joint gravitational and electromagnetic observation provides insight into astrophysics, dense matter, gravitation, and cosmology. Published by the American Physical Society 2017

Journal Article The Lagrange Multiplier Test and its Applications to Model Specification in Econometrics Get access T. S. Breusch, T. S. Breusch University of Southampton Search for other works by this author on: Oxford Academic Google Scholar A. R. Pagan A. R. Pagan Australian National University Search for other works by this author on: Oxford Academic Google Scholar The Review of Economic Studies, Volume 47, Issue 1, 1980, Pages 239–253, https://doi.org/10.2307/2297111 Published: 01 January 1980

Context. We present the second Gaia data release, Gaia DR2, consisting of astrometry, photometry, radial velocities, and information on astrophysical parameters and variability, for sources brighter than magnitude 21. In addition epoch astrometry and photometry are provided for a modest sample of minor planets in the solar system. Aims. A summary of the contents of Gaia DR2 is presented, accompanied by a discussion on the differences with respect to Gaia DR1 and an overview of the main limitations which are still present in the survey. Recommendations are made on the responsible use of Gaia DR2 results. Methods. The raw data collected with the Gaia instruments during the first 22 months of the mission have been processed by the Gaia Data Processing and Analysis Consortium (DPAC) and turned into this second data release, which represents a major advance with respect to Gaia DR1 in terms of completeness, performance, and richness of the data products. Results. Gaia DR2 contains celestial positions and the apparent brightness in G for approximately 1.7 billion sources. For 1.3 billion of those sources, parallaxes and proper motions are in addition available. The sample of sources for which variability information is provided is expanded to 0.5 million stars. This data release contains four new elements: broad-band colour information in the form of the apparent brightness in the G BP (330–680 nm) and G RP (630–1050 nm) bands is available for 1.4 billion sources; median radial velocities for some 7 million sources are presented; for between 77 and 161 million sources estimates are provided of the stellar effective temperature, extinction, reddening, and radius and luminosity; and for a pre-selected list of 14 000 minor planets in the solar system epoch astrometry and photometry are presented. Finally, Gaia DR2 also represents a new materialisation of the celestial reference frame in the optical, the Gaia -CRF2, which is the first optical reference frame based solely on extragalactic sources. There are notable changes in the photometric system and the catalogue source list with respect to Gaia DR1, and we stress the need to consider the two data releases as independent. Conclusions. Gaia DR2 represents a major achievement for the Gaia mission, delivering on the long standing promise to provide parallaxes and proper motions for over 1 billion stars, and representing a first step in the availability of complementary radial velocity and source astrophysical information for a sample of stars in the Gaia survey which covers a very substantial fraction of the volume of our galaxy.

Rockström, J., W. Steffen, K. Noone, Å. Persson, F. S. Chapin, III, E. Lambin, T. M. Lenton, M. Scheffer, C. Folke, H. Schellnhuber, B. Nykvist, C. A. De Wit, T. Hughes, S. van der Leeuw, H. Rodhe, S. Sörlin, P. K. Snyder, R. Costanza, U. Svedin, M. Falkenmark, L. Karlberg, R. W. Corell, V. J. Fabry, J. Hansen, B. Walker, D. Liverman, K. Richardson, P. Crutzen, and J. Foley. 2009. Planetary boundaries:exploring the safe operating space for humanity. Ecology and Society 14(2): 32. https://doi.org/10.5751/ES-03180-140232

My first exposure to Support Vector Machines came this spring when heard Sue Dumais present impressive results on text categorization using this analysis technique. This issue's collection of essays should help familiarize our readers with this interesting new racehorse in the Machine Learning stable. Bernhard Scholkopf, in an introductory overview, points out that a particular advantage of SVMs over other learning algorithms is that it can be analyzed theoretically using concepts from computational learning theory, and at the same time can achieve good performance when applied to real problems. Examples of these real-world applications are provided by Sue Dumais, who describes the aforementioned text-categorization problem, yielding the best results to date on the Reuters collection, and Edgar Osuna, who presents strong results on application to face detection. Our fourth author, John Platt, gives us a practical guide and a new technique for implementing the algorithm efficiently.

Scaling factors for obtaining fundamental vibrational frequencies, low-frequency vibrations, zero-point vibrational energies (ZPVE), and thermal contributions to enthalpy and entropy from harmonic frequencies determined at 19 levels of theory have been derived through a least-squares approach. Semiempirical methods (AM1 and PM3), conventional uncorrelated and correlated ab initio molecular orbital procedures [Hartree−Fock (HF), Møller−Plesset (MP2), and quadratic configuration interaction including single and double substitutions (QCISD)], and several variants of density functional theory (DFT: B-LYP, B-P86, B3-LYP, B3-P86, and B3-PW91) have been examined in conjunction with the 3-21G, 6-31G(d), 6-31+G(d), 6-31G(d,p), 6-311G(d,p), and 6-311G(df,p) basis sets. The scaling factors for the theoretical harmonic vibrational frequencies were determined by a comparison with the corresponding experimental fundamentals utilizing a total of 1066 individual vibrations. Scaling factors suitable for low-frequency vibrations were obtained from least-squares fits of inverse frequencies. ZPVE scaling factors were obtained from a comparison of the computed ZPVEs (derived from theoretically determined harmonic vibrational frequencies) with ZPVEs determined from experimental harmonic frequencies and anharmonicity corrections for a set of 39 molecules. Finally, scaling factors for theoretical frequencies that are applicable for the computation of thermal contributions to enthalpy and entropy have been derived. A complete set of recommended scale factors is presented. The most successful procedures overall are B3-PW91/6-31G(d), B3-LYP/6-31G(d), and HF/6-31G(d).

Ecological extinction caused by overfishing precedes all other pervasive human disturbance to coastal ecosystems, including pollution, degradation of water quality, and anthropogenic climate change. Historical abundances of large consumer species were fantastically large in comparison with recent observations. Paleoecological, archaeological, and historical data show that time lags of decades to centuries occurred between the onset of overfishing and consequent changes in ecological communities, because unfished species of similar trophic level assumed the ecological roles of overfished species until they too were overfished or died of epidemic diseases related to overcrowding. Retrospective data not only help to clarify underlying causes and rates of ecological change, but they also demonstrate achievable goals for restoration and management of coastal ecosystems that could not even be contemplated based on the limited perspective of recent observations alone.

BACKGROUND: Taxonomic classification of marker-gene sequences is an important step in microbiome analysis. RESULTS: We present q2-feature-classifier ( https://github.com/qiime2/q2-feature-classifier ), a QIIME 2 plugin containing several novel machine-learning and alignment-based methods for taxonomy classification. We evaluated and optimized several commonly used classification methods implemented in QIIME 1 (RDP, BLAST, UCLUST, and SortMeRNA) and several new methods implemented in QIIME 2 (a scikit-learn naive Bayes machine-learning classifier, and alignment-based taxonomy consensus methods based on VSEARCH, and BLAST+) for classification of bacterial 16S rRNA and fungal ITS marker-gene amplicon sequence data. The naive-Bayes, BLAST+-based, and VSEARCH-based classifiers implemented in QIIME 2 meet or exceed the species-level accuracy of other commonly used methods designed for classification of marker gene sequences that were evaluated in this work. These evaluations, based on 19 mock communities and error-free sequence simulations, including classification of simulated "novel" marker-gene sequences, are available in our extensible benchmarking framework, tax-credit ( https://github.com/caporaso-lab/tax-credit-data ). CONCLUSIONS: Our results illustrate the importance of parameter tuning for optimizing classifier performance, and we make recommendations regarding parameter choices for these classifiers under a range of standard operating conditions. q2-feature-classifier and tax-credit are both free, open-source, BSD-licensed packages available on GitHub.

Suppose you are given some data set drawn from an underlying probability distribution P and you want to estimate a "simple" subset S of input space such that the probability that a test point drawn from P lies outside of S equals some a priori specified value between 0 and 1. We propose a method to approach this problem by trying to estimate a function f that is positive on S and negative on the complement. The functional form of f is given by a kernel expansion in terms of a potentially small subset of the training data; it is regularized by controlling the length of the weight vector in an associated feature space. The expansion coefficients are found by solving a quadratic programming problem, which we do by carrying out sequential optimization over pairs of input patterns. We also provide a theoretical analysis of the statistical performance of our algorithm. The algorithm is a natural extension of the support vector algorithm to the case of unlabeled data.

IMPORTANCE: Limited information exists about the epidemiology, recognition, management, and outcomes of patients with the acute respiratory distress syndrome (ARDS). OBJECTIVES: To evaluate intensive care unit (ICU) incidence and outcome of ARDS and to assess clinician recognition, ventilation management, and use of adjuncts-for example prone positioning-in routine clinical practice for patients fulfilling the ARDS Berlin Definition. DESIGN, SETTING, AND PARTICIPANTS: The Large Observational Study to Understand the Global Impact of Severe Acute Respiratory Failure (LUNG SAFE) was an international, multicenter, prospective cohort study of patients undergoing invasive or noninvasive ventilation, conducted during 4 consecutive weeks in the winter of 2014 in a convenience sample of 459 ICUs from 50 countries across 5 continents. EXPOSURES: Acute respiratory distress syndrome. MAIN OUTCOMES AND MEASURES: The primary outcome was ICU incidence of ARDS. Secondary outcomes included assessment of clinician recognition of ARDS, the application of ventilatory management, the use of adjunctive interventions in routine clinical practice, and clinical outcomes from ARDS. RESULTS: Of 29,144 patients admitted to participating ICUs, 3022 (10.4%) fulfilled ARDS criteria. Of these, 2377 patients developed ARDS in the first 48 hours and whose respiratory failure was managed with invasive mechanical ventilation. The period prevalence of mild ARDS was 30.0% (95% CI, 28.2%-31.9%); of moderate ARDS, 46.6% (95% CI, 44.5%-48.6%); and of severe ARDS, 23.4% (95% CI, 21.7%-25.2%). ARDS represented 0.42 cases per ICU bed over 4 weeks and represented 10.4% (95% CI, 10.0%-10.7%) of ICU admissions and 23.4% of patients requiring mechanical ventilation. Clinical recognition of ARDS ranged from 51.3% (95% CI, 47.5%-55.0%) in mild to 78.5% (95% CI, 74.8%-81.8%) in severe ARDS. Less than two-thirds of patients with ARDS received a tidal volume 8 of mL/kg or less of predicted body weight. Plateau pressure was measured in 40.1% (95% CI, 38.2-42.1), whereas 82.6% (95% CI, 81.0%-84.1%) received a positive end-expository pressure (PEEP) of less than 12 cm H2O. Prone positioning was used in 16.3% (95% CI, 13.7%-19.2%) of patients with severe ARDS. Clinician recognition of ARDS was associated with higher PEEP, greater use of neuromuscular blockade, and prone positioning. Hospital mortality was 34.9% (95% CI, 31.4%-38.5%) for those with mild, 40.3% (95% CI, 37.4%-43.3%) for those with moderate, and 46.1% (95% CI, 41.9%-50.4%) for those with severe ARDS. CONCLUSIONS AND RELEVANCE: Among ICUs in 50 countries, the period prevalence of ARDS was 10.4% of ICU admissions. This syndrome appeared to be underrecognized and undertreated and associated with a high mortality rate. These findings indicate the potential for improvement in the management of patients with ARDS. TRIAL REGISTRATION: clinicaltrials.gov Identifier: NCT02010073.

In phylogenetic analyses of molecular sequence data, partitioning involves estimating independent models of molecular evolution for different sets of sites in a sequence alignment. Choosing an appropriate partitioning scheme is an important step in most analyses because it can affect the accuracy of phylogenetic reconstruction. Despite this, partitioning schemes are often chosen without explicit statistical justification. Here, we describe two new objective methods for the combined selection of best-fit partitioning schemes and nucleotide substitution models. These methods allow millions of partitioning schemes to be compared in realistic time frames and so permit the objective selection of partitioning schemes even for large multilocus DNA data sets. We demonstrate that these methods significantly outperform previous approaches, including both the ad hoc selection of partitioning schemes (e.g., partitioning by gene or codon position) and a recently proposed hierarchical clustering method. We have implemented these methods in an open-source program, PartitionFinder. This program allows users to select partitioning schemes and substitution models using a range of information-theoretic metrics (e.g., the Bayesian information criterion, akaike information criterion [AIC], and corrected AIC). We hope that PartitionFinder will encourage the objective selection of partitioning schemes and thus lead to improvements in phylogenetic analyses. PartitionFinder is written in Python and runs under Mac OSX 10.4 and above. The program, source code, and a detailed manual are freely available from www.robertlanfear.com/partitionfinder.