SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
EDGE-pro into DESeq, DESeq run error? epistatic Bioinformatics 5 09-02-2015 04:32 AM
interpretation of fastQC result liaojinyue Bioinformatics 1 07-16-2014 12:30 PM
strange plotDispEsts result from Deseq sadiexiaoyu Bioinformatics 2 06-03-2013 12:25 PM
Strange Result in Deseq Output cburke11 Bioinformatics 2 12-04-2012 08:35 AM
interpreting DESeq result, padj - values Azazel Bioinformatics 0 10-06-2010 05:49 PM

Reply
 
Thread Tools
Old 08-07-2014, 04:35 AM   #21
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

Normally you create a character vector containing filenames that you would use (it's usually best to just make a CSV file with all of this and then load it). Otherwise, yes, you have to type it all by hand.
dpryan is offline   Reply With Quote
Old 08-07-2014, 07:43 AM   #22
huma Asif
Member
 
Location: Japan

Join Date: Oct 2010
Posts: 53
Default

while its going let me explain you a problem about DESeq
i tried to use input file with genenames but it was giving me this error
Error in round(countData) : non-numeric argument to mathematical function

so what i did in analysis i cut the column with gene name
now problem is the interpretation of result as in output i have id and i am confused how to combine this info with my gene
i want to attach my count file but here is no option to attach



output (total)
id baseMean baseMeanA baseMeanB foldChange log2FoldChange pval padj
18 32.80688778 4.323701761 61.2900738 14.17537036 3.825314524 0.034361637 0.836923637

output to see just significant
id baseMean baseMeanA baseMeanB foldChange log2FoldChange pval padj
15959 3564.773938 16.43006669 7113.11781 432.9329846 8.757999912 8.02E-20 1.39E-15

problem here is what this id is telling me and how i am going to combine with gene names

please do reply atleast this DESeq output will give me an idea to compare with others
huma Asif is offline   Reply With Quote
Old 08-07-2014, 12:37 PM   #23
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

That ID corresponds to whatever you counted on in your annotation file. If the corresponding names are in the file then just load it (see the GenomicFeatures package) and use that. Alternatively, there are annotation packages for many organisms, so just load the appropriate one in R and use it.
dpryan is offline   Reply With Quote
Old 08-07-2014, 02:04 PM   #24
huma Asif
Member
 
Location: Japan

Join Date: Oct 2010
Posts: 53
Default

thank you dpryan
huma Asif is offline   Reply With Quote
Old 08-08-2014, 05:11 AM   #25
huma Asif
Member
 
Location: Japan

Join Date: Oct 2010
Posts: 53
Default

can you help me how i can get rid of 1,2,3, in first column
> head(rawcounts)
Gene.ID A B
1 WASH7P 1 0
2 FAM138A 0 0
3 FAM138F 0 0
4 OR4F5 0 0
5 LOC100132062 0 0
6 LOC100132287 0 0
i tried
rawcounts <- read.delim("counts.txt", header=TRUE, sep="/t", row.names="Gene.ID")
but it doesnot work
huma Asif is offline   Reply With Quote
Old 08-08-2014, 08:12 AM   #26
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

You might just try:
Code:
rawcounts <- read.delim("counts.txt", header=T, row.names=1)
dpryan is offline   Reply With Quote
Reply

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 01:48 PM.


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