Entering edit mode
Diya v
▴
20
@diya-v-3664
Last seen 9.5 years ago
Hi,
I have a question.I have 2control samples(v7,v8) and 2 treatment samples(g1,g4).I have the data generated from aroma.affymetrix and I have applied log2 on it and load the data for doing statistical analysis by limma.Is my design matrix correct?
dat.m<-read.csv("datam.csv",skip=0)
genes<-rownames(dat.m)
groups<-c(C,C,T,T)
groups<-as.factor(groups)
design<-model.matrix(~groups)
design
fit<-lmFit(dat.m, design)
fit<-eBayes(fit)
toptable(fit, coef=2)
Thanks, Diya
please visit my site Sarang188