intspan 0.7.0

Command line tools for IntSpan related bioinformatics operations
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90

use clap::*;
use intspan::*;
use serde_yaml::Value;
use std::collections::BTreeMap;

// Create clap subcommand arguments
pub fn make_subcommand<'a>() -> Command<'a> {
    Command::new("combine")
        .about("Combine multiple sets of runlists in a yaml file")
        .after_help(
            r###"
It's expected that the YAML file contains multiple sets of runlists,
otherwise this command will make no effects

"###,
        )
        .arg(
            Arg::new("infile")
                .help("Sets the input file to use")
                .required(true)
                .index(1),
        )
        .arg(
            Arg::new("op")
                .long("op")
                .takes_value(true)
                .default_value("union")
                .forbid_empty_values(true)
                .help("Operations: intersect, union, diff or xor"),
        )
        .arg(
            Arg::new("outfile")
                .short('o')
                .long("outfile")
                .takes_value(true)
                .default_value("stdout")
                .forbid_empty_values(true)
                .help("Output filename. [stdout] for screen"),
        )
}

// command implementation
pub fn execute(args: &ArgMatches) -> std::result::Result<(), Box<dyn std::error::Error>> {
    //----------------------------
    // Loading
    //----------------------------
    let yaml: BTreeMap<String, Value> = read_yaml(args.value_of("infile").unwrap());
    let s_of = yaml2set_m(&yaml);
    let chrs = chrs_in_sets(&s_of);

    let op = args.value_of("op").unwrap();

    //----------------------------
    // Operating
    //----------------------------
    let mut res: BTreeMap<String, IntSpan> = BTreeMap::new();
    fill_up_s(&mut res, &chrs);

    let names: Vec<_> = s_of.keys().cloned().collect();
    let first = names[0].clone();

    for name in names {
        let set = s_of.get(name.as_str()).unwrap();
        for chr in set.keys() {
            if name == first {
                let intspan = set.get(chr).unwrap();
                res.entry(chr.to_string()).and_modify(|e| e.merge(intspan));
            } else {
                let mut intspan_op = res.get(chr).unwrap().copy();
                intspan_op = match op {
                    "intersect" => intspan_op.intersect(set.get(chr).unwrap()),
                    "diff" => intspan_op.diff(set.get(chr).unwrap()),
                    "union" => intspan_op.union(set.get(chr).unwrap()),
                    "xor" => intspan_op.xor(set.get(chr).unwrap()),
                    _ => panic!("Invalid IntSpan Op"),
                };
                //                eprintln!("Op {}: {}", op, intspan_op.to_string());
                res.insert(chr.into(), intspan_op);
            }
        }
    }

    //----------------------------
    // Output
    //----------------------------
    let out_yaml = set2yaml(&res);
    write_yaml(args.value_of("outfile").unwrap(), &out_yaml)?;

    Ok(())
}