Thanks Likun,
I was confused by the q-value. It all makes sense now.
- Nick
I was confused by the q-value. It all makes sense now.
- Nick
You are currently viewing the SEQanswers forums as a guest, which limits your access. Click here to register now, and join the discussion
Control Gene Count RPKM total no. reads gene length ABAT 811 1.57049 87082676 5930 Mutant ABAT 4264 10.3046 69779755 5930
library(DEGseq) sample1 <- "D:/data/sample1.txt" sample2 <- "D:/data/sample2.txt" refFlat <- "D:/data/refFlat.txt" mapResultBatch1 <- c(sample1) mapResultBatch2 <- c(sample2) outputDir <- "D:/data/DEGseqExample" DEGseq(mapResultBatch1, mapResultBatch2, fileFormat = "eland",refFlat = refFlat, outputDir = outputDir, method = "LRT")
hist(LogVal(Sample1),main=label1,xlab="log2(Number of reads mapped to a gene)",col=4,breaks=100,freq=FALSE,ylim=c(0,0.5))
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, Today, 11:49 AM
|
0 responses
13 views
0 likes
|
Last Post
by seqadmin
Today, 11:49 AM
|
||
Started by seqadmin, Yesterday, 08:47 AM
|
0 responses
16 views
0 likes
|
Last Post
by seqadmin
Yesterday, 08:47 AM
|
||
Started by seqadmin, 04-11-2024, 12:08 PM
|
0 responses
61 views
0 likes
|
Last Post
by seqadmin
04-11-2024, 12:08 PM
|
||
Started by seqadmin, 04-10-2024, 10:19 PM
|
0 responses
60 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 10:19 PM
|
Comment