TODO

[zhenglei-gao]

• error message popup / Error message X
• Gerald error issue ??
• 90% CI for formation fractions? ??
• NA cases examination, p-value to >0.5? X
• Fitting error with MCMC from Johannes X
• formation fraction transformation.

[zhenglei-gao]

• param_org and param_transformed

[zhenglei-gao]

[zhenglei-gao]

compare_multi_kinmod(Scan01,newparms=rbind(Fit1$par,Fit$par))
$SSR
                         ALL     ParentW   ParentS       MetW
Initial Model    38783.35022 11705.51241 207.67849 8754.60534
With New Parms 1    64.26150    13.67971  12.32477   21.68911
With New Parms 2    64.25599    13.67184  12.34968   21.67516
                        MetS
Initial Model    18115.55398
With New Parms 1    16.56791
With New Parms 2    16.55930

$SSR.unweighted
                         ALL     ParentW   ParentS       MetW
Initial Model    11705.51241 11705.51241 207.67849 8754.60534
With New Parms 1    13.67971    13.67971  12.32477   21.68911
With New Parms 2    13.67184    13.67184  12.34968   21.67516
                        MetS
Initial Model    18115.55398
With New Parms 1    16.56791
With New Parms 2    16.55930

[zhenglei-gao]

> cbind(s1$par[,5],summary(Fit)$par[,5])
                              [,1]          [,2]
M0_ParentW            0.000000e+00 1.170668e-132
k_ParentW_to_sink     4.846992e-01  1.408129e-03
k_ParentS_to_sink     3.466559e-01  2.617248e-01
k_MetW_to_sink        5.000000e-01  4.999967e-01
k_MetS_to_sink        4.592472e-01  2.340146e-01
k_ParentW_to_ParentS  0.000000e+00  1.367091e-26
k_ParentW_to_MetW     3.216226e-34  5.121430e-24
k_ParentS_to_ParentW  1.214202e-15  1.917490e-03
k_ParentS_to_MetS     1.056144e-41  1.352622e-02
k_MetW_to_MetS        2.833840e-09  7.145216e-07
k_MetS_to_MetW       1.307252e-315  1.067206e-06

[zhenglei-gao]

👍

[zhenglei-gao]

http://orion.uwaterloo.ca/~hwolkowi/henry/reports/talks.d/t09talks.d/09waterloomatlab.d/optimTipsWebinar/html/optimTipsTricksWalkthrough.html

Matlab optim

[zhenglei-gao]

Wishlist:

1. Multiple Starting Point.
2. Different Parameter Transformations.
3. Comparison of the optimization routines/Benchmark Test?