library(pheatmap)
library(argparse)
parser = ArgumentParser()
parser$add_argument("--fpkm", help = "fpkm file with sep='\t', genes as rows and samples as cols; gene_id in 1st col", required = TRUE)
parser$add_argument("--type", help = "the scale type", default = 'row')
parser$add_argument("--show", help = "show rownames", default = 'F')
parser$add_argument("--crow", help = "cluster row", default = 'TRUE')
parser$add_argument("--ccol", help = "cluster column", default = 'TRUE')
args <- parser$parse_args()
data = args$fpkm
type = args$type
## 不能对绘图参数传参？
#show = args$show
#crow = args$crow
#ccol = args$ccol
data <- read.table(data, header = T, sep = "\t", row.names = 1, check.names = FALSE)
## 对数转换
data <- log10(data + 1)
data
## 去除全部是0的行
data <- data[which(rowSums(data) > 0),]
pheatmap(data, scale = type, col = colorRampPalette(c("green4","black","red4"))(75), show_rownames = T, cluster_row = F, cluster_cols = F, filename = "heatmap.pdf")
pheatmap(data, scale = type, col = colorRampPalette(c("green4","black","red4"))(75), show_rownames = T, cluster_row = F, cluster_cols = F, filename = "heatmap.png")
R15
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着，像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上，一...
- 文/苍兰香墨 我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
推荐阅读更多精彩内容
- 手机行业步入4月份，OPPO依旧热得烫手——从3月31日晚上OPPO的超视野异形屏旗舰手机R15发布，到4月1日0...
- NDK r13升级到 r15编译问题，报以下错误： 可以看到提示 ECmedia 库依赖于一些未定义的module...
- 一：问题 今天在 OPPO R15 上调试代码，安装后发现无法启动。报错如下 报错的重点在 这里Caused b...
- 好事多磨！上个月OPPO放出有关ColorOS 6的消息震动了手机圈，也让众多手机爱好者对这个以无边界美学理念加身...
