require(GenABEL.data)
data(srdta)
snp.names(srdta, begin = 50000, end = 100000)
snp.names(srdta, begin = 50000, end = 100000, chromosome = "1")
# does not make sense with these data:
snp.names(srdta, begin = 50000, end = 100000, chromosome = "X")
# again makes sense:
snp.names(srdta, end = 100000)
snp.names(srdta, begin = 2200000)
# show summary for SNPs in region between 50,000 and 100,000
a <- snp.names(srdta, begin = 50000, end = 100000)
summary(srdta@gtdata[,a])
Run the code above in your browser using DataLab