Software Description and Links
The link http://www.math.missouri.edu/~calvin/entropyfiles/index.html
contains Matlab program files for three methods of estimating the minimal
entropy distance d(a,b).
The minimal entropy distance d(a,b), assuming a rich probability path
minimzing arc, is computed by the Matlab program e02daug.m via e02daug(a,b)
which calls functions e02faug, e02fdyn, e02h, e02HN, e02Hofh, e02init,
e02Jacsq, e02phi and e02regr.
The function e02init initializes variables ``hatsig'' and ``hatz'' for use
of the differential system determined by e02fdyn which implements
Proposition 12(Explicit f) for f.
The functions e02Hofh and e02HN use e02h (which is the code given in
Proposition 12(Explicit f)) to evaluate the entropy function of a vector.
The program e02phi gives the result of following the solution of the initial
value problem across the s interval. Newton's method uses e02Jacsq, the
numerical Jacobian matrix. Linear regression to improve the arc length
estimate is accomplished via e02regr. Finally, the differential system is
augmented by e02faug which adds step 13 of Section (deflated) in order
that d can be computed without quadrature methods by solving the augmented
system.
The program e02dlinear uses an input argument N, an integer of at least
100, to estimate the value of the entropy integral over the linear path
containing N equally spaced points on the closed linear path joining
a to b. It calls e02Hofh and e02h to generate the functional values
followed by e02quad which calls e02regr to carry out a composite
Simpson's rule. Note that N can be any integer as opposed to for the
usual Simpson's rule which requires and even number of subintervals,
which would have made N odd.
The elementary metric d_G(a,b) is computed by e02d_G(a,b) which
calls the function e02G.
Programs e02dlinear and e02d_G will run on any pair of probablility
vectors a and b of the same length N.
We used a package called Weighbor to create the tree structures found in the examples section.
In addition C programs implementing to generate the dendronic
tree encapsulated post-script were written see the link Will's code.
web page:
http://www.math.missouri.edu/~calvin/entropyfiles/index.html.
|