Skip to content

Commit

Permalink
Update plink_gtex.sh
Browse files Browse the repository at this point in the history
  • Loading branch information
zhangdb authored May 20, 2022
1 parent 492cba9 commit 78ce662
Showing 1 changed file with 2 additions and 2 deletions.
4 changes: 2 additions & 2 deletions script/geno_prep/plink_gtex.sh
Original file line number Diff line number Diff line change
Expand Up @@ -2,7 +2,7 @@

cd $SIGNET_TMP_ROOT/tmpg

echo -e "Preprocessing using Plink ...\n"
echo -e "Preprocessing genotype data with PLINK...\n"

plink --vcf GTEx_snp.vcf --double-id --make-bed --missing --freq --recode --out Geno_GTEx_$tissue
plink --bfile Geno_GTEx_$tissue --recodeA --mac 5 --out GTEx_$tissue
Expand All @@ -15,4 +15,4 @@ awk {'print $5'} GTEx_$tissue.frq | tail -n+2 > ${resg}_snps.maf
tail -n+2 GTEx_$tissue.raw |cut -d " " -f 7- > clean_Genotype.data
sed -e 's/NA/0/g' clean_Genotype.data > ${resg}_clean_Genotype_repNA.data

echo -e "Plink preprocessing finished\n"
echo -e "Preprocessing genotype data ... completed!\n"

0 comments on commit 78ce662

Please sign in to comment.