Calculate linkage disequilibrium from Hudson's ms simulated binary data
Hello, I have simulated genetic data from Hudson's ms (and COSI) software in this form:
001010
110101
010100
001010
010100
where 0 and 1 represent the ancestral and derived allele, respectively. Each line represents a chromosome from an individual, whereas each column represents a SNP. Is there any way of calculating LD-related measures from this data, such as D' or r2?. I know how to do it from real genetic data (A,C,G,T) but I have not found how to do it from binary (0,1) data. Any suggestion is very welcome !
• 2,564 views
•
link
updated
by
Biostar
•
written
by
Mr Locuace
•
written
by
Mr Locuace
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Blog: The Complete Guide to Metatranscriptome Sequencing
written by Novogene •
-
Linkage disequilibrium from VCF file of variants in bacterial metagenomes
written by matthias.hotzinger •
-
Convert vcf files with phased genotypes to standard haplotype format
written by Mr Locuace
-
Tutorial: It's about the journey, and not the destination
written by theobroma22
-
Getting derived alleles from alternate/reference allele data
written by Mr Locuace
For the number of recombination events I use RecMin. The input file format is quite similar to the data you have. http://www.stats.ox.ac.uk/~myers/RecMin.html
Thanks very much @microfuge .I will try it !. Can it calculate r2 or D'?