Entering edit mode
Xiaolian Yuan
▴
30
@xiaolian-yuan-5974
Last seen 10.3 years ago
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/diver
si
ty/data_download/hapmap_plink_files/diversity_rice32.bgross.chr12.hapm
ap
"
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]]