scMEGA
scMEGA copied to clipboard
NetCentPlot error
Hi! thanks for the jobs I've got the object to build the network on
head(df.grn2) tf gene weights 77 ATF4 GABARAPL1 0.8221601 85 ATF4 GRIK3 0.8021985 100 ATF4 KIT 0.8130254 102 ATF4 KRAS 0.8125377 166 ATF4 SEC22B 0.8518962 223 ATF6 ACAP3 0.8103851 V(netobj)$type <- ifelse(V(V(netobj)$type <- ifelse(V(netobj)$name %in% dfgrn$tf,"TF/Gene","Gene"))$name %in% df.grn2$tf,"TF/Gene","Gene") netobj <- graph_from_data_frame(df.grn2,directed = TRUE)
Looking at the object, it looks normal
netobj IGRAPH 54cba5f DN-B 355 788 --
- attr: name (v/c), type (v/c), weights (e/n)
- edges from 54cba5f (vertex names): [1] ATF4 ->GABARAPL1 ATF4 ->GRIK3 ATF4 ->KIT ATF4 ->KRAS [5] ATF4 ->SEC22B ATF6 ->ACAP3 BACH2->ALX1 BACH2->ANKH [9] BACH2->CDKN1C BACH2->CRISPLD1 BACH2->EBF2 BACH2->ETV6 [13] BACH2->GLIPR1 BACH2->HECTD2 BACH2->IGFBP7 BACH2->ITGB8
- ... omitted several edges However, when I run
p <- NetCentPlot(netobj, "RUNX1") Error in layout_with_focus(graph, v = focus, weights = weights, iter = niter, : g must be a connected graph.
I checked netobj
is_connected(netobj, mode = "weak") [1] FALSE
I don't know how do I get the final NetCentPlot, Can you offer any help? thanks a lot!