Ld Analysis

Question: Ld Analysis

8.3 years ago by

mary • 210

Bologna university

mary • 210 wrote:

Dear Friend

I have a population with 1066 bulls that just sire was genotyped by Illumina Bovine SNP50K BeadChip . my data is not trios, just in some animal, father and son was genotyped, but I have pedigree for all animals. I used plink for calculate r2 chromosome by chromosome in first step. as I am new in LD analysis, i confused in analysis , I want to know:

1- is pedigree important in this data analysis?

2- how can i calculate D' for my data?

3- when I used plink for haplotype reconstruction (phasing) as a command "plink --bfile mydata --blocks " the file that it gave me was empty??!?! i don't know its correct or not?

is there any one that can help me?

linkage ld • 3.8k views

ADD COMMENT • link •

modified 8.3 years ago • written 8.3 years ago by mary • 210

8.3 years ago by

Zev.Kronenberg ♦ 11k

United States

Zev.Kronenberg ♦ 11k wrote:

This article is a great place to start:

http://www.ncbi.nlm.nih.gov/pmc/articles/PMC1226024/

You don't need family / related individuals to calculate LD. Using related individuals will inflate your LD estimate for a population.

ADD COMMENT • link modified 8.3 years ago • written 8.3 years ago by Zev.Kronenberg ♦ 11k

8.3 years ago by

Bioch'Ti • 1.0k

France (Avignon)

Bioch'Ti • 1.0k wrote:

Hi Maryam,

1 - Pedigree is not important in LD analysis, but will, as Zev told you, inflate LD estimates. Recently, an [R] package LDcorsSV (Mangin et al., 2011, Heredity) taking into account Kinship and/or population structure was released. You may try it. A bit long computing time but it works.

2 - You can use different softwares to compute LD estimation but check for which estimates is the most appropriate one (D', r, r2). One of the most used software is 'Plink!' but you can also use 'Tassel' for example.

3 - The bfile option in Plink needs a special format if I remember well.

ADD COMMENT • link written 8.3 years ago by Bioch'Ti • 1.0k

8.3 years ago by

mary • 210

Bologna university

mary • 210 wrote:

Thanks a lot

1- I prefer don't include the pedigree, because I want compare two breed that for second breed i don't have pedigree.

2- I think the plink just calculate r and r2. i can do haploview, but i don't know I need do haplotype reconstruction before input data to haploview or not? (i can do haplotype Reconstruction by fastPHASE), but plink is easy and fast

3- I did plink --file mydata --make-bed --out mydata, so i had special format for do blocks command. but i didn't understand why plink doesn't do that?

thank you very much for help

ADD COMMENT • link written 8.3 years ago by mary • 210

Please log in to add an answer.

Similar posts • Search »