data(SNPInfo)
Grove ML, Yu B, Cochran BJ, Haritunians T, Bis JC, Taylor KD, Hansen M, Borecki I, Cupples LA, Fornage M, Gudnason V, Harris T, Katherisan S, Kraaij R, Launer LJ, Levy D, Liu Y, Mosley T, Peloso GM, Psaty BM, Rich SS, Rivadeneira F, Siscovick DS, Smith AV, Uitterlinden A, van Duijn CM, Wilson JG, O'Donnell CJ, Rotter JI, Boerwinkle E. Best practices and joint calling of the Illumina HumanExome BeadChip: the CHARGE consortium. PLoS One [submitted]
data(SNPInfo)
##summary of the data set:
summary(as.numeric(table(SNPInfo$SKATgene)))
hist(table(SNPInfo$SKATgene),nclass = 300,xlim=c(0,50),
main = "SNPs per gene", xlab ="#SNPs", ylab = "#Genes")
Run the code above in your browser using DataLab