Function to remove centromere columns and rows from a full Hi-C contact matrix

remove_centromere(mat)

Arguments

mat	A full Hi-C matrix

Value

A list of (1) the column/row numbers of the centromere and (2) the Hi-c matrix with the centromere removed

Examples

m <- matrix(rpois(100, 5), 10, 10)
m[5,] <- 0
m[,5] <- 0
remove_centromere(m)
#> $centromere_loc
#> [1] 5
#> 
#> $mat
#>       [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9]
#>  [1,]    6    3    6    2    7    4    2    5    3
#>  [2,]    5    4    2    8    5    9    1    5    5
#>  [3,]    3    5    7    5    6    6    7    2    5
#>  [4,]    3    0    6    4    3    4    6    5    5
#>  [5,]    3    3    6    4    4    2    7    4    6
#>  [6,]    2    9    4    5    9    4    1    4    4
#>  [7,]    5    5    7   10    6    6    7    0    5
#>  [8,]    4    6    6    4    4    5    4    4    6
#>  [9,]    5    3    4    6    4    3    5    5    9
#>