首页
学习
活动
专区
圈层
工具
发布
社区首页 >问答首页 >使用heatmaply将行名添加到热图

使用heatmaply将行名添加到热图
EN

Stack Overflow用户
提问于 2017-11-27 10:46:05
回答 1查看 674关注 0票数 0

我在Rstudio中绘制了一些基因表达的热图,但是行名(应该是样本或基因名称)被数字取代。我已经检查过some posts了,但是没有机会。抱歉,如果这是一个简单的问题或重复的问题,我是新来R的。

这是我的R脚本:

代码语言:javascript
复制
###############################################################

library(heatmaply)

filename <- "TFzebra-TMM-name.matrix"

head(filename)

my_data <- read.table(filename, sep='\t', quote='', stringsAsFactors=FALSE, header=TRUE)
head(my_data)
dim(my_data)
nrow(my_data)
ncol(my_data)
row.names(my_data) <- my_data$samples
head (my_data)
my_data <- my_data[, -1]
head(my_data)
data <- my_data/rowSums(my_data)
my_matrix <- as.matrix(data) 

heatmaply(data, xlab = "Features", ylab = "Transcription Factors", 
          scale = "column",
          main = "Data transformation using 'scale'",
          margins = c(60,100,40,20))

注意:

代码语言:javascript
复制
> head (my_data)
    samples   J1 J2 J3    H1    H2    H3
1   zf-C2H2 0.00  0  0 0.507 0.232 0.540
2     baz2a 0.00  0  0 0.355 0.342 0.132
3    baz2ab 0.00  0  0 0.091 0.070 0.163
4    kdm5ba 0.00  0  0 1.835 0.856 1.527
5      dbpa 0.00  0  0 5.344 5.355 0.000
6 LOC555983 0.02  0  0 0.294 0.634 0.835

NOTE2:

代码语言:javascript
复制
> my_data <- my_data[, -1]
> head(my_data)
    J1 J2 J3    H1    H2    H3
1 0.00  0  0 0.507 0.232 0.540
2 0.00  0  0 0.355 0.342 0.132
3 0.00  0  0 0.091 0.070 0.163
4 0.00  0  0 1.835 0.856 1.527
5 0.00  0  0 5.344 5.355 0.000
6 0.02  0  0 0.294 0.634 0.835
EN

回答 1

Stack Overflow用户

发布于 2017-11-28 16:29:41

在删除第一列之前,您应该将第一列添加到row.names属性(使用rownames函数,是的,我知道,它不太一致)。这就是:

代码语言:javascript
复制
rownames(my_data) <- my_data[, 1]
my_data <- my_data[, -1]
library(heatmaply)
heatmaply(my_data) # should now work

对于未来-请记住,既然你是在要求别人从他们的时间(免费)给你帮助,你必须首先向他们表明你已经尽了最大努力最大限度地利用他们的时间。因此,在你发帖之前,请确保你的问题包括了所有需要的信息,这样人们就会回答它(特别是让它包括一个自包含的例子)。请阅读以下两个关于如何编写好问题的指南,这将有助于确保您在未来获得所需的帮助:http://stackoverflow.com/help/how-to-ask

票数 1
EN
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接:

https://stackoverflow.com/questions/47503276

复制
相关文章

相似问题

领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档