Limit Findsimilar() To A Subset Of Genes In The Whole Cufflinks Data Set [Cummerbund] · Biostar

This is a test version of Biostars. For the public version, visit https://www.biostars.org.

Limit Findsimilar() To A Subset Of Genes In The Whole Cufflinks Data Set [Cummerbund]

Hello folks,

I am running into what I hope is not a fundamental limitation of cummeRbund's findSimilar() function.

findSimilar is supposedly an object method of class CuffSet (itself being basically an interface to the SQLite DB tables). So I can't really figure out a straightforward way to run findSimilar() on a set of genes that is LESS than the entire cufflinks experiment. However, I am a Python guy. R is a necessary tool that I use, but I am not fluent. Basically, at the conversational level, probably. So maybe (I hope) I am missing something.

A quick description of what I would like to do:

I have a list of genes that are 1-to-1 orthologs in multiple species
One of these genes has an RNA abundance profile that I am keen on
I would like to pull out the top 'n' genes with similar abundance profiles THAT ARE ALSO 1-to-1 ORTHOLOGS like geneX
I have a list of these 1-to-1 ortholog genes and a CuffGeneSet object initialized with them

Can anyone help me figure this out?

Thanks,

Gus

cummerbund rna-seq

• 2,429 views • link

updated 12.3 years ago by Biostar

• written 12.9 years ago by wadunn83 •

0 answers

No answers yet.

Log in to answer this question.

More posts like this

SNP calling pipeline - testing different programs
written 11.2 years ago by Eric Normandeau

<p>I am evaluating different possible pipelines to perform the following:</p> <ul> <li>Map Ion Proton reads on a reference genome</li> <li>Call SNPs</li> <li>Filter the SNPs</li> </ul> …
Cummerbund Help
written 12.5 years ago by madkitty

<p>I have 3 samples in 3 folders under the working directory of R,</p> <ul> <li>~/Sample1 </li> <li>~/Sample2 </li> <li>~/Sample3</li> </ul> <p>How can I point to …
Cogent.Db.Ensembl Cookbook Example Returns 'Programmingerror', New Bug?
written 12.6 years ago by a1ultima

<p>I get an error using PyCogent to query EnsEMBl's database - I suspect this is a new bug resulting from either a new update in …
Codon Alignment Via Python?
written 12.6 years ago by a1ultima

<p>I have pairs of coding DNA sequences which I wish to perform pairwise <strong>codon alignments</strong> via <strong>Python</strong>, I have "half completed" the process.</p> <p><strong>So far..</strong></p> …
Python 2 Vs Python 3 In Bioinformatics
written 12.7 years ago by Eric Normandeau

<p>Hi,</p> <p>I have postponed switching to Python 3 for a few years now for the reasons listed below.</p> <p><strong>Pros of Python 2</strong></p> <ul> <li>Used by …
Tool: Bedops V2.2 Released
written 13.3 years ago by Alex Reynolds

<p><strong><a href='http://code.google.com/p/bedops/'>BEDOPS</a></strong> is a suite of tools to address common questions raised in genomic studies — mostly with regard to overlap and proximity relationships between …
Convert Bam To Fastq With Pe & Se Reads
written 14.4 years ago by Pasta •

<p>Hi there,</p> <p><strong>EDIT</strong></p> <p>I filtered my PE reads based on quality. Filtering generated 3 kind of files: </p> <ul> <li>1 file with PE of good …
New To Limma And Linear Modeling, Is Limma Appropriate For My Data?
written 14.4 years ago by Ric •

<p>I have a fairly complex microarray experiment for which I think limma might be a good approach, but I’m having some difficultly figuring out how …
Well-Annotated Functional Pathway
written 14.6 years ago by Andrew Su

<p><strong>The dilemna</strong>: I am submitting a grant proposal <em>tomorrow morning</em>. However, inspiration hit this morning and I need to make a new figure on short …
Are There Any Substitutes For Sequencher By Gene Codes?
written 15.5 years ago by Blunders •

<p>If suggesting a substitute, please included the following information if possible.</p> <ul> <li>URL</li> <li>Platforms</li> <li>Version Number</li> <li>Owned by</li> <li>License type</li> <li>Source Code</li> <li>Date of last …