next up previous contents
Next:  Examples Up:  Methods, algorithms and a Previous:  Methods, algorithms and a   Contents


 The method

The question is : Given a set of incomplete and noisy data (say, Fo,h with its sigma(Fh) and PHI_h), which map (of a large number of maps consistent with the observed data) is the one that will minimise the probability of misinterpreting it ? Stating the same problem in a different way, we could ask (i) which map (of the set of admissible maps) will only show features for which there is evidence in the data, or, (ii) which map makes the least assumptions about the data (especially the missing data, but also the distribution of errors in the observed).

Clearly, the Fo,hexp(iPHI_h) synthesis is not the map we want : Not only we assume that all missing data have F = 0 (a rather improbable event), but also that Fh = Fo,h,forallh. Gull, S.F. & Daniell, G.J.4, suggested that the map we need is the one for which the configurational entropy - SUM(j)mjlogmj, where mj is the density at the grid point j of the map, reaches a maximum. It is easy to see that - SUM(j)mjlogmj reaches a maximum when mj = e-1,forallj, that is, when the map has a uniform density, and thus, contain no information. Maximising - SUM(j)mjlogmj subject to the constraint that the map is consistent with the observed data, gives the MAXENT map.

The consistency with the observed data is described in terms of the difference between the observed data and those calculated from a trial map, weighted by the standard deviation of the measurement. If Fc,h is the calculated value of the datum h, Fo,h its observed value and sigma(Fh) the standard deviation of the observation, then the statistic

SUM(h)(F_c,h - F_o,h)^2 / sigma(F_h)^2

possesses a Chi^2 distribution with an expected value equal to the number of data points. Maximising - SUM(j)mjlogmj subject to the constraint SUM(h) | Fc,h - Fo,h^2/sigma(Fh)2 = n, where n is the number of data points, gives the basic iteration formula :

mj = exp{ - 1 + lambdaSUM(h)(Fo,h - Fc,h) / sigma(Fh)^2exp(2piijh)}

Given Fo,h, sigma(Fh) and an positive multiplier lambda, this equation can determine the densities mj on a map. The program GraphEnt applies this formula iteratively (starting from a uniform map) until convergence (as judged by the value of Chi^2) is achieved. Although this algorithm is neither the most efficient nor the most stable, it is relatively easy to code and it leads (at least in the case of Patterson functions), to the same results as other, more complex algorithms5.



Footnotes

... G.J.4
Gull, S.F. & Daniell, G.J., (1978), Nature, 272, 686-690.
... algorithms5
Skilling, J. & Bryan, R.K., (1984), Mon. Not. R. astr. Soc., 211, 111-124.

next up previous contents
Next:  Examples Up:  Methods, algorithms and a Previous:  Methods, algorithms and a   Contents
NMG, Nov 2002