The invention relates to a method for determining haplotypes. The method is capable of providing phase information from a set of unphased genotype observations in an efficient robust manner. The method is capable of handling large datasets that include missing data.