Я учусь, как использовать R и для упражнения я использую экспериментальную ExpresionSet из Bioconductor: http://bioconductor.org/packages/release/data/experiment/html/leukemiasEset.htmlНи цвета, ни Key Main Title появляется (правильно) на карте тепла (тепло map.2 в г)
Так просто: Моя карта тепла не отображает цветную клавишу, и главный заголовок отображается на левом сайте. Я отказался от проблемы с цветовым ключом после поиска в Интернете, как ее решить, и не нашел решения. Но затем я добавил заголовок и так как он не появляется должным образом либо, я подумал, что, может быть, и то и другое в связи с тем же (неизвестное для меня) причине ...
Вот мой код:
difexp <- exprs(leukemiasEset)[c(which(adjPval < 0.01)),c(1:12, 25:36)]
heatmap.2(difexp,
trace = "none",
cexCol = 0.6,
ColSideColors = as.character(as.numeric(factor(leukemiasEset$LeukemiaType[c(1:12, 25:36)]))),
main = "Differentially expressed genes in ALL and CLL samples",
cex.main = 1.5,
key = TRUE,
keysize = 1.5,
density.info = "histogram")
Где:
pvalues <- c()
for(i in 1:nrow(exprs(leukemiasEset))) {
R <- t.test(exprs(leukemiasEset)[i, leukemiasEset$LeukemiaType == "ALL"],
exprs(leukemiasEset)[i, leukemiasEset$LeukemiaType == "CLL"],
var.equal = TRUE)
pvalues <- c(pvalues, R$p.value)
}
adjPval <- p.adjust(pvalues, method = "fdr")
И вот как это выглядит:
Так как я новичок, я думаю, что это, вероятно, «легкий peasy» для экспертов ... большое вам спасибо!