Question: Linkage disequilibrium analysis of EUR populations from 1000 Genomes phase3 data: Considerations of singletons vs family trios
0
gravatar for Scott
4.6 years ago by
Scott80
Canada
Scott80 wrote:

I am interested in calculating LD and R^2 values for a region of interest using Haploview. I am interested in the data from phase 3 of the 1000 Genomes, so I have been loading the data in linkage format from .ped and .info files as downloaded from 1000 Genomes. My question is how the structure of the population (singletons versus family trios) influences LD analysis. When I load files with data from family trios, it treats them as singletons by default. Is it correct to conduct LD analysis with family trios, or should this be avoided? 

ADD COMMENTlink modified 24 months ago by Biostar ♦♦ 20 • written 4.6 years ago by Scott80
3
gravatar for chrchang523
4.6 years ago by
chrchang5234.9k
United States
chrchang5234.9k wrote:

The usual approach is to exclude all "nonfounders"--everyone who has a recorded ancestor also in the dataset, or who has e.g. a shared parent with another person in the dataset (even if said parent wasn't genotyped).  (This is plink's default when calculating LD/r^2.)

ADD COMMENTlink written 4.6 years ago by chrchang5234.9k

Thanks a lot. That was my suspicion. Do you know if there is any easy way to exclude certain individuals from either the generation of the .ped file, or the haploview analysis once loaded? I know you can view the individuals used in Haploview, but I am afraid I will have to edit the .ped file manually. I have all of the ancestry information from the 1000 Genomes site, so that's not a problem. I am just wondering if I have to filter the individuals myself. 

ADD REPLYlink written 4.6 years ago by Scott80
1

If you can get your data into plink .ped/.map format,

plink --file [...] --filter-founders --recode --out [...]

ADD REPLYlink written 4.6 years ago by chrchang5234.9k

Thanks again. I'll check out plink

ADD REPLYlink written 4.6 years ago by Scott80
2
gravatar for Adam
4.5 years ago by
Adam990
United States
Adam990 wrote:

The final Phase 3 release has had 31 known related individuals removed. Details of the removed individuals are provided in the release folder on the FTP site (/vol1/ftp/release/20130502/). The genotypes for these individuals can be found in the supporting subfolders. 

ADD COMMENTlink written 4.5 years ago by Adam990
0
gravatar for Scott
4.5 years ago by
Scott80
Canada
Scott80 wrote:

After checking this manually it looks like all of the VCF files from Phase 3 of the 1000 Genomes Project obtained using the "data slicer" tool contain only non-related individuals. 

ADD COMMENTlink written 4.5 years ago by Scott80
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1008 users visited in the last hour