[R] constrast in edge R

[Jeff Newmiller]

You seem to be writing in a jargon that is not typical on this list, and I cannot identify what package you are using. I suspect that you need to be asking your question in the Bioconductor support area.
---------------------------------------------------------------------------
Jeff Newmiller                        The     .....       .....  Go Live...
DCN:<jdnewmil at dcn.davis.ca.us>        Basics: ##.#.       ##.#.  Live Go...
                                      Live:   OO#.. Dead: OO#..  Playing
Research Engineer (Solar/Batteries            O.O#.       #.O#.  with
/Software/Embedded Controllers)               .OO#.       .OO#.  rocks...1k
--------------------------------------------------------------------------- 
Sent from my phone. Please excuse my brevity.

On January 14, 2015 11:41:21 AM PST, Noha Osman <nmo_138 at usc.edu> wrote:
>Hey everyone,
>
>
>I have data frame called subtest as following:
>RNA.LATER.MEN2B_S5 RNA.LATER.ROSA_S4 RNA.MEN2B.1_S2 RNA.MEN2B.2_S3
>RNA.ROSA_S1
>1 13707 13866 12193 12671 10178
>2 0 0 0 0 1
>3 7165 5002 1256 1341 2087
>6 8537 16679 9042 9620 19168
>10 19438 25234 15563 16419 16582
>16 3 3 11 3 5
>I would like to analysis the MEN samples to ROSA samples I did script
>like that .
>> group=c("LMS5","LRS4","MS2","MS3","RS1")
>> y=DGEList(counts=data.matrix(subtest),group=group ,genes=genes)
>> indices=which(rowSums(cpm(y)>1)<3)
>> y=y[-indices,]
>> y=calcNormFactors(y)
>> design=model.matrix(~0+group,data=subtest)
>> cont.matrix=
>makeContrasts(groupLMS5-groupMS2-groupMS3-(groupLRS4-groupRS1),levels=design)
>> fit <- glmFit(y, design,dispersion=y$trended.dispersion)
>> lrt.all <- glmLRT(fit, contrast=cont.matrix)
>> topTags(lrt.all)
>summary(dexp<-decideTestsDGE(lrt.all,p=0.05,adjust="BH"))
>-1 1
>0 36
>1 14248
>I got one gene as down regulated and 14248 as upregulated so I think I
>did something wrong.please Iam a new user in that package and I want to
>get a appropriate number of up and down regulated genes to downstream
>analysis
>
>	[[alternative HTML version deleted]]
>
>______________________________________________
>R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
>https://stat.ethz.ch/mailman/listinfo/r-help
>PLEASE do read the posting guide
>http://www.R-project.org/posting-guide.html
>and provide commented, minimal, self-contained, reproducible code.