【发布时间】:2018-08-24 09:40:31
【问题描述】:
我的图表是
nodes = data.frame(GENES = c("GENE_A", "GENE_B", "GENE_C", "GENE_D", "GENE_E","GENE_F", "GENE_G", "GENE_H", "GENE_I", "GENE_J"),TYPES = c("LIGAND", "RECEPTOR", "PROTEIN", "RECEPTOR", "LIGAND", "PROTEIN", "LIGAND","RECEPTOR", "PROTEIN", "LIGAND"))
edges = data.frame(FROM=c("GENE_A","GENE_B","GENE_C","GENE_C","GENE_D","GENE_E","GENE_F","GENE_A","GENE_H","GENE_H","GENE_I","GENE_J","GENE_J"),TO=c("GENE_B","GENE_D","GENE_I","GENE_G","GENE_F","GENE_I","GENE_B","GENE_G","GENE_E","GENE_F","GENE_C","GENE_D","GENE_J"),STATUS=c('ACTIVE','INACTIVE','UNKNOWN','INACTIVE','INACTIVE','UNKNOWN','ACTIVE', 'ACTIVE','ACTIVE','INACTIVE','UNKNOWN',"ACTIVE","ACTIVE"))
g <- graph_from_data_frame(edges, directed = TRUE, vertices = nodes)
我想子图化并简化下一个节点d = data.frame(NODES = "GENE_B", "GENE_D", "GENE_H"),但每次我simplify我都会丢失边缘的所有属性。
作为输出,我想获得d 中具有直接边 (order = 1) 的节点图。我试过了
subg <- induced.subgraph(graph = g, vids = unlist(neighborhood(graph = g, order = 1, nodes = d, mode = 'out'))) 但简化后我失去了边缘属性,无法追踪E(g)$STATUS
【问题讨论】: