if (FALSE) {
## global intrinsic dimension estimate
x <- gendata(DataName='SwissRoll',n=300)
x <- gendata(DataName='SwissRoll',n=300,p=3,q=2)
estcorint <- corint(x=x,k1=5,k2=10)
print(estcorint)
estmle <- lbmle(x=x,k1=3,k2=5) ## estimation by mle
print(estmle)
estnii <- nni(x=x) ## estimation by nearest neighbor information
print(estnni)
estconvU <- convU(x=x) ## estimation by convergence property of U-stats
print(estconvU)
estpackG <- pack(x=x,greedy=TRUE) ## estimation by the packing number method with greedy algorithm
print(estpackG)
estpackC <- pack(x=x,greedy=FALSE) ## estimation by the packing number method by clutering
print(estpackC)
## local intrinsic dimension estimate
tmp <- gendata(DataName='ldbl',n=300)
x <- tmp$x
estmada <- mada(x=x,local=TRUE)
head(estmada) ## estimated local intrinsic dimensions by mada
head(tmp$tDim) ## true local intrinsic dimensions
estside <- side(x=x,local=TRUE)
head(estside) ## estimated local intrinsic dimensions by side
}
Run the code above in your browser using DataLab