Question: How to set the control group in microarray analysis from illumina data?
gravatar for dodausp
5.0 years ago by
dodausp140 wrote:

Hi, there

I am sorry for for the sounded simple question, but how do I set up my control group when processing my data to be used for analysis?

Say I have 3 groups in triplicates (all extracted from Genome Studio), being:

control (C1, C2, C3); Treatment A (A1, A2 ,A3); and Treatment B (B1, B2, B3).

Now, I am able to load the files to R (using Beadarray package), together with the control (probes) table - it generates a EListRaw object. I am also able to run the normalization (neqc), fit linear model (lmFit), and empirical Bayes analysis (eBayes) with no problems. And as expected, you have 3 possible contrast analysis in there: C vs A, C vs B and A vs B.

However, here is my question, how do I compare A vs B where both are already normalized to the control, or put it more simple, how do I compare log2(A/C) vs log2(B/C)? Is there any function in R where I can do this?

Please, any help will be much appreciated.



ADD COMMENTlink modified 5.0 years ago by vassialk190 • written 5.0 years ago by dodausp140
gravatar for vassialk
5.0 years ago by
vassialk190 wrote:

Use Expander, MeV and Genesis and JMP, Cytoscape softwares then go to R and Python. Writing code is a big skill and science. You are not obliged to do it for your task.

When you say "compare", you perhaps say about running ttest or ANOVA on groups or test like LIMMA.

ADD COMMENTlink written 5.0 years ago by vassialk190
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1583 users visited in the last hour