Search
Question: using custom BSgenome with MEDIPS
0
gravatar for Vining, Kelly
5.0 years ago by
Vining, Kelly220
Vining, Kelly220 wrote:
________________________________ From: Vining, Kelly Sent: Tuesday, December 17, 2013 10:58 AM To: 'Lukas Chavez' Cc: bioconductor@r-project.org Subject: RE: [BioC] Question about MEDIPS extend parameter Hello BioC group, have created a BSgenome object for my genome of interest, and it is in my current working directory. I now want to use it with MEDIPS. How do I make it appear under >available.genomes()? Is it necessary for my custom genome to appear in this list for me to be able to use it with MEDIPS? Or is there a different method for accessing this object? Thanks much, --Kelly V. [[alternative HTML version deleted]]
ADD COMMENTlink modified 5.0 years ago by Lukas Chavez480 • written 5.0 years ago by Vining, Kelly220
0
gravatar for Lukas Chavez
5.0 years ago by
Lukas Chavez480
USA/La Jolla/UCSD
Lukas Chavez480 wrote:
Hi Kelly, let us assume that the package name of your successfully created custom BSgenome for e.g. 'oryza sativa' is 'BSgenome.OryzaSativa' and that it is installed in your R environment using R CMD INSTALL <package>. First load your custom BSgenome object into your R session using library(BSgenome.OryzaSativa) Note, BSgenome.OryzaSativa will not appear in the list returned by available.genomes(). However, in case the library is loaded you should be able to immediately use it as usual in MEDIPS like library(MEDIPS) BSgenome= "BSgenome.OryzaSativa" MSet = MEDIPS.createSet(file="oryza.bam", BSgenome=BSgenome) Lukas On Wed, Dec 18, 2013 at 6:57 AM, Vining, Kelly <kelly.vining@oregonstate.edu> wrote: > > ________________________________ > From: Vining, Kelly > Sent: Tuesday, December 17, 2013 10:58 AM > To: 'Lukas Chavez' > Cc: bioconductor@r-project.org > Subject: RE: [BioC] Question about MEDIPS extend parameter > > Hello BioC group, > have created a BSgenome object for my genome of interest, and it is in my > current working directory. I now want to use it with MEDIPS. How do I make > it appear under > >available.genomes()? Is it necessary for my custom genome to appear in > this list for me to be able to use it with MEDIPS? Or is there a different > method for accessing this object? > > Thanks much, > --Kelly V. > > > > [[alternative HTML version deleted]] > > _______________________________________________ > Bioconductor mailing list > Bioconductor@r-project.org > https://stat.ethz.ch/mailman/listinfo/bioconductor > Search the archives: > http://news.gmane.org/gmane.science.biology.informatics.conductor > [[alternative HTML version deleted]]
ADD COMMENTlink written 5.0 years ago by Lukas Chavez480
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 160 users visited in the last hour