ntitation data of a DNA sample from the signal, G is a matrix that represents the genotypes in the DNA sample, with a column j giving the alleles for individual j, w is a weight column vector that represents relative proportions of template DNA in the sample and e is an error vector, wherein the solution includes calculation of data variance sigma**(2) from the linear model d = G.w+e together with