Question: How to use 1000 Genomes data for LDheatmap package in R
gravatar for mqzhu
5.2 years ago by
United States
mqzhu0 wrote:

I am trying to visualize LD blocks within 1Mb flanking a SNP. And I don't want to use Haploview because it uses Hap Map 3 (build 17 assembly) which is quite outdated. So I downloaded SNP data from 1000 Genomes phase 3, using the online tool "VCF to PED converter". I got .ped and .info files. Then I used an R package ‘LDheatmap’ (which can calculate the LD in r^2 and can visualize LD in heatmap). But the files (.ped, .info files) from 1000 Genomes are not compatible input files for LDheatmap.

The example data set for LDheatmap, "CEUData", contains a data frame and a vector. The format is like this:

  • CEUSNP: A dataframe of SNP genotypes. Each row represents an individual. Each column represents a SNP. SNP IDs are headers of each column.
  • CEUDist: A vector of integers, representing SNP physical map locations on the chromosome.

Does any one know how to convert .ped and .info files from 1000 Genomes into compatible input files (dataframe and vector) for LDheatmap package in R?

heatmap snp ldheatmap lingkage R • 2.8k views
ADD COMMENTlink modified 5.2 years ago • written 5.2 years ago by mqzhu0

You could use Haploview with the downloaded 1000 data.

ADD REPLYlink written 5.2 years ago by Maxime Lamontagne2.2k

Did you ever fix this?

ADD REPLYlink written 2.1 years ago by s.w.vanderlaan40
Please log in to add an answer.


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