Hello to all,

I have started agilent microarray data analysis (one color), I have used R bioconductor for analysis steps i used in are as follows now i am not able to get heatmap of differentially expressed genes . Please help me.

library(limma)

T <-readTargets("Targets.txt")

x <-read.maimages(targets, source="agilent",green.only=TRUE)
y <-backgroundCorrect(x, "normexp", offset=50)
y <-normalizeBetweenArrays(y, method="quantile")

y.ave <- avereps(y, ID=y$genes$ProbeName)

f <- factor(targets$Condition, levels = unique(targets$Condition))

design <- model.matrix(~0 + f)
colnames(design) <- levels(f)

fit <- lmFit(y.ave, design)

contrast.matrix <- makeContrasts("control-drought", "control-heat", "drought-control","control-combined","heat-drought", levels=design)

fit2 <- contrasts.fit(fit, contrast.matrix)
fit2 <- eBayes(fit2)

C.top<-as.numeric(rownames(topTable(fit2, coef=1,n=50)))

heatmap(y.ave$M[C.top,])
Error in heatmap(y.ave$M[C.top, ]) : 'x' must be a numeric matrix