The hardware and bandwidth for this mirror is donated by dogado GmbH, the Webhosting and Full Service-Cloud Provider. Check out our Wordpress Tutorial.
If you wish to report a bug, or if you are interested in having us mirror your free-software or open-source project, please feel free to contact us at mirror[@]dogado.de.
HDCD contains efficient implementations of several multiple change-point detection algorithms, including Efficient Sparsity Adaptive Change-point estimator (ESAC) and Informative sparse projection for estimating change-points (Inspect).
You can install the development version of HDCD from GitHub with:
# install.packages("devtools")
devtools::install_github("peraugustmoen/HDCD")This is a basic example which shows you how to run ESAC:
library(HDCD)
n = 50
p = 50
set.seed(100)
# Generating data
X = matrix(rnorm(n*p), ncol = n, nrow=p)
# Adding a single sparse change-point (at location \eta = 25):
X[1:5, 26:n] = X[1:5, 26:n] +2
# Vanilla ESAC:
res = ESAC(X)
res$changepoints
#> [1] 25
# Manually setting leading constants for \lambda(t) and \gamma(t)
res = ESAC(X,
threshold_d = 2, threshold_s = 2, #leading constants for \lambda(t)
threshold_d_test = 2, threshold_s_test = 2 #leading constants for \gamma(t)
)
res$changepoints #estimated change-point locations
#> [1] 25
# Empirical choice of thresholds:
res = ESAC(X, empirical = TRUE, N = 100, tol = 1/100)
res$changepoints
#> [1] 25
# Manual empirical choice of thresholds (equivalent to the above)
thresholds_emp = ESAC_calibrate(n,p, N=100, tol=1/100)
res = ESAC(X, thresholds_test = thresholds_emp[[1]])
res$changepoints
#> [1] 25These binaries (installable software) and packages are in development.
They may not be fully stable and should be used with caution. We make no claims about them.
Health stats visible at Monitor.