### R code from vignette source 'vignettes/GGBase/inst/doc/newSNPloc.Rnw' ################################################### ### code chunk number 1: lkcl (eval = FALSE) ################################################### ## getClass("snpLocs") ## data(hsSnpLocs) ## hsSnpLocs ################################################### ### code chunk number 2: lkc (eval = FALSE) ################################################### ## snpLocs.Hs(chrnum(20), rsid("rs6060535")) ################################################### ### code chunk number 3: doco (eval = FALSE) ################################################### ## humanSNPlocs = list() ## library(SNPlocs.Hsapiens.dbSNP.20090506) ## 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) ## } ## } ################################################### ### code chunk number 4: doco2 (eval = FALSE) ################################################### ## require(org.Hs.eg.db) ## chrl = org.Hs.egCHRLENGTHS ## offs = c(0, cumsum(as.double(chrl[1:22]))) ## ################################################### ### code chunk number 5: junk (eval = FALSE) ################################################### ## #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) ## #}) ################################################### ### code chunk number 6: docont (eval = FALSE) ################################################### ## 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.20090506 = getv("SNPlocs.Hsapiens.dbSNP.20090506")))