# NOT RUN {
## create 'swiss' roll dataset
X = aux.gensamples(dname="swiss")
## try different k values
out1 = est.clustering(X, kmin=5)
out2 = est.clustering(X, kmin=25)
out3 = est.clustering(X, kmin=50)
## print the results
line1 = paste0("* est.clustering : kmin=5 gives ",round(out1$estdim,2))
line2 = paste0("* est.clustering : kmin=25 gives ",round(out2$estdim,2))
line3 = paste0("* est.clustering : kmin=50 gives ",round(out3$estdim,2))
cat(paste0(line1,"\n",line2,"\n",line3))
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab