Thank you,
Teja
Code:
countsTable <- read.delim("countstable.txt", header=TRUE, stringsAsFactors=TRUE) rownames( countsTable ) <- countsTable$gene countsTable <- countsTable[ , -1 ] conds <- factor(c("tumor","tumor","tumor","tumor","tumor","tumor","tumor","tumor")) cds <- newCountDataSet( countsTable, conds ) cds <- estimateSizeFactors( cds ) cds <- counts( cds, normalized=TRUE ) cdsBlind <- estimateDispersions( cds, method="blind" ) vsd <- getVarianceStabilizedData( cdsBlind ) select <- order(res$pval)[1:250] colors <- colorRampPalette(c("white","darkblue"))(100) heatmap( vsd[select,], col = colors, scale = "none")
Leave a comment: