lbmle {ider}

Maximum Likelihood Estimation of Intrinsic Dimension.

Description

lbmle estimate the intrinsic dimension of a given dataset.

Usage

lbmle(x = NULL, k1 = NULL, k2 = NULL, BC = TRUE, DM = FALSE, p = NULL)

Arguments

x

data matrix or distance matrix given by as.matrix(dist(x)).

k1

first k-NN parameter.

k2

second k-NN parameter.

BC

whether bias is corrected or not. logical.

DM

whether 'x' is distance matrix or not. logical.

p

ambient dimension used for automatically define 'k1' and 'k2'.

Details

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.

Value

Estimated global intrinsic dimension.

Author(s)

Hideitsu Hino [e-mail]

References

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/

Examples

x <- gendata(DataName='SwissRoll',n=300)
estmle <- lbmle(x=x,k1=3,k2=5)
print(estmle)