Hello, I'm having trouble constructing a proper design matrix for statistical analyses.
I am working with 2-color Agilent array data for multiple arrays. I've bg-corrected and normalized the expression values.
However, the design of the arrays are not uniform.
The goal of the experiment is to see if there's differential expression between pre-treatment and post-treatment with a drug.
However, the microarrays were formatted as follows:
pre = pre- drug treatment
post = post- drug treatment
file Cy3 Cy5
Array 1 file1.txt pre-sample1 pre-sample2
Array 2 file2.txt pre-sample3 post-sample3
Array 3 file3.txt pre-sample4 pre-sample5
Array 4 file4.txt post-sample6 pre-sample6
Array 5 file5.txt post-sample7 pre-sample7
Array 6 file6.txt pre-sample8 post-sample8
I've been following this guide http://koti.mbnet.fi/tuimala/oppaat/r2.pdf and I'm a bit lost to what the design matrix would look entail.
Is there any way to design a matrix with the pre to pre arrays in there or do I need to repeat the analysis twice, one analysis with just Array 1 & 3, a second analysis with the remaining groups?
Any guidance would be appreciated.
edit: is design <- modelMatrix(targets, ref="pre") an accurate design matrix?
0 answers
No answers yet.
Log in to answer this question.
It's not clear exactly what you're trying to compare, but have a look at section 8.3, 8.4, 10.2, 11.2.2 here: http://www.bioconductor.org/packages/release/bioc/vignettes/limma/inst/doc/usersguide.pdf
if you set Cy3 and Cy5, to contain only the factors "pre" and "post", then you should be able to make the design matrix as you suggest.
you probably want your Cy3 and Cy5 columns to be "pre" or "post".
Yep, that's what I did to use the
modelMatrixfunction, but still not sure whether this is the correct design. I just joined this lab, and the pre-pre pairs are a QC measure while the pre-post pairs are to test for differential expression. If that's the case, I would leave the pre-pre pairs out, correct?