R/topology.R

In mccallm/crgnet: Gene network modeling via TopNet reveals functional dependencies between diverse tumor-critical mediator genes

topology <- function(fits){
  ngene <- nrow(fits[[1]]$parents)
  nexp <- length(fits[[1]]$trajectories)
  grph <- matrix(0, nrow=ngene, ncol=nexp)
  for(k in 1:length(fits)){
    gtmp <- t(fits[[k]]$parents + 1)
    for(j in 1:ncol(gtmp)){
      grph[j, gtmp[,j]] <- grph[j, gtmp[,j]] + 1
    }
  }
  cgrph = grph / length(fits)
  return(cgrph)
}