File: de.swift

package info (click to toggle)
libgoby-java 3.3.1%2Bdfsg2-11
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid, trixie
  • size: 58,108 kB
  • sloc: java: 78,105; cpp: 5,011; xml: 3,170; python: 2,108; sh: 1,575; ansic: 277; makefile: 114
file content (75 lines) | stat: -rw-r--r-- 2,203 bytes parent folder | download | duplicates (2)
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

type alignment {
  string basename;
}
type tsv;
type textfile;

app (textfile t) version_goby() {
   goby "1g" "version" stdout=@filename(t);
}

(string result) tr(string text, string from, string to) {
    app {
        tr from to stdin=text stdout=result;  
    }
}


app (tsv stats) alignment_to_annotation_counts(
string groupId1,
string group1_basenames,
string groupId2,
string group2_basenames,
string spaceSeparatedBasenames,
string annotationFile,
string useWeights,
string adjustGCBias) {


   goby "3g" "alignment-to-annotation-counts" spaceSeparatedBasenames
   "--groups " @strcat(groupId1,"=",group1_basenames,
   "/",
   groupId2,"=",group2_basenames )
   "--compare " @strcat(groupId1,"/",groupId2)
   "--annotation" annotationFile
   "--use-weights" useWeights
   "--adjust-gc-bias" adjustGCBias
   "--parallel --normalization-methods aligned-count"
   "--include-annotation-types gene --write-annotation-counts false"
   "--eval group-averages"
   "--stats" @stats ;

}



(tsv t) call_de(     string groupId1,
                     string group1_basenames,
                     string groupId2,
                     string group2_basenames,
                     string annotationFile,
                     string useWeights,
                     string adjustGCBias) {

  // string spaceSeparatedBasename1 = @regexp(group1_basenames, "[,]"," ");
  string spaceSeparatedBasename1 = tr(group1_basenames, "[,]","[ ]");
  trace(spaceSeparatedBasename1);
  string spaceSeparatedBasename2 = tr(group2_basenames, "[,]","[ ]");
  trace(spaceSeparatedBasename2);
  string allBasenames = @strcat(spaceSeparatedBasename1," ",spaceSeparatedBasename2);

  string statsFilename=@strcat(groupId1 ,"-",groupId2 ,"-",useWeights ,"-",adjustGCBias ,".tsv");
  tsv stats <single_file_mapper;file=statsFilename>;

  trace(tr("AAABBBCCC","[B]","[Z]"));
 
 /* t=  alignment_to_annotation_counts(groupId1="A", group1_basenames,
                                   groupId2="B", group2_basenames=group2_basenames,
                                   spaceSeparatedBasenames=allBasenames,
                                   annotationFile,
                                   useWeights, adjustGCBias);
   */
}