Skip to content

Commit

Permalink
Update plink.sh
Browse files Browse the repository at this point in the history
  • Loading branch information
zhangdb authored May 20, 2022
1 parent 89cc69c commit 492cba9
Showing 1 changed file with 2 additions and 2 deletions.
4 changes: 2 additions & 2 deletions script/geno_prep/plink.sh
Original file line number Diff line number Diff line change
@@ -1,9 +1,9 @@
## This sciprt will apply plink commands to the files

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

plink --silent --noweb --ped $pedfile --map $mapfile --geno $geno --hwe $hwe --recode --snps-only --list-duplicate-vars suppress-first --out $SIGNET_TMP_ROOT/tmpg/clean_Genotype &&
plink --silent --noweb --file $SIGNET_TMP_ROOT/tmpg/clean_Genotype --mind $mind --exclude $SIGNET_TMP_ROOT/tmpg/clean_Genotype.dupvar --recode --out $SIGNET_TMP_ROOT/tmpg/clean_Genotype --set-hh-missing &&
plink --silent --noweb --file $SIGNET_TMP_ROOT/tmpg/clean_Genotype --recodeA --out $SIGNET_TMP_ROOT/tmpg/clean_Genotype --set-hh-missing &&
tail -n+2 $SIGNET_TMP_ROOT/tmpg/clean_Genotype.raw > $SIGNET_TMP_ROOT/tmpg/clean_Genotype.data &&
echo -e "Plink preprocessing finished\n"
echo -e "Preprocessing genotype data... Completed!\n"

0 comments on commit 492cba9

Please sign in to comment.