lbmle
estimate the intrinsic dimension of a given dataset.
lbmle(x = NULL, k1 = NULL, k2 = NULL, BC = TRUE, DM = FALSE, p = NULL)
Estimated global intrinsic dimension.
data matrix or distance matrix given by as.matrix(dist(x)).
first k-NN parameter.
second k-NN parameter.
whether bias is corrected or not. logical.
whether 'x'
is distance matrix or not. logical.
ambient dimension used for automatically define 'k1'
and 'k2'
.
Hideitsu Hino hideitsu.hino@gmail.com
The likelihood of the rate parameter of the Poisson process, which characterize the behaviour of the distance from a point to another point in the given dataset, is considered, and the maximum likelihood estimator (MLE) for the intrinsic dimension is derived. The original method proposed by Levina and Bickel contains a known bias, and it is corrected by Mackay and Ghahramani. This function implements both, with the default the bias corrected estimate.
E. Levina and P. J. Bickel. Maximum likelihood estimation of intrinsic dimension. Advances in Neural Information Processing Systems 17, 2005.
D. MacKay and Z. Ghahramani. http://www.inference.org.uk/mackay/dimension/
x <- gendata(DataName='SwissRoll',n=300)
estmle <- lbmle(x=x,k1=3,k2=5)
print(estmle)
Run the code above in your browser using DataLab