Pith. sign in

REVIEW 3 major objections 4 minor 34 references

Estimation of Expected Euler Characteristic Curves of Nonstationary Smooth Gaussian Random Fields

T0 review · 3 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash

Pith's one-line read The paper shows that a Hermite projection of pinned Euler characteristic curves unbiasedly estimates the Lipschitz-Killing curvatures of nonstationary smooth Gaussian fields, and that a multiplier bootstrap extends this to the limiting…

desk verdict Solid new estimator for LKCs via Hermite projection, with a bootstrap extension that is promising but has a real proof gap in the consistency theorem that needs an extra assumption. read the letter →

arxiv 1908.02493 v2 pith:475AWMYL submitted 2019-08-07 math.ST stat.TH

classification math.STstat.TH MSC 62M4060G6062G09
keywords randomfieldtheoryLipschitz-KillingcurvaturesexpectedEulercharacteristiccurveGaussianrelatedfieldsHermiteprojectionestimatormultiplierbootstrapnonstationaryexcursionsets
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The expected Euler characteristic curve of a smooth Gaussian field is a finite linear combination of known density functions, with coefficients given by the Lipschitz-Killing curvatures of the domain under the metric induced by the field. The paper proves that those coefficients can be estimated by projecting the observed, pinned Euler characteristic curve onto the density functions using a weighted Hermite inner product. Under smoothness and finite-moment conditions on the field's critical points, the resulting Hermite projection estimator is unbiased, has finite variance, and averages over repeated observations to a consistent, asymptotically normal estimate. A multiplier bootstrap version applies the same projection to Gaussian multiplier fields built from standardized residuals, consistently estimating the curvatures of the Gaussian limiting field of non-Gaussian or non-iid data. This matters because LKC estimation no longer needs stationarity or isotropy, so EEC-based thresholding can be used on arbitrary curved domains and on data such as brain images and cosmic microwave background maps.

What carries the argument

The workhorse is the Hermite projector on a weighted $L^2$ space with inner product $\langle g,h\rangle=\int g(u)h(u)e^{u^2/2}\,du$: the EC densities $\rho_d(u)=(2\pi)^{-(d+1)/2}H_{d-1}(u)e^{-u^2/2}$ are orthogonal in this inner product, so each LKC is the projection coefficient of the pinned EEC curve onto $\rho_d$. The estimator evaluates this indefinite integral as a finite sum over the field's critical values, because the EC curve is constant between critical levels. The bootstrap extension constructs a Gaussian multiplier field $R_g^{(N)}=\sum_{n=1}^N g_n R_n$ from standardized residuals $R_n$, which conditional on the data is a mean-zero Gaussian field with covariance equal to the empirical correlation; applying the projector to this field and letting the number of multiplier samples grow yields the bHPE.

What would settle it

Simulate data from a non-Gaussian field with a known Gaussian limit but a spatial correlation whose empirical second derivatives converge non-uniformly, such as spatially varying smoothing; run the bHPE for increasing sample sizes and check whether the estimates approach the true limiting LKCs. If the estimates fail to approach despite marginal correlations looking correct, the uniform-convergence condition (R2) is doing the load-bearing work.

Watch

Extended reading notes

Core claim

The central claim is that the Lipschitz-Killing curvatures are recoverable from a single observed Euler characteristic curve by linear projection: after pinning the curve by subtracting the known baseline $L_0\Phi^+$, the $d$-th LKC is estimated by $\hat L_d = (2\pi)^{d/2}/(d-1)! \int_{-\infty}^{\infty} H_{d-1}(u)\,\chi_f^\circ(u)\,du$. For Gaussian fields satisfying assumptions (G1)-(G4), this estimate is unbiased and has finite variance, and its sample average is consistent and asymptotically normal. For non-Gaussian fields satisfying a functional CLT, the bootstrap HPE applies the same projector to Gaussian multiplier fields generated from standardized residuals; the estimator equals the LKCs of the Gaussian field with the empirical correlation, and converges to the LKCs of the limiting field as the sample size grows. Plugging either LKC estimate into the Gaussian kinematic formula gives a smooth estimator of the EEC curve that inherits a functional CLT and yields confidence bands and consistent threshold estimates.

Load-bearing premise

The bootstrap estimator rests on the empirical correlation of the standardized residuals, and its derivatives up to order two, converging uniformly to the true correlation; if that convergence fails, the constructed Gaussian multiplier field does not have the right geometry and the LKC estimates no longer converge to the limiting field's LKCs.

Editorial extensions

If this is right

  • LKC estimation no longer requires stationarity or isotropy: any domain where the Euler characteristic of excursion sets can be computed, including curved surfaces like the sphere or the cortical surface, becomes usable.
  • The smooth HPE of the EEC curve is an orthogonal projection of the nonparametric average and has lower variance; the paper's functional CLT gives both pointwise and simultaneous confidence bands for the true EEC curve.
  • The threshold solving $\widehat{\mathrm{EEC}}(u)=\alpha$ is a consistent estimator of the excursion threshold and is asymptotically Gaussian, so significance levels for fMRI activation maps and the expected number of false clusters can be reported with standard errors.
  • For non-Gaussian observations, the bHPE estimates the LKCs of the Gaussian limiting field without knowing the field's mean or variance, extending the method to general linear models and standardized residual fields.
  • In the cosmic microwave background example, the observed field's second LKC lies about 2.93 standard deviations from the mean of the simulations, giving a quantitative topological comparison between observed and simulated sky maps.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The same projection structure should extend to Gaussian-related fields covered by the Gaussian kinematic formula, such as chi-squared and t fields, by replacing the EC densities with the corresponding Gaussian-related densities; the paper notes the extension for plug-in EEC estimators but does not develop the HPE there.
  • The bHPE's consistency leans on uniform convergence of the empirical correlation and its derivatives, so in short samples the estimator is effectively estimating the LKCs of the empirical correlation; a practical diagnostic would be to compare bHPE outputs under different residual constructions on the same data.
  • The critical-value representation of the HPE suggests an online or streaming implementation: an algorithm that passes through the field once and records critical values and topology changes yields the estimator without storing the entire field, connecting naturally to streaming Euler-curve algorithms.
  • The CMB comparison points to a reusable diagnostic: comparing observed LKC estimates to simulation means can flag physical model mismatch, and the same logic could be applied to other cosmological maps to probe non-Gaussianity in the early universe.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. The paper proposes Hermite projection estimators (HPE) for the Lipschitz-Killing curvatures (LKCs) of nonstationary smooth Gaussian random fields, based on projecting the observed Euler characteristic (EC) curve onto the EC densities in a weighted L2 space. A multiplier bootstrap variant (bHPE) is introduced to estimate the LKCs of the Gaussian limiting field of non-Gaussian or non-iid data. The paper proves unbiasedness, finite variance, consistency, and CLTs for the HPE and the resulting plug-in EEC estimator, derives asymptotic confidence bands, and analyzes the estimated detection threshold. The methods are evaluated on simulated 2D fields and applied to CMB and fMRI data. The main results are Theorems 1-5, with proofs in the appendix and code provided online.

Significance. If the results hold, this is a valuable contribution to random field theory and its applications in neuroimaging and cosmology. It offers a computationally simple, theoretically grounded alternative to stationarity-based LKC estimators, and the multiplier bootstrap extension to non-Gaussian limiting fields is novel and practically important. The paper ships reproducible code, includes extensive simulations comparing against existing methods, and gives concrete, falsifiable predictions (e.g., unbiasedness of the HPE under (G1)-(G4a), consistency of the bHPE under (R1)-(R2)). The proofs of Theorems 1, 2, 4, and 5 are largely rigorous, and the simulation results convincingly support the finite-sample claims. However, the proof of the central bHPE consistency theorem (Theorem 3(ii)) has a gap, and the theorem relies on an unverified regularity assumption about the Gaussian multiplier field, so the full consistency claim is not yet established.

major comments (3)
  1. [Section 3, Theorem 3(ii); Supplementary B.4] The proof of the convergence L(ˆr^{(N)}) → L(r) in Theorem 3(ii) requires convergence of the fourth-order mixed partial derivatives of the empirical correlation, because the Riemannian curvature tensor is expressed in B.4 in terms of E[∂d∂d' R_g^{(N)} ∂d''∂d''' R_g^{(N)}] = ∂^4_{s,s'} ˆr^{(N)}(s,s')|_{s=s'}. However, assumption (R2) only assumes uniform almost sure convergence of ˆr^{(N)} and its partial derivatives up to order 2. The sentence in B.4 asserting that this fourth-order quantity 'converges by (R2)' is therefore not justified. Since the LKC map is not, in general, a continuous function of only the second derivatives of the correlation, the convergence L(ˆr^{(N)}) → L(r) is an additional assumption rather than a consequence of (R2). This is a load-bearing gap for the main bHPE consistency claim. The authors should strengthen (R2) to include uniform convergence of derivatives up to order 4 (with appropriate domination) or prove continuity of the LKC functional under a weaker metric that only involves second derivatives.
  2. [Section 3, Theorem 3] The theorem assumes that, for almost all residuals R1,...,RN and all sufficiently large N, the Gaussian multiplier field R_g^{(N)} satisfies (G1)-(G4a). No sufficient conditions on the observed fields or residuals are provided under which this assumption holds. Since (G4) is itself admitted in Section 2 to be an active research topic, the bHPE consistency result is conditional on an unverified regularity property. The paper should either supply checkable conditions under which the GMF inherits (G1)-(G4a) from the residuals, or explicitly state Theorem 3 as a conditional result. As written, the abstract's claim that the bHPE 'allows estimation of LKCs ... of the limiting field' is stronger than what the theorem establishes.
  3. [Section 2, Corollary 1; Supplementary B.3] The proof of Corollary 1 infers superpolynomial decay of Cov[χf(u), χf(v)] from the finiteness of the double integrals in (10) for all d,d'. Finiteness of all polynomial moments of an integrable function does not by itself imply pointwise decay faster than every polynomial unless additional regularity (e.g., monotonicity or continuity of the covariance kernel in its tails) is assumed. This does not affect the main theorems, but as written the corollary's proof is incomplete. The authors should either prove the decay under explicit additional conditions or reformulate the corollary as a statement about polynomial moments of the covariance kernel.
minor comments (4)
  1. [Section 3, Theorem 3 statement] The phrase 'for almost all N > N′' is awkward because N is the fixed sample size in the theorem. It would be clearer to state 'for all N ≥ N0 almost surely' or 'for all sufficiently large N, almost surely over the residuals.'
  2. [Section 2, Theorem 2(ii)] In the statement 'the covariances σdd′ are finite for all d,d′ = 1,...,N', the upper limit N appears to be a typo and should be D, the dimension of the manifold.
  3. [Section 3, assumption (R2)] Assumption (R2) says that ˆr^{(N)}(s,s') and its partial derivatives up to order 2 converge uniformly, but the proof of Theorem 3(ii) uses fourth-order mixed partial derivatives. The notation should clarify whether derivatives are taken with respect to both arguments s and s', and the assumption should be stated in a form that makes the required order explicit.
  4. [Section 4, Theorem 4] Theorem 4 is stated 'under the assumptions (G1)-(G4)', but part (i) only requires that the arbitrary estimator ˆL^{(N)} be consistent and part (ii) only requires the CLT for ˆL^{(N)}. The hypotheses should be stated directly in terms of the estimator ˆL^{(N)} rather than referencing (G1)-(G4), which are conditions on the underlying field and not on the estimator.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the HPE and bHPE are plug-in/projection estimators whose theory rests on the externally cited Gaussian kinematic formula; the proof gap in Theorem 3 is a regularity issue, not a circular reduction.

full rationale

The central derivation is not circular. The HPE defines the LKCs as Hermite-projection coefficients of the EEC curve via the Gaussian kinematic formula (Eq. 6), and estimates them by applying the same linear functional to an observed EC curve (Eq. 9). This is a plug-in/method-of-moments estimator, not a fit: no parameter is adjusted to reproduce a target LKC, and unbiasedness follows from linearity plus the external GKF. The bHPE likewise uses Gaussian multiplier fields whose conditional correlation is the empirical correlation (Eq. 17), so the estimator is the plug-in L(rhat^(N)); consistency of plug-in requires continuity of the LKC map, which the paper attempts to show in B.4. The skeptical concern that B.4 implicitly needs fourth-order convergence of the empirical correlation while (R2) only supplies derivatives up to order 2 is a genuine correctness or assumption gap, not a circularity: the conclusion is not definitionally identical to the input. The self-citations in the paper (Telschow and Schwartzman 2019, Telschow et al. 2019) concern sufficient smoothness conditions and confidence-band constructions; they are not load-bearing for the main estimator identities and are not used to forbid alternatives. No equation in the paper reduces to its own input, and no fitted parameter is renamed as a prediction. Therefore the appropriate finding is no significant circularity, with the caveat that Theorem 3's consistency proof should be strengthened with an explicit fourth-order condition.

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

No free parameters are fitted; the estimators are linear projections with known basis functions. The paper relies on the GKF and on regularity assumptions about the fields and residuals that are stated but not derived.

assumptions (4)
  • standard math Gaussian kinematic formula (GKF) expressing EEC as finite linear combination of EC densities with LKC coefficients.
    Used in Eq (1) and throughout; the entire projection approach assumes this theorem.
  • domain assumption The random field satisfies (G1)-(G4): mean zero, variance one, C^2 paths, nondegenerate derivative distribution, local Holder condition on second derivatives, and finite moments of the number of critical points.
    Assumptions stated in Section 2; needed for unbiasedness and finite variance in Theorem 2.
  • domain assumption Standardized residuals satisfy (R1) and (R2): vanishing sum, unit sum of squares, and uniform convergence of empirical correlation and its derivatives to the true correlation.
    Used in Theorem 3 for consistency of the bHPE.
  • ad hoc to paper The Gaussian multiplier field conditioned on residuals satisfies (G1)-(G4a) for almost all N.
    Assumed in Theorem 3; no sufficient conditions are provided.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Estimation of Expected Euler Characteristic Curves of Nonstationary Smooth Gaussian Random Fields." pith.science (2026). https://pith.science/paper/475AWMYL

@misc{pith2026190802493,
  author       = {Pith},
  title        = {Pith review of: Estimation of Expected Euler Characteristic Curves of Nonstationary Smooth Gaussian Random Fields},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/475AWMYL}},
  note         = {Machine review of arXiv:1908.02493}
}
read the original abstract

The expected Euler characteristic (EEC) curve of excursion sets of a Gaussian random field is used to approximate the distribution of its supremum for high thresholds. Viewed as a function of the excursion threshold, the EEC is expressed by the Gaussian kinematic formula (GKF) as a linear function of the Lipschitz-Killing curvatures (LKCs) of the field, which solely depend on the domain and covariance function of the field. So far its use for non-stationary Gaussian fields over non-trivial domains has been limited because in this case the LKCs are difficult to estimate. In this paper, consistent estimators of the LKCs are proposed as linear projections of "pinned" observed Euler characteristic curves and a linear parametric estimator of the EEC curve is obtained, which is more efficient than its nonparametric counterpart for repeated observations. A multiplier bootstrap modification reduces the variance of the estimator, and allows estimation of LKCs and EEC of the limiting field of non-Gaussian fields satisfying a functional CLT. The proposed methods are evaluated using simulations of 2D fields and illustrated in thresholding of 3D fMRI brain activation maps and cosmological simulations on the 2-sphere.

Figures

Figures reproduced from arXiv: 1908.02493 by the authors.

Figure 1
Figure 1. (left) Single realization of the isotropic Gaussian random field (32). (middle) [PITH_FULL_IMAGE:figures/full_fig_p007_1.png] view at source ↗
Figure 2
Figure 2. Isotropic Gaussian field (ν = 5 and L = 50): comparison of mean and standard deviation of different LKC estimators. Black lines represent the theoretical LKCs. LKC estimation In the Gaussian case, [PITH_FULL_IMAGE:figures/full_fig_p022_2.png] view at source ↗
Figure 3
Figure 3. Isotropic non-Gaussian field (ν = 5 and L = 50): comparison of mean and standard deviation of different LKC estimators. Black lines represent the theoretical LKCs. 23 [PITH_FULL_IMAGE:figures/full_fig_p023_3.png] view at source ↗
Figures from the paper (8 more)
Figure 4
Figure 4. Figure 4: Isotropic Gaussian field: Covariance functions of the sample average EC curve [PITH_FULL_IMAGE:figures/full_fig_p024_4.png]
Figure 5
Figure 5. Figure 5: Isotropic Gaussian field: pointwise coverage of the EEC curve under respec [PITH_FULL_IMAGE:figures/full_fig_p025_5.png]
Figure 6
Figure 6. Figure 6: Scale space field: comparison of mean and standard deviation of different LKC [PITH_FULL_IMAGE:figures/full_fig_p026_6.png]
Figure 7
Figure 7. Figure 7: Scale space: covariance functions of the sample average EC curve (left) and the [PITH_FULL_IMAGE:figures/full_fig_p027_7.png]
Figure 8
Figure 8. Figure 8: Scale space: pointwise coverage of the EEC curve under respectively “true” [PITH_FULL_IMAGE:figures/full_fig_p028_8.png]
Figure 9
Figure 9. Figure 9: CMB example: (top) cleaned and smoothed observed CMB field. (left/right) [PITH_FULL_IMAGE:figures/full_fig_p029_9.png]
Figure 10
Figure 10. Figure 10: fMRI example: Observed EC curves and estimates with pointwise confidence [PITH_FULL_IMAGE:figures/full_fig_p043_10.png]
Figure 11
Figure 11. Figure 11: fMRI activation at a significance level of FWER=0.05 (left) and CER=1 (right), [PITH_FULL_IMAGE:figures/full_fig_p044_11.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

34 extracted references · 33 canonical work pages

  1. [1]

    P. A. Ade, N. Aghanim, C. Armitage-Caplan, M. Arnaud, M. Ashdown, F. Atrio-Barandela, J. Aumont, C. Baccigalupi, A. J. Banday, R. Barreiro, et al. Planck 2013 results. xxiii. isotropy and statistics of the cmb. Astronomy & Astrophysics, 571: 0 A23, 2014

  2. [2]

    P. A. Ade, N. Aghanim, M. Arnaud, M. Ashdown, J. Aumont, C. Baccigalupi, A. Banday, R. Barreiro, J. Bartlett, N. Bartolo, et al. Planck 2015 results-xii. full focal plane simulations. Astronomy & Astrophysics, 594: 0 A12, 2016

  3. [3]

    R. J. Adler and J. E. Taylor. Random fields and geometry. Springer Science & Business Media, 2009

  4. [4]

    R. J. Adler, K. Bartz, S. C. Kou, and A. Monod. Estimating thresholding levels for random fields via euler characteristics. arXiv preprint arXiv:1704.08562, 2017

  5. [5]

    Aza \" s and M

    J.-M. Aza \" s and M. Wschebor. Level sets and extrema of random processes and fields. John Wiley & Sons, 2009

  6. [6]

    Bartlett, M

    J. Bartlett, M. Bucher, J. Cardoso, G. Castex, J. Delabrouille, K. Ganga, Y. Giraud-H \'e raud, M. Le Jeune, G. Patanchon, M. Piat, et al. Planck 2015 results: Ix. diffuse component separation: Cmb maps. Astronomy and Astrophysics, 594: 0 A9--A9, 2016

  7. [7]

    E. T. Bullmore, J. Suckling, S. Overmeyer, S. Rabe-Hesketh, E. Taylor, and M. J. Brammer. Global, voxel, and cluster tests, by theory and permutation, for a difference between two groups of structural mr images of the brain. IEEE transactions on medical imaging, 18 0 (1): 0 32--42, 1999

  8. [8]

    Cheng and A

    D. Cheng and A. Schwartzman. Multiple testing of local maxima for detection of peaks in random fields. Annals of statistics, 45 0 (2): 0 529, 2017

Show all 34 references
  1. [9]

    Cheng, V

    D. Cheng, V. Cammarota, Y. Fantaye, D. Marinucci, A. Schwartzman, et al. Multiple testing of local maxima for detection of peaks on the (celestial) sphere. Bernoulli, 26 0 (1): 0 31--60, 2020

  2. [10]

    F. Eicker. Asymptotic normality and consistency of the least squares estimators for families of linear regressions. The Annals of Mathematical Statistics, pages 447--456, 1963

  3. [11]

    Eklund, T

    A. Eklund, T. E. Nichols, and H. Knutsson. Cluster failure: Why fmri inferences for spatial extent have inflated false-positive rates. Proceedings of the national academy of sciences, 113 0 (28): 0 7900--7905, 2016

  4. [12]

    Estrade and J

    A. Estrade and J. Fournier. Number of critical points of a gaussian random field: Condition for a finite variance. Statistics & Probability Letters, 118: 0 94--99, 2016

  5. [13]

    K. M. Gorski, E. Hivon, A. J. Banday, B. D. Wandelt, F. K. Hansen, M. Reinecke, and M. Bartelmann. Healpix: A framework for high-resolution discretization and fast analysis of data distributed on the sphere. The Astrophysical Journal, 622 0 (2): 0 759, 2005

  6. [14]

    J. R. Gott, A. L. Melott, and M. Dickinson. The sponge-like topology of large-scale structure in the universe. The Astrophysical Journal, 306: 0 341--357, 1986

  7. [15]

    Heiss and H

    T. Heiss and H. Wagner. Streaming algorithm for euler characteristic curves of multidimensional images. In International Conference on Computer Analysis of Images and Patterns, pages 397--409. Springer, 2017

  8. [16]

    Hikage, Y

    C. Hikage, Y. Suto, I. Kayo, A. Taruya, T. Matsubara, M. S. Vogeley, F. Hoyle, J. R. Gott III, J. Brinkmann, and S. collaboration. Three-dimensional genus statistics of galaxies in the sdss early data release. Publications of the Astronomical Society of Japan, 54 0 (5): 0 707-...

  9. [17]

    S. J. Kiebel, J.-B. Poline, K. J. Friston, A. P. Holmes, and K. J. Worsley. Robust smoothness estimation in statistical parametric maps using standardized residuals from the general linear model. Neuroimage, 10 0 (6): 0 756--766, 1999

  10. [18]

    Liebl and M

    D. Liebl and M. Reimherr. Fast and fair simultaneous confidence bands for functional parameters. arXiv preprint arXiv:1910.00131, 2019

  11. [19]

    J. M. Moran, E. Jolly, and J. P. Mitchell. Social-cognitive deficits in normal aging. Journal of neuroscience, 32 0 (16): 0 5553--5561, 2012

  12. [20]

    T. E. Nichols. Multiple testing corrections, nonparametric methods, and random field theory. Neuroimage, 62 0 (2): 0 811--815, 2012

  13. [21]

    Pranav, H

    P. Pranav, H. Edelsbrunner, R. Van de Weygaert, G. Vegter, M. Kerber, B. J. Jones, and M. Wintraecken. The topology of the cosmic web in terms of persistent betti numbers. Monthly Notices of the Royal Astronomical Society, 465 0 (4): 0 4281--4310, 2017

  14. [22]

    Pranav, R

    P. Pranav, R. J. Adler, T. Buchert, H. Edelsbrunner, B. J. Jones, A. Schwartzman, H. Wagner, and R. Van de Weygaert. Unexpected topology of the temperature fluctuations in the cosmic microwave background. Astronomy & Astrophysics, 627: 0 A163, 2019

  15. [23]

    D. O. Siegmund and K. J. Worsley. Testing for a signal with unknown location and scale in a stationary gaussian random field. The Annals of Statistics, pages 608--639, 1995

  16. [24]

    Sommerfeld, S

    M. Sommerfeld, S. Sain, and A. Schwartzman. Confidence regions for spatial excursion sets from repeated random field observations, with an application to climate. Journal of the American Statistical Association, 113 0 (523): 0 1327--1340, 2018

  17. [25]

    Taylor, A

    J. Taylor, A. Takemura, R. J. Adler, et al. Validity of the expected euler characteristic heuristic. The Annals of Probability, 33 0 (4): 0 1362--1396, 2005

  18. [26]

    J. E. Taylor and K. J. Worsley. Detecting sparse signals in random fields, with an application to brain mapping. Journal of the American Statistical Association, 102 0 (479): 0 913--928, 2007

  19. [27]

    J. E. Taylor et al. A gaussian kinematic formula. The Annals of Probability, 34 0 (1): 0 122--158, 2006

  20. [28]

    F. J. Telschow and A. Schwartzman. Simultaneous confidence bands for functional data using the gaussian kinematic formula. arXiv preprint arXiv:1901.06386, 2019

  21. [29]

    F. J. Telschow, M. R. Pierrynowski, and S. F. Huckemann. Confidence tubes for curves on so (3) and identification of subject-specific gait change after kneeling. arXiv preprint arXiv:1909.06583, 2019

  22. [30]

    A. W. Van der Vaart. Asymptotic statistics, volume 3. Cambridge university press, 2000

  23. [31]

    K. J. Worsley. An unbiased estimator for the roughness of a multivariate Gaussian random field. Department of Mathematics and Statistics, McGill University Montr \'e al, Canada, 1996

  24. [32]

    K. J. Worsley, A. C. Evans, S. Marrett, and P. Neelin. A three-dimensional statistical analysis for cbf activation studies in human brain. Journal of Cerebral Blood Flow & Metabolism, 12 0 (6): 0 900--918, 1992

  25. [33]

    K. J. Worsley, S. Marrett, P. Neelin, A. C. Vandal, K. J. Friston, and A. C. Evans. A unified statistical approach for determining significant signals in images of cerebral activation. Human brain mapping, 4 0 (1): 0 58--73, 1996

  26. [34]

    K. J. Worsley, J. E. Taylor, F. Tomaiuolo, and J. Lerch. Unified univariate and multivariate random field theory. Neuroimage, 23: 0 S189--S195, 2004

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.