您好老兄,我在用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()
代码如上所示,只能出来这个折线图,热图出不来,保存下来的也是这个折线图。请问这是什么原因?

您好老兄,我在用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()
代码如上所示,只能出来这个折线图,热图出不来,保存下来的也是这个折线图。请问这是什么原因?