Bayesian inference involves two main computational challenges. First, in estimating the parameters of some model for the data, the posterior distribution may well be highly multi-modal: a regime in which the convergence to stationarity of traditional Markov Chain Monte Carlo (MCMC) techniques becomes incredibly slow. Second, in selecting between a set of competing models the necessary estimation of the Bayesian evidence for each is, by definition, a (possibly high-dimensional) integration over the entire parameter space; again this can be a daunting computational task, although new Monte Carlo (MC) integration algorithms offer solutions of ever increasing efficiency. Nested sampling (NS) is one such contemporary MC strategy targeted at calculation of the Bayesian evidence, but which also enables posterior inference as a by-product, thereby allowing simultaneous parameter estimation and model selection. The widely-used MultiNest algorithm presents a particularly efficient implementation of the NS technique for multi-modal posteriors. In this paper we discuss importance nested sampling (INS), an alternative summation of the MultiNest draws, which can calculate the Bayesian evidence at up to an order of magnitude higher accuracy than `vanilla' NS with no change in the way MultiNest explores the parameter space. This is accomplished by treating as a (pseudo-)importance sample the totality of points collected by MultiNest, including those previously discarded under the constrained likelihood sampling of the NS algorithm. We apply this technique to several challenging test problems and compare the accuracy of Bayesian evidences obtained with INS against those from vanilla NS.

The (supermassive black hole mass, $M_\text{BH}$)-(bulge stellar mass, $M_{\rm*,sph}$) relation is, obviously, derived using two quantities. We endeavor to provide accurate values for the latter via detailed multicomponent galaxy decompositions for the current full sample of 43 spiral galaxies having directly measured $M_\text{BH}$ values; 35 of these galaxies have been alleged to contain pseudobulges, 21 have water maser measurements, and three appear bulgeless. This more than doubles the previous sample size of spiral galaxies with a finessed image analysis. We have analyzed near-infrared images, accounting for not only the bulge, disk (exponential, truncated, or inclined), and bar but also for spiral arms and rings and additional central components (active galactic nuclei (AGNs), etc.). A symmetric Bayesian analysis finds $\log\left(M_\text{BH}/M_{\odot}\right)=\left(2.44_{-0.31}^{+0.35}\right)\log\left\{M_{\rm*,sph}/[\upsilon(1.15\times10^{10}\,M_{\odot})]\right\}+(7.24\pm0.12)$, with $\upsilon$ a stellar mass-to-light ratio term. The level of scatter equals that about the $M_{\rm BH}$-$\sigma_*$ relation. The nonlinear slope rules out the idea that many mergers, coupled with the central limit theorem, produced this scaling relation, and it corroborates previous observational studies and simulations, which have reported a near-quadratic slope at the low-mass end of the $M_\text{BH}$-$M_{\rm*,sph}$ diagram. Furthermore, bulges with AGNs follow this relation; they are not offset by an order of magnitude, and models that have invoked AGN feedback to establish a linear $M_{\rm BH}$-$M_{\rm*,sph}$ relation need revisiting. We additionally present an updated $M_\text{BH}$-(S\'ersic index, $n_\text{sph}$) relation for spiral galaxy bulges with a comparable level of scatter and a new $M_{\rm*,sph}$-(spiral-arm pitch angle, $\phi$) relation.

Black hole mass ($M_{BH}$) scaling relations are typically derived using the properties of a galaxy's bulge and samples dominated by (high-mass) early-type galaxies. Studying late-type galaxies should provide greater insight into the mutual growth of black holes and galaxies in more gas-rich environments. We have used 40 spiral galaxies to establish how $M_{BH}$ scales with both the total stellar mass ($M_{*,tot}$) and the disk's stellar mass, having measured the spheroid (bulge) stellar mass ($M_{*,sph}$) and presented the $M_{BH}$-$M_{*,sph}$ relation in Paper I. The relation involving $M_{*,tot}$ may be beneficial for estimating $M_{BH}$ either from pipeline data or at higher redshift, conditions that are not ideal for the accurate isolation of the bulge. A symmetric Bayesian analysis finds $\log\left(M_{BH}/M_{\odot}\right)=\left(3.05_{-0.49}^{+0.57}\right)\log\left\{M_{*,tot}/[\upsilon(6.37\times10^{10}\,M_{\odot})]\right\}+(7.25_{-0.14}^{+0.13})$. The scatter from the regression of $M_{BH}$ on $M_{*,tot}$ is 0.66 dex; compare 0.56 dex for $M_{BH}$ on $M_{*,sph}$ and $0.57$ dex for $M_{BH}$ on $\sigma_*$. The slope is $>2$ times that obtained using core-S\'ersic early-type galaxies, echoing a similar result involving $M_{*,sph}$, and supporting a varied growth mechanism among different morphological types. This steeper relation has consequences for galaxy/black hole formation theories, simulations, and predicting black hole masses. We caution that (i) an $M_{BH}$-$M_{*,tot}$ relation built from a mixture of early- and late-type galaxies will find an arbitrary slope of approximately 1-3, with no physical meaning beyond one's sample selection, and (ii) evolutionary studies of the $M_{BH}$-$M_{*,tot}$ relation need to be mindful of the galaxy types included at each epoch. We additionally update the $M_{*,tot}$-($\textit{face-on}$ spiral arm pitch angle) relation.

Machine learning techniques offer a precious tool box for use within astronomy to solve problems involving so-called big data. They provide a means to make accurate predictions about a particular system without prior knowledge of the underlying physical processes of the data. In this article, and the companion papers of this series, we present the set of Generalized Linear Models (GLMs) as a fast alternative method for tackling general astronomical problems, including the ones related to the machine learning paradigm. To demonstrate the applicability of GLMs to inherently positive and continuous physical observables, we explore their use in estimating the photometric redshifts of galaxies from their multi-wavelength photometry. Using the gamma family with a log link function we predict redshifts from the PHoto-z Accuracy Testing simulated catalogue and a subset of the Sloan Digital Sky Survey from Data Release 10. We obtain fits that result in catastrophic outlier rates as low as ~1% for simulated and ~2% for real data. Moreover, we can easily obtain such levels of precision within a matter of seconds on a normal desktop computer and with training sets that contain merely thousands of galaxies. Our software is made publicly available as an user-friendly package developed in Python, R and via an interactive web application (https://cosmostatisticsinitiative.shinyapps.io/CosmoPhotoz). This software allows users to apply a set of GLMs to their own photometric catalogues and generates publication quality plots with minimum effort from the user. By facilitating their ease of use to the astronomical community, this paper series aims to make GLMs widely known and to encourage their implementation in future large-scale projects, such as the Large Synoptic Survey Telescope.

Astronomical observations of extended sources, such as cubes of integral field spectroscopy (IFS), encode auto-correlated spatial structures that cannot be optimally exploited by standard methodologies. This work introduces a novel technique to model IFS datasets, which treats the observed galaxy properties as realizations of an unobserved Gaussian Markov random field. The method is computationally efficient, resilient to the presence of low-signal-to-noise regions, and uses an alternative to Markov Chain Monte Carlo for fast Bayesian inference, the Integrated Nested Laplace Approximation (INLA). As a case study, we analyse 721 IFS data cubes of nearby galaxies from the CALIFA and PISCO surveys, for which we retrieve the maps of the following physical properties: age, metallicity, mass and extinction. The proposed Bayesian approach, built on a generative representation of the galaxy properties, enables the creation of synthetic images, recovery of areas with bad pixels, and an increased power to detect structures in datasets subject to substantial noise and/or sparsity of sampling. A snippet code to reproduce the analysis of this paper is available in the COIN toolbox, together with the field reconstructions of the CALIFA and PISCO samples.

Approximate Bayesian Computation (ABC) enables parameter inference for complex physical systems in cases where the true likelihood function is unknown, unavailable, or computationally too expensive. It relies on the forward simulation of mock data and comparison between observed and synthetic catalogues. Here we present cosmoabc, a Python ABC sampler featuring a Population Monte Carlo (PMC) variation of the original ABC algorithm, which uses an adaptive importance sampling scheme. The code is very flexible and can be easily coupled to an external simulator, while allowing to incorporate arbitrary distance and prior functions. As an example of practical application, we coupled cosmoabc with the numcosmo library and demonstrate how it can be used to estimate posterior probability distributions over cosmological parameters based on measurements of galaxy clusters number counts without computing the likelihood function. cosmoabc is published under the GPLv3 license on PyPI and GitHub and documentation is available at http://goo.gl/SmB8EX

We developed a hierarchical Bayesian model (HBM) to investigate how the presence of Seyfert activity relates to their environment, herein represented by the galaxy cluster mass, $M_{200}$, and the normalized cluster-centric distance, $r/r_{200}$. We achieved this by constructing an unbiased sample of galaxies from the Sloan Digital Sky Survey, with morphological classifications provided by the Galaxy Zoo Project. A propensity score matching approach is introduced to control for the effects of confounding variables: stellar mass, galaxy colour, and star formation rate. The connection between Seyfert-activity and environmental properties in the de-biased sample is modelled within an HBM framework using the so-called logistic regression technique, suitable for the analysis of binary data (e.g., whether or not a galaxy hosts an AGN). Unlike standard ordinary least square fitting methods, our methodology naturally allows modelling the probability of Seyfert-AGN activity in galaxies on their natural scale, i.e. as a binary variable. Furthermore, we demonstrate how an HBM can incorporate information of each particular galaxy morphological type in a unified framework. In elliptical galaxies, our analysis indicates a strong correlation of Seyfert-AGN activity with $r/r_{200}$, and a weaker correlation with the mass of the host. In spiral galaxies these trends do not appear, suggesting that the link between Seyfert activity and the properties of spiral galaxies are independent of the environment.

Whether or not the initial star cluster mass function is established through a universal, galactocentric-distance-independent stochastic process, on the scales of individual galaxies, remains an unsolved problem. This debate has recently gained new impetus through the publication of a study that concluded that the maximum cluster mass in a given population is not solely determined by size-of-sample effects. Here, we revisit the evidence in favor and against stochastic cluster formation by examining the young ($\lesssim$ a few $\times 10^8$ yr-old) star cluster mass--galactocentric radius relation in M33, M51, M83, and the Large Magellanic Cloud. To eliminate size-of-sample effects, we first adopt radial bin sizes containing constant numbers of clusters, which we use to quantify the radial distribution of the first- to fifth-ranked most massive clusters using ordinary least-squares fitting. We supplement this analysis with an application of quantile regression, a binless approach to rank-based regression taking an absolute-value-distance penalty. Both methods yield, within the $1\sigma$ to $3\sigma$ uncertainties, near-zero slopes in the diagnostic plane, largely irrespective of the maximum age or minimum mass imposed on our sample selection, or of the radial bin size adopted. We conclude that, at least in our four well-studied sample galaxies, star cluster formation does not necessarily require an environment-dependent cluster formation scenario, which thus supports the notion of stochastic star cluster formation as the dominant star cluster-formation process within a given galaxy.

In this paper, the third in a series illustrating the power of generalized linear models (GLMs) for the astronomical community, we elucidate the potential of the class of GLMs which handles count data. The size of a galaxy's globular cluster population $N_{\rm GC}$ is a prolonged puzzle in the astronomical literature. It falls in the category of count data analysis, yet it is usually modelled as if it were a continuous response variable. We have developed a Bayesian negative binomial regression model to study the connection between $N_{\rm GC}$ and the following galaxy properties: central black hole mass, dynamical bulge mass, bulge velocity dispersion, and absolute visual magnitude. The methodology introduced herein naturally accounts for heteroscedasticity, intrinsic scatter, errors in measurements in both axes (either discrete or continuous), and allows modelling the population of globular clusters on their natural scale as a non-negative integer variable. Prediction intervals of 99% around the trend for expected $N_{\rm GC}$comfortably envelope the data, notably including the Milky Way, which has hitherto been considered a problematic outlier. Finally, we demonstrate how random intercept models can incorporate information of each particular galaxy morphological type. Bayesian variable selection methodology allows for automatically identifying galaxy types with different productions of GCs, suggesting that on average S0 galaxies have a GC population 35% smaller than other types with similar brightness.

Machine learning techniques offer a plethora of opportunities in tackling big data within the astronomical community. We present the set of Generalized Linear Models as a fast alternative for determining photometric redshifts of galaxies, a set of tools not commonly applied within astronomy, despite being widely used in other professions. With this technique, we achieve catastrophic outlier rates of the order of ~1%, that can be achieved in a matter of seconds on large datasets of size ~1,000,000. To make these techniques easily accessible to the astronomical community, we developed a set of libraries and tools that are publicly available.

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.

We investigate the utility to computational Bayesian analyses of a particular family of recursive marginal likelihood estimators characterized by the (equivalent) algorithms known as "biased sampling" or "reverse logistic regression" in the statistics literature and "the density of states" in physics. Through a pair of numerical examples (including mixture modeling of the well-known galaxy data set) we highlight the remarkable diversity of sampling schemes amenable to such recursive normalization, as well as the notable efficiency of the resulting pseudo-mixture distributions for gauging prior sensitivity in the Bayesian model selection context. Our key theoretical contributions are to introduce a novel heuristic ("thermodynamic integration via importance sampling") for qualifying the role of the bridging sequence in this procedure and to reveal various connections between these recursive estimators and the nested sampling technique.

In astronomical and cosmological studies one often wishes to infer some properties of an infinite-dimensional field indexed within a finite-dimensional metric space given only a finite collection of noisy observational data. Bayesian inference offers an increasingly-popular strategy to overcome the inherent ill-posedness of this signal reconstruction challenge. However, there remains a great deal of confusion within the astronomical community regarding the appropriate mathematical devices for framing such analyses and the diversity of available computational procedures for recovering posterior functionals. In this brief research note I will attempt to clarify both these issues from an "applied statistics" perpective, with insights garnered from my post-astronomy experiences as a computational Bayesian / epidemiological geostatistician.

We have analysed the growth of Brightest Group Galaxies and Brightest Cluster Galaxies (BGGs/BCGs) over the last 3 billion years using a large sample of 883 galaxies from the Galaxy And Mass Assembly Survey. By comparing the stellar mass of BGGs and BCGs in groups and clusters of similar dynamical masses, we find no significant growth between redshift $z=0.27$ and $z=0.09$. We also examine the number of BGGs/BCGs that have line emission, finding that approximately 65 per cent of BGGs/BCGs show H$\alpha$ in emission. From the galaxies where the necessary spectroscopic lines were accurately recovered (54 per cent of the sample), we find that half of this (i.e. 27 per cent of the sample) harbour on-going star formation with rates up to $10\,$M$_{\odot}$yr$^{-1}$, and the other half (i.e. 27 per cent of the sample) have an active nucleus (AGN) at the centre. BGGs are more likely to have ongoing star formation, while BCGs show a higher fraction of AGN activity. By examining the position of the BGGs/BCGs with respect to their host dark matter halo we find that around 13 per cent of them do not lie at the centre of the dark matter halo. This could be an indicator of recent cluster-cluster mergers. We conclude that BGGs and BCGs acquired their stellar mass rapidly at higher redshifts as predicted by semi-analytic models, mildly slowing down at low redshifts.

In this brief research note I present a generalized version of the Savage-Dickey Density Ratio for representation of the Bayes factor (or marginal likelihood ratio) of nested statistical models; the new version takes the form of a Radon-Nikodym derivative and is thus applicable to a wider family of probability spaces than the original (restricted to those admitting an ordinary Lebesgue density). A derivation is given following the measure-theoretic construction of Marin & Robert (2010), and the equivalent estimator is demonstrated in application to a distributional modeling problem.

We present photometric measurements for the galaxies - and when possible their bulges and disks - in the 0.05<z<0.0585 groups of the Zurich Environmental Study (ZENS); these measurements include (B-I) colors, color gradients and maps, color dispersions, as well as stellar masses and star-formation rates. The ZENS galaxies are classified into quenched, moderately star-forming, and strongly star-forming using a combination of spectral features and FUV-to-optical colors; this approach optimally distinguishes quenched systems from dust-reddened star-forming galaxies. The latter contribute up to 50% to the (B-I) "red sequence" at ~10^10Msun. At fixed morphological or spectral type, we find that galaxy stellar masses are largely independent of environment, and especially of halo mass. As a first utilization of our photometric database, we study, at fixed stellar mass and Hubble type, how (B-I) colors, color gradients and color dispersion of disk satellites depend on group mass (M_GROUP), group-centric distance (R/R_200) and large-scale structure overdensity. The strongest environmental trend is found for disk-dominated satellites with M_GROUP and R/R_200. At M<10^10 Msun, disk-dominated satellites are redder in the inner regions of the groups than in the outer parts. At M>10^10 Msun, these satellites have shallower color gradients in higher mass groups and in the cores of groups compared with lower mass groups and the outskirts of groups. Stellar population analyses and semi-analytic models suggest that disk-dominated satellites undergo quenching of star formation in their outer disks, on timescales ~2 Gyr, as they progressively move inside the group potential.

In the second paper of this series we extend our Bayesian reanalysis of the evidence for a cosmic variation of the fine structure constant to the semi-parametric modelling regime. By adopting a mixture of Dirichlet processes prior for the unexplained errors in each instrumental subgroup of the benchmark quasar dataset we go some way towards freeing our model selection procedure from the apparent subjectivity of a fixed distributional form. Despite the infinite-dimensional domain of the error hierarchy so constructed we are able to demonstrate a recursive scheme for marginal likelihood estimation with prior-sensitivity analysis directly analogous to that presented in Paper I, thereby allowing the robustness of our posterior Bayes factors to hyper-parameter choice and model specification to be readily verified. In the course of this work we elucidate various similarities between unexplained error problems in the seemingly disparate fields of astronomy and clinical meta-analysis, and we highlight a number of sophisticated techniques for handling such problems made available by past research in the latter. It is our hope that the novel approach to semi-parametric model selection demonstrated herein may serve as a useful reference for others exploring this potentially difficult class of error model.

We review the evidence behind recent claims of spatial variation in the fine structure constant deriving from observations of ionic absorption lines in the light from distant quasars. To this end we expand upon previous non-Bayesian analyses limited by the assumptions of an unbiased and strictly Normal distribution for the "unexplained errors" of the benchmark quasar dataset. Through the technique of reverse logistic regression we estimate and compare marginal likelihoods for three competing hypotheses---(i) the null hypothesis (no cosmic variation), (ii) the monopole hypothesis (a constant Earth-to-quasar offset), and (iii) the monopole+dipole hypothesis (a cosmic variation manifest to the Earth-bound observer as a North-South divergence)---under a variety of candidate parametric forms for the unexplained error term. Our analysis reveals weak support for a skeptical interpretation in which the apparent dipole effect is driven solely by systematic errors of opposing sign inherent in measurements from the two telescopes employed to obtain these observations. Throughout we seek to exemplify a 'best practice' approach to Bayesian model selection with prior-sensitivity analysis; in a companion paper we extend this methodology to a semi-parametric framework using the infinite-dimensional Dirichlet process.

We present structural measurements for the galaxies in the 0.05<z<0.0585 groups of the Zurich Environmental Study, aimed at establishing how galaxy properties depend on four environmental parameters: group halo mass M_GROUP, group-centric distance R/R_200, ranking into central or satellite, and large-scale structure density delta_LSS. Global galaxy structure is quantified both parametrically and non-parametrically. We correct all these measurements for observational biases due to PSF blurring and surface brightness effects as a function of galaxy size, magnitude, steepness of light profile and ellipticity. Structural parameters are derived also for bulges, disks and bars. We use the galaxy bulge-to-total ratios (B/T), together with the calibrated non-parametric structural estimators, to implement a quantitative morphological classification that maximizes purity in the resulting morphological samples. We investigate how the concentration C of satellite galaxies depends on galaxy mass for each Hubble type, and on M_GROUP, R/R_200 and delta_LSS. At galaxy masses M>10^10 M_sun, the concentration of disk satellites increases with increasing stellar mass, separately within each morphological bin of B/T. The known increase in concentration with stellar mass for disk satellites is thus due, at least in part, to an increase in galaxy central stellar density at constant B/T. The correlation between concentration and galaxy stellar mass becomes progressively steeper for later morphological types. The concentration of disk satellites shows a barely significant dependence on delta_LSS or R/R_200. The strongest environmental effect is found with group mass for M>10^10 M_sun disk-dominated satellites, which are ~10% more concentrated in high mass groups than in lower mass groups.

The Zurich Environmental Study (ZENS) is based on a sample of ~1500 galaxy members of 141 groups in the mass range ~10^12.5-14.5 M_sun within the narrow redshift range 0.05<z<0.0585. ZENS adopts novel approaches, here described, to quantify four different galactic environments, namely: (1) the mass of the host group halo; (2) the projected halo-centric distance; (3) the rank of galaxies as central or satellites within their group halos; and (4) the filamentary large-scale structure (LSS) density. No self-consistent identification of a central galaxy is found in ~40% of <10^13.5 M_sun groups, from which we estimate that ~15% of groups at these masses are dynamically unrelaxed systems. Central galaxies in relaxed and unrelaxed groups have in general similar properties, suggesting that centrals are regulated by their mass and not by their environment. Centrals in relaxed groups have however ~30% larger sizes than in unrelaxed groups, possibly due accretion of small satellites in virialized group halos. At M>10^10 M_sun, satellite galaxies in relaxed and unrelaxed groups have similar size, color and (specific) star formation rate distributions; at lower galaxy masses, satellites are marginally redder in relaxed relative to unrelaxed groups, suggesting quenching of star formation in low-mass satellites by physical processes active in relaxed halos. Finally, relaxed and unrelated groups show similar stellar mass conversion efficiencies, peaking at halo masses around 10^12.5 M_sun. In the enclosed ZENS catalogue we publish all environmental diagnostics as well as the galaxy structural and photometric measurements described in companion ZENS papers II and III.

Abridged. We use COSMOS to study in a self-consistent way the change in the number densities of quenched early-type galaxies (Q-ETGs) of a given size over the interval 0.2 < z < 1.0 to study the claimed size evolution of these galaxies. At 10^10.5<Mgalaxy<10^11 Msun, we see no change in the number density of compact Q-ETGs, while at >10^11 Msun we find a decrease by 30%. In both mass bins, the increase of the median sizes of Q-ETGs with time is primarily caused by the addition to the size function of larger and more diffuse Q-ETGs. At all masses, compact Q-ETGs become systematically redder towards later epochs, with a (U-V) difference consistent with passive evolution of their stellar populations, indicating that they are a population that does not appreciably evolve in size. At all epochs, the larger Q-ETGs (at least in the lower mass bin) have average rest-frame colors systematically bluer than those of the more compact Q-ETGs, suggesting that the former are younger than the latter. The idea that new, large, Q-ETGs are responsible for the observed growth in the median size of the population at a given mass is supported by the sizes and number of the star-forming galaxies that are expected to be progenitors of the new Q-ETGs over the same period. In the low mass bin, the new Q-ETG have 30% smaller sizes than their star-forming progenitors. This is likely due to the fading of their disks after they cease star-formation. Comparison with higher z shows that the median size of newly-quenched galaxies roughly scales, at constant mass, as (1+z)^-1. The dominant cause of the size evolution seen in the Q-ETG population is thus that the average sizes of individual Q-ETGs scale with the average density of the Universe at the time when they were quenched, with subsequent size changes in individual objects through eg merging of secondary importance, especially at masses <10^11 Msun.

The Galaxy And Mass Assembly (GAMA) survey is a multiwavelength photometric and spectroscopic survey, using the AAOmega spectrograph on the Anglo-Australian Telescope to obtain spectra for up to ~300000 galaxies over 280 square degrees, to a limiting magnitude of r_pet < 19.8 mag. The target galaxies are distributed over 0<z<0.5 with a median redshift of z~0.2, although the redshift distribution includes a small number of systems, primarily quasars, at higher redshifts, up to and beyond z=1. The redshift accuracy ranges from sigma_v~50km/s to sigma_v~100km/s depending on the signal-to-noise of the spectrum. Here we describe the GAMA spectroscopic reduction and analysis pipeline. We present the steps involved in taking the raw two-dimensional spectroscopic images through to flux-calibrated one-dimensional spectra. The resulting GAMA spectra cover an observed wavelength range of 3750<lambda<8850 A at a resolution of R~1300. The final flux calibration is typically accurate to 10-20%, although the reliability is worse at the extreme wavelength ends, and poorer in the blue than the red. We present details of the measurement of emission and absorption features in the GAMA spectra. These measurements are characterised through a variety of quality control analyses detailing the robustness and reliability of the measurements. We illustrate the quality of the measurements with a brief exploration of elementary emission line properties of the galaxies in the GAMA sample. We demonstrate the luminosity dependence of the Balmer decrement, consistent with previously published results, and explore further how Balmer decrement varies with galaxy mass and redshift. We also investigate the mass and redshift dependencies of the [NII]/Halpha vs [OIII]/Hbeta spectral diagnostic diagram, commonly used to discriminate between star forming and nuclear activity in galaxies.

We have measured the clustering properties of low-redshift (z < 0.3) sub-mm galaxies detected at 250 micron in the Herschel-ATLAS Science Demonstration Phase (SDP) field. We selected a sample for which we have high-quality spectroscopic redshifts, obtained from reliably matching the 250-micron sources to a complete (for r < 19.4) sample of galaxies from the GAMA database. Both the angular and spatial clustering strength are measured for all z < 0.3 sources as well as for five redshift slices with thickness delta z=0.05 in the range 0.05 < z < 0.3. Our measured spatial clustering length r_0 is comparable to that of optically-selected, moderately star-forming (blue) galaxies: we find values around 5 Mpc. One of the redshift bins contains an interesting structure, at z = 0.164.

Context: The mass-metallicity relationship (MMR) of star-forming galaxies is well-established, however there is still some disagreement with respect to its exact shape and its possible dependence on other observables. Aims: We measure the MMR in the Galaxy And Mass Assembly (GAMA) survey. We compare our measured MMR to that measured in the Sloan Digital Sky Survey (SDSS) and study the dependence of the MMR on various selection criteria to identify potential causes for disparities seen in the literature. Methods: We use strong emission line ratio diagnostics to derive oxygen abundances. We then apply a range of selection criteria for the minimum signal-to-noise in various emission lines, as well as the apparent and absolute magnitude to study variations in the inferred MMR. Results: The shape and position of the MMR can differ significantly depending on the metallicity calibration and selection used. After selecting a robust metallicity calibration amongst those tested, we find that the mass-metallicity relation for redshifts 0.061< z<0.35 in GAMA is in reasonable agreement with that found in the SDSS despite the difference in the luminosity range probed. Conclusions: In view of the significant variations of the MMR brought about by reasonable changes in the sample selection criteria and method, we recommend that care be taken when comparing the MMR from different surveys and studies directly. We also conclude that there could be a modest level of evolution over 0.06<z<0.35 within the GAMA sample.

We use the GAMA I dataset combined with GALEX, SDSS and UKIDSS imaging to construct the low-redshift (z<0.1) galaxy luminosity functions in FUV, NUV, ugriz, and YJHK bands from within a single well constrained volume of 3.4 x 10^5 (Mpc/h)^{3}. The derived luminosity distributions are normalised to the SDSS DR7 main survey to reduce the estimated cosmic variance to the 5 per cent level. The data are used to construct the cosmic spectral energy distribution (CSED) from 0.1 to 2.1 \mum free from any wavelength dependent cosmic variance for both the elliptical and non-elliptical populations. The two populations exhibit dramatically different CSEDs as expected for a predominantly old and young population respectively. Using the Driver et al. (2008) prescription for the azimuthally averaged photon escape fraction, the non-ellipticals are corrected for the impact of dust attenuation and the combined CSED constructed. The final results show that the Universe is currently generating (1.8 +/- 0.3) x 10^{35} h W Mpc^{-3} of which (1.2 +/- 0.1) x 10^{35} h W Mpc^{-3} is directly released into the inter-galactic medium and (0.6 +/- 0.1) x 10^{35} h W Mpc^{-3} is reprocessed and reradiated by dust in the far-IR. Using the GAMA data and our dust model we predict the mid and far-IR emission which agrees remarkably well with available data. We therefore provide a robust description of the pre- and post dust attenuated energy output of the nearby Universe from 0.1micron to 0.6mm. The largest uncertainty in this measurement lies in the mid and far-IR bands stemming from the dust attenuation correction and its currently poorly constrained dependence on environment, stellar mass, and morphology.

We measure the two-point angular correlation function of a sample of 4,289,223 galaxies with r < 19.4 mag from the Sloan Digital Sky Survey as a function of photometric redshift, absolute magnitude and colour down to M_r - 5log h = -14 mag. Photometric redshifts are estimated from ugriz model magnitudes and two Petrosian radii using the artificial neural network package ANNz, taking advantage of the Galaxy and Mass Assembly (GAMA) spectroscopic sample as our training set. The photometric redshifts are then used to determine absolute magnitudes and colours. For all our samples, we estimate the underlying redshift and absolute magnitude distributions using Monte-Carlo resampling. These redshift distributions are used in Limber's equation to obtain spatial correlation function parameters from power law fits to the angular correlation function. We confirm an increase in clustering strength for sub-L* red galaxies compared with ~L* red galaxies at small scales in all redshift bins, whereas for the blue population the correlation length is almost independent of luminosity for ~L* galaxies and fainter. A linear relation between relative bias and log luminosity is found to hold down to luminosities L~0.03L*. We find that the redshift dependence of the bias of the L* population can be described by the passive evolution model of Tegmark & Peebles (1998). A visual inspection of a random sample of our r < 19.4 sample of SDSS galaxies reveals that about 10 per cent are spurious, with a higher contamination rate towards very faint absolute magnitudes due to over-deblended nearby galaxies. We correct for this contamination in our clustering analysis.

"Approximate Bayesian Computation" (ABC) represents a powerful methodology for the analysis of complex stochastic systems for which the likelihood of the observed data under an arbitrary set of input parameters may be entirely intractable-the latter condition rendering useless the standard machinery of tractable likelihood-based, Bayesian statistical inference (e.g. conventional Markov Chain Monte Carlo simulation; MCMC). In this article we demonstrate the potential of ABC for astronomical model analysis by application to a case study in the morphological transformation of high redshift galaxies. To this end we develop, first, a stochastic model for the competing processes of merging and secular evolution in the early Universe; and second, through an ABC-based comparison against the observed demographics of massive (M_gal > 10^11 M_sun) galaxies (at 1.5 < z < 3) in the CANDELS/EGS dataset we derive posterior probability densities for the key parameters of this model. The "Sequential Monte Carlo" (SMC) implementation of ABC exhibited herein, featuring both a self-generating target sequence and self-refining MCMC kernel, is amongst the most efficient of contemporary approaches to this important statistical algorithm. We highlight as well through our chosen case study the value of careful summary statistic selection, and demonstrate two modern strategies for assessment and optimisation in this regard. Ultimately, our ABC analysis of the high redshift morphological mix returns tight constraints on the evolving merger rate in the early Universe and favours major merging (with disc survival or rapid reformation) over secular evolution as the mechanism most responsible for building up the first generation of bulges in early-type disks.

We present a detailed investigation into the effects of galaxy environment on their star formation rates (SFR) using galaxies observed in the Galaxy and Mass Assembly Survey (GAMA). We use three independent volume-limited samples of galaxies within z < 0.2 and Mr < -17.8. We investigate the known SFR-density relationship and explore in detail the dependence of SFR on stellar mass and density. We show that the SFR-density trend is only visible when we include the passive galaxy population along with the star-forming population. This SFR-density relation is absent when we consider only the star-forming population of galaxies, consistent with previous work. While there is a strong dependence of the EWH?a on density we find, as in previous studies, that these trends are largely due to the passive galaxy population and this relationship is absent when considering a "star-forming" sample of galaxies. We find that stellar mass has the strongest influence on SFR and EWH?a with the environment having no significant effect on the star-formation properties of the star forming population. We also show that the SFR-density relationship is absent for both early and late-type star-forming galaxies. We conclude that the stellar mass has the largest impact on the current SFR of a galaxy, and any environmental effect is not detectable. The observation that the trends with density are due to the changing morphology fraction with density implies that the timescales must be very short for any quenching of the SFR in infalling galaxies. Alternatively galaxies may in fact undergo predominantly in-situ evolution where the infall and quenching of galaxies from the field into dense environments is not the dominant evolutionary mode.

Galaxy and Mass Assembly (GAMA) is a project to study galaxy formation and evolution, combining imaging data from ultraviolet to radio with spectroscopic data from the AAOmega spectrograph on the Anglo-Australian Telescope. Using data from phase 1 of GAMA, taken over three observing seasons, and correcting for various minor sources of incompleteness, we calculate galaxy luminosity functions (LFs) and their evolution in the ugriz passbands. At low redshift, z < 0.1, we find that blue galaxies, defined according to a magnitude-dependent but non-evolving colour cut, are reasonably well fit over a range of more than ten magnitudes by simple Schechter functions in all bands. Red galaxies, and the combined blue-plus-red sample, require double power-law Schechter functions to fit a dip in their LF faintward of the characteristic magnitude M* before a steepening faint end. This upturn is at least partly due to dust-reddened disk galaxies. We measure evolution of the galaxy LF over the redshift range 0.002 < z < 0.5 using both a parametric fit and by measuring binned LFs in redshift slices. The characteristic luminosity L* is found to increase with redshift in all bands, with red galaxies showing stronger luminosity evolution than blue galaxies. The comoving number density of blue galaxies increases with redshift, while that of red galaxies decreases, consistent with prevailing movement from blue cloud to red sequence. As well as being more numerous at higher redshift, blue galaxies also dominate the overall luminosity density beyond redshifts z = 0.2. At lower redshifts, the luminosity density is dominated by red galaxies in the riz bands, by blue galaxies in u and g.

We present near-IR surface photometry (2D-profiling) for a sample of 29 nearby galaxies for which super-massive black hole (SMBH) masses are constrained. The data is derived from the UKIDSS-LASS survey representing a significant improvement in image quality and depth over previous studies based on 2MASS data. We derive the spheroid luminosity and spheroid S\'ersic index for each galaxy with GALFIT3 and use these data to construct SMBH mass -bulge luminosity ($M_{\rm bh}$--$L$) and SMBH - S\'ersic index ($M_{\rm bh}$--$n$) relations. The best fit K-band relation for elliptical and disk galaxies is $\log(M_{\rm bh}/M_{\odot})= -0.36(\pm 0.03) (M_{\rm K} + 18) + 6.17(\pm 0.16)$ with an intrinsic scatter of 0.4$^{+0.09}_{-0.06}$dex whilst for elliptical galaxies we find $\log(M_{\rm bh}/M_{\odot})= -0.42(\pm 0.06) (M_{\rm K} + 22) + 7.5(\pm 0.15)$ with an intrinsic scatter of 0.31$^{+0.087}_{-0.047}$dex. Our revised $M_{\rm bh}$--$L$ relation agrees closely with the previous near-IR constraint by \citet{tex:G07}. The lack of improvement in the intrinsic scatter in moving to higher quality near-IR data suggests that the SMBH relations are not currently limited by the quality of the imaging data but is either intrinsic or a result of uncertainty in the precise number of required components required in the profiling process. Contrary to expectation (see \citealt{tex:GD07a}) a relation between SMBH mass and the S\'ersic index was not found at near-IR wavelengths. This latter outcome is believed to be explained by the generic inconsistencies between 1D and 2D galaxy profiling which are currently under further investigation.

This paper describes the first catalogue of photometrically-derived stellar mass estimates for intermediate-redshift (z < 0.65) galaxies in the Galaxy And Mass Assembly (GAMA) spectroscopic redshift survey. These masses, as well as the full set of ancillary stellar population parameters, will be made public as part of GAMA data release 2. Although the GAMA database does include NIR photometry, we show that the quality of our stellar population synthesis fits is significantly poorer when these NIR data are included. Further, for a large fraction of galaxies, the stellar population parameters inferred from the optical-plus-NIR photometry are formally inconsistent with those inferred from the optical data alone. This may indicate problems in our stellar population library, or NIR data issues, or both; these issues will be addressed for future versions of the catalogue. For now, we have chosen to base our stellar mass estimates on optical photometry only. In light of our decision to ignore the available NIR data, we examine how well stellar mass can be constrained based on optical data alone. We use generic properties of stellar population synthesis models to demonstrate that restframe colour alone is in principle a very good estimator of stellar mass-to-light ratio, M*/Li. Further, we use the observed relation between restframe (g-i) and M*/Li for real GAMA galaxies to argue that, modulo uncertainties in the stellar evolution models themselves, (g-i) colour can in practice be used to estimate M*/Li to an accuracy of < ~0.1 dex. This 'empirically calibrated' (g-i)-M*/Li relation offers a simple and transparent means for estimating galaxies' stellar masses based on minimal data, and so provides a solid basis for other surveys to compare their results to z < ~0.4 measurements from GAMA.

We use the high angular resolution in the near-infrared of the WFC3 on HST to determine YHVz color-color selection criteria to identify and characterize 1.5<z<3.5 galaxies in the HUDF09 and ERS (GOODS-South) fields. The WFC3 NIR images reveal galaxies at these redshifts that were undetected in the rest-frame UV HUDF/GOODS images, as well as true centers and regular disks in galaxies classified as highly irregular in rest-frame UV light. Across the 1.5<z<2.15 redshift range, regular disks are unveiled in the WFC3 images of ~25% of both intermediate and high mass galaxies, i.e., above 10^10 Msun. Meanwhile, galaxies maintaining diffuse and/or irregular morphologies in the rest-frame optical light---i.e., not yet dynamically settled---at these epochs are almost entirely restricted to masses below 10^11 Msun. In contrast at 2.25 < z < 3.5 these diffuse and/or irregular structures overwhelmingly dominate the morphological mix in both the intermediate and high mass regimes, while no regular disks, and only a small fraction (25%) of smooth spheroids, are evident above 10^11 Msun. Strikingly, by 1.5 < z < 2.25 roughly 2 out of every 3 galaxies at the highest masses are spheroids. In our small sample, the fraction of star-forming galaxies at these mass scales decreases concurrently from ~60% to ~5%. If confirmed, this indicates that z~2 is the epoch of both the morphological transformation and quenching of star-formation which assemble the first substantial population of massive ellipticals.

We investigate the properties of satellite galaxies that surround isolated hosts within the redshift range 0.01 < z < 0.15, using data taken as part of the Galaxy And Mass Assembly survey. Making use of isolation and satellite criteria that take into account stellar mass estimates, we find 3514 isolated galaxies of which 1426 host a total of 2998 satellites. Separating the red and blue populations of satellites and hosts, using colour-mass diagrams, we investigate the radial distribution of satellite galaxies and determine how the red fraction of satellites varies as a function of satellite mass, host mass and the projected distance from their host. Comparing the red fraction of satellites to a control sample of small neighbours at greater projected radii, we show that the increase in red fraction is primarily a function of host mass. The satellite red fraction is about 0.2 higher than the control sample for hosts with 11.0 < log M_* < 11.5, while the red fractions show no difference for hosts with 10.0 < log M_* < 10.5. For the satellites of more massive hosts the red fraction also increases as a function of decreasing projected distance. Our results suggest that the likely main mechanism for the quenching of star formation in satellites hosted by isolated galaxies is strangulation.

Using the complete GAMA-I survey covering ~142 sq. deg. to r=19.4, of which ~47 sq. deg. is to r=19.8, we create the GAMA-I galaxy group catalogue (G3Cv1), generated using a friends-of-friends (FoF) based grouping algorithm. Our algorithm has been tested extensively on one family of mock GAMA lightcones, constructed from Lambda-CDM N-body simulations populated with semi-analytic galaxies. Recovered group properties are robust to the effects of interlopers and are median unbiased in the most important respects. G3Cv1 contains 14,388 galaxy groups (with multiplicity >= 2$), including 44,186 galaxies out of a possible 110,192 galaxies, implying ~40% of all galaxies are assigned to a group. The similarities of the mock group catalogues and G3Cv1 are multiple: global characteristics are in general well recovered. However, we do find a noticeable deficit in the number of high multiplicity groups in GAMA compared to the mocks. Additionally, despite exceptionally good local spatial completeness, G3Cv1 contains significantly fewer compact groups with 5 or more members, this effect becoming most evident for high multiplicity systems. These two differences are most likely due to limitations in the physics included of the current GAMA lightcone mock. Further studies using a variety of galaxy formation models are required to confirm their exact origin.

The stellar initial mass function (IMF) describes the distribution in stellar masses produced from a burst of star formation. For more than fifty years, the implicit assumption underpinning most areas of research involving the IMF has been that it is universal, regardless of time and environment. We measure the high-mass IMF slope for a sample of low-to-moderate redshift galaxies from the Galaxy And Mass Assembly survey. The large range in luminosities and galaxy masses of the sample permits the exploration of underlying IMF dependencies. A strong IMF-star formation rate dependency is discovered, which shows that highly star forming galaxies form proportionally more massive stars (they have IMFs with flatter power-law slopes) than galaxies with low star formation rates. This has a significant impact on a wide variety of galaxy evolution studies, all of which rely on assumptions about the slope of the IMF. Our result is supported by, and provides an explanation for, the results of numerous recent explorations suggesting a variation of or evolution in the IMF.

We use multiwavelength data from the Galaxy And Mass Assembly (GAMA) and Herschel ATLAS (H-ATLAS) surveys to compare the relationship between various dust obscuration measures in galaxies. We explore the connections between the ultraviolet (UV) spectral slope, $\beta$, the Balmer decrement, and the far infrared (IR) to $150\,$nm far ultraviolet (FUV) luminosity ratio. We explore trends with galaxy mass, star formation rate (SFR) and redshift in order to identify possible systematics in these various measures. We reiterate the finding of other authors that there is a large scatter between the Balmer decrement and the $\beta$ parameter, and that $\beta$ may be poorly constrained when derived from only two broad passbands in the UV. We also emphasise that FUV derived SFRs, corrected for dust obscuration using $\beta$, will be overestimated unless a modified relation between $\beta$ and the attenuation factor is used. Even in the optimum case, the resulting SFRs have a significant scatter, well over an order of magnitude. While there is a stronger correlation between the IR to FUV luminosity ratio and $\beta$ parameter than with the Balmer decrement, neither of these correlations are particularly tight, and dust corrections based on $\beta$ for high redshift galaxy SFRs must be treated with caution. We conclude with a description of the extent to which the different obscuration measures are consistent with each other as well as the effects of including other galactic properties on these correlations.

I present a critical review of techniques for estimating confidence intervals on binomial population proportions inferred from success counts in small-to-intermediate samples. Population proportions arise frequently as quantities of interest in astronomical research; for instance, in studies aiming to constrain the bar fraction, AGN fraction, SMBH fraction, merger fraction, or red sequence fraction from counts of galaxies exhibiting distinct morphological features or stellar populations. However, two of the most widely-used techniques for estimating binomial confidence intervals--the 'normal approximation' and the Clopper & Pearson approach--are liable to misrepresent the degree of statistical uncertainty present under sampling conditions routinely encountered in astronomical surveys, leading to an ineffective use of the experimental data (and, worse, an inefficient use of the resources expended in obtaining that data). Hence, I provide here an overview of the fundamentals of binomial statistics with two principal aims: (i) to reveal the ease with which (Bayesian) binomial confidence intervals with more satisfactory behaviour may be estimated from the quantiles of the beta distribution using modern mathematical software packages (e.g. R, matlab, mathematica, IDL, python); and (ii) to demonstrate convincingly the major flaws of both the 'normal approximation' and the Clopper & Pearson approach for error estimation.

We present an analysis of the properties of the lowest Halpha-luminosity galaxies (L_Halpha<4x10^32 W; SFR<0.02 Msun/yr) in the Galaxy And Mass Assembly (GAMA) survey. These galaxies make up the the rise above a Schechter function in the number density of systems seen at the faint end of the Halpha luminosity function. Above our flux limit we find that these galaxies are principally composed of intrinsically low stellar mass systems (median stellar mass =2.5x10^8 Msun) with only 5/90 having stellar masses M>10^10 Msun. The low SFR systems are found to exist predominantly in the lowest density environments (median density ~0.02 galaxy Mpc^-2 with none in environments more dense than ~1.5 galaxy Mpc^-2). Their current specific star formation rates (SSFR; -8.5 < log(SSFR[yr^-1])<-12.) are consistent with their having had a variety of star formation histories. The low density environments of these galaxies demonstrates that such low-mass, star-forming systems can only remain as low-mass and forming stars if they reside sufficiently far from other galaxies to avoid being accreted, dispersed through tidal effects or having their gas reservoirs rendered ineffective through external processes.

The Galaxy And Mass Assembly (GAMA) survey has been operating since February 2008 on the 3.9-m Anglo-Australian Telescope using the AAOmega fibre-fed spectrograph facility to acquire spectra with a resolution of R~1300 for 120,862 SDSS selected galaxies. The target catalogue constitutes three contiguous equatorial regions centred at 9h (G09), 12h (G12) and 14.5h (G15) each of 12 x 4 sq.deg to limiting fluxes of r < 19.4, r < 19.8, and r < 19.4 mag respectively (and additional limits at other wavelengths). Spectra and reliable redshifts have been acquired for over 98 per cent of the galaxies within these limits. Here we present the survey footprint, progression, data reduction, redshifting, re-redshifting, an assessment of data quality after 3 years, additional image analysis products (including ugrizYJHK photometry, Sersic profiles and photometric redshifts), observing mask, and construction of our core survey catalogue (GamaCore). From this we create three science ready catalogues: GamaCoreDR1 for public release, which includes data acquired during year 1 of operations within specified magnitude limits (February 2008 to April 2008); GamaCoreMainSurvey containing all data above our survey limits for use by the GAMA team and collaborators; and GamaCoreAtlasSv containing year 1, 2 and 3 data matched to Herschel-ATLAS Science Demonstration data. These catalogues along with the associated spectra, stamps and profiles can be accessed via the GAMA website: http://www.gama-survey.org/

In order to generate credible 0.1-2 {\mu}m SEDs, the GAMA project requires many Gigabytes of imaging data from a number of instruments to be re-processed into a standard format. In this paper we discuss the software infrastructure we use, and create self-consistent ugrizYJHK photometry for all sources within the GAMA sample. Using UKIDSS and SDSS archive data, we outline the pre-processing necessary to standardise all images to a common zeropoint, the steps taken to correct for seeing bias across the dataset, and the creation of Gigapixel-scale mosaics of the three 4x12 deg GAMA regions in each filter. From these mosaics, we extract source catalogues for the GAMA regions using elliptical Kron and Petrosian matched apertures. We also calculate S\'ersic magnitudes for all galaxies within the GAMA sample using SIGMA, a galaxy component modelling wrapper for GALFIT 3. We compare the resultant photometry directly, and also calculate the r band galaxy LF for all photometric datasets to highlight the uncertainty introduced by the photometric method. We find that (1) Changing the object detection threshold has a minor effect on the best-fitting Schechter parameters of the overall population (M* +/- 0.055mag, {\alpha} +/- 0.014, {\Phi}* +/- 0.0005 h^3 Mpc^{-3}). (2) An offset between datasets that use Kron or Petrosian photometry regardless of the filter. (3) The decision to use circular or elliptical apertures causes an offset in M* of 0.20mag. (4) The best-fitting Schechter parameters from total-magnitude photometric systems (such as SDSS modelmag or S\'ersic magnitudes) have a steeper faint-end slope than photometry dependent on Kron or Petrosian magnitudes. (5) Our Universe's total luminosity density, when calculated using Kron or Petrosian r-band photometry, is underestimated by at least 15%.

We present self-consistent star formation rates derived through pan-spectral analysis of galaxies drawn from the Galaxy and Mass Assembly (GAMA) survey. We determine the most appropriate form of dust obscuration correction via application of a range of extinction laws drawn from the literature as applied to Halpha, [O{II}] and UV luminosities. These corrections are applied to a sample of 31,508 galaxies from the GAMA survey at z < 0.35. We consider several different obscuration curves, including those of Milky Way, Calzetti (2001) and Fischera and Dopita (2005) curves and their effects on the observed luminosities. At the core of this technique is the observed Balmer decrement, and we provide a prescription to apply optimal obscuration corrections using the Balmer decrement. We carry out an analysis of the star formation history (SFH) using stellar population synthesis tools to investigate the evolutionary history of our sample of galaxies as well as to understand the effects of variation in the Initial Mass Function (IMF) and the effects this has on the evolutionary history of galaxies. We find that the Fischera and Dopita (2005) obscuration curve with an R_{v} value of 4.5 gives the best agreement between the different SFR indicators. The 2200A feature needed to be removed from this curve to obtain complete consistency between all SFR indicators suggesting that this feature may not be common in the average integrated attenuation of galaxy emission. We also find that the UV dust obscuration is strongly dependent on the SFR.

We investigate the (large-scale) bar fraction in a mass-complete sample of M > 10^10.5 Msun disk galaxies at 0.2 < z < 0.6 in the COSMOS field. The fraction of barred disks strongly depends on mass, disk morphology, and specific star formation rate (SSFR). At intermediate stellar mass (10^10.5 < M < 10^11 Msun) the bar fraction in early-type disks is much higher, at all redshifts, by a factor ~2, than that in late-type disks. This trend is reversed at higher stellar mass (M > 10^11 Msun), where the fraction of bars in early-type disks becomes significantly lower, at all redshifts, than that in late-type disks. The bar fractions for galaxies with low and high SSFRs closely follow those of the morphologically-selected early-type and late-type populations, respectively. This indicates a close correspondence between morphology and SSFR in disk galaxies at these earlier epochs. Interestingly, the total bar fraction in 10^10.5 < M < 10^11 Msun disks is built up by a factor of ~2 over the redshift interval explored, while for M > 10^11 Msun disks it remains roughly constant. This indicates that, already by z ~ 0.6, spectral and morphological transformations in the most massive disk galaxies have largely converged to the familiar Hubble sequence that we observe in the local Universe, while for intermediate mass disks this convergence is ongoing until at least z ~ 0.2. Moreover, these results highlight the importance of employing mass-limited samples for quantifying the evolution of barred galaxies. Finally, the evolution of the barred galaxy populations investigated does not depend on the large-scale environmental density (at least, on the scales which can be probed with the available photometric redshifts).

We have determined the luminosity function of 250um-selected galaxies detected in the ~14 sq.deg science demonstration region of the Herschel-ATLAS project out to a redshift of z=0.5. Our findings very clearly show that the luminosity function evolves steadily out to this redshift. By selecting a sub-group of sources within a fixed luminosity interval where incompleteness effects are minimal, we have measured a smooth increase in the comoving 250um luminosity density out to z=0.2 where it is 3.6+1.4-0.9 times higher than the local value.

We present colour-colour diagrams of detected sources in the Herschel-ATLAS Science Demonstration Field from 100 to 500 microns using both PACS and SPIRE. We fit isothermal modified black bodies to the spectral energy distribution (SED) to extract the dust temperature of sources with counterparts in Galaxy And Mass Assembly (GAMA) or SDSS surveys with either a spectroscopic or a photometric redshift. For a subsample of 330 sources detected in at least three FIR bands with a significance greater than 3 $\sigma$, we find an average dust temperature of $(28 \pm 8)$K. For sources with no known redshift, we populate the colour-colour diagram with a large number of SEDs generated with a broad range of dust temperatures and emissivity parameters, and compare to colours of observed sources to establish the redshift distribution of this sample. For another subsample of 1686 sources with fluxes above 35 mJy at 350 microns and detected at 250 and 500 microns with a significance greater than 3$\sigma$, we find an average redshift of $2.2 \pm 0.6$.

We combine data from the MGC, SDSS and UKIDSS LAS surveys to produce ugrizYJHK luminosity functions and densities from within a common, low redshift volume (z<0.1, ~71,000 h_1^-3 Mpc^3 for L* systems) with 100 per cent spectroscopic completeness. In the optical the fitted Schechter functions are comparable in shape to those previously reported values but with higher normalisations (typically 0, 30, 20, 15, 5 per cent higher phi*-values in u, g, r, i, z respectively over those reported by the SDSS team). We attribute these to differences in the redshift ranges probed, incompleteness, and adopted normalisation methods. In the NIR we find significantly different Schechter function parameters (mainly in the M* values) to those previously reported and attribute this to the improvement in the quality of the imaging data over previous studies. This is the first homogeneous measurement of the extragalactic luminosity density which fully samples both the optical and near-IR regimes. Unlike previous compilations that have noted a discontinuity between the optical and near-IR regimes our homogeneous dataset shows a smooth cosmic spectral energy distribution (CSED). After correcting for dust attenuation we compare our CSED to the expected values based on recent constraints on the cosmic star-formation history and the initial mass function.

We describe the spectroscopic target selection for the Galaxy And Mass Assembly (GAMA) survey. The input catalogue is drawn from the Sloan Digital Sky Survey (SDSS) and UKIRT Infrared Deep Sky Survey (UKIDSS). The aim is to measure redshifts for galaxies in three 4x12 degree regions at 9h, 12h and 14.5h, on the celestial equator, with magnitude selections r < 19.4, z < 18.2 and K(AB) < 17.6 over all three regions, and r < 19.8 in the 12-h region. The target density is 1080 per sq. deg. in the 12-h region and 720 per sq. deg. in the other regions. The average GAMA target density and area are compared with completed and ongoing galaxy redshift surveys. The GAMA survey implements a highly complete star-galaxy separation that jointly uses an intensity-profile separator (delta(sg) = r-band psf mag - model mag) as per the SDSS) and a colour separator. The colour separator is defined as delta(sg,jk) = J - K - f(g-i), where f(g-i) is a quadratic fit to the J-K colour of the stellar locus over the range 0.3 < g-i < 2.3. All galaxy populations investigated are well separated with delta(sg,jk) > 0.2. From two years out of a three-year AAOmega program on the Anglo-Australian Telescope, we have obtained 79599 unique galaxy redshifts. Previously known redshifts in the GAMA region bring the total up to 98497. The median galaxy redshift is 0.2 with 99% at z < 0.5. We present some of the global statistical properties of the survey, including K-band galaxy counts, colour-redshift relations and preliminary n(z).

A heuristic greedy algorithm is developed for efficiently tiling spatially dense redshift surveys. In its first application to the Galaxy and Mass Assembly (GAMA) redshift survey we find it rapidly improves the spatial uniformity of our data, and naturally corrects for any spatial bias introduced by the 2dF multi object spectrograph. We make conservative predictions for the final state of the GAMA redshift survey after our final allocation of time, and can be confident that even if worse than typical weather affects our observations, all of our main survey requirements will be met.

The Galaxy And Mass Assembly (GAMA) project is the latest in a tradition of large galaxy redshift surveys, and is now underway on the 3.9m Anglo-Australian Telescope at Siding Spring Observatory. GAMA is designed to map extragalactic structures on scales of 1kpc - 1Mpc in complete detail to a redshift of z~0.2, and to trace the distribution of luminous galaxies out to z~0.5. The principal science aim is to test the standard hierarchical structure formation paradigm of Cold Dark Matter (CDM) on scales of galaxy groups, pairs, discs, bulges and bars. We will measure (1) the Dark Matter Halo Mass Function (as inferred from galaxy group velocity dispersions); (2) baryonic processes, such as star formation and galaxy formation efficiency (as derived from Galaxy Stellar Mass Functions); and (3) the evolution of galaxy merger rates (via galaxy close pairs and galaxy asymmetries). Additionally, GAMA will form the central part of a new galaxy database, which aims to contain 275,000 galaxies with multi-wavelength coverage from coordinated observations with the latest international ground- and space-based facilities: GALEX, VST, VISTA, WISE, HERSCHEL, GMRT and ASKAP. Together, these data will provide increased depth (over 2 magnitudes), doubled spatial resolution (0.7"), and significantly extended wavelength coverage (UV through Far-IR to radio) over the main SDSS spectroscopic survey for five regions, each of around 50 deg^2. This database will permit detailed investigations of the structural, chemical, and dynamical properties of all galaxy types, across all environments, and over a 5 billion year timeline.

We investigate the origin of the galaxy color-concentration bimodality at the bright-end of the luminosity function (M(B) - 5 log h < -18 mag) with regard to the bulge-disc nature of galaxies. Via (2D) surface brightness profile modeling with GIM2D, we subdivide the local galaxy population in the Millennium Galaxy Catalogue into one-component and two-component systems. We reveal that one-component (elliptical and disc-only) systems define the two peaks of the galaxy color-concentration distribution (with total stellar mass densities of 0.7 +/- 0.1 and 1.3 +/- 0.1 x 10^8 h Msol Mpc^-3 respectively), while two-component systems contribute to both a bridging population and the red, concentrated peak (with total stellar mass densities of 1.1 +/- 0.1 and 1.8 +/- 0.2 x 10^8 h Msun Mpc^-3 respectively). Moreover, luminous, `bulge-less, red discs' and `disc-less, blue bulges' (blue ellipticals) are exceptionally rare (with volume-densities of 1.7 +/- 0.3 and 1.1 \pm 0.1 x 10^-4 h^3 Mpc^-3 respectively). Finally, within the two-component population we confirm a previously-reported correlation between bulge and disc color (with a mean offset of only <(u-r)_bulge-(u-r)_disc> = 0.22 +/- 0.02 mag).

