BPFrequency<- | R Documentation |
BPFrequency
slot in a
genomicProfiles
object.
Setter method for BPFrequency
slot in a
genomicProfiles
object.
If base pair frequency is unknown, BPFrequency
will compute base pair
frequency from a DNA sequence.
BPFrequency(object)<-value
object |
|
value |
A vector of length 4 containing the probability of occurrence of each base pair (A,C,G,T) in order. Default value is c(0.25,0.25,0.25,0.25). A A |
Default value is c(0.25,0.25,0.25,0.25)
When generating a Postion Weight Matrix from a Position Frequency Matrix,
the probability of occurrence of each base pair (Base Pair Frequency) is
necessary (as originally described by Gary Stormo). It is possible to
set custom values for BPFrequency
with a vector of length 4
containing the probability of occurrence of each base pair (A,C,G,T) in order.
If Base pair frequency is unknown, BPFrequency
will compute base pair
frequency from a DNA sequence when building a
genomicProfiles
object.
The nature of this sequence can be aBSgenome
object or a
DNAStringSet
. In order to decrease run time,
it is advised to use DNAStringSet
.
Returns a genomicProfiles
object with an updated
value for BPFrequency
.
Patrick C.N. Martin <pcnmartin@gmail.com>
Zabet NR, Adryan B (2015) Estimating binding properties of transcription factors from genome-wide binding profiles. Nucleic Acids Res., 43, 84–94. Patrick C.N. Martin and Nicolae Radu Zabe (2020) Dissecting the binding mechanisms of transcription factors to DNA using a statistical thermodynamics framework. CSBJ, 18, 3590-3605.
data(ChIPanalyserData) # path to Position Frequency Matrix PFM <- file.path(system.file("extdata",package="ChIPanalyser"),"BEAF-32.pfm") #As an example of genome, this example will run on the Drosophila genome if(!require("BSgenome.Dmelanogaster.UCSC.dm6", character.only = TRUE)){ if (!requireNamespace("BiocManager", quietly=TRUE)) install.packages("BiocManager") BiocManager::install("BSgenome.Dmelanogaster.UCSC.dm6") } library(BSgenome.Dmelanogaster.UCSC.dm6) DNASequenceSet <- getSeq(BSgenome.Dmelanogaster.UCSC.dm6) # Building genomicProfiles object GPP <- genomicProfiles(PFM=PFM,PFMFormat = "JASPAR", BPFrequency=DNASequenceSet) # Updating BPFrequency ## !! Note!! BPFrequency is used to compute PWM from PFM ## IF updated after building GPP, then it will not influence PWM ## Advised to build with BPFrequency directly BPFrequency(GPP) <- DNASequenceSet BPFrequency(GPP) <- c(0.25,0.25,0.25,0.25)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.