# NOT RUN {
set.seed(1)
arr <- array(rnorm(12000),dim=c(40,25,30))
arr[25:33,12:20, 20:23] <- 10
# Getting events
out <- get_clusters_3d(arr, thres=0.985)
out <- tune_cpdbee_3D(arr, out$data[ ,1:3])
out$best
# }
Run the code above in your browser using DataLab