We explore the use of the color-log(n) plane (where n is the global Sersic index) as a tool for subdividing the high redshift galaxy population in a physically-motivated manner. Using a series of volume-limited samples out to z=1.5 in the Hubble Ultra Deep Field (UDF) we confirm the correlation between color-log(n) plane position and visual morphology observed locally and in other high redshift studies in the color and/or structure domain. Via comparison to a low redshift sample from the Millennium Galaxy Catalogue we quantify evolution by color-log(n) type, accounting separately for the specific selection and measurement biases against each. Specifically, we measure decreases in B-band surface brightness of 1.57 +/- 0.22 mag/sq.arcsec and 1.65 +/- 0.22 mag/sq.arcsec for `blue, diffuse' and `red, compact' galaxies respectively between redshift unity and the present day.

We use the Hubble Ultra Deep Field to study the galaxy luminosity-size (M-Re) distribution. With a careful analysis of selection effects due to both detection completeness and measurement reliability we identify bias-free regions in the M-Re plane for a series of volume-limited samples. By comparison to a nearby survey also having well defined selection limits, namely the Millennium Galaxy Catalogue, we present clear evidence for evolution in surface brightness since z ~ 0.7. Specifically, we demonstrate that the mean, rest-frame B-band surface brightness for galaxies in a sample spanning 8 magnitudes in luminosity between M = -22 and -14 mag increases by ~1.0 mag/arcsec^2 from z ~ 0.1 to z ~ 0.7. We also highlight the importance of considering surface brightness dependent measurement biases in addition to incompleteness biases. In particular, the increasing, systematic under-estimation of Kron fluxes towards low surface brightnesses may cause diffuse, yet luminous, systems to be mistaken for faint, compact objects.

