Deleted:Creating Heatmaps Of Diffrentially Expressed Genes In Angilent Single-Color Array Analysis
0
0
Entering edit mode
3.0 years ago
Aki ▴ 20

I am using limma to analyze Agilent's single color arrays to find differentially expressed genes, using the following procedure.

targets <- limma::readTargets("targets.txt")
project <- limma::read.maimages(targets,source="agilent", green.only=TRUE)
project.bgc <- backgroundCorrect(project, method="normexp", offset=50)
project.NormData <-normalizeBetweenArrays(project.bgc,method="quantile")
design <- paste(targets$Target, sep="")
design <- factor(design)
comparisonmodel <- model.matrix(~0+design)
colnames(comparisonmodel) <- levels(design)
project.fit <- lmFit(project.NormData, comparisonmodel)
project.fit.eBayes <- eBayes(project.fit)
CaseControl <- makeContrasts(TR-Control, levels=comparisonmodel)
CaseControl.fitmodel <- contrasts.fit(project.fit.eBayes, CaseControl)
CaseControl.fitmodel.eBayes <- eBayes(CaseControl.fitmodel)

The problem is that I don't know how to make a heatmap.

Does anyone tell a way to do this?

heatmap limma • 323 views
ADD COMMENT
This thread is not open. No new answers may be added
Traffic: 865 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6