I am new to microarray data analysis. I am analyzing microarray experiment using PROC MIXED in SAS and use genotype, replication, time point, genotypetimepoint as fixed factors and replicationtimepoint as random factor. I want to find differentially expressed genes for each genotype*timepoint combination. Do I look at lsmeans file and associated p-values or diffs file and associated p-values.
Thanks, Priyanka
•
written
by
Priyanka Surana •
3 answers
Maybe you can try Anovarray. AnovArray is a set of SAS macros for the Analysis of Variance of Gene expression. It can be applied to the quantification of biological and technological variation and to the identification of differentially expressed genes in two or more conditions.http://www.bioinformatics.org/project/?group_id=465
Thank you everyone. They were very helpful comments and answers. I have a good starting point to go from. :) Priyanka
Log in to answer this question.
-
Extracting the interaction term in bulk RNA-Seq data
written by james.zhang20 •
-
Estimating power in RNA-seq differential expression
written by ams686 •
-
Mixed model analysis using lme4
written by b.g.tamang •
-
linear mixed model for gene expression data in R
written by evafinegan •
-
Why kinship is a random effect while population structure is a fixed effect in GWAS studies?
written by rimgubaev
-
genomic inflation factor and how to use it to estimate adjusted P-values in GWAS?
written by Ana
-
Requesting help with RNA-Seq differential gene analysis
written by sagnik •
-
limma 2 factor contrast matrix design
written by Ming Tommy Tang
-
Deseq Standard Comparison Or Multifactorial / Model Fitting?
written by gaelgarcia05
-
Measure Significance Of Genes Associated With Multiple Phenotypes
written by Khader Shameer
You might want to consider not to use SAS... It is of course possible to do microarray statistics in SAS but the mainstream approach really is to do it in R/Bioconductor and many complete packages exist to assist you in doing so. If you are new to the field, like you say you are, I would just switch to Bioconductor and learn that, even if you have prior knowledge of SAS.
Since your question is essentially about how to interpret a generalized mixed linear model and is not inherently tied to microarray data, you might consider asking it on stats.stackexchange.com, where the readership is mostly statisticians.