Mercurial > repos > elixir-it > vinyl_optimizer
diff old_version/optimizer.pl @ 3:221db2eb3c8e draft default tip
Uploaded
| author | elixir-it |
|---|---|
| date | Wed, 22 Jul 2020 19:23:14 +0000 |
| parents | 6e4eb4856874 |
| children |
line wrap: on
line diff
--- a/old_version/optimizer.pl Wed Jul 22 19:20:30 2020 +0000 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,163 +0,0 @@ -#!/usr/bin/perl -w -use Cwd; - -%arguments= -( -"fileR"=>"", #file: vcf file of affected individuals -"fileC"=>"", #file: vcf file of unaffected individuals -"ofile"=>"", #name: name of the output files -"AD"=>"T", -"XL"=>"F", # -"disease_clinvar"=>[4,6], #numeric mandadory, multiple values -"score_AF"=>[2,4], #numeric mandatory, multiple values -"score_functional"=>[4,6], #numeric mandatory, multiple values -"score_NS"=>[2,4], #numeric mandatory, multiple values -"score_nIND"=>[2,4], #numeric mandatory, multiple values -"scoreeQTL"=>[2,4], #numeric mandatory, multiple values -"scoreG"=>[3,5], #numeric mandatory, multiple values -"scoreI"=>[1..3], -"disease"=>"", #name optional -"similarD"=>"", #file optional -"lgenes"=>"", #file optional -"leQTL"=>"qfile", #file mandatory, but default value -"keywords"=>"kfile", #file mandatory, but default value -"effects"=>"efile", #file mandatory, but default value -"ifile"=>"inter_Hs.file" -"nind"=>5, -"AF"=>0.0001 -); - -@arguments=@ARGV; -for ($i=0;$i<=$#ARGV;$i+=2) -{ - $act=$ARGV[$i]; - $act=~s/-//g; - $val=$ARGV[$i+1]; - if (exists $arguments{$act}) - { - unless ($val=~/\:/) - { - $arguments{$act}=$val; - }else{ - ($s,$e)=(split(/\:/,$val)); - @vls=(); - print "$s $e\n"; - for ($I=$s;$I<=$e;$I++) - { - push(@vls,$I); - } - print "@vls\n"; - $arguments{$act}=\@vls; - } - }else{ - warn("$act: unknown argument\n"); - @valid=keys %arguments; - warn("Valid arguments are @valid\n"); - die("All those moments will be lost in time, like tears in rain.\n Time to die!\n"); - } - #print "$act $val\n"; -} - -%compose_ARG= -( - "disease"=>1, #name optional - "similarD"=>1, #file optional - "lgenes"=>1, #file optional - "leQTL"=>1, #file mandatory, but default value - "keywords"=>1, #file mandatory, but default value - "effects"=>1, #file - "AF"=>1, - "nind"=>1, - "AD"=>1, - "XL"=>1, - "ifile"=>1 - -); - -$dir=getcwd(); -print "$dir\n"; -$add_string=" "; - -foreach $a (keys %compose_ARG) -{ - $add_string.="-$a $arguments{$a} " if $arguments{$a} ne ""; -} - -#print "$add_string"; - - -@DC=@{$arguments{"disease_clinvar"}}; #=>[2..8], #numeric mandadory, but default value -@AF=@{$arguments{"score_AF"}}; #numeric mandatory, but default value -@FUN=@{$arguments{"score_functional"}}; #=>[2..8], #numeric mandatory, but default value -@NS=@{$arguments{"score_NS"}}; #=>[2..8], #numeric mandatory, but default value -@NI=@{$arguments{"score_nIND"}}; #=>[2..8], #numeric mandatory, but default value -@Q=@{$arguments{"scoreeQTL"}}; #numeric mandatory, but default value -@G=@{$arguments{"scoreG"}}; #=>[1..3], #numeric mandatory, but default value -@I=@{$arguments{"scoreI"}}; - -$Rfile=$arguments{"fileR"}; -$Tfile=$arguments{"fileC"}; -$Ofile=$arguments{"ofile"}; -open(O,">$Ofile"); -print O "disease_clinvar\tscore_AF\tscore_functional\tscore_NS\tscore_nIND\tscoreeQTL\tscoreG\tCut_off\tPosD\tPosH\tPvalue\tFold_change\n"; -$commandR="perl $dir/score_complete_alt_M.pl -vcf $dir/$Rfile -ofile $dir/R.csv -ovcfile $dir/R.vcf $add_string"; -$commandT="perl $dir/score_complete_alt_M.pl -vcf $dir/$Tfile -ofile $dir/T.csv -ovcfile $dir/T.vcf $add_string"; - - -foreach $dc (@DC) -{ - - foreach $af (@AF) - { - - foreach $fun (@FUN) - { - - foreach $ns (@NS) - { - - foreach $ni (@NI) - { - - - foreach $Q (@Q) - { - foreach $G (@G) - { - foreach $I (@I) - { - $commandR.="-score_AF $af "; - $commandT.="-score_AF $af "; - $commandR.="-score_functional $fun "; - $commandT.="-score_functional $fun "; - $commandR.="-score_NS $ns "; - $commandT.="-score_NS $ns "; - $commandR.="-score_nIND $ni "; - $commandT.="-score_nIND $ni "; - $commandR.="-scoreeQTL $Q "; - $commandT.="-scoreeQTL $Q "; - $commandR.="-disease_clinvar $dc "; - $commandT.="-disease_clinvar $dc "; - $commandR.="-scoreG $G "; - $commandT.="-scoreG $G "; - $commandR.="-scoreI $I "; - $commandT.="-scoreI $I "; - - #print "$commandR\n$commandT\n"; - system($commandR)==0||die("Non giro Rfile"); - system($commandT)==0||die("Non giro Tfile"); - $commandR="perl $dir/score_complete_alt_M.pl -vcf $dir/$Rfile -ofile $dir/R.csv -ovcfile $dir/R.vcf $add_string"; - $commandT="perl $dir/score_complete_alt_M.pl -vcf $dir/$Tfile -ofile $dir/T.csv -ovcfile $dir/T.vcf $add_string"; - - $r=`Rscript --vanilla $dir/survival_optM.R $dir/R.csv $dir/T.csv`; - - print O "$dc\t$af\t$fun\t$ns\t$ni\t$Q\t$G\t$r\n"; - } - #die(); - } - } - } - } - } - } -}
