Hi, SAMMate uses edgeR Bioconductor package to detect differentially expressed genes. After you get read counts, please load edgeR and follow the commands below. It represents two groups with two replicates in each group.
Best,
Dongxiao
> x <- read.delim("fileofcounts.txt",row.names="Symbol")
> group <- factor(c(1,1,2,2))
> y <- DGEList(counts=x,group=group)
> y <- calcNormFactors(y)
> y <- estimateCommonDisp(y)
> y <- estimateTagwiseDisp(y)
> et <- exactTest(y)
> topTags(et)
--
You received this message because you are subscribed to the Google Groups "asammate" group.
To unsubscribe from this group and stop receiving emails from it, send an email to asammate+u...@googlegroups.com.
For more options, visit https://groups.google.com/groups/opt_out.