
We present a method of calibrating the properties of photometric redshift
bins as part of a larger Markov Chain Monte Carlo (MCMC) analysis for the
inference of cosmological parameters. The redshift bins are characterised by
their mean and variance, which are varied as free parameters and marginalised
over when obtaining the cosmological parameters. We demonstrate that the
likelihood function for crosscorrelations in an angular power spectrum
framework tightly constrains the properties of bins such that they may be well
determined, reducing their influence on cosmological parameters and avoiding
the bias from poorly estimated redshift distributions. We demonstrate that even
with only three photometric and three spectroscopic bins, we can recover
accurate estimates of the mean redshift of a bin to within $\Delta\mu \approx
34 \times10^{3}$ and the width of the bin to $\Delta\sigma \approx
1\times10^{3}$ for galaxies near $z = 1$. This indicates that we may be able
to bring down the photometric redshift errors to a level which is in line with
the requirements for the next generation of cosmological experiments.

A new model independent method is presented for the analysis of pulsar timing
data and the estimation of the spectral properties of an isotropic
gravitational wave background (GWB). We show that by rephrasing the likelihood
we are able to eliminate the most costly aspects of computation normally
associated with this type of data analysis. When applied to the International
Pulsar Timing Array Mock Data Challenge data sets this results in speedups of
approximately 2 to 3 orders of magnitude compared to established methods. We
present three applications of the new likelihood. In the low signal to noise
regime we sample directly from the power spectrum coefficients of the GWB
signal realization. In the high signal to noise regime, where the data can
support a large number of coefficients, we sample from the joint probability
density of the power spectrum coefficients for the individual pulsars and the
GWB signal realization. Critically in both these cases we need make no
assumptions about the form of the power spectrum of the GWB, or the individual
pulsars. Finally we present a method for characterizing the spatial correlation
between pulsars on the sky, making no assumptions about the form of that
correlation, and therefore providing the only truly general Bayesian method of
confirming a GWB detection from pulsar timing data.

A new catalogue of extrasolar planets is presented by reanalysing a
selection of published radial velocity data sets using EXOFIT (Balan & Lahav
2009). All objects are treated on an equal footing within a Bayesian framework,
to give orbital parameters for 94 exoplanetary systems. Model selection
(between one and twoplanet solutions) is then performed, using both a visual
flagging method and a standard chisquare analysis, with agreement between the
two methods for 99% of the systems considered. The catalogue is to be made
available online, and this 'proof of concept' study may be maintained and
extended in the future to incorporate all systems with publicly available
radial velocity data, as well as transit and microlensing data.

We present the results of the GREAT08 Challenge, a blind analysis challenge
to infer weak gravitational lensing shear distortions from images. The primary
goal was to stimulate new ideas by presenting the problem to researchers
outside the shear measurement community. Six GREAT08 Team methods were
presented at the launch of the Challenge and five additional groups submitted
results during the 6 month competition. Participants analyzed 30 million
simulated galaxies with a range in signal to noise ratio, pointspread function
ellipticity, galaxy size, and galaxy type. The large quantity of simulations
allowed shear measurement methods to be assessed at a level of accuracy
suitable for currently planned future cosmic shear observations for the first
time. Different methods perform well in different parts of simulation parameter
space and come close to the target level of accuracy in several of these. A
number of fresh ideas have emerged as a result of the Challenge including a
reexamination of the process of combining information from different galaxies,
which reduces the dependence on realistic galaxy modelling. The image
simulations will become increasingly sophisticated in future GREAT challenges,
meanwhile the GREAT08 simulations remain as a benchmark for additional
developments in shear measurement algorithms.

The GRavitational lEnsing Accuracy Testing 2008 (GREAT08) Challenge focuses
on a problem that is of crucial importance for future observations in
cosmology. The shapes of distant galaxies can be used to determine the
properties of dark energy and the nature of gravity, because light from those
galaxies is bent by gravity from the intervening dark matter. The observed
galaxy images appear distorted, although only slightly, and their shapes must
be precisely disentangled from the effects of pixelisation, convolution and
noise. The worldwide gravitational lensing community has made significant
progress in techniques to measure these distortions via the Shear TEsting
Program (STEP). Via STEP, we have run challenges within our own community, and
come to recognise that this particular image analysis problem is ideally
matched to experts in statistical inference, inverse problems and computational
learning. Thus, in order to continue the progress seen in recent years, we are
seeking an infusion of new ideas from these communities. This document details
the GREAT08 Challenge for potential participants. Please visit
http://www.great08challenge.info for the latest information.

Retrieval of orbital parameters of extrasolar planets poses considerable
statistical challenges.Due to sparse sampling, measurement errors, parameters
degeneracy and modelling limitations, there are no unique values of basic
parameters, such as period and eccentricity. Here, we estimate the orbital
parameters from radial velocity data in a Bayesian framework by utilizing
Markov Chain Monte Carlo (MCMC) simulations with the MetropolisHastings
algorithm. We follow a methodology recently proposed by Gregory and Ford. Our
implementation of MCMC is based on the objectoriented approach outlined by
Graves. We make our resulting code, ExoFit, publicly available with this paper.
It can search for either one or two planets as illustrated on mock data. As an
example we reanalysed the orbital solution of companions to HD 187085 and HD
159868 from the published radial velocity data. We confirm the degeneracy
reported for orbital parameters of the companion to HD 187085, and show that a
loweccentricity orbit is more probable for this planet. For HD 159868, we
obtained slightly different orbital solution and a relatively high 'noise'
factor indicating the presence of an unaccounted signal in the radial velocity
data. ExoFit is designed in such a way that it can be extended for a variety of
probability models, including different Bayesian priors.