MuChPoint | R Documentation |
Produce a block-wise estimation of a symmetric matrix.
MuChPoint(Y, Lmax = nrow(Y)/2, N = NULL, cores = 1, verbose = TRUE)
Y |
symmetric matrix of observations. |
Lmax |
a positive integer less than number of columns (and number of rows).
By default, |
N |
a positive integer vector less than number of columns (and number of rows).
N is used when the break-points are known.
By default, |
cores |
a positive integer giving the number of cores used. If you use windows, the parallelization is impossible. By default, 1. |
verbose |
logical. To display the progression bars. By default TRUE. |
Article: BRAULT V., OUADAH S., SANSONNET L. and LEVY-LEDUC C. Nonparametric homogeneity tests and multiple change-point estimation for analyzing large Hi-C data matrices. Journal of Multivariate Analysis, 2017
require(MuChPoint) mu=c(rep(c(rep(1,25),rep(0,25)),3))%*%t(rep(c(rep(0,25),rep(1,25)),3)) Y=matrix(rnorm(150^2,0,5),150)+mu+t(mu) Y=as.matrix(Matrix::forceSymmetric(Y)) res=MuChPoint(Y) plot(res,Y,L=5,shiny=FALSE) plot(res,Y,L=1:10,shiny=FALSE,ask=FALSE)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.