Skip to content

绘制热图代码输出折线图 #158

@Bio-info-ai

Description

@Bio-info-ai

你好,我用绘制热图的代码,折线图,两者一起的代码,均出现折线图:

library(ClusterGVis)
library(ComplexHeatmap) # 用于扩展热图功能

读取你的水稻表达矩阵,假设名为 'rice_exps.csv'

格式示例:第一列为基因ID,后续列为5个时期的表达量

rice_exps <- read.csv("RD15_Sy1_tpmave.csv", row.names = 1, header = TRUE, check.names = FALSE)

3. (可选但推荐) 数据标准化

对于转录组TPM/FPKM数据,通常进行Z-score标准化以使不同基因间表达量可比

rice_expr_zscore <- scale(rice_exps)

绘制不同聚类数量下的组内平方和变化图

getClusters(obj = rice_expr_zscore)

设置随机种子以保证结果可重复

set.seed(123)
clust_res <- clusterData(obj = rice_expr_zscore,
clusterMethod="mfuzz",
clusterNum=8) # 将6替换为你确定的k值

绘制热图

visCluster(object = clust_res, plot.type = "heatmap")

绘制表达趋势折线图

visCluster(object = clust_res, plot.type = "line")

尝试绘制组合图

visCluster(object = clust_res, plot.type = "both")

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions