if(require("PWMEnrich.Dmelanogaster.background")){
data(MotifDb.Dmel.PFM)
PFMtoPWM(MotifDb.Dmel.PFM) # convert to PWM with uniform background
prior = getBackgroundFrequencies("dm3", quick=TRUE) # get background for drosophila (quick mode on a reduced dataset)
PFMtoPWM(MotifDb.Dmel.PFM, prior.params=prior) # convert with genomic background
}
Run the code above in your browser using DataLab