R: nel2igraph и PN.amalgamation - неправильно созданный играф - PullRequest
Новая
       10

R: nel2igraph и PN.amalgamation - неправильно созданный играф

0 голосов
/ 06 июня 2018

Я столкнулся с проблемой с пакетом shp2graph.Я хочу использовать функцию PN.amalgamation, которая прекрасно работает (см. Ниже).После этого я хотел бы создать объект igraph.Здесь код не может этого сделать.Я могу прекрасно создавать объекты igraph для каждого неамальгамированного объекта shp2graph.

Вот мой пример кода, который в основном является копией вставки из описания пакета shp2graph: 

library(igraph)
library(shp2graph)

data(ORN)
rtNEL<-readshpnw(ORN.nt, ELComputed=TRUE)
res.sl<-SL.extraction(rtNEL[[2]],rtNEL[[3]])
res.me<-ME.simplification(res.sl[[1]],res.sl[[2]],DegreeL=res.sl[[4]]) 
res.pn<-PN.amalgamation(res.me[[1]],res.me[[2]],DegreeL=res.me[[4]])
ptcoords<-Nodes.coordinates(res.pn[[1]])
plot(ORN.nt)
points(ptcoords, col="green")
plot(ORN.nt)
points(Nodes.coordinates(rtNEL[[2]]), col="red")

# igraph created from amalgamation is wrong
test <- nel2igraph(nodelist= res.pn[[1]], edgelist=res.pn[[2]], Directed = TRUE)
plot(test,vertex.size=1,edge.width=1,edge.arrow.size=0,vertex.label=NA)

# res.me is one step before amalgamation
test <- nel2igraph(nodelist= res.me[[1]], edgelist=res.me[[2]], Directed = TRUE)
plot(test,vertex.size=1,edge.width=1,edge.arrow.size=0,vertex.label=NA)

Любая помощь приветствуется.

1 Ответ

0 голосов
/ 06 июня 2018

Я обнаружил, что ошибка как-то связана с взаимодействием с пакетом igraph.Проблема в том, что метки узлов, создаваемые PN.amalgamation, больше не являются непрерывными;некоторые отсутствуют, так как мы их удалили.Тем не менее, igraph почему-то все еще пытается их создать и выдает следующее предупреждение:

Для тех, у кого возникла такая же проблема, это обходной путь, который заново устанавливает метки.

СоздатьВаша собственная nel2igraph функция: 

nel2igraph_corr <- function (nodelist, edgelist, weight = NULL, eadf = NULL, Directed = FALSE) 
{


  nodes <- nodelist[, 1]
  Ne <- length(edgelist[, 1])
  Nn <- length(nodes)

    for (i in 1:Nn) {
      kk <- nodelist[i,][[1]]
      edgelist[which(edgelist[,c(2)]==kk),2] <- i
      edgelist[which(edgelist[,c(3)]==kk),3] <- i

      nodelist[i,][[1]] <- i
  }

    if (!is.null(weight)) {
    if (length(weight) != Ne && is.numeric(weight)) 
      stop("Please give right edge weight, which must be numeric and the same length as edges elment")
  }
  if (!is.null(eadf)) {
    if (length(eadf[, 1]) != Ne) 
      stop("The eadf must be numeric and the same length as edges elment")
  }



  gr <- graph.edgelist(unique(edgelist[, c(2, 3)]), directed = T)
  gr <- set.vertex.attribute(gr, "x", V(gr), Nodes.coordinates(nodelist)[,1])
  gr <- set.vertex.attribute(gr, "y", V(gr), Nodes.coordinates(nodelist)[, 
                                                                         2])
  gr.es <- E(gr)
  if (!is.null(weight)) 
    gr <- set.edge.attribute(gr, "weight", gr.es, weight)
  if (!is.null(eadf)) {
    eanms <- colnames(eadf)
    n <- length(eanms)
    for (i in 1:n) gr <- set.edge.attribute(gr, eanms[i], 
                                            gr.es, eadf[, i])
  }
  gr
}
...