Question: Heatmap of Agilent two colors Dye Swap experiment
0
gravatar for antgomo
2.9 years ago by
antgomo20
Spain
antgomo20 wrote:

Hi, I have an Agilent Dye Swap experiment, two colors, where I have my samples Resistant and Sensible

Distributed like this (this is my targets.txt file)

SampleNumber    FileName    Cy3    Cy5
1    A.txt    S    R
2    B.txt    R    S
3    C.txt    S    R
4    D.txt    R    S
5    E.txt    S    R
6    F.txt    R    S
7    B1.txt    R    S
8    C1.txt    S    R
9    D1.txt    R    S
10    E1.txt    S    R
11    F1.txt    R    S

I analyzed the data using the nex code

library(limma)

targets <- readTargets("targets_all.txt")
RG <- read.maimages(targets,source="agilent.median")
RG <- backgroundCorrect(RG, method="normexp", offset=16)
MA <- normalizeWithinArrays(RG, method="loess")
MA.avg <- avereps(MA, ID=MA$genes$ProbeName)
   
design <- modelMatrix(targets, ref="S")
fit2 <- eBayes(fit)
output <- topTable(fit2, adjust="fdr", coef="R", number=nrow(fit2), genelist=MA.avg$genes)
genes_diff_expr.005<- subset(output,output$adj.P.Val < 0.05)
pdf("heatmap_diff_exprs_R_S.pdf")
heatmap.2(MA$M[rownames(MA$M) %in% genes_diff_expr.005$ProbeName,], trace="none", col=greenred(10),cexRow=0.2,labRow=NA,cexCol=0.7)
dev.off()

Now, I want to see how my samples cluster for this signature. I've been doing so many heatmaps in years, mainly using Affymetrix or Agilent one color experiments, but I don't know how to use to plot all the samples distributed in the two channels (R). Using MA$M, I got the log2ratios, the point is that the name of my columns are as in targets A,B,C and so on, and I want R or S, but distributed in two channels I don't know how to specify which is sensible or resistant

any ideas?

Thanks in advance

ADD COMMENTlink modified 2.6 years ago by Biostar ♦♦ 20 • written 2.9 years ago by antgomo20
0
gravatar for Manvendra Singh
2.9 years ago by
Manvendra Singh1.9k
Berlin, Germany
Manvendra Singh1.9k wrote:

draw a boxplot to ensure that "MA" is normalized

make your rownames as gene_ID or probe ID of the data "MA"

############################ calculating top differentially expressed genes PS: you can choose the number in percentage,

#depends what fraction of genes you want to consider in R, please have this libraries if you don't have

## analysis starts from here

library(genefilter)
percentage<-c(0.990)
sds<-rowSds(MA)
sel<-(sds>quantile(sds,percentage))
set<-MA[sel, ]
distmeth<-c("euclidian")
D<-dist(t(set), method=distmeth)
treemeth<-c("average")
hc<-hclust(D, method=treemeth)

############### clustering
plot(hc)

######## for heatmap 
library(gplots)
library(RColorBrewer)
hr <- hclust(as.dist(1-cor(t(set), method="pearson")), method="centroid"); 
hc <- hclust(as.dist(1-cor(set, method="spearman")), method="centroid")  

# Cuts the tree and creates color vector for clusters.
mycl <- cutree(hr,k=4, h=max(hr$height)/1.5); 
mycolhc <- rainbow(length(unique(mycl)), start=0.1, end=0.9); mycolhc <- mycolhc[as.vector(mycl)] ; myheatcol <- bluered(75)
heatmap.2(as.matrix(set), Rowv=as.dendrogram(hr), Colv=as.dendrogram(hc), col=myheatcol, scale="row", density.info="none", trace="none", RowSideColors=mycolhc,cexRow=1.5, cexCol=1.5, keysize=1,margins=c(12,10))
ADD COMMENTlink modified 2.9 years ago by Ram12k • written 2.9 years ago by Manvendra Singh1.9k

OK. Thanks Manu, but this is not the point. I got my diff expressed genes, and I know to build the heatmap. The problem is that I have my R and S samples distributed in different channels and what I want to see is my heatmap of R and S, maybe I was no clear at this point, not see how A,B and C (arrays) are clustering

ADD REPLYlink written 2.9 years ago by antgomo20

How did you solve your problem? Thanks!

ADD REPLYlink written 6 months ago by thyagoleal0
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 684 users visited in the last hour