[BioC] 2-way Anova analysis

[James W. MacDonald]

Sharon wrote:
>  Hi All,
> 
> I am trying to analyze affy data for strain, treatment and their interaction
> effects (RMA normalized, log2 expression values). The design looks like
> 
>    strain A --> treatment 1       strain B --> treatment 1
>                     treatment 2                         treatment 2
> 
> so, I used anova(aov(expra ~ strain + treatment + strain*treatment, m)),
> where m <- data.frame(strain, treatment, expr)
> 
> I have got the corresponding pvalues for strain, treatment & interaction
> effects( this for all the probsets).  Now, I would like to separate the
> number of significant probesets for strain A and strain B from strain effect
> (I would like to do the same for treatment effect).  How can I do this? Any
> help is very much appreciated.

I am not sure exactly what you are asking for here, but you are probably 
better off using the limma package. After fitting the model you can then 
use e.g., decideTests() and vennDiagram() to get the number of probesets 
significant in contrasts of interest.

There is an excellent User's Guide that accompanies this package that 
will answer about 90% of the questions you may have.

HTH,

Jim


> 
> Thanks in Advance.
> Regards,
>  Sharon
> 
> 	[[alternative HTML version deleted]]
> 
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/bioconductor


-- 
James W. MacDonald
University of Michigan
Affymetrix and cDNA Microarray Core
1500 E Medical Center Drive
Ann Arbor MI 48109
734-647-5623



**********************************************************
Electronic Mail is not secure, may not be read every day, and should not be used for urgent or sensitive issues.