################################################### ### chunk number 1: lkcl eval=FALSE ################################################### ## #line 65 "vignettes/GGBase/inst/doc/newSNPloc.Rnw" ## getClass("snpLocs") ## data(hsSnpLocs) ## hsSnpLocs ################################################### ### chunk number 2: lkc eval=FALSE ################################################### ## #line 73 "vignettes/GGBase/inst/doc/newSNPloc.Rnw" ## snpLocs.Hs(chrnum(20), rsid("rs6060535")) ################################################### ### chunk number 3: doco eval=FALSE ################################################### ## #line 82 "vignettes/GGBase/inst/doc/newSNPloc.Rnw" ## humanSNPlocs = list() ## library(SNPlocs.Hsapiens.dbSNP.20071016) ## if (file.exists("humanSNPlocs.rda")) ## load("humanSNPlocs.rda") else { ## for (i in c(as.character(1:22), "X", "Y")) { ## curc = getSNPlocs(paste("chr", i, sep="")) ## rsid.int = as.integer(curc[,1]) ## loc.int = as.integer(curc[,3]) ## humanSNPlocs[[i]] = rbind(rsid=rsid.int, loc=loc.int) ## # cat(i) ## } ## } ################################################### ### chunk number 4: doco2 eval=FALSE ################################################### ## #line 99 "vignettes/GGBase/inst/doc/newSNPloc.Rnw" ## require(org.Hs.eg.db) ## chrl = org.Hs.egCHRLENGTHS ## offs = c(0, cumsum(as.double(chrl[1:22]))) ## ################################################### ### chunk number 5: junk eval=FALSE ################################################### ## #line 105 "vignettes/GGBase/inst/doc/newSNPloc.Rnw" ## #setClass("snpLocs", representation(locEnv="environment", ## # offsets="numeric", organism="character", versions="character")) ## # ## #setMethod("show", "snpLocs", function(object) { ## # cat("snpLocs instance, organism ", object@organism, "\n") ## # cat("based on:\n") ## # print(object@versions) ## #}) ################################################### ### chunk number 6: docont eval=FALSE ################################################### ## #line 117 "vignettes/GGBase/inst/doc/newSNPloc.Rnw" ## el = new.env() ## getv = function(x) installed.packages()[x, "Version"] ## for (i in names(humanSNPlocs)) ## assign(i, humanSNPlocs[[i]], el) ## hsSnpLocs = new("snpLocs", locEnv=el, offsets=offs, ## organism="Hs", versions=c( ## org.Hs.eg.db=getv("org.Hs.eg.db"), ## SNPlocs.Hsapiens.dbSNP.20071016 = getv("SNPlocs.Hsapiens.dbSNP.20071016")))