Skip to content

Math related to branch length estimation for Brownian motion continuous character evolution on unrooted phylogenies.

Notifications You must be signed in to change notification settings

argriffing/fels-reml-jan-2014

Repository files navigation

sample output:

incidence matrix:
[[ 1.  0.  0.  0. -1.  0.]
 [ 0.  1.  0.  0. -1.  0.]
 [ 0.  0.  1.  0.  0. -1.]
 [ 0.  0.  0.  1.  0. -1.]
 [ 0.  0.  0.  0.  1. -1.]]

reference branch lengths:
[ 1.10047901  0.18570314  0.48115285  1.05849769  2.37489145]

centered covariance matrix for the reference branch lengths:
[[ 1.32057691  0.44879187 -0.81251629 -0.9568525 ]
 [ 0.44879187  0.86318897 -0.58382232 -0.72815853]
 [-0.81251629 -0.58382232  1.01091383  0.38542477]
 [-0.9568525  -0.72815853  0.38542477  1.29958625]]

test branch lengths:
[ 0.8664449   0.27201223  0.28129483  1.08485675  5.94596636]

cross entropy:
4.56607962074

sampling a bunch of data...

sample data covariance matrix:
[[ 1.31769524  0.44591047 -0.80975888 -0.95384683]
 [ 0.44591047  0.86340588 -0.58334058 -0.72597577]
 [-0.80975888 -0.58334058  1.01089956  0.3821999 ]
 [-0.95384683 -0.72597577  0.3821999   1.2976227 ]]

average log likelihoods:
-4.56869723934

gradient of cross entropy:
[-0.05854639 -0.01943114 -0.03443427 -0.04078485 0.0428304 ]

hessian of cross entropy:
[[  4.97582035e-01   4.63804250e-01   4.88878348e-04  -1.08723293e-05  5.62837524e-04]
 [  4.63804250e-01   5.07361276e-01  -3.16744791e-04  -4.65100077e-04 -2.65770425e-03]
 [  4.88878348e-04  -3.16744791e-04   4.02400603e-01   3.32841627e-01  1.62992650e-03]
 [ -1.08723293e-05  -4.65100077e-04   3.32841627e-01   3.33038731e-01 -6.56618008e-04]
 [  5.62837524e-04  -2.65770425e-03   1.62992650e-03  -6.56618008e-04 -1.13641266e-03]]

eigenvalues of hessian of cross entropy:
[-0.00134004  0.03313779  0.03877968  0.70236462  0.96630417]

minimizing cross entropy...

branch lengths at minimum:
[ 1.10047688  0.18569188  0.48108829  1.05858838  2.3742567 ]

minimum cross entropy:
4.43149459225

gradient at minimum:
[ -4.52529508e-06  -3.58763508e-05  -2.49821793e-05   6.70903445e-06 -4.05207841e-05]

hessian at minimum:
[[  3.07947920e-01   2.69648831e-01   6.00973212e-04   1.24067924e-04 1.27116018e-03]
 [  2.69648831e-01   5.33753118e-01   2.11096999e-02   4.35799255e-03 4.46505931e-02]
 [  6.00973212e-04   2.11096999e-02   3.31782426e-01   1.64997661e-01 2.88342576e-02]
 [  1.24067924e-04   4.35799255e-03   1.64997661e-01   2.33629845e-01 5.95268903e-03]
 [  1.27116018e-03   4.46505931e-02   2.88342576e-02   5.95268903e-03 6.09893418e-02]]

eigenvalues of hessian at minimum:
[ 0.04869748  0.10964552  0.13732916  0.45567661  0.71675389]

About

Math related to branch length estimation for Brownian motion continuous character evolution on unrooted phylogenies.

Resources

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published

Languages