Hi,
Trying to do the following and getting an error of Model matrix not full rank:
sampleNames<-c("MMV1", "MMV2", "MMV3", "FMV1", "FMV2", "FMV3", "MME7", "MME8", "MME9", "FME1", "FME2", "FME3")
sampleSex<-c("male", "male", "male", "female", "female", "female", "male", "male", "male", "female", "female", "female")
sampleTreatment<-c("vehicle", "vehicle", "vehicle", "vehicle", "vehicle", "vehicle", "BB", "BB", "BB", "BB", "BB", "BB")
sampleBatch<-c("1","2","3","1","2","3","4","5","6","4","5","6")
colData<-data.frame(sampleName=sampleNames, sex=sampleSex, treatment=sampleTreatment, batch=sampleBatch)
se<-SummarizedExperiment(assays=list(counts=counts),rowRanges=rowRanges, colData=colData)
table(sampleSex, sampleTreatment)
dds<-DESeqDataSet(se, design= ~ batch + sex + treatment)
converting counts to integer mode
Error in checkFullRank(modelMatrix) :
the model matrix is not full rank, so the model cannot be fit as specified.
One or more variables or interaction terms in the design formula are linear
combinations of the others and must be removed.
Please read the vignette section 'Model matrix not full rank':
> sessionInfo()
R version 3.4.1 (2017-06-30)
Platform: x86_64-apple-darwin15.6.0 (64-bit)
Running under: macOS Sierra 10.12.5
Matrix products: default
BLAS: /System/Library/Frameworks/Accelerate.framework/Versions/A/Frameworks/vecLib.framework/Versions/A/libBLAS.dylib
LAPACK: /Library/Frameworks/R.framework/Versions/3.4/Resources/lib/libRlapack.dylib
locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
attached base packages:
[1] parallel stats4 stats graphics grDevices utils datasets methods base
other attached packages:
[1] DiffBind_2.4.8 ggplot2_2.2.1 DESeq2_1.16.1
[4] chipenrich_2.0.1 rtracklayer_1.36.6 TxDb.Mmusculus.UCSC.mm10.knownGene_3.4.0
[7] org.Mm.eg.db_3.4.1 TxDb.Mmusculus.UCSC.mm10.ensGene_3.4.0 GenomicFeatures_1.28.5
[10] AnnotationDbi_1.38.2 csaw_1.10.0 BiocParallel_1.10.1
[13] SummarizedExperiment_1.6.5 DelayedArray_0.2.7 matrixStats_0.52.2
[16] Biobase_2.36.2 GenomicRanges_1.28.6 GenomeInfoDb_1.12.3
[19] IRanges_2.10.5 S4Vectors_0.14.7 BiocGenerics_0.22.1
[22] BiocInstaller_1.26.1
loaded via a namespace (and not attached):
[1] amap_0.8-14 TH.data_1.0-8 colorspace_1.3-2 rjson_0.2.15 hwriter_1.3.2 htmlTable_1.11.1
[7] XVector_0.16.0 base64enc_0.1-3 rstudioapi_0.7 MatrixModels_0.4-1 ggrepel_0.7.0 bit64_0.9-7
[13] mvtnorm_1.0-6 org.Dr.eg.db_3.4.1 codetools_0.2-15 splines_3.4.1 chipenrich.data_2.0.0 geneplotter_1.54.0
[19] knitr_1.17 Formula_1.2-2 Rsamtools_1.28.0 annotate_1.54.0 cluster_2.0.6 GO.db_3.4.1
[25] pheatmap_1.0.8 graph_1.54.0 readr_1.1.1 compiler_3.4.1 GOstats_2.42.0 backports_1.1.2
[31] assertthat_0.2.0 Matrix_1.2-12 lazyeval_0.2.1 limma_3.32.10 org.Rn.eg.db_3.4.1 quantreg_5.34
[37] acepack_1.4.1 htmltools_0.3.6 tools_3.4.1 bindrcpp_0.2 gtable_0.2.0 glue_1.2.0
[43] GenomeInfoDbData_0.99.0 Category_2.42.1 systemPipeR_1.10.2 dplyr_0.7.4 ShortRead_1.34.2 Rcpp_0.12.14
[49] Biostrings_2.44.2 nlme_3.1-131 gdata_2.18.0 stringr_1.2.0 gtools_3.5.0 XML_3.98-1.9
[55] polspline_1.1.12 org.Hs.eg.db_3.4.1 edgeR_3.18.1 MASS_7.3-47 zoo_1.8-0 zlibbioc_1.22.0
[61] scales_0.5.0 hms_0.4.0 sandwich_2.4-0 RBGL_1.52.0 SparseM_1.77 RColorBrewer_1.1-2
[67] BBmisc_1.11 memoise_1.1.0 gridExtra_2.3 rms_5.1-1 biomaRt_2.32.1 rpart_4.1-11
[73] latticeExtra_0.6-28 stringi_1.1.6 RSQLite_2.0 Rhtslib_1.8.0 genefilter_1.58.1 checkmate_1.8.5
[79] caTools_1.17.1 rlang_0.1.6 pkgconfig_2.0.1 BatchJobs_1.7 bitops_1.0-6 lattice_0.20-35
[85] bindr_0.1 GenomicAlignments_1.12.2 htmlwidgets_0.9 bit_1.1-12 GSEABase_1.38.2 AnnotationForge_1.18.2
[91] plyr_1.8.4 magrittr_1.5 sendmailR_1.2-1 R6_2.2.2 gplots_3.0.1 Hmisc_4.0-3
[97] multcomp_1.4-8 DBI_0.7 mgcv_1.8-22 pillar_1.1.0 foreign_0.8-69 survival_2.41-3
[103] RCurl_1.95-4.8 nnet_7.3-12 tibble_1.4.1 org.Dm.eg.db_3.4.1 KernSmooth_2.23-15 locfit_1.5-9.1
[109] grid_3.4.1 data.table_1.10.4-3 blob_1.1.0 digest_0.6.14 xtable_1.8-2 brew_1.0-6
[115] munsell_0.4.3
vignette('DESeq2')

In addition this should be insightful: