R/mat_d_col.R

In jeitziner/TTMap: Two-Tier Mapper: a clustering tool based on topological data analysis

mat_d_col <-
function(mat){
    for(i in seq_len(dim(mat)[1])){
        s <- grep(1, mat[,i])
            for(j in s){
                if(j == i) mat <- mat
                else{
                    w <- grep(1, mat[j,])
                    mat[i, w] <- rep(1, length(w))
                    mat[j, ] <- rep(0, dim(mat)[1])
                }
            }
    }
    return(mat)}