annotate rDiff/src/process_command_line_args.m @ 2:233c30f91d66

updated python based GFF parsing module which will handle GTF/GFF/GFF3 file types
author vipints <vipin@cbio.mskcc.org>
date Tue, 08 Oct 2013 07:15:44 -0400
parents 0f80a5141704
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
1 function [CFG]=process_command_line_args(CFG,ARGS)
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
2 % [CFG]=process_command_line_args(CFG,ARGS)
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
3
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
4 %This function adds the command line arguments to the CFG config variable
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
5
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
6 %Parse the ARGS
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
7
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
8 %Check wether octave or matlab is used
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
9 if size(ver('Octave'),1)
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
10 INTERPR = 1;
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
11 else
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
12 INTERPR = 0;
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
13 end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
14
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
15 %turn of warning
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
16 if INTERPR
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
17 warning('off', 'Octave:precedence-change');
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
18 warning('off', 'Octave:function-name-clash');
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
19 warning('off', '');
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
20 warning('off', 'Octave:num-to-str');
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
21 warning('off', 'Octave:function-name-clash');
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
22 warning('off', 'Octave:divide-by-zero');
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
23 warning('off', 'Octave:future-time-stamp');
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
24 warning('off', 'solve_qp:constraints');
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
25 warning('off', 'Octave:assign-as-truth-value');
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
26 warning('off', 'Octave:matlab-incompatible');
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
27 else
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
28 warning('off', 'MATLAB:typeaheadBufferOverflow');
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
29 end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
30
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
31 %Seperate the Variable Field
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
32 if INTERPR
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
33 ARGS = strsplit(ARGS,';');
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
34 else
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
35 ARGS=regexp(ARGS,';','split');
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
36 end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
37
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
38 %Assign the local variables
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
39
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
40 %iterate over Fields
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
41 for i=1:length(ARGS)
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
42 if isempty(ARGS{i})
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
43 continue
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
44 end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
45 if INTERPR
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
46 VALS = strsplit(ARGS{i},':');
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
47 else
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
48 VALS=regexp(ARGS{i},':','split');
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
49 end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
50
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
51
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
52 if length(VALS)>2
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
53 error([" more than one field for variable: " VALS{1} ":" VALS{:} "\n Maybe there are colons in the input argument?"])
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
54 end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
55
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
56 if strcmp(VALS{1},"RDIFF_RES_DIR"),RDIFF_RES_DIR=VALS{2};continue,end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
57 if strcmp(VALS{1},"RDIFF_INPUT_DIR"),RDIFF_INPUT_DIR=VALS{2};continue,end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
58 if strcmp(VALS{1},"BAM_INPUT1"),BAM_INPUT1=VALS{2};continue,end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
59 if strcmp(VALS{1},"BAM_INPUT2"),BAM_INPUT2=VALS{2};continue,end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
60 if strcmp(VALS{1},"GFF_INPUT"),GFF_INPUT=VALS{2};continue,end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
61 if strcmp(VALS{1},"READ_LENGTH"),READ_LENGTH=str2num(VALS{2});continue,end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
62 if strcmp(VALS{1},"MIN_READ_LENGTH"),MIN_READ_LENGTH=str2num(VALS{2});continue,end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
63 if strcmp(VALS{1},"EST_GENE_EXPR"),EST_GENE_EXPR=str2num(VALS{2});continue,end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
64 if strcmp(VALS{1},"ONLY_GENE_EXPR"),ONLY_GENE_EXPR=str2num(VALS{2});continue,end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
65 if strcmp(VALS{1},"VAR_PATH1"),VAR_PATH1=VALS{2};continue,end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
66 if strcmp(VALS{1},"VAR_PATH2"),VAR_PATH2=VALS{2};continue,end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
67 if strcmp(VALS{1},"SAVE_VAR1"),SAVE_VAR1=VALS{2};continue,end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
68 if strcmp(VALS{1},"SAVE_VAR2"),SAVE_VAR2=VALS{2};continue,end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
69 if strcmp(VALS{1},"PRED_VAR1"),PRED_VAR1=VALS{2};continue,end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
70 if strcmp(VALS{1},"PRED_VAR2"),PRED_VAR2=VALS{2};continue,end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
71 if strcmp(VALS{1},"ONLY_GENE_START"),ONLY_GENE_START=str2num(VALS{2});continue,end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
72 if strcmp(VALS{1},"SUBSAMPLE"),SUBSAMPLE=str2num(VALS{2});continue,end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
73 if strcmp(VALS{1},"CLIP"),CLIP=str2num(VALS{2});continue,end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
74 if strcmp(VALS{1},"BOOTSTRAP"),BOOTSTRAP=str2num(VALS{2});continue,end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
75 if strcmp(VALS{1},"TEST_METH_NAME"),TEST_METH_NAME=VALS{2};continue,end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
76 if strcmp(VALS{1},"MERGE_SAMPLE"),MERGE_SAMPLE=str2num(VALS{2});continue,end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
77 end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
78
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
79 %Process Bamfiles
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
80 if INTERPR
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
81 BAMS1 = strsplit(BAM_INPUT1,',');
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
82 BAMS2 = strsplit(BAM_INPUT2,',');
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
83 else
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
84 BAMS1=regexp(BAM_INPUT1,',','split');
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
85 BAMS2=regexp(BAM_INPUT2,',','split');
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
86 end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
87
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
88 CFG.BAM_FILES={BAMS1{:},BAMS2{:}};
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
89
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
90 %Name of the experiment. Use the FILENAMES if the entries are empty.
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
91 CFG.NAMES=CFG.BAM_FILES;
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
92 for i=1:length(CFG.NAMES)
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
93 CFG.NAMES{i}=strrep(CFG.NAMES{i},"/","_");
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
94 end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
95
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
96 % Give the directory where the input-files are
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
97 CFG.data_dir = [RDIFF_INPUT_DIR '/'];
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
98
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
99 % Indicate to which sample the bam-files belong
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
100 CFG.SAMPLES=[repmat(1,1,size(BAMS1,2)),repmat(2,1,size(BAMS2,2))];
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
101
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
102 %Process directories
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
103
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
104
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
105
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
106 % Location of the gene structure
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
107 CFG.genes_path=GFF_INPUT;
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
108
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
109 % Output directory
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
110
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
111 CFG.out_base = [RDIFF_RES_DIR '/'];
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
112
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
113 % Output directory for temporary files
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
114 CFG.out_base_temp = [CFG.out_base '/temp/'];
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
115 mkdir(CFG.out_base_temp);
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
116
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
117
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
118 %Check which method to perform
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
119 if strcmp(TEST_METH_NAME,'poisson')
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
120 CFG.perform_poisson=1;
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
121 end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
122 if strcmp(TEST_METH_NAME,'param')
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
123 CFG.perform_parametric=1;
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
124 end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
125 if strcmp(TEST_METH_NAME,'nonparam')
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
126 CFG.perform_nonparametric=1;
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
127 end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
128 if strcmp(TEST_METH_NAME,'mmd')
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
129 CFG.perform_mmd=1;
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
130 end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
131
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
132
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
133 %Process arguments for gene expression estimation
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
134 CFG.estimate_gene_expression=EST_GENE_EXPR;
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
135 CFG.only_gene_expression=ONLY_GENE_EXPR;
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
136
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
137
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
138 %Set options for the variance function
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
139 CFG.merge_sample1=MERGE_SAMPLE;
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
140 CFG.merge_sample2=MERGE_SAMPLE;
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
141
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
142 %If samples contains leass than one sample, merge replicates
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
143 if length(BAMS1)<2
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
144 CFG.merge_sample1=1;
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
145 end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
146 if length(BAMS2)<2
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
147 CFG.merge_sample2=1;
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
148 end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
149
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
150 %Use predefined parameters
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
151 CFG.predefined_variance_function1=[];
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
152 if not(isempty(PRED_VAR1))
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
153 if INTERPR
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
154 VALS = strsplit(PRED_VAR1,',');
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
155 else
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
156 VALS=regexp(PRED_VAR1,',','split');
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
157 end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
158 for i=1:length(VALS)
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
159 CFG.predefined_variance_function1(end+1)=str2num(VALS{i});
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
160 end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
161 end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
162 CFG.predefined_variance_function2=[];
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
163 if not(isempty(PRED_VAR2))
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
164 if INTERPR
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
165 VALS = strsplit(PRED_VAR2,',');
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
166 else
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
167 VALS=regexp(PRED_VAR2,',','split');
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
168 end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
169 for i=1:length(VALS)
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
170 CFG.predefined_variance_function2(end+1)=str2num(VALS{i});
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
171 end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
172 end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
173
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
174
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
175 if not(isempty(SAVE_VAR1))
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
176 CFG.save_variance_function_1=SAVE_VAR1;
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
177 else
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
178 CFG.save_variance_function_1='variance_function_1.mat';
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
179 end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
180 if not(isempty(SAVE_VAR2))
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
181 CFG.save_variance_function_2=SAVE_VAR2;
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
182 else
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
183 CFG.save_variance_function_2='variance_function_2.mat';
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
184 end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
185
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
186 if not(isempty(VAR_PATH1))
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
187 CFG.variance_function_1=VAR_PATH1;
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
188 end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
189 if not(isempty(VAR_PATH2))
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
190 CFG.variance_function_2=VAR_PATH2;
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
191 end
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
192
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
193 %Option not implemented yet
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
194 CFG.compute_variance_function_1=1;
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
195 CFG.compute_variance_function_2=1;
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
196
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
197 %use only gene starts and stops for rDiff.nonparametric variance
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
198 %function esitmation
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
199 CFG.only_gene_start=ONLY_GENE_START;
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
200
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
201 %Process read arguments
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
202 CFG.sequenced_length=READ_LENGTH;
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
203 CFG.min_read_length=min(CFG.sequenced_length,MIN_READ_LENGTH);
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
204
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
205 CFG.rDiff_subsample=SUBSAMPLE;
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
206 CFG.rDiff_nonparametric_subsample_variance_estimation=CFG.rDiff_subsample;
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
207
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
208 CFG.bases_to_clip=CLIP;
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
209
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
210 %Process arguments for rDiff.nonparametric
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
211 CFG.bootstraps=BOOTSTRAP;
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
212
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
213 return