################################################### ### chunk number 1: loadlibs ################################################### #line 43 "vignettes/annotate/inst/doc/useProbeInfo.Rnw" library("annotate") library("rae230a.db") library("rae230aprobe") ################################################### ### chunk number 2: selprobe ################################################### #line 55 "vignettes/annotate/inst/doc/useProbeInfo.Rnw" ps = names(as.list(rae230aACCNUM)) myp = ps[1001] myA = get(myp, rae230aACCNUM) wp = rae230aprobe$Probe.Set.Name == myp myPr = rae230aprobe[wp,] ################################################### ### chunk number 3: getACC ################################################### #line 87 "vignettes/annotate/inst/doc/useProbeInfo.Rnw" myseq = getSEQ(myA) nchar(myseq) library("Biostrings") mybs = DNAString(myseq) match1 = matchPattern(as.character(myPr[1,1]), mybs) match1 as.matrix(match1) myPr[1,5] ################################################### ### chunk number 4: getRev ################################################### #line 108 "vignettes/annotate/inst/doc/useProbeInfo.Rnw" myp = ps[100] myA = get(myp, rae230aACCNUM) wp = rae230aprobe$Probe.Set.Name == myp myPr = rae230aprobe[wp,] myseq = getSEQ(myA) mybs = DNAString(myseq) Prstr = as.character(myPr[1,1]) match2 = matchPattern(Prstr, mybs) ## expecting 0 (no match) length(match2) match2 = matchPattern(reverseComplement(DNAString(Prstr)), mybs) nchar(match2) nchar(myseq) - as.matrix(match2) myPr[1,5] ################################################### ### chunk number 5: ################################################### #line 160 "vignettes/annotate/inst/doc/useProbeInfo.Rnw" sessionInfo()