How to seperate the groupA and B in Linux

11 views
Skip to first unread message

Bo Zhu

unread,
Aug 19, 2013, 9:49:39 AM8/19/13
to asam...@googlegroups.com
Hi,
  Due to the large Bam file, I am using linux version of SAMMate. I am wondering how it recognize the different groups as in windows version. Because from command line it only mentioned to use the comma.


Best
Bo

Dongxiao Zhu

unread,
Aug 19, 2013, 11:45:25 AM8/19/13
to asam...@googlegroups.com

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.

Tin Nguyen Chi

unread,
Aug 19, 2013, 3:07:14 PM8/19/13
to asam...@googlegroups.com
Hi Bo,

The current linux version of SAMMate does not have the module to group
the input as the Windows version. Sorry for the inconvenience, we hope
to upgrade it soon.

Thanks

Bo Zhu

unread,
Aug 20, 2013, 11:04:47 AM8/20/13
to asam...@googlegroups.com
Thanks for you guys reply. It will be really cool to add this function like Trinity. Meanwhile, I am wondering if next version can add a memory usage on Windows GUI version as it always occupy all my memories...

Best
Bo

在 2013年8月19日星期一UTC+2下午9时07分14秒,Tin Nguyen写道:
Reply all
Reply to author
Forward
0 new messages