R/meanNet.R

In ECGen/coNet: Model and analyze co-occurrence based networks

#' Calculate the mean of a set of networks.
#' 
#' Given a set of networks in matrix format, this function will
#' calculate the cell-wise (Hadamard) mean.
#' 
#' @param x A list of networks to be used for calculating distances.
#' @param zero.na LOGICAL: Should NA values be treated as zeros?
#' @return Returns a matrix of means.
#' @author Matthew K. Lau
#' @seealso \code{\link{coNet}} %% ~~objects to See Also as \code{\link{help}}, ~~~
#' @examples
#' 
#' net.l <- lapply(1:10, function(x) matrix(runif(100), nrow = 10))
#' meanNet(net.l)
#' 
#' @export meanNet
meanNet <-function(x, zero.na = TRUE){
    x <- Reduce("+", x) / sum(unlist(x))
    if (zero.na){x[is.na(x)] <- 0}
    return(x)
}