Hello guys,
if anyone knows, could you please tell me why is this happening:
i ran cufflinks on galaxy with default parameters and had satisfactory results.
I then ran the same samples with same parameters except changing max intron length from 300000 to 600000
in the second run have the exact number of transcripts but the FPKM values are much much lower..
Any suggestion?
Thanks,
ibseq
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
how many replicates in each condition do you have?
you could also use SAMseq (samr v2 R-package). this package works with many kinds of designs: paired, quantitative, right censored (like overall survival).
in my hands, SAMseq produced most significant genes (followed by edgeR, baySeq, DESeq, NOIseq, and far far behind cuffdiff) , which were rather robust in bootstrap validations.
my design: 12 normal vs 12 cancer (paired, means from the same patient).
Leave a comment:
-
Cuffdiff supports replicates but does not handle paired replicates to my knowledge.
Btw, I would recommend using DESeq instead of DEGseq, the spelling is similar but the internal statistical modelling is very different.Last edited by Thomas Doktor; 02-02-2012, 04:14 AM.
Leave a comment:
-
I agree that we'd better stick to cuffdiff for differentially expressed gene analysis. Doe cuffdiff have "paired"-analysis feature for the data with replicates? The paired-analysis feature is the main reason I want to use edgeR.
Leave a comment:
-
edgeR and DEGseq take raw counts. They then do their own normalizations. Taking results from cufflinks and trying to use this in any of these programs is not a good approach, even though a lot of people try it for some reason. If you want to use the output of Cufflinks for differential expression, then I would stick to the Cufflinks pipeline and use Cuffdiff.
Otherwise, extract read counts for each gene from your bam/sam/bed file and use this as input for edgeR/DEGseq.
Leave a comment:
-
Cufflinks, differentially expressed genes
Hi,
I am trying to run edgeR or DEGseq using the output from cufflinks.
I usually use mapped reads count as an input to edgeR or DEGseq. What cufflinks output do I need to use for an input to edgeR or DEGseq? I am thinking about adding "coverage" of each isoform for a gene from isoforms.fpkm_tracking file. Does this make sense?
Thank you!Tags: None
Latest Articles
Collapse
-
by seqadmin
The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...-
Channel: Articles
07-08-2024, 03:19 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 07-25-2024, 06:46 AM
|
0 responses
9 views
0 likes
|
Last Post
by seqadmin
07-25-2024, 06:46 AM
|
||
Started by seqadmin, 07-24-2024, 11:09 AM
|
0 responses
26 views
0 likes
|
Last Post
by seqadmin
07-24-2024, 11:09 AM
|
||
Started by seqadmin, 07-19-2024, 07:20 AM
|
0 responses
160 views
0 likes
|
Last Post
by seqadmin
07-19-2024, 07:20 AM
|
||
Started by seqadmin, 07-16-2024, 05:49 AM
|
0 responses
127 views
0 likes
|
Last Post
by seqadmin
07-16-2024, 05:49 AM
|
Leave a comment: