:- use_module(library(debug)).
:- use_module(library(lib)).

:- lib(gbn).
:- lib(debug_call).

:- debug(mye).

/** mye.

Run myeloma experiments as reported on paper.

Data copied from myeloma/rea.18/wgs/data/rea_min20.dat.

*/

mye :-
        debug_call( mye, start, true ),
    DatF = 'data/mye/mye_min20.dat',
    E = 2,
    GBNOpts = [ data(DatF),
                setting(edge_penalty,E),
                dir(OsOdir)
                % debug(true)
              ],
    gbn( GBNOpts ),
        debug_call( mye, start, fisher_nets ),
    gbn_fisher_nets( dir(OsOdir) ),    
        debug_call( mye, start, fam_hmaps ),
    gbn_fam_hmaps( [dir(OsOdir),not(true)] ),   % only for binaries
        debug_call( mye, start, gates_nets ),
    gbn_gates_nets( dir(OsOdir) ),
        debug_call( mye, start, svg_legend ),
    gbn_svg_legend( dir(OsOdir) ),
        debug_call( mye, end, true ).   

display_var_as( hyper, "HRD" ).
display_var_as( t4_14, "t(4;14)" ).
display_var_as( t11_14, "t(11;14)" ).
display_var_as( t14_16, "t(14;16)" ).