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



      Love Cricket? Check out live scores, photos, video highlights and more. Click here http://cricket.yahoo.com
	[[alternative HTML version deleted]]

