[BioC] DE in EdgeR

R [guest] guest at bioconductor.org
Thu Nov 1 11:38:53 CET 2012


I am testing DE between two groups that each contain 10 time points. 

I am creating my design matrix like this:

Condition <- factor(c("KO","KO","KO","KO","KO","KO","KO","KO","KO","WT","WT","WT","WT","WT","WT","WT","WT","WT"))
Condition <- relevel(Condition, ref="WT")
Time <- factor(c("0H","3H","6H","9H","12H","18H","21H","24H","27H","0H","3H","6H","9H","12H","18H","21H","24H","27H"))

design <- model.matrix(~Time+Time:Condition)


How do I proceed if i want to compare the whole group (WT) against the other group (KO) considering the time points as replicates. 



 -- output of sessionInfo(): 

R version 2.14.1 (2011-12-22)
Platform: x86_64-pc-linux-gnu (64-bit)

locale:
 [1] LC_CTYPE=nb_NO.UTF-8       LC_NUMERIC=C               LC_TIME=nb_NO.UTF-8       
 [4] LC_COLLATE=nb_NO.UTF-8     LC_MONETARY=nb_NO.UTF-8    LC_MESSAGES=nb_NO.UTF-8   
 [7] LC_PAPER=C                 LC_NAME=C                  LC_ADDRESS=C              
[10] LC_TELEPHONE=C             LC_MEASUREMENT=nb_NO.UTF-8 LC_IDENTIFICATION=C       

attached base packages:
[1] splines   stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
 [1] BiocInstaller_1.2.1   DESeq_1.6.1           locfit_1.5-8          org.Mm.eg.db_2.6.4   
 [5] RSQLite_0.11.2        DBI_0.2-5             AnnotationDbi_1.16.19 Biobase_2.14.0       
 [9] edgeR_2.4.6           limma_3.10.3         

loaded via a namespace (and not attached):
 [1] annotate_1.32.3    genefilter_1.36.0  geneplotter_1.32.1 grid_2.14.1       
 [5] IRanges_1.12.6     lattice_0.20-0     RColorBrewer_1.0-5 survival_2.36-10  
 [9] tools_2.14.1       xtable_1.7-0

--
Sent via the guest posting facility at bioconductor.org.



More information about the Bioconductor mailing list