Skip to content

visCluster画不出热图 #169

@PanSX-Dr

Description

@PanSX-Dr

您好老兄,我在用clusterGVIS画热图的时候,前边所有的代码都是照搬的教程,但是只能出折线图,最后出不来人途。
pbmc.markers <- all_markers %>%
dplyr::group_by(cluster) %>%
dplyr::top_n(n = 20, wt = avg_log2FC)
st.data <- prepareDataFromscRNA(object = hcc_sce,
diffData = pbmc.markers,
showAverage = TRUE)
st.data

enrich for clusters

enrich <- enrichCluster(object = st.data,
OrgDb = org.Hs.eg.db,
type = "KEGG",
organism = "hsa",
pvalueCutoff = 0.05,
topn = 5,
seed = 5201314)
head(enrich)
markGenes = unique(pbmc.markers$gene)[sample(1:length(unique(pbmc.markers$gene)),40,
replace = F)]

line plot

visCluster(object = st.data,
plot.type = "line")

heatmap plot

pdf('sc1.pdf',height = 10,width = 6,onefile = F)
visCluster(object = st.data,
plot.type = "heatmap",
column_names_rot = 45,
markGenes = markGenes,
cluster.order = c(1:9))
dev.off()

add bar plot

pdf('sc2.pdf',height = 10,width = 14,onefile = F)
visCluster(object = st.data,
plot.type = "both",
column_names_rot = 45,
show_row_dend = F,
markGenes = markGenes,
markGenes.side = "left",
annoTerm.data = enrich,
line.side = "left",
cluster.order = c(1:9),
go.col = rep(jjAnno::useMyCol("stallion",n = 9),each = 5),
add.bar = T)
dev.off()
代码如上所示,只能出来这个折线图,热图出不来,保存下来的也是这个折线图。请问这是什么原因?

Image

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