Hello Bioconductor,



I have SNP data in HapMap format and would like to read it in R as
snp.Matrix for further analysis. But I have problem to read my hapmap
data using function read.HapMap.data {snpMatrix}




library(chopsticks)

library(hexbin)

snp_data="ftp://ftp.gramene.org/pub/gramene/CURRENT_RELEASE/data/diversi
ty/data_download/hapmap_plink_files/diversity_rice32.bgross.chr12.hapmap
"

result <- read.HapMap.data(snp_data)

error:Unexpected header: rs#	alleles	chrom	pos	strand	assembly#.....

conversion failed.



I have tried to clean up the header with simple names, but it did not
work either.



Any idea how to solve it?



Thanks!



Xiaolian
	[[alternative HTML version deleted]]