This is a test version of Biostars. For the public version, visit https://www.biostars.org.
New To Limma Package And Need Help With Design Matrix.

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?

r limma microarray agilent bioconductor

Yep, that's what I did to use the modelMatrix function, 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?

0 answers

No answers yet.

Log in to answer this question.