Skip to content
Snippets Groups Projects
Commit 70090f1e authored by Sebastien Moretti's avatar Sebastien Moretti
Browse files

Fix matrix_prob name, path and option

parent 7deeeec5
No related branches found
No related tags found
No related merge requests found
......@@ -138,13 +138,13 @@ if ($motifLib == "promoter"){
}
$command='awk \'BEGIN{print ">log-odds matrix TATA-box"}$1 == "WM"{print $3*100, $4*100, $5*100, $6*100}\' '.$tf_file.' > '.$mat_temp;
exec($command);
$command = '/home/local/perl/matrix_prob.pl -e '.$coEval.' -bg '.$bg.' '.$mat_temp.' 2> /dev/null | grep -o "[0-9]*.[0-9]*%";';
$command = 'matrix_prob -e '.$coEval.' -b '.$bg.' '.$mat_temp.' 2> /dev/null | grep -o "[0-9]*.[0-9]*%";';
$tf_cutoff = shell_exec($command);
$tf_cutoff = str_replace("\n", '', $tf_cutoff);
}else{
$tf_file = "./wwwtmp/wmxMat_".$unique.".wmx";
exec("awk 'BEGIN{RS=\">\";}\$0 ~ \"$tf\"{print \">\"\$0}' $pwmFile > $mat_temp");
$command = '/home/local/perl/matrix_prob.pl -e '.$coEval.' -bg '.$bg.' '.$mat_temp.' 2> /dev/null | grep -o "[0-9]*.[0-9]*%";';
$command = 'matrix_prob -e '.$coEval.' -b '.$bg.' '.$mat_temp.' 2> /dev/null | grep -o "[0-9]*.[0-9]*%";';
#$command1 = 'CO="97.2%";';
#print "$command\n";
$tf_cutoff = shell_exec($command);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment