Intrinsic dimension estimation for locally undersampled data

Identifying the minimal number of parameters needed to describe a dataset is a challenging problem known in the literature as intrinsic dimension estimation. All the existing intrinsic dimension estimators are not reliable whenever the dataset is locally undersampled, and this is at the core of the so called curse of dimensionality. Here we introduce a new intrinsic dimension estimator that leverages on simple properties of the tangent space of a manifold and extends the usual correlation integral estimator to alleviate the extreme undersampling problem. Based on this insight, we explore a multiscale generalization of the algorithm that is capable of (i) identifying multiple dimensionalities in a dataset, and (ii) providing accurate estimates of the intrinsic dimension of extremely curved manifolds. We test the method on manifolds generated from global transformations of high-contrast images, relevant for invariant object recognition and considered a challenge for state-of-the-art intrinsic dimension estimators.

Processing, analyzing and extracting information from high dimensional data is at the core of the modern research in machine learning and pattern recognition.One of the main challenges in this field is to decompose and compress, without losing information, the redundant representations of complex data that are produced across diverse scientific disciplines, including computer vision, signal processing and bioinformatics.Manifold learning and dimensional reduction [1][2][3] are the main techniques employed to perform this task.Several of these approaches work under the reasonable assumption that the points (or samples) of a dataset, represented as vectors of real numbers lying in a space of large embedding dimension D, actually belong to a manifold M, whose intrinsic dimension (ID) d is much lower than D. The problem of providing accurate estimates for this number has been recognized multiple times and in different contexts: in psychometry by Shepard [4,5], in computer science by Trunk [6], Fukunaga and Olsen [7], in physics by Grassberger, Procaccia [8,9], and Takens [10].More recently, IDE has been reconsidered with the advent of big data analysis, artificial intelligence and demanding molecular dynamics simulations, and several estimators to measure the intrinsic dimension have been proposed [11][12][13][14][15][16][17][18][19][20][21][22].

Standard algorithms for IDE and extreme locally undersampled regime
Algorithms for intrinsic dimension estimation can be roughly classified in two groups [12].Projective methods compute the eigenvalues of the D ×D covariance matrix C X of the data X, defined as (C X ) ij = 1/N N µ=1 x µ i x µ j , where x µ i is the i-th component of the µ-th sample vector of the dataset x µ (µ = 1, . . ., N ).The ID is then estimated by looking for jumps in the magnitude of the sorted eigenvalues of C X (see top left panel in Fig. 1).Principal component analysis (PCA) is the main representative of this class of algorithms.Both a global (gPCA) and a multiscale version (mPCA) of the algorithm are used [15,23].In the former one evaluates the covariance matrix on the whole dataset X, whereas in the latter one performs the spectral analysis on local subsets X(x 0 , r c ) of X, obtained by selecting one particular point x 0 and including in the local covariance matrix only those points that lie inside a cutoff radius r c , which is then varied.
The main limitation of the global PCA is that it can only detect the correct ID of linearly embedded manifolds (i.e.linear manifolds R d ∈ R D embedded trivially via rotations and translations), systematically overestimating the ID of curved/non-linearly embedded datasets.The mPCA could in principle fix this issue.However, PCA only works if the number of samples N d log d, otherwise being inconclusive (see top left panel of Fig. 1).This is a major drawback in the multiscale case, since to measure the correct ID the cutoff radius r c needs to be small enough, implying that the sampling of the manifold must be dense enough to guarantee that a sufficient number of samples lies inside the subsets X(x 0 , r c ).Another technical issue that makes mPCA difficult to employ for IDE is the fact that the amplitude of the "jump" in the magnitude of the sorted eigenvalues depends on the data, and the choice of a threshold size is somewhat arbitrary.
The other group of estimators belongs to the so-called geometric (or fractal) methods.Their common ancestor is the correlation dimension (CorrDim) introduced by Grassberger and Procaccia [9] to measure the fractal dimension of strange attractors in the context of dynamical systems.This estimator is based on the observation that the density of neighbours (also known as correlation integral in the literature) with a given cutoff distance r scales as ρ X (r) ∼ r d for r → 0 and therefore one can extract the ID by measuring the slope of the linear part of ρ as a function of r in log-log scale, since the relation d = lim r→0 log ρ X (r)/ log r holds (see bottom left panel of Fig. 1).
CorrDim is very effective for the estimation of low IDs (d 10), whereas it systematically underestimates in the case of manifolds with larger IDs.This drawback is well known in the literature [18] and is only partially mitigated by more recent and advanced generalizations of CorrDim based on k-nearest-neighbors distances [12].The reason why all these algorithms systematically fail for d 10 is due to a fundamental limitation of most geometric methods: indeed it is possible to prove [24] that the accurate estimation of the ID requires a number of samples N which grows exponentially in the intrinsic dimension d.As a consequence, one observes a systematic undersampling of the small radius region of the density of neighbors ρ X (r), as shown for the CorrDim estimator in the bottom left panel of Fig. 1.
Both mPCA and CorrDim, as well as their more recent generalizations such as DANCo [18], are based on the fundamental fact that, locally, samples in a dataset are effectively drawn uniformly from a d-dimensional disk {x x x ∈ R d : ||x|| ≤ 1} linearly embedded in R D .This is the informal way to state rigorous results based on the tangent space approximation to smooth manifolds and embeddings [18,25].On one side, local neighbourhoods of large ID datasets (d 10) need a number of points N exponential in the ID to be sufficiently sampled; on the other hand, the tangent space approximation needs dense sampling of small patches of the manifold to be used in practice.This incompatibility of requirements is the so called curse of dimensionality, and defines a theoretical limit for all multiscale projective and geometric ID estimators.
In order to try and break the curse of dimensionality, additional information about the probability distribution of the data must be assumed.The main ingredient of our spell is the assumption that data are locally isotropic.This suggests to consider the average correlation integral for hyperspheres {x x x ∈ R d : ||x|| = 1}.We will leverage on this observation to develop a novel geometric ID estimator for linearly embedded manifolds, which departs from the small radius limit of the This method works whenever the number of samples N in the dataset is sufficiently large (N d log d), and the jump occurs between the d-th and (d + 1)-th eigenvalues (green line in the plot corresponding to a dataset of N = 1000 samples drawn uniformly from a hypercube in dimension d = 100).In the opposite undersampled regime N < d, PCA is inconclusive (orange line in the plot, same dataset with N = 50 samples).(Bottom left) Geometric methods, such as the CorrDim or the more general k-nearest-neighbours estimators, are based on the scaling of the density of neighbours ρ(r/r 0 ) at a small cutoff distance r with respect to the average diameter r 0 of the dataset.In particular, ρ(r/r 0 ) ∼ (r/r 0 ) d for r/r 0 → 0 independently on the details of the datasets, so that a log-log linear fit estimates the ID as the slope of the fitting line.However, the small r regime is exponentially (in the ID, [24]) difficult to sample.This effect is at the origin of the so-called curse of dimensionality, and it induces a systematic underestimation of the ID.As N increases, though one is able to sample smaller r regions as shown in the plot (orange line N = 50, green line N = 1000 as above), convergence to the true dimension (d = 100 in the plot) is not possible.(Right) The density of neighbours ρ is more generally defined at any cutoff distance r.The r r 0 regime is the one used to compute CorrDim.In the remaining region r r 0 , the density ρ increases and eventually approaches one, indicating that the underlying proximity graph (see insets) gets more and more connected.We observe that this region is easily sampled at any fixed N (plot at N = 50), but the functional form of ρ in this informative regime is in principle dependent on the details of the dataset.), for d = (5,15,30).We are able to efficiently extract the correct ID even though this is a highly non-uniformly sampled dataset, whose ρ displays manifold-dependent features (in this case, step-like patterns).Moreover, we observe that, as we increase d, the density of neighbours of this dataset quickly converges to our functional form.It is worth noticing that the whole functional form ( 2 We observe that at N ∼ 100 we have an error of the order of 1% almost independently on the ID, and that ID estimation is possible also in the extreme undersampled N < d regime. (Bottom right) The FCI method estimates the correct ID even when the data are corrupted by noise.Here we consider a linearly embedded hypercube dataset H 40,60 and add on the top of that a 60-dimensional gaussian noise of standard deviation σ.We observe a sharp transition in the estimated ID between the regime in which the noise is a perturbation (σ 0.1 and dest = 40) and the regime in which the noise covers the signal (σ 0.2 and dest = 60).
density of neighbors ρ X (r) and considers this quantity at finite radius r.We show in the following that this method overcomes the extreme undersampling issue and displays a remarkable robustness to non-uniform sampling and noise.Based on the intriguing features reported above, we propose a multiscale generalization of the estimator, capable of providing the correct ID for datasets extracted from highly curved or multidimensional manifolds.

Full Correlation Integral (FCI) Estimator
The tangent space approximation suggests that a special role in the IDE problem is played by uniformly sampled d-dimensional disks linearly embedded.The average correlation integral for the boundary of this manifold, which is the (d − 1)-dimensional sphere of radius r s , can be analytically evaluated as (see Materials and methods) where 2 F 1 is the (2, 1)-hypergeometric function, Ω d is the d-dimensional solid angle and r = r/r s is the adimensional cutoff radius.We take (2) as the definition of the full correlation integral, to stress that we work away from the small radius limit employed for CorrDim.
It is worth noticing that the FCI has a sigmoidal shape which is steeper as the ID grows (see for instance the black lines in the top left panel of Fig. 2).This observation translates into a simple exact algorithm to determine the ID d of linearly embedded spherical datasets, by performing a non-linear regression of the empirical density of neighbours using the FCI in (2).More in general, this protocol is exact for linearly embedded Euclidean spaces sampled with a rotational invariant probability distribution, by projecting onto the unit sphere and adding one to the ID estimated on this new dataset.We summarize our FCI estimator in two steps: 1. compute the center of mass b of the empirical data as b = 1/N µ x µ and translate each datapoint by this quantity, so that the resulting dataset is centered at the origin.Then normalize each sample; 2. measure the empirical correlation integral of the dataset as a function of the radius r and perform a non-linear regression of this empirical density of neighbours using the FCI in (2) as the non-linear model and d and r s as the free parameters; as the normalization step artificially removes one degree of freedom, increase the estimated ID by one.
We discuss the technical details regarding the fitting protocol in the Methods.

Robustness of the FCI estimator
We now provide strong numerical evidence that the FCI estimator goes well beyond the exact results summarized above by testing it on multiple synthetic non-spherical datasets.First, we notice that manifold-dependent features tend to disappear from the empirical correlation integral as the ID grows, quickly converging to the FCI prediction for angularly uniform data.In the top left panel of Fig. 2, we highlight this effect by showing the empirical correlation integral for three single instances (d = 5, 15, 30) of the dataset uniformly drawn from {0, 1} d , linearly embedded in D = 60 dimensions (D d,60 ).Among the many datasets breaking rotational invariance that we have investigated, this is one of those where manifold-dependent features are more pronounced: we observe a ladder-type pattern that quickly disappears as the ID grows from 5 to 30.Nonetheless, even in the low-dimensional case, where the steps may affect the ID estimation, the FCI method works.We stress that in many relevant cases, including linearly embedded and uniformly sampled hypercubes (H d,D ), deviations from (2) are negligible even in low dimension.In the bottom left panel of Fig. 2, we substantiate this point by comparing the empirical FCI the three datasets D 15,60 , H 15,60 , and a rotationally invariant dataset sampled with radial Gaussian distribution G 15,60 .
The FCI estimator shares some similarities with the one recently introduced in [19] by Granata and Carnevale; here the authors use a derivative of our empirical correlation integral as the nonlinear model to fit the mid-height section of the curve.Compared to the method of Granata and Carnevale, the FCI estimator has two additional major strengths.First, the normalization procedure sets a common typical scale for all datasets, making the comparison with (2) straightforward.Second, and perhaps more importantly, our non-linear fit is performed by taking into account the whole functional form in (2), and not only the mid-height local portion of the empirical FCI; as the top left panel of Fig. 2 shows, manifold-dependent features are hardly avoided if the fit is performed locally.
Our estimator is robust to extreme undersampling and to noise, as highlighted in the center and right panels of Fig. 2 for the H d,D dataset.In a broad ID range (up to d = 10 3 ), we provide accurate estimates even in the regime N < d with a relative average error that decays quickly with the number of samples N (almost independently from the ID); as a matter of fact, for N = 100 the error is already below one percent.The method is also particularly robust to Gaussian noise, showing a sharp crossover between the regime where the ID is correctly retrieved to a phase where the noise covers the signal.
So far, we have applied and verified the performance of our algorithm only on linearly embedded manifolds.More in general, we have verified that the FCI estimator correctly identifies the ID even in the case of simple non-linear polynomial embeddings or of slightly curved manifolds.However, the method presented above is global, thus it is expected to fail on manifolds with high intrinsic curvature.As in the case of the global PCA [15,23], we can overcome this issue by providing a suitable multiscale generalization of the FCI estimator.The multiscale FCI method selects single points in the dataset and their neighbours at a fixed maximum distance r cutoff , which is then varied.The FCI estimator is then used on each neighbourhood, giving an estimation for a local ID dest(r cutoff ).Crucially, the robustness to extreme local undersampling of the FCI estimator allows to shrink the radius of the neighbourhoods r cutoff , giving a reliable estimate of local IDs, that appear as pronounced plateaux in the dest vs r cutoff plot.Alternatively, one can look at dest(n), where n is the number of nearest neighbours used in the estimation, keeping more control on the number of points used in the local estimation.(Left) We present an illustrative application of the multiscale FCI method in the case of the Swiss Roll dataset (N = 2000) for two particular samples.Both samples hint to the correct ID estimation dest = 2.We observe that the sample extracted from the highly curved inner region of the Swiss Roll provides an overestimated ID whereas the sample from the outer and flatter region exhibits a plateau at the correct ID.This suggests to use the minimum dest reached as an estimator for the true ID, but a more careful analysis is needed.(Top center) We apply the multiscale FCI method to a highly curved manifold C 6,12 (N = 10000) introduced in [13], challenging to all ID estimators (see also Materials and methods for its parametric defintion).The local dest spans the range between the true ID dest ∼ 6 to the embedding dimension D ∼ 12. (Bottom center) We show the multiscale FCI analysis on a multidimensional manifold which is built as the union of two intersecting hypercubes datasets H 20,50 and H 30,50 , each one consisting of N = 1000 samples.At small number of neighbours, we observe spurious effects due to vary sparse sampling of the neighbourhoods; quick convergence to the true IDs is then observed.(Right) As a last validation test, we generate an artifical dataset of bitmap images with multiple "blobs" with five degrees of freedom each (see Materials and methods), that we use as a proxy for curved manifolds of transformations of high contrast images.The multiscale analysis works nicely in either the one (d = 5) and three (d = 15) blob cases, although for multiple blobs we observe that the high curvature reflects in band that spans more than ten dimensions.

Multiscale FCI estimator
We perform a multiscale analysis of the FCI estimator by selecting a random sample x 0 in the dataset and a cutoff radius r c .We then apply the FCI estimation method to the set of points whose distance from x 0 is less then r c .In this way, we obtain a local ID estimate d x0 (r c ) that depends on the cutoff radius.Varying r c and x 0 , we obtain a family of curves that describes the local ID of the dataset at different scales (see also left panel in Fig. 3).Another possibility to perform the multiscale analysis is to control the number of neighbors used for the local ID estimation; this parameter is clearly in one-to-one correspondence with the cutoff radius.First, we look at a paradigmatic curved manifold studied in the literature, the Swiss roll SR 2,3 (left panel in Fig. 3).In general, we observe three regimes: for very small cutoff radius, the local estimation by the FCI is not reliable due to the extreme scarcity of neighbouring points (N 20).For very large cutoff radius, the local ID converges to the global FCI estimation, as more and more samples lie inside the cutoff radius.In between, if the manifold is sampled densely enough, we observe a plateau.In the central bottom panel of Fig. 3, we show two representative samples of the SR 2,3 dataset (N = 2000), chosen from regions with very different curvatures.When the curvature is large, the height of the plateau identifies an overestimated ID, which drifts towards the embedding dimension as expected.On the other hand, when the curvature is small, the correct ID is identified.This observation suggests that the best estimator for the multiscale analysis could be the minimum ID identified in the plateau region.use this heuristic in the tests below, but we leave a more detailed investigation of the multiscale method for future work.Now we move to multidimensional datasets and non-trivially curved manifolds, which are considered challenges for state-of-the-art ID estimators.In the bottom center panel of Fig. 3, we show the multiscale analysis for an instance of the dataset H 20,50 ∪ H 30,50 (N = 1000 + 1000), representing two intersecting hypercubes with different IDs.Two plateaux at d = 20 and d = 30 are clearly visible and allow to infer the multidimensionality of the dataset (in this case the measured local dimension is displayed as a function of the number of neighbors used for the estimation).In the top center panel of Fig. 3, we display the same analysis for an instance of the dataset C 6,12 (N = 2500), first introduced in [13] and considered a challenging dataset for ID estimation for its high intrinsic curvature.Although a thorough comparison of the performance of our algorithm with state-of-art ID estimators is out of the scope of our investigation (see [18] or [26] for nice recent meta-analysis), we observe that our prediction d 5.9 is pretty accurate (state-of-the-art estimators such as DANCo find d 6.9 on the C 6,12 dataset sampled in the same conditions).
This excellent accuracy on highly curved manifolds, combined -at the same time-with the removal of the well known underestimation issue common to all geometric methods, provides two major advantages of our algorithm over other traditional schemes used for ID estimation and suggests to test and validate it on the manifolds of global transfomations (e.g.translations, rotations, dilations) generated by high contrast images, relevant in the context of invariant object recognition.These manifolds often feature high local curvature, since even infinitesimal transformations produce almost orthogonal tangent spaces [27,28].
As a more demanding test for our multiscale FCI estimator, we report a preliminary investigation of a manifold (artificially generated) that belongs to this class, where we can keep under control the intrinsic dimensionality (see right panels of fig.3).We consider bitmap images with multiple blobs (possibly overlapping) with five degrees of freedom each (two for translations, one for rotation, two for asphericity and dilation, see also Materials and methods).Even here the multiscale analysis provides a reliable indicator, as we can convince by looking at the one and three blob cases (d = 5 and d = 15 respectively).It is worth noticing that whereas single points are affected by the high curvature of the manifold (resulting in higher ID estimations), the minimum of the plateaux works again nicely as the estimator of the correct ID.

Discussion
In this manuscript we introduced the FCI estimator for the ID of spherically sampled, linearly embedded datasets, showing that it is robust to noise and non-idealities and more importantly that it works effectively in the extreme undersampled regime (N < d).
We performed a multiscale analysis of the FCI estimator on challenging datasets, featuring high curvature and multidimensionality, showing that we can extract the correct ID as the minimum local d est .Further work will be needed to fully explore this observation and to construct a proper multiscale estimator.
We leave open to future investigations the analysis of high dimensional manifolds of high contrast images taken from the Machine Learning literature, as well as the possibility of combining our estimator with state-of-the-art techniques for dimensional reduction and manifold learning, or, even more ambitiously, to elaborate on it in order to propose a novel more effective toolbox for these tasks.

Average correlation integral for uniformly sampled hyperspheres
Here we derive the average correlation integral ((2)) for a dataset uniformly sampled from the hypersphere S d or radius r s , i.e.

Empirical full correlation integral and fitting procedure
The empirical full correlation integral is easily computed in Mathematica 12 using the following one-liner: The non-linear fit for the FCI estimator was performed using the default function FindFit, with d and r s as free parameters.The empirical full correlation integral was preprocessed before the application of FindFit by extracting a RandomSample of min 1000, N (N −1) 2 of its points to speed-up the fitting procedure.

Description of the datasets
In this section we briefly describe the datasets used in the presented numerical simulations.In the following, a linear embedding is the map ι : R d → R D that appends D − d zeros to its argument, and then rotates it in R D by a randomly chosen rotation matrix.

High-contrast images datasets
The high-contrast images shown in Figure 3 (right panel) were generated by assigning to each pixel of a l × l bitmap the following value v i,j : a i,j = cos θ(j − dx) + sin θ(i + dy) b i,j = − sin θ(j − dx) + cos θ(i + dy) v i,j = 1 − a 2 i,j + e 2 b 2 i,j (1 + e 2 )s 2 Morover, any pixel of value less then 0.01 was manually set to 0 to increase the contrast of the image.

Figure 1 :
Figure 1: Standard projective [Principal Component Analysis (PCA)] and geometric [Correlation Dimension (CorrDim) and its generalizations] methods for intrinsic dimension estimation fail due to undersampling.(Top left) PCA estimates the intrinsic dimension d of linearly embedded datasets by detecting abrupt jumps in the magnitude of the sorted eigenvalues of the corresponding correlation matrix.This method works whenever the number of samples N in the dataset is sufficiently large (N d log d), and the jump occurs between the d-th and (d + 1)-th eigenvalues (green line in the plot corresponding to a dataset of N = 1000 samples drawn uniformly from a hypercube in dimension d = 100).In the opposite undersampled regime N < d, PCA is inconclusive (orange line in the plot, same dataset with N = 50 samples).(Bottom left) Geometric methods, such as the CorrDim or the more general k-nearest-neighbours estimators, are based on the scaling of the density of neighbours ρ(r/r 0 ) at a small cutoff distance r with respect to the average diameter r 0 of the dataset.In particular, ρ(r/r 0 ) ∼ (r/r 0 ) d for r/r 0 → 0 independently on the details of the datasets, so that a log-log linear fit estimates the ID as the slope of the fitting line.However, the small r regime is exponentially (in the ID,[24]) difficult to sample.This effect is at the origin of the so-called curse of dimensionality, and it induces a systematic underestimation of the ID.As N increases, though one is able to sample smaller r regions as shown in the plot (orange line N = 50, green line N = 1000 as above), convergence to the true dimension (d = 100 in the plot) is not possible.(Right) The density of neighbours ρ is more generally defined at any cutoff distance r.The r r 0 regime is the one used to compute CorrDim.In the remaining region r r 0 , the density ρ increases and eventually approaches one, indicating that the underlying proximity graph (see insets) gets more and more connected.We observe that this region is easily sampled at any fixed N (plot at N = 50), but the functional form of ρ in this informative regime is in principle dependent on the details of the dataset.

Figure 2 :
Figure 2: Intrinsic Dimension (ID) estimation is possible in the extreme undersampled regime for arbitrarly large ID with the Full Correlation Integral (FCI) estimator, in the case of linearly embedded and slightly curved manifolds, possibily non uniformly sampled and with noise.(Top left) We show the density of neighbours ρ of preprocessed (centered and normalized) data (number of samples N = 500) extracted from {0, 1} d linearly embedded in D = 60 dimensions (D d,60), for d =(5,15, 30).We are able to efficiently extract the correct ID even though this is a highly non-uniformly sampled dataset, whose ρ displays manifold-dependent features (in this case, step-like patterns).Moreover, we observe that, as we increase d, the density of neighbours of this dataset quickly converges to our functional form.It is worth noticing that the whole functional form (2) is needed for the fit; in fact, a local fit of the slope of ρ at half-height would result in an incorrect ID estimation.(Bottom left) We show the density of neighbours ρ of preprocessed data (N = 500) extracted uniformly from {0, 1} d , [0, 1] d and from R d with multivariate gaussian distribution, for d = 15 and linearly embedded in D = 60 dimensions.All plot lines are compatible with the same functional form (2), pointing to an intriguing manifestation of "universality" for high-dimensional data.(Center) To highlight the predictive power of the FCI method for a broad spectrum of dimensionalities (ranging from d = 4 to d = 200), we exhibit the estimated ID versus the number of sample points N for the linearly embedded hypercube H d,500 .Error bars are computed by averaging over 10 samples for each pair (N, d).(Top right) We asses quantitatively the predictive power of the FCI method by computing the average relative error |(dest − d)/d| (over 20 random instances) of the estimated ID in the range 5 ≤ d ≤ 1000, 5 ≤ N ≤ 1000.We observe that at N ∼ 100 we have an error of the order of 1% almost independently on the ID, and that ID estimation is possible also in the extreme undersampled N < d regime.(Bottomright) The FCI method estimates the correct ID even when the data are corrupted by noise.Here we consider a linearly embedded hypercube dataset H 40,60 and add on the top of that a 60-dimensional gaussian noise of standard deviation σ.We observe a sharp transition in the estimated ID between the regime in which the noise is a perturbation (σ 0.1 and dest = 40) and the regime in which the noise covers the signal (σ 0.2 and dest = 60).
) is needed for the fit; in fact, a local fit of the slope of ρ at half-height would result in an incorrect ID estimation.(Bottom left) We show the density of neighbours ρ of preprocessed data (N = 500) extracted uniformly from {0, 1} d , [0, 1] d and from R d with multivariate gaussian distribution, for d = 15 and linearly embedded in D = 60 dimensions.All plot lines are compatible with the same functional form (2), pointing to an intriguing manifestation of "universality" for high-dimensional data.(Center) To highlight the predictive power of the FCI method for a broad spectrum of dimensionalities (ranging from d = 4 to d = 200), we exhibit the estimated ID versus the number of sample points N for the linearly embedded hypercube H d,500 .Error bars are computed by averaging over 10 samples for each pair (N, d).(Top right) We asses quantitatively the predictive power of the FCI method by computing the average relative error |(dest − d)/d| (over 20 random instances) of the estimated ID in the range 5 ≤ d ≤ 1000, 5 ≤ N ≤ 1000.

d = 5 , D = 81 2 d = 15 , D = 81 2 N = 4000 Figure 3 :
Figure 3:The multiscale generalization of the FCI provides a state-of-the-art tool to tackle the ID estimation for complex datasets with multidimensional features and high intrinsic curvature.The multiscale FCI method selects single points in the dataset and their neighbours at a fixed maximum distance r cutoff , which is then varied.The FCI estimator is then used on each neighbourhood, giving an estimation for a local ID dest(r cutoff ).Crucially, the robustness to extreme local undersampling of the FCI estimator allows to shrink the radius of the neighbourhoods r cutoff , giving a reliable estimate of local IDs, that appear as pronounced plateaux in the dest vs r cutoff plot.Alternatively, one can look at dest(n), where n is the number of nearest neighbours used in the estimation, keeping more control on the number of points used in the local estimation.(Left) We present an illustrative application of the multiscale FCI method in the case of the Swiss Roll dataset (N = 2000) for two particular samples.Both samples hint to the correct ID estimation dest = 2.We observe that the sample extracted from the highly curved inner region of the Swiss Roll provides an overestimated ID whereas the sample from the outer and flatter region exhibits a plateau at the correct ID.This suggests to use the minimum dest reached as an estimator for the true ID, but a more careful analysis is needed.(Top center) We apply the multiscale FCI method to a highly curved manifold C 6,12 (N = 10000) introduced in[13], challenging to all ID estimators (see also Materials and methods for its parametric defintion).The local dest spans the range between the true ID dest ∼ 6 to the embedding dimension D ∼ 12. (Bottom center) We show the multiscale FCI analysis on a multidimensional manifold which is built as the union of two intersecting hypercubes datasets H 20,50 and H 30,50 , each one consisting of N = 1000 samples.At small number of neighbours, we observe spurious effects due to vary sparse sampling of the neighbourhoods; quick convergence to the true IDs is then observed.(Right) As a last validation test, we generate an artifical dataset of bitmap images with multiple "blobs" with five degrees of freedom each (see Materials and methods), that we use as a proxy for curved manifolds of transformations of high contrast images.The multiscale analysis works nicely in either the one (d = 5) and three (d = 15) blob cases, although for multiple blobs we observe that the high curvature reflects in band that spans more than ten dimensions.

l
side of the bitmap fixed at 81 ∆x horizontal translation uniform in (−, ) ∆y vertical translation uniform in (−, ) s size uniform in () e eccentricity uniform in () θ angle of the major axis uniform in (−π/2, π/2)