LDheatmap包使用

绘制染色体特定片段的LDheatmap
使用LDheatmap包
genetics转换数据格式

snp<-read.csv("片段截取的hmp文件",sep=",",header=T)
library(genetics)
num<-ncol(snp) 
 for(i in 1:num){
   snp[,i]<-as.genotype(snp[,i]) 
# convert the columns of the original data frame into  genotype objects
 }