The variance of three-locus linkage disequilibria for an equilibrium infinite alleles model is solved numerically on a computer, using identity coefficients. It is shown that the variance of three-locus linkage disequilibrium created by random drift, although smaller than the variance of two-locus linkage disequilibrium, is of the same order of magnitude. Hence third-order disequilibria are not necessarily good indications of selection. The formula for the variance of linkage disequilibrium is given when there is no recombination between the genes. This model can also be interpreted as intragenic recombination between three sites within a gene.