R/matrixtoGNAR.R

In nunesmatt/GNAR: Methods for Fitting Network Time Series Models

matrixtoGNAR <- function(input.mat){
  stopifnot(nrow(input.mat)==ncol(input.mat))
  if(min(input.mat) < 0){
    cat("WARNING: negative entries present in original matrix")
  }
  edges <- dist <- vector(mode="list", length=nrow(input.mat))

  for(ii in 1:nrow(input.mat)){
    nz <- which(input.mat[ii,] != 0)
    edges[[ii]] <- nz
    dist[[ii]] <- 1/input.mat[ii,nz]
  }
  out <- list(edges=edges, dist=dist)
  class(out) <- "GNARnet"
  return(out)
}