convergence_detection.ml 1.28 KB
Newer Older
Carine Rey's avatar
Carine Rey committed
1 2 3 4 5 6 7 8 9 10 11 12 13
open Bistro.EDSL
open Bistro.Std
open Core
open File_formats
open Bistro_bioinfo.Std

type pcoc_out
type diffsel_out

type det_out =
  | Pcoc_out
  | Diffsel_out

14
let pcoc ?plot_complete ?gamma  ~(faa:aminoacid_fasta workflow) ~(tree:_ workflow) :  (*`pcoc*) det_out directory workflow =
Carine Rey's avatar
Carine Rey committed
15 16 17 18 19 20 21 22
  let env = docker_image ~account:"carinerey" ~name:"pcoc" ~tag:"06212018" () in
  workflow ~descr:"convergence_detection.pcoc" [
    cmd "pcoc_det.py" ~env [
        opt "-t" dep tree;
        opt "-m" string "-";
        opt "-aa"  dep faa ;
        opt "-o"  ident dest ;
        option ( flag string "--gamma" ) gamma;
23
        option ( flag string "--plot --plot_complete_ali" ) plot_complete;
Carine Rey's avatar
Carine Rey committed
24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40
    ]
  ]


let diffsel ~(phy_n:nucleotide_phylip workflow) ~(tree: _ workflow) ~(w_every:int) ~(n_cycles: int) :  (*`diffsel*) det_out directory workflow =
  let env = docker_image ~account:"vlanore" ~name:"diffsel" ~tag:"v1.0" () in
  (*_build/diffsel -t data/samhd1.tree -d data/samhd1.ali -ncond 3 -x 1 10000 myrun*)
  workflow ~descr:"convergence_detection.diffsel" [
    mkdir_p dest;
    cmd "_build/diffsel" ~env [
        opt "-t" dep tree;
        opt "-d"  dep phy_n ;
        opt "-ncond"  int 2 ;
        opt "-x" seq [ int w_every; string " "; int n_cycles];
        seq [ ident dest ; string "/myrun"] ;
    ]
  ]