We use the Hubble Ultra Deep Field (UDF) to study the galaxy luminosity-size (M-Re) distribution. With a careful analysis of selection effects due to both detection completeness and measurement reliability we identify bias-free regions in the M-Re plane for a series of volume-limited samples. We also investigate the colour-log(n) distribution of these galaxies and further subdivide our data by structural type to separately study compact and diffuse objects. By comparison to the nearby Millennium Galaxy Catalogue, we present tenative evidence for evolution of diffuse, disk-like galaxies with redshift -- both in mean surface brightness and the slope of the M-Re relation. In contrast we find no evidence of structural evolution in the compact galaxy M-Re relation over this redshift range, although there is a suggestion of colour evolution. We also highlight the importance of considering surface brightness dependent measurement biases in addition to incompleteness biases. In particular, the increasing, systematic under-estimation of Kron fluxes towards low surface brightnesses may cause diffuse, yet luminous, systems to be mistaken for faint, compact objects.

We have modelled the light distribution in 10095 galaxies from the Millennium Galaxy Catalogue (MGC), providing publically available structural catalogues for a large, representative sample of galaxies in the local Universe. Three different models were used: (1) a single Sersic function for the whole galaxy, (2) a bulge-disc decomposition model using a de Vaucouleurs (R^{1/4}) bulge plus exponential disc, (3) a bulge-disc decomposition model using a Sersic (R^{1/n}) bulge plus exponential disc. Repeat observations for 700 galaxies demonstrate that stable measurements can be obtained for object components with a half-light radius comparable to, or larger than, the seeing half-width at half maximum. We show that with careful quality control, robust measurements can be obtained for large samples such as the MGC. We use the catalogues to show that the galaxy colour bimodality is due to the two-component nature of galaxies (i.e. bulges and discs) and not to two distinct galaxy populations. We conclude that understanding galaxy evolution demands the routine bulge-disc decomposition of the giant galaxy population at all redshifts.

Using 10 095 galaxies (B < 20 mag) from the Millennium Galaxy Catalogue, we derive B-band luminosity distributions and selected bivariate brightness distributions for the galaxy population. All subdivisions extract highly correlated sub-sets of the galaxy population which consistently point towards two overlapping distributions. A clear bimodality in the observed distribution is seen in both the rest-(u-r) colour and log(n) distributions. The rest-(u-r) colour bimodality becomes more pronounced when using the core colour as opposed to global colour. The two populations are extremely well separated in the colour-log(n) plane. Using our sample of 3 314 (B < 19 mag) eyeball classified galaxies, we show that the bulge-dominated, early-type galaxies populate one peak and the bulge-less, late-type galaxies occupy the second. The early- and mid-type spirals sprawl across and between the peaks. This constitutes extremely strong evidence that the fundamental way to divide the luminous galaxy population is into bulges and discs and that the galaxy bimodality reflects the two component nature of galaxies and not two distinct galaxy classes. We argue that these two-components require two independent formation mechanisms/processes and advocate early bulge formation through initial collapse and ongoing disc formation through splashback, infall and merging/accretion. We calculate the B-band luminosity-densities and stellar-mass densities within each subdivision and estimate that the z ~ 0 stellar mass content in spheroids, bulges and discs is 35 +/- 2 per cent, 18 +/- 7 and 47 +/- 7 per cent respectively. [Abridged]