Source author record

Dan Cheng

Dan Cheng appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

12works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 published item(s)

preprint2026arXiv

Multiple Testing of Local Extrema for Detection of Structural Breaks in Piecewise Linear Models

In this paper, we propose a new generic method for detecting the number and locations of structural breaks or change points in piecewise linear models under stationary Gaussian noise. Our method transforms the change point detection problem into identifying local extrema (local maxima and local minima) through kernel smoothing and differentiation of the data sequence. By computing p-values for all local extrema based on peak height distributions of smooth Gaussian processes, we utilize the Benjamini-Hochberg procedure to identify significant local extrema as the detected change points. Our method can distinguish between two types of change points: continuous breaks (Type I) and jumps (Type II). We study three scenarios of piecewise linear signals, namely pure Type I, pure Type II and a mixture of Type I and Type II change points. The results demonstrate that our proposed method ensures asymptotic control of the False Discover Rate (FDR) and power consistency, as sequence length, slope changes, and jump size increase. Furthermore, compared to traditional change point detection methods based on recursive segmentation, our approach only requires a single test for all candidate local extrema, thereby achieving the smallest computational complexity proportionate to the data sequence length. Additionally, numerical studies illustrate that our method maintains FDR control and power consistency, even in non-asymptotic cases when the size of slope changes or jumps is not large. We have implemented our method in the R package "dSTEM" (available from https://cran.r-project.org/web/packages/dSTEM).

preprint2023arXiv

An approximation to peak detection power using Gaussian random field theory

We study power approximation formulas for peak detection using Gaussian random field theory. The approximation, based on the expected number of local maxima above the threshold $u$, $\mathbb{E}[M_u]$, is proved to work well under three asymptotic scenarios: small domain, large threshold, and sharp signal. An adjusted version of $\mathbb{E}[M_u]$ is also proposed to improve accuracy when the expected number of local maxima $\mathbb{E}[M_{-\infty}]$ exceeds 1. Cheng and Schwartzman (2018) developed explicit formulas for $\mathbb{E}[M_u]$ of smooth isotropic Gaussian random fields with zero mean. In this paper, these formulas are extended to allow for rotational symmetric mean functions, so that they are suitable for power calculations. We also apply our formulas to 2D and 3D simulated datasets, and the 3D data is induced by a group analysis of fMRI data from the Human Connectome Project to measure performance in a realistic setting.

preprint2023arXiv

The expected Euler characteristic approximation to excursion probabilities of Gaussian vector fields

Let $\{(X(t), Y(s)): t\in T, s\in S\}$ be an $\mathbb{R}^2$-valued, centered, unit-variance smooth Gaussian vector field, where $T$ and $S$ are compact rectangles in $\mathbb{R}^N$. It is shown that, as $u\to \infty$, the joint excursion probability $\mathbb{P} \{\sup_{t\in T} X(t) \geq u, \sup_{s\in S} Y(s) \geq u \}$ can be approximated by $\mathbb{E}\{χ(A_u)\}$, the expected Euler characteristic of the excursion set $A_u=\{(t,s)\in T\times S: X(t) \ge u, Y(s) \ge u\}$, such that the error is super-exponentially small. This verifies the expected Euler characteristic heuristic (cf. Taylor, Takemura and Alder (2005), Alder and Taylor (2007)) for a large class of smooth Gaussian vector fields.

preprint2022arXiv

Ultrashort Pulse Generation in Modeless Laser Cavity

We demonstrate experimentally that random phase modulation of an erbium-doped fiber ring-laser by an intra-cavity electro-optic phase modulator did not inhibit ultrashort-pulse operation. Stable and self-starting ultrashort-pulse operation with a single pulse circulating in the cavity was achieved even when the phase modulator was driven with random sequences sufficiently fast and strong to render the laser cavity modeless, in the sense that heterodyning of the laser output did not show any spectral lines corresponding to a mode spectrum. No significant change in measured pulse characteristics was observed, compared to conventional mode-locking in the unmodulated cavity. The insensitivity to the random phase modulation is expected, given the lack of phase-sensitive elements in the cavity.

preprint2016arXiv

Excursion probability of Gaussian random fields on sphere

Let $X=\{X(x): x\in\mathbb{S}^N\}$ be a real-valued, centered Gaussian random field indexed on the $N$-dimensional unit sphere $\mathbb{S}^N$. Approximations to the excursion probability ${\mathbb{P}}\{\sup_{x\in\mathbb{S}^N}X(x)\ge u\}$, as $u\to\infty$, are obtained for two cases: (i) $X$ is locally isotropic and its sample functions are non-smooth and; (ii) $X$ is isotropic and its sample functions are twice differentiable. For case (i), the excursion probability can be studied by applying the results in Piterbarg (Asymptotic Methods in the Theory of Gaussian Processes and Fields (1996) Amer. Math. Soc.), Mikhaleva and Piterbarg (Theory Probab. Appl. 41 (1997) 367--379) and Chan and Lai (Ann. Probab. 34 (2006) 80--121). It is shown that the asymptotics of ${\mathbb{P}}\{\sup_{x\in\mathbb {S}^N}X(x)\ge u\}$ is similar to Pickands' approximation on the Euclidean space which involves Pickands' constant. For case (ii), we apply the expected Euler characteristic method to obtain a more precise approximation such that the error is super-exponentially small.

preprint2016arXiv

Expected Number and Height Distribution of Critical Points of Smooth Isotropic Gaussian Random Fields

We obtain formulae for the expected number and height distribution of critical points of smooth isotropic Gaussian random fields parameterized on Euclidean space or spheres of arbitrary dimension. The results hold in general in the sense that there are no restrictions on the covariance function of the field except for smoothness and isotropy. The results are based on a characterization of the distribution of the Hessian of the Gaussian field by means of the family of Gaussian orthogonally invariant (GOI) matrices, of which the Gaussian orthogonal ensemble (GOE) is a special case. The obtained formulae depend on the covariance function only through a single parameter (Euclidean space) or two parameters (sphere), and include the special boundary case of random Laplacian eigenfunctions.

preprint2016arXiv

The mean Euler characteristic and excursion probability of Gaussian random fields with stationary increments

Let $X=\{X(t),t\in {\mathbb{R}}^N\}$ be a centered Gaussian random field with stationary increments and $X(0)=0$. For any compact rectangle $T\subset {\mathbb{R}}^N$ and $u\in {\mathbb{R}}$, denote by $A_u=\{t\in T:X(t)\geq u\}$ the excursion set. Under $X(\cdot)\in C^2({\mathbb{R}}^N)$ and certain regularity conditions, the mean Euler characteristic of $A_u$, denoted by ${\mathbb{E}}\{φ(A_u)\}$, is derived. By applying the Rice method, it is shown that, as $u\to\infty$, the excursion probability ${\mathbb{P}}\{\sup_{t\in T}X(t)\geq u\}$ can be approximated by ${\mathbb{E}}\{φ(A_u)\}$ such that the error is exponentially smaller than ${\mathbb{E}}\{φ(A_u)\}$. This verifies the expected Euler characteristic heuristic for a large class of Gaussian random fields with stationary increments.

preprint2015arXiv

Excursion Probabilities of Isotropic and Locally Isotropic Gaussian Random Fields on Manifolds

Let $X= \{X(p), p\in M\}$ be a centered Gaussian random field, where $M$ is a smooth Riemannian manifold. For a suitable compact subset $D\subset M$, we obtain the approximations to excursion probability $\mathbb{P}\{\sup_{p\in D} X(p) \ge u \}$, as $u\to \infty$, for two cases: (i) $X$ is smooth and isotropic; (ii) $X$ is non-smooth and locally isotropic. For case (i), the expected Euler characteristic approximation is formulated explicitly; while for case (ii), it is shown that the asymptotics is similar to Pickands' approximation on Euclidean space which involves Pickands' constant and the volume of $D$. These extend the results in \citep{Cheng:2014} from sphere to general Riemannian manifolds.

preprint2015arXiv

Excursion Probability of Certain Non-centered Smooth Gaussian Random Fields

Let $X = \{X(t): t\in T \}$ be a non-centered, unit-variance, smooth Gaussian random field indexed on some parameter space $T$, and let $A_u(X,T) = \{t\in T: X(t)\geq u\}$ be the excursion set of $X$ exceeding level $u$. Under certain smoothness and regularity conditions, it is shown that, as $u\to \infty$, the excursion probability $\mathbb{P}\{\sup_{t\in T} X(t)\ge u \}$ can be approximated by the expected Euler characteristic of $A_u(X,T)$, denoted by $\mathbb{E}\{χ(A_u(X,T))\}$, such that the error is super-exponentially small. This verifies the expected Euler characteristic heuristic for a large class of non-centered smooth Gaussian random fields and provides a much more accurate approximation compared with those existing results by the double sum method. The explicit formulae for $\mathbb{E}\{χ(A_u(X,T))\}$ are also derived for two cases: (i) $T$ is a rectangle and $X-\mathbb{E} X$ is stationary; (ii) $T$ is an $N$-dimensional sphere and $X-\mathbb{E} X$ is isotropic.

preprint2014arXiv

Distribution of the Height of Local Maxima of Gaussian Random Fields

Let $\{f(t): t\in T\}$ be a smooth Gaussian random field over a parameter space $T$, where $T$ may be a subset of Euclidean space or, more generally, a Riemannian manifold. For any local maximum of $f(t)$ located at $t_0$ in the interior of $T$, we provide general formulae and asymptotic approximations for both the tail distribution of the height of a local maximum $\mathbb{P}\{f(t_0)>u | t_0 \text{is a local maximum of} f(t) \}$ and the overshoot distribution of a local maximum $\mathbb{P}\{f(t_0)>u+v | t_0 \text{is a local maximum of} f(t) \text{and} f(t_0)>v\}$. Assuming further that $f$ is isotropic, we apply techniques from random matrix theory related to the Gaussian orthogonal ensemble to compute such conditional probabilities explicitly when $T$ is Euclidean or a sphere of arbitrary dimension. Such calculations are motivated by the statistical problem of detecting peaks in the presence of smooth Gaussian noise.

preprint2014arXiv

Multiple Testing of Local Maxima for Detection of Peaks in Random Fields

A topological multiple testing scheme is presented for detecting peaks in images under stationary ergodic Gaussian noise, where tests are performed at local maxima of the smoothed observed signals. The procedure generalizes the one-dimensional scheme of Schwartzman et al. (2011) to Euclidean domains of arbitrary dimension. Two methods are developed according to two different ways of computing p-values: (i) using the exact distribution of the height of local maxima (Cheng and Schwartzman, 2014), available explicitly when the noise field is isotropic; (ii) using an approximation to the overshoot distribution of local maxima above a pre-threshold (Cheng and Schwartzman, 2014), applicable when the exact distribution is unknown, such as when the stationary noise field is non-isotropic. The algorithms, combined with the Benjamini-Hochberg procedure for thresholding p-values, provide asymptotic strong control of the False Discovery Rate (FDR) and power consistency, with specific rates, as the search space and signal strength get large. The optimal smoothing bandwidth and optimal pre-threshold are obtained to achieve maximum power. Simulations show that FDR levels are maintained in non-asymptotic conditions. The methods are illustrated in a nanoscopy image analysis problem of detecting fluorescent molecules against the image background.