excludeDoubleMatchingProbes: Exclude ProbeLevelSet probes that match more than once in genome

Description

Function that will remove probes from ProbeLevelSet if they have more than one match in a given genome.

Usage

excludeDoubleMatchingProbes(object, genome="BSgenome.Hsapiens.UCSC.hg18", verbose=TRUE, directions=c("matchForwardSense", "matchForwardAntisense", "matchReverseSense", "matchReverseAntisense"), previousData = NULL)

Arguments

object

A ProbeLevelSet class object.

genome

character string with the name of the BSGenome in which sequences should be found. Defaults to the human genome.

verbose

TRUE or FALSE.

directions

character string with elements from c("matchForwardSense", "matchForwardAntisense", "matchReverseSense", "matchReverseAntisense"). Defines which directions (complementary and reverse mirrorings) that should be scanned. Defaults to all directions.

previousData

Optional: The output from a call to findSequenceInGenome. If given the scanning will be skipped, and the probes will be omitted directly. Useful in cases were datasets from the same region needs to be processed.

Value

Details

This function will take quite a while to run, so if you have many sequences, overnight runs are recommended. BSgenome contains some alternative versions of chromosomes. They are marked with an underscore. This function automatically disregards chromosome names with an underscore, and this is known to work for the human genome. Nevertheless, check the output printed to terminal if all chromosomes are included. The function is a wrapper around findSequenceInGenome, which can be used for purposes that are more flexible (although that function really is just following the example in the BSgenome package)

At present, there is no functionality to check matches with known SNP or known splice forms taken into account.

Examples

Run this code


	## Not run: 
# 	#you can run this, but it takes a lot of time
# 	probelevelsetwithnodoubles<-excludeDoubleMatchingProbes(probelevelset)
# 	## End(Not run)

Run the code above in your browser using DataLab

Last chance! 50% off unlimited learning