SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Error in CummeRbund while creating database Richard Barker RNA Sequencing 4 08-05-2013 04:59 PM
Creating a vcf for one gene from scatch Nino Bioinformatics 0 05-20-2013 08:56 AM
gene tracking in cummeRbund stephenhart Bioinformatics 1 03-19-2013 12:09 PM
Creating Gene Sets - CummeRbund SHeaph Bioinformatics 0 01-17-2013 07:56 AM
cufflinks cummeRbund set condition Jetse Bioinformatics 1 12-11-2012 04:48 AM

Reply
 
Thread Tools
Old 06-29-2013, 02:35 AM   #1
uqfgaiti
Member
 
Location: Australia, Brisbane

Join Date: Nov 2012
Posts: 13
Default cummeRbund creating gene set

Hi all,

I am using cummeRbund, an R visualization package for Cufflinks high-throughput sequencing data.

I want to create a gene sets to be able to analyse the differentially expressed genes I'm interested in.
THis is what the manual suggests:

> data(sampleData)
> myGeneIds<-sampleIDs
> myGeneIds

[1] "XLOC_001363" "XLOC_001297" "XLOC_001339" "XLOC_000132"
[5] "XLOC_001265" "XLOC_000151" "XLOC_001359" "XLOC_000069"
[9] "XLOC_000170" "XLOC_000105" "XLOC_001262" "XLOC_001348"
[13] "XLOC_001411" "XLOC_001369" "XLOC_000158" "XLOC_001370"
[17] "XLOC_001263" "XLOC_000115" "XLOC_000089" "XLOC_001240"

> myGenes<-getGenes(cuff,myGeneIds)

where sampleIDs is a vector of 'gene_id' or 'gene_short_name' values to identify the genes I wish to select.

I made a normal .txt file with my genes, which looks like this:

"XLOC_000012"
"XLOC_000033"
"XLOC_000071"
"XLOC_000080"
"XLOC_000081"
"XLOC_000102"
"XLOC_000105"
"XLOC_000113"
"XLOC_000156"
"XLOC_000191"
"XLOC_000212"
"XLOC_000229"
"XLOC_000230"
"XLOC_000241"
"XLOC_000242"
"XLOC_000293"
"XLOC_000294"
"XLOC_000328"
"XLOC_000356"
"XLOC_000398"
"XLOC_000417"
"XLOC_000434"
"XLOC_000442"
etc...

Then, in R, I run these commands:

myGeneIds<-read.table("my txt file")
myGenes<-getGenes(cuff,myGeneIds)
> myGenes
CuffGeneSet instance for 1 genes

and what I got is just the analysis of the first gene in the txt file. But I have 577 genes in the txt file.

How can I make the vector with my gene_id in order to make cummeRbund work correctly?

THanks for help
Federico
uqfgaiti is offline   Reply With Quote
Reply

Tags
cuffdiff, cummerbund

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 04:13 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO