Skip to content

Commit

Permalink
Update gexp_prep_gtex.sh
Browse files Browse the repository at this point in the history
  • Loading branch information
zhangdb authored May 20, 2022
1 parent d76aef4 commit 82fba7d
Showing 1 changed file with 3 additions and 4 deletions.
7 changes: 3 additions & 4 deletions script/gexp_prep/gexp_prep_gtex.sh
Original file line number Diff line number Diff line change
Expand Up @@ -7,13 +7,13 @@ cp $SIGNET_TMP_ROOT/tmpg/subjid_wb_common.txt subjid_wb_common.txt

awk '{print $2}' $reads | tail -n+4 > GTEx_genename

echo -e "Looking for the samples for $tissue tissue \n"
echo -e "Looking for the samples for $tissue tissue... \n"

Rscript $SIGNET_SCRIPT_ROOT/gexp_prep/gexp_prep_gtex.R "reads='$reads'" "tpm='$tpm'" "tissue='$tissue'" "anno='$anno'" &&

tabix -l $SIGNET_TMP_ROOT/tmpg/GTEx_snp.vcf.gz > vcf_chr_list

echo -e "Begin normalizing \n"
echo -e "Normalizing gene expressions...\n"

python3 $SIGNET_SCRIPT_ROOT/gexp_prep/eqtl_prepare_expression_igt2log.py $tpm \
$reads \
Expand All @@ -25,5 +25,4 @@ python3 $SIGNET_SCRIPT_ROOT/gexp_prep/eqtl_prepare_expression_igt2log.py $tpm \
--sample_id_list sample_ids_tissue.txt


echo -e "Normalizing finished \n"

echo -e "Normalizing gene expressions... Completed!\n"

0 comments on commit 82fba7d

Please sign in to comment.