Is there a way to get a list of all homogygous sites for each sample with plink (regardless if they are in a run)
2
0
Entering edit mode
13 months ago
curious ▴ 810

With plink I understand that you get runs of homogygosity with --homozyg, but I would like to a list of all sites that are homozygous for each sample. Is there any way to do this?

plink • 1.4k views
ADD COMMENT
0
Entering edit mode
13 months ago
bk11 ★ 3.0k

You can run like this-

plink --file your_data --homozyg --homozyg-group

Please check the Runs of homozygosity in the link below-

https://zzz.bwh.harvard.edu/plink/ibdibs.shtml#homo

ADD COMMENT
0
Entering edit mode

I am a bit confused though, I am looking for a list of all sitse that are homozygous for A1 regardless of whether they are in a run of homozygosity. If I understand correctly that command gives me a file plink.hom.overlap that provides information on overlapping ROHs. Is there some way I can parse plink.hom.overlap to get what I need, I can't quite see how.

ADD REPLY
0
Entering edit mode

Please keep reading down the documentation listed above.

ADD REPLY
0
Entering edit mode
12 months ago

For one sample, a simple way to do this is to filter for minor-allele-count <= 0:

plink2 --bfile ... --indv <sample ID> --max-mac 0 --write-snplist

For many samples, you can use a shell-script for-loop, or plink2's --loop-cats flag. (Given a choice, I recommend learning about shell-script for-loops first, since they have many other applications.)

ADD COMMENT

Login before adding your answer.

Traffic: 2136 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6