-
Spectral Methods for Matrix Product Factorization
Authors:
Saieed Akbari,
Yi-Zheng Fan,
Fu-Tao Hu,
Babak Miraftab,
Yi Wang
Abstract:
A graph $G$ is factored into graphs $H$ and $K$ via a matrix product if there exist adjacency matrices $A$, $B$, and $C$ of $G$, $H$, and $K$, respectively, such that $A = BC$. In this paper, we study the spectral aspects of the matrix product of graphs, including regularity, bipartiteness, and connectivity. We show that if a graph $G$ is factored into a connected graph $H$ and a graph $K$ with no…
▽ More
A graph $G$ is factored into graphs $H$ and $K$ via a matrix product if there exist adjacency matrices $A$, $B$, and $C$ of $G$, $H$, and $K$, respectively, such that $A = BC$. In this paper, we study the spectral aspects of the matrix product of graphs, including regularity, bipartiteness, and connectivity. We show that if a graph $G$ is factored into a connected graph $H$ and a graph $K$ with no isolated vertices, then certain properties hold. If $H$ is non-bipartite, then $G$ is connected. If $H$ is bipartite and $G$ is not connected, then $K$ is a regular bipartite graph, and consequently, $n$ is even. Furthermore, we show that trees are not factorizable, which answers a question posed by Maghsoudi et al.
△ Less
Submitted 4 July, 2024;
originally announced July 2024.
-
The largest Laplacian eigenvalue and the balancedness of simplicial complexes
Authors:
Yi-Zheng Fan,
Hui-Feng Wu,
Yi Wang
Abstract:
Let $K$ be a simplical complex, and let $\mathcal{L}_i^{up}(K), \mathcal{Q}_i^{up}(K)$ be the $i$-th up Laplacian and signless Laplacian of $K$, respectively. In this paper we proved that the largest eigenvalue of $\mathcal{L}_i^{up}(K)$ is not greater than the largest eigenvalue of $\mathcal{Q}_i^{up}(K)$; furthermore, if $K$ is $(i+1)$-path connected, then the equality holds if and only if the…
▽ More
Let $K$ be a simplical complex, and let $\mathcal{L}_i^{up}(K), \mathcal{Q}_i^{up}(K)$ be the $i$-th up Laplacian and signless Laplacian of $K$, respectively. In this paper we proved that the largest eigenvalue of $\mathcal{L}_i^{up}(K)$ is not greater than the largest eigenvalue of $\mathcal{Q}_i^{up}(K)$; furthermore, if $K$ is $(i+1)$-path connected, then the equality holds if and only if the $i$-th incidence signed graph $B_i(K)$ of $K$ is balanced. As an application we provided an upper bound for the largest eigenvalue of the $i$-th up Laplacian of $K$, which improves the bound given by Horak and Jost and generalizes the result of Anderson and Morley on graphs.We characterized the balancedness of simplicial complexes under operations such as wedge sum, join, Cartesian product and duplication of motifs. For each $i \ge 0$, by using wedge sum or duplication of motifs, we can construct an infinitely many $(i+1)$-path connected simplicial complexes $K$ with $B_i(K)$ being balanced.
△ Less
Submitted 29 May, 2024;
originally announced May 2024.
-
Zero-one Grothendieck Polynomials
Authors:
Yiming Chen,
Neil J. Y. Fan,
Zelin Ye
Abstract:
Fink, Mészáros and St.Dizier showed that the Schubert polynomial $\mathfrak{S}_w(x)$ is zero-one if and only if $w$ avoids twelve permutation patterns. In this paper, we prove that the Grothendieck polynomial $\mathfrak{G}_w(x)$ is zero-one, i.e., with coefficients either 0 or $\pm$1, if and only if $w$ avoids six patterns. As applications, we show that zero-one homogeneous Grothendieck polynomial…
▽ More
Fink, Mészáros and St.Dizier showed that the Schubert polynomial $\mathfrak{S}_w(x)$ is zero-one if and only if $w$ avoids twelve permutation patterns. In this paper, we prove that the Grothendieck polynomial $\mathfrak{G}_w(x)$ is zero-one, i.e., with coefficients either 0 or $\pm$1, if and only if $w$ avoids six patterns. As applications, we show that zero-one homogeneous Grothendieck polynomials are Lorentzian, partially confirming three conjectures of Huh, Matherne, Mészáros and St.Dizier. Moreover, we verify several conjectures on the support and coefficients of Grothendieck polynomials posed by Mészáros, Setiabrata and St.Dizier for the case of zero-one Grothendieck polynomials.
△ Less
Submitted 8 May, 2024;
originally announced May 2024.
-
Bivariate first-order random coefficient integer-valued autoregressive processes based on modified negative binomial operator
Authors:
Yixuan Fan,
Dehui Wang
Abstract:
In this paper, a new bivariate random coefficient integer-valued autoregressive process based on modified negative binomial operator with dependent innovations is proposed. Basic probabilistic and statistical properties of this model are derived. To estimate unknown parameters, Yule-Walker, conditional least squares and conditional maximum likelihood methods are considered and evaluated by Monte C…
▽ More
In this paper, a new bivariate random coefficient integer-valued autoregressive process based on modified negative binomial operator with dependent innovations is proposed. Basic probabilistic and statistical properties of this model are derived. To estimate unknown parameters, Yule-Walker, conditional least squares and conditional maximum likelihood methods are considered and evaluated by Monte Carlo simulations. Asymptotic properties of the estimators are derived. Moreover, coherent forecasting and possible extension of the proposed model is provided. Finally, the proposed model is applied to the monthly crime datasets and compared with other models.
△ Less
Submitted 27 April, 2024;
originally announced April 2024.
-
Height Filtrations and Base Loci on Flag Bundles over a Curve
Authors:
Yangyu Fan,
Wenbin Luo,
Binggang Qu
Abstract:
Let $k$ be an algebraically closed field of characteristic zero. Let $C/k$ be a projective smooth curve with function field $K=k(C)$ and $G/k$ be a connected reductive group. Let $F$ be a principal $G$-bundle on $C$. Let $P \subseteq G$ be a parabolic subgroup and $λ: P \longrightarrow G$ be a strictly anti-dominant character. Then $F/P \longrightarrow C$ is a flag bundle and…
▽ More
Let $k$ be an algebraically closed field of characteristic zero. Let $C/k$ be a projective smooth curve with function field $K=k(C)$ and $G/k$ be a connected reductive group. Let $F$ be a principal $G$-bundle on $C$. Let $P \subseteq G$ be a parabolic subgroup and $λ: P \longrightarrow G$ be a strictly anti-dominant character. Then $F/P \longrightarrow C$ is a flag bundle and $\mathcal{L}_λ=F \times_P k_λ$ on $F/P$ is a relatively ample line bundle.
Let $h_{\mathcal{L}_λ}$ be the height function on $X=(F/P)_K$ induced by $\mathcal{L}_λ$. We compute its height filtration and successive minima: the height filtration is given by a Bruhat decomposition and successive minima are given by slopes.
An interesting application is that, in this case, Zhang's inequality on successive minima can be enhanced into an equality. This is proved from the convex geometry point of view.
Let $f \in N^1(F/P)$ be the numerical class of a vertical fiber. We also compute the augmented base loci $\mathrm{B}_+(\mathcal{L}_λ-tf)$ for any $t \in \mathbb{R}$ and it turns out that they are almost the same as the height filtration, which accords with the philosophy in Arakelov geometry that the positivity of the line bundle $\mathcal{L}_λ$ is equivalent to the positivity of its induced height function $h_{\mathcal{L}_λ}$. As a corollary, we compute the $k$-th movable cones of flag bundles over curves for all $k$.
△ Less
Submitted 11 March, 2024;
originally announced March 2024.
-
Spectral bipartite Turan problems on linear hypergraphs
Authors:
Chuan-Ming She,
Yi-Zheng Fan,
Liying Kang
Abstract:
Let $F$ be a graph and let $\mathcal{B}_r(F)$ be the class of $r$-uniform Berge-$F$ hypergraphs. In this paper, by establishing a relationship between the spectral radius of the adjacency tensor of a uniform hypergraph and its local structure via walks, we give a spectral asymptotic bound for $\mathcal{B}_{r}(C_3)$-free linear $r$-uniform hypergraphs and upper bounds for the spectral radii of…
▽ More
Let $F$ be a graph and let $\mathcal{B}_r(F)$ be the class of $r$-uniform Berge-$F$ hypergraphs. In this paper, by establishing a relationship between the spectral radius of the adjacency tensor of a uniform hypergraph and its local structure via walks, we give a spectral asymptotic bound for $\mathcal{B}_{r}(C_3)$-free linear $r$-uniform hypergraphs and upper bounds for the spectral radii of $\mathcal{B}_{r}(K_{2,t})$-free or $\{\mathcal{B}_{r}(K_{s,t}),\mathcal{B}_{r}(C_{3})\}$-free linear $r$-uniform hypergraphs, where $C_{3}$ and $K_{s,t}$ are respectively the triangle and the complete bipartite graph with one part having $s$ vertices and the other part having $t$ vertices. Our work implies an upper bound for the number of edges of $\{\mathcal{B}_{r}(K_{s,t}),\mathcal{B}_{r}(C_{3})\}$-free linear $r$-uniform hypergraphs, and extends some known work on (spectral) extreme problems of hypergraphs.
△ Less
Submitted 4 March, 2024;
originally announced March 2024.
-
Analysis of the particle relaxation method for generating uniform particle distributions in smoothed particle hydrodynamics
Authors:
Yu Fan,
Xiaoliang Li,
Shuoguo Zhang,
Xiangyu Hu,
Nikolaus A. Adams
Abstract:
We establish a theoretical framework of the particle relaxation method for uniform particle generation of Smoothed Particle Hydrodynamics. We achieve this by reformulating the particle relaxation as an optimization problem. The objective function is an integral difference between discrete particle-based and smoothed-analytical volume fractions. The analysis demonstrates that the particle relaxatio…
▽ More
We establish a theoretical framework of the particle relaxation method for uniform particle generation of Smoothed Particle Hydrodynamics. We achieve this by reformulating the particle relaxation as an optimization problem. The objective function is an integral difference between discrete particle-based and smoothed-analytical volume fractions. The analysis demonstrates that the particle relaxation method in the domain interior is essentially equivalent to employing a gradient descent approach to solve this optimization problem, and we can extend such an equivalence to the bounded domain by introducing a proper boundary term. Additionally, each periodic particle distribution has a spatially uniform particle volume, denoted as characteristic volume. The relaxed particle distribution has the largest characteristic volume, and the kernel cut-off radius determines this volume. This insight enables us to control the relaxed particle distribution by selecting the target kernel cut-off radius for a given kernel function.
△ Less
Submitted 1 March, 2024;
originally announced March 2024.
-
A Pieri type formula for motivic Chern classes of Schubert cells in Grassmannians
Authors:
Neil J. Y. Fan,
Peter L. Guo,
Changjian Su,
Rui Xiong
Abstract:
We prove a Pieri formula for motivic Chern classes of Schubert cells in the equivariant K-theory of Grassmannians, which is described in terms of ribbon operators on partitions. Our approach is to transform the Schubert calculus over Grassmannians to the calculation in a certain affine Hecke algebra. As a consequence, we derive a Pieri formula for Segre motivic classes of Schubert cells in Grassma…
▽ More
We prove a Pieri formula for motivic Chern classes of Schubert cells in the equivariant K-theory of Grassmannians, which is described in terms of ribbon operators on partitions. Our approach is to transform the Schubert calculus over Grassmannians to the calculation in a certain affine Hecke algebra. As a consequence, we derive a Pieri formula for Segre motivic classes of Schubert cells in Grassmannians. We apply the Pieri formulas to establish a relation between motivic Chern classes and Segre motivic classes, extending a well-known relation between the classes of structure sheaves and ideal sheaves. As another application, we find a symmetric power series representative for the class of the dualizing sheaf of a Schubert variety.
△ Less
Submitted 6 February, 2024;
originally announced February 2024.
-
Hitting probability for Reflected Brownian Motion at Small Target
Authors:
Yuchen Fan
Abstract:
We derive the asymptotic behavior of hitting probability at small target of size $O(ε)$ for reflected Brownian motion in domains with suitable smooth boundary conditions, where the boundary of domain contains both reflecting part, absorbing part and target. In this case the domain could be localized near the target and explicit computations are possible. The asymptotic behavior is only related to…
▽ More
We derive the asymptotic behavior of hitting probability at small target of size $O(ε)$ for reflected Brownian motion in domains with suitable smooth boundary conditions, where the boundary of domain contains both reflecting part, absorbing part and target. In this case the domain could be localized near the target and explicit computations are possible. The asymptotic behavior is only related to $ε$ up to some multiplicative constants that depends on the domain and boundary conditions.
△ Less
Submitted 1 February, 2024;
originally announced February 2024.
-
Automatic Implementation of Neural Networks through Reaction Networks--Part II: Error Analysis
Authors:
Yuzhen Fan,
Xiaoyu Zhang,
Chuanhou Gao,
Denis Dochain
Abstract:
This paired article aims to develop an automated and programmable biochemical fully connected neural network (BFCNN) with solid theoretical support. In Part I, a concrete design for BFCNN is presented, along with the validation of the effectiveness and exponential convergence of computational modules. In this article, we establish the framework for specifying the realization errors by monitoring t…
▽ More
This paired article aims to develop an automated and programmable biochemical fully connected neural network (BFCNN) with solid theoretical support. In Part I, a concrete design for BFCNN is presented, along with the validation of the effectiveness and exponential convergence of computational modules. In this article, we establish the framework for specifying the realization errors by monitoring the errors generated from approaching equilibrium points in individual modules, as well as their vertical propagation from upstream modules and horizontal accumulation from previous iterations. Ultimately, we derive the general error upper bound formula for any iteration and illustrate its exponential convergence order with respect to the phase length of the utilized chemical oscillator. The numerical experiments, based on the classification examples, reveal the tendency of total errors related to both the phase length and iteration number.
△ Less
Submitted 13 January, 2024;
originally announced January 2024.
-
The spectra of Laplace operators on covering simplicial complexes
Authors:
Yi-Zheng Fan,
Yi-Min Song,
Yi Wang
Abstract:
In this paper, by the representation theory of symmetric group, we give a decomposition of the Laplace operator (in matrix form) of a covering simplicial complex into the direct sum of some matrices, including the Laplace operator of the underlying simplicial complex. So, the spectrum of a covering simplicial complex can be expressed into a union of the spectrum of the underlying simplicial comple…
▽ More
In this paper, by the representation theory of symmetric group, we give a decomposition of the Laplace operator (in matrix form) of a covering simplicial complex into the direct sum of some matrices, including the Laplace operator of the underlying simplicial complex. So, the spectrum of a covering simplicial complex can be expressed into a union of the spectrum of the underlying simplicial complex and the spectra of some other matrices, which implies a result of Horak and Jost. In particular, we show that the spectrum of a $2$-fold covering simplicial complex is the union the that of the underlying simplicial complex and that of an incidence-signed simplicial complex, which is an analog of Bilu and Linial's result on graphs. Finally we show that the dimension of the cohomology group of a covering simplicial complex is greater than or equal to that of the underlying simplicial complex.
△ Less
Submitted 19 December, 2023;
originally announced December 2023.
-
Harnessing the Power of Neural Operators with Automatically Encoded Conservation Laws
Authors:
Ning Liu,
Yiming Fan,
Xianyi Zeng,
Milan Klöwer,
Lu Zhang,
Yue Yu
Abstract:
Neural operators (NOs) have emerged as effective tools for modeling complex physical systems in scientific machine learning. In NOs, a central characteristic is to learn the governing physical laws directly from data. In contrast to other machine learning applications, partial knowledge is often known a priori about the physical system at hand whereby quantities such as mass, energy and momentum a…
▽ More
Neural operators (NOs) have emerged as effective tools for modeling complex physical systems in scientific machine learning. In NOs, a central characteristic is to learn the governing physical laws directly from data. In contrast to other machine learning applications, partial knowledge is often known a priori about the physical system at hand whereby quantities such as mass, energy and momentum are exactly conserved. Currently, NOs have to learn these conservation laws from data and can only approximately satisfy them due to finite training data and random noise. In this work, we introduce conservation law-encoded neural operators (clawNOs), a suite of NOs that endow inference with automatic satisfaction of such conservation laws. ClawNOs are built with a divergence-free prediction of the solution field, with which the continuity equation is automatically guaranteed. As a consequence, clawNOs are compliant with the most fundamental and ubiquitous conservation laws essential for correct physical consistency. As demonstrations, we consider a wide variety of scientific applications ranging from constitutive modeling of material deformation, incompressible fluid dynamics, to atmospheric simulation. ClawNOs significantly outperform the state-of-the-art NOs in learning efficacy, especially in small-data regimes.
△ Less
Submitted 4 June, 2024; v1 submitted 18 December, 2023;
originally announced December 2023.
-
Automatic Implementation of Neural Networks through Reaction Networks -- Part I: Circuit Design and Convergence Analysis
Authors:
Yuzhen Fan,
Xiaoyu Zhang,
Chuanhou Gao,
Denis Dochain
Abstract:
Information processing relying on biochemical interactions in the cellular environment is essential for biological organisms. The implementation of molecular computational systems holds significant interest and potential in the fields of synthetic biology and molecular computation. This two-part article aims to introduce a programmable biochemical reaction network (BCRN) system endowed with mass a…
▽ More
Information processing relying on biochemical interactions in the cellular environment is essential for biological organisms. The implementation of molecular computational systems holds significant interest and potential in the fields of synthetic biology and molecular computation. This two-part article aims to introduce a programmable biochemical reaction network (BCRN) system endowed with mass action kinetics that realizes the fully connected neural network (FCNN) and has the potential to act automatically in vivo. In part I, the feedforward propagation computation, the backpropagation component, and all bridging processes of FCNN are ingeniously designed as specific BCRN modules based on their dynamics. This approach addresses a design gap in the biochemical assignment module and judgment termination module and provides a novel precise and robust realization of bi-molecular reactions for the learning process. Through equilibrium approaching, we demonstrate that the designed BCRN system achieves FCNN functionality with exponential convergence to target computational results, thereby enhancing the theoretical support for such work. Finally, the performance of this construction is further evaluated on two typical logic classification problems.
△ Less
Submitted 30 November, 2023;
originally announced November 2023.
-
On the adaptation of causal forests to manifold data
Authors:
Yiyi Huo,
Yingying Fan,
Fang Han
Abstract:
Researchers often hold the belief that random forests are "the cure to the world's ills" (Bickel, 2010). But how exactly do they achieve this? Focused on the recently introduced causal forests (Athey and Imbens, 2016; Wager and Athey, 2018), this manuscript aims to contribute to an ongoing research trend towards answering this question, proving that causal forests can adapt to the unknown covariat…
▽ More
Researchers often hold the belief that random forests are "the cure to the world's ills" (Bickel, 2010). But how exactly do they achieve this? Focused on the recently introduced causal forests (Athey and Imbens, 2016; Wager and Athey, 2018), this manuscript aims to contribute to an ongoing research trend towards answering this question, proving that causal forests can adapt to the unknown covariate manifold structure. In particular, our analysis shows that a causal forest estimator can achieve the optimal rate of convergence for estimating the conditional average treatment effect, with the covariate dimension automatically replaced by the manifold dimension. These findings align with analogous observations in the realm of deep learning and resonate with the insights presented in Peter Bickel's 2004 Rietz lecture.
△ Less
Submitted 26 December, 2023; v1 submitted 27 November, 2023;
originally announced November 2023.
-
Dimensionality reduction of networked systems with separable coupling-dynamics: theory and applications
Authors:
Chengyi Tu,
Ying Fan,
Tianyu Shi
Abstract:
Complex dynamical systems are prevalent in various domains, but their analysis and prediction are hindered by their high dimensionality and nonlinearity. Dimensionality reduction techniques can simplify the system dynamics by reducing the number of variables, but most existing methods do not account for networked systems with separable coupling-dynamics, where the interaction between nodes can be…
▽ More
Complex dynamical systems are prevalent in various domains, but their analysis and prediction are hindered by their high dimensionality and nonlinearity. Dimensionality reduction techniques can simplify the system dynamics by reducing the number of variables, but most existing methods do not account for networked systems with separable coupling-dynamics, where the interaction between nodes can be decomposed into a function of the node state and a function of the neighbor state. Here, we present a novel dimensionality reduction framework that can effectively capture the global dynamics of these networks by projecting them onto a low-dimensional system. We derive the reduced system's equation and stability conditions, and propose an error metric to quantify the reduction accuracy. We demonstrate our framework on two examples of networked systems with separable coupling-dynamics: a modified susceptible-infected-susceptible model with direct infection and a modified Michaelis-Menten model with activation and inhibition. We conduct numerical experiments on synthetic and empirical networks to validate and evaluate our framework, and find a good agreement between the original and reduced systems. We also investigate the effects of different network structures and parameters on the system dynamics and the reduction error. Our framework offers a general and powerful tool for studying complex dynamical networks with separable coupling-dynamics.
△ Less
Submitted 27 November, 2023;
originally announced November 2023.
-
Perfect codes in 2-valent Cayley digraphs on abelian groups
Authors:
Shilong Yu,
Yuefeng Yang,
Yushuang Fan,
Xuanlong Ma
Abstract:
For a digraph $Γ$, a subset $C$ of $V(Γ)$ is a perfect code if $C$ is a dominating set such that every vertex of $Γ$ is dominated by exactly one vertex in $C$. In this paper, we classify strongly connected 2-valent Cayley digraphs on abelian groups admitting a perfect code, and determine completely all perfect codes of such digraphs.
For a digraph $Γ$, a subset $C$ of $V(Γ)$ is a perfect code if $C$ is a dominating set such that every vertex of $Γ$ is dominated by exactly one vertex in $C$. In this paper, we classify strongly connected 2-valent Cayley digraphs on abelian groups admitting a perfect code, and determine completely all perfect codes of such digraphs.
△ Less
Submitted 16 June, 2024; v1 submitted 29 October, 2023;
originally announced October 2023.
-
Ichino period for CM forms
Authors:
Li Cai,
Yangyu Fan,
Yong Jiang
Abstract:
In both local and global settings, we establish explicit relations between Ichino triple product period and Waldspurger toric periods for CM forms via the theta lifting and the see-saw principle.
In both local and global settings, we establish explicit relations between Ichino triple product period and Waldspurger toric periods for CM forms via the theta lifting and the see-saw principle.
△ Less
Submitted 4 October, 2023;
originally announced October 2023.
-
Hypergraph coverings and Ramanujan Hypergraphs
Authors:
Yi-Min Song,
Yi-Zheng Fan,
Zhengke Miao
Abstract:
In this paper we investigate Ramanujan hypergraphs by using hypergraph coverings. We first show that the spectrum of a $k$-fold covering $\bar{H}$ of a connected hypergraph $H$ contains the spectrum of $H$, and that it is the union of the spectrum of $H$ and the spectrum of an incidence-signed hypergraph with $H$ as underlying hypergraph if $k=2$, which generalizes Bilu-Linial result on graph cove…
▽ More
In this paper we investigate Ramanujan hypergraphs by using hypergraph coverings. We first show that the spectrum of a $k$-fold covering $\bar{H}$ of a connected hypergraph $H$ contains the spectrum of $H$, and that it is the union of the spectrum of $H$ and the spectrum of an incidence-signed hypergraph with $H$ as underlying hypergraph if $k=2$, which generalizes Bilu-Linial result on graph coverings. We give a lower bound for the second largest eigenvalue of a $d$-regular hypergraph by universal cover, which generalizes Alon-Boppana bound on $d$-regular graphs and Feng-Li bound on $(d,r)$-regular hypergraphs. By using interlacing family of polynomials, we prove that every $(d,r)$-regular hypergraph has a right-sided Ramanujan $2$-covering, and has a left-sided Ramanujan $2$-covering if the roots of the matching polynomial of its incident graph satisfy some condition. By Ramanujan $2$-coverings, we prove the existence of some families of infinite many left-sided or right-sided $(d,r)$-regular Ramanujan hypergraphs under certain conditions on $d$ and $r$.
△ Less
Submitted 4 March, 2024; v1 submitted 2 October, 2023;
originally announced October 2023.
-
SOFARI: High-Dimensional Manifold-Based Inference
Authors:
Zemin Zheng,
Xin Zhou,
Yingying Fan,
Jinchi Lv
Abstract:
Multi-task learning is a widely used technique for harnessing information from various tasks. Recently, the sparse orthogonal factor regression (SOFAR) framework, based on the sparse singular value decomposition (SVD) within the coefficient matrix, was introduced for interpretable multi-task learning, enabling the discovery of meaningful latent feature-response association networks across differen…
▽ More
Multi-task learning is a widely used technique for harnessing information from various tasks. Recently, the sparse orthogonal factor regression (SOFAR) framework, based on the sparse singular value decomposition (SVD) within the coefficient matrix, was introduced for interpretable multi-task learning, enabling the discovery of meaningful latent feature-response association networks across different layers. However, conducting precise inference on the latent factor matrices has remained challenging due to orthogonality constraints inherited from the sparse SVD constraint. In this paper, we suggest a novel approach called high-dimensional manifold-based SOFAR inference (SOFARI), drawing on the Neyman near-orthogonality inference while incorporating the Stiefel manifold structure imposed by the SVD constraints. By leveraging the underlying Stiefel manifold structure, SOFARI provides bias-corrected estimators for both latent left factor vectors and singular values, for which we show to enjoy the asymptotic mean-zero normal distributions with estimable variances. We introduce two SOFARI variants to handle strongly and weakly orthogonal latent factors, where the latter covers a broader range of applications. We illustrate the effectiveness of SOFARI and justify our theoretical results through simulation examples and a real data application in economic forecasting.
△ Less
Submitted 26 September, 2023;
originally announced September 2023.
-
Bumpless pipe dreams meet Puzzles
Authors:
Neil J. Y. Fan,
Peter L. Guo,
Rui Xiong
Abstract:
Knutson and Zinn-Justin recently found a puzzle rule for the expansion of the product $\mathfrak{G}_{u}(x,t)\cdot \mathfrak{G}_{v}(x,t)$ of two double Grothendieck polynomials indexed by permutations with separated descents. We establish its triple Schubert calculus version in the sense of Knutson and Tao, namely, a formula for expanding $\mathfrak{G}_{u}(x,y)\cdot \mathfrak{G}_{v}(x,t)$ in differ…
▽ More
Knutson and Zinn-Justin recently found a puzzle rule for the expansion of the product $\mathfrak{G}_{u}(x,t)\cdot \mathfrak{G}_{v}(x,t)$ of two double Grothendieck polynomials indexed by permutations with separated descents. We establish its triple Schubert calculus version in the sense of Knutson and Tao, namely, a formula for expanding $\mathfrak{G}_{u}(x,y)\cdot \mathfrak{G}_{v}(x,t)$ in different secondary variables. Our rule is formulated in terms of pipe puzzles, incorporating both the structures of bumpless pipe dreams and classical puzzles. As direct applications, we recover the separated-descent puzzle formula by Knutson and Zinn-Justin (by setting $y=t$) and the bumpless pipe dream model of double Grothendieck polynomials by Weigandt (by setting $v=\operatorname{id}$ and $x=t$). Moreover, we utilize the formula to partially confirm a positivity conjecture of Kirillov about applying a skew operator to a Schubert polynomial.
△ Less
Submitted 1 September, 2023;
originally announced September 2023.
-
Fourier-Mukai numbers of K3 categories of very general special cubic fourfolds
Authors:
Yu-Wei Fan,
Kuan-Wen Lai
Abstract:
We give counting formulas for the number of Fourier-Mukai partners of the K3 category of a very general special cubic fourfold.
We give counting formulas for the number of Fourier-Mukai partners of the K3 category of a very general special cubic fourfold.
△ Less
Submitted 26 July, 2023;
originally announced July 2023.
-
Weak Galerkin methods for the Stokes eigenvalue problem
Authors:
Yunying Fan,
Qilong Zhai
Abstract:
In this paper, we rewrite the Stokes eigenvalue problem as an Elliptic eigenvalue problem restricted to subspace, and introduce an abstract framework of solving abstract elliptic eigenvalue problem to give the WG scheme, error estimates and asymptotic lower bounds. Besides, we introduce a new stabilizer and several inequalities to prove GLB properties. Some numerical examples are provided to valid…
▽ More
In this paper, we rewrite the Stokes eigenvalue problem as an Elliptic eigenvalue problem restricted to subspace, and introduce an abstract framework of solving abstract elliptic eigenvalue problem to give the WG scheme, error estimates and asymptotic lower bounds. Besides, we introduce a new stabilizer and several inequalities to prove GLB properties. Some numerical examples are provided to validate our theoretical analysis.
△ Less
Submitted 18 July, 2023;
originally announced July 2023.
-
ARK: Robust Knockoffs Inference with Coupling
Authors:
Yingying Fan,
Lan Gao,
Jinchi Lv
Abstract:
We investigate the robustness of the model-X knockoffs framework with respect to the misspecified or estimated feature distribution. We achieve such a goal by theoretically studying the feature selection performance of a practically implemented knockoffs algorithm, which we name as the approximate knockoffs (ARK) procedure, under the measures of the false discovery rate (FDR) and $k$-familywise er…
▽ More
We investigate the robustness of the model-X knockoffs framework with respect to the misspecified or estimated feature distribution. We achieve such a goal by theoretically studying the feature selection performance of a practically implemented knockoffs algorithm, which we name as the approximate knockoffs (ARK) procedure, under the measures of the false discovery rate (FDR) and $k$-familywise error rate ($k$-FWER). The approximate knockoffs procedure differs from the model-X knockoffs procedure only in that the former uses the misspecified or estimated feature distribution. A key technique in our theoretical analyses is to couple the approximate knockoffs procedure with the model-X knockoffs procedure so that random variables in these two procedures can be close in realizations. We prove that if such coupled model-X knockoffs procedure exists, the approximate knockoffs procedure can achieve the asymptotic FDR or $k$-FWER control at the target level. We showcase three specific constructions of such coupled model-X knockoff variables, verifying their existence and justifying the robustness of the model-X knockoffs framework. Additionally, we formally connect our concept of knockoff variable coupling to a type of Wasserstein distance.
△ Less
Submitted 4 June, 2024; v1 submitted 10 July, 2023;
originally announced July 2023.
-
Quantifying Distributional Model Risk in Marginal Problems via Optimal Transport
Authors:
Yanqin Fan,
Hyeonseok Park,
Gaoqian Xu
Abstract:
This paper studies distributional model risk in marginal problems, where each marginal measure is assumed to lie in a Wasserstein ball centered at a fixed reference measure with a given radius. Theoretically, we establish several fundamental results including strong duality, finiteness of the proposed Wasserstein distributional model risk, and the existence of an optimizer at each radius. In addit…
▽ More
This paper studies distributional model risk in marginal problems, where each marginal measure is assumed to lie in a Wasserstein ball centered at a fixed reference measure with a given radius. Theoretically, we establish several fundamental results including strong duality, finiteness of the proposed Wasserstein distributional model risk, and the existence of an optimizer at each radius. In addition, we show continuity of the Wasserstein distributional model risk as a function of the radius. Using strong duality, we extend the well-known Makarov bounds for the distribution function of the sum of two random variables with given marginals to Wasserstein distributionally robust Markarov bounds. Practically, we illustrate our results on four distinct applications when the sample information comes from multiple data sources and only some marginal reference measures are identified. They are: partial identification of treatment effects; externally valid treatment choice via robust welfare functions; Wasserstein distributionally robust estimation under data combination; and evaluation of the worst aggregate risk measures.
△ Less
Submitted 3 July, 2023;
originally announced July 2023.
-
Spectral extremal graphs for edge blow-up of star forests
Authors:
Jing Wang,
Zhenyu Ni,
Liying Kang,
Yi-zheng Fan
Abstract:
The edge blow-up of a graph $G$, denoted by $G^{p+1}$, is obtained by replacing each edge of $G$ with a clique of order $p+1$, where the new vertices of the cliques are all distinct. Yuan [J. Comb. Theory, Ser. B, 152 (2022) 379-398] determined the range of the Turán numbers for edge blow-up of all bipartite graphs and the exact Turán numbers for edge blow-up of all non-bipartite graphs. In this p…
▽ More
The edge blow-up of a graph $G$, denoted by $G^{p+1}$, is obtained by replacing each edge of $G$ with a clique of order $p+1$, where the new vertices of the cliques are all distinct. Yuan [J. Comb. Theory, Ser. B, 152 (2022) 379-398] determined the range of the Turán numbers for edge blow-up of all bipartite graphs and the exact Turán numbers for edge blow-up of all non-bipartite graphs. In this paper we prove that the graphs with the maximum spectral radius in an $n$-vertex graph without any copy of edge blow-up of star forests are the extremal graphs for edge blow-up of star forests when $n$ is sufficiently large.
△ Less
Submitted 29 June, 2023;
originally announced June 2023.
-
Shifting numbers of abelian varieties via bounded t-structures
Authors:
Yu-Wei Fan
Abstract:
The shifting numbers measure the asymptotic amount by which an endofunctor of a triangulated category translates inside the category, and are analogous to Poincare translation numbers that are widely used in dynamical systems. Motivated by this analogy, Fan-Filip raised the following question: ``Do the shifting numbers define a quasimorphism on the group of autoequivalences of a triangulated categ…
▽ More
The shifting numbers measure the asymptotic amount by which an endofunctor of a triangulated category translates inside the category, and are analogous to Poincare translation numbers that are widely used in dynamical systems. Motivated by this analogy, Fan-Filip raised the following question: ``Do the shifting numbers define a quasimorphism on the group of autoequivalences of a triangulated category?" An affirmative answer was given by Fan-Filip for the bounded derived category of coherent sheaves on an elliptic curve or an abelian surface, via properties of the spaces of Bridgeland stability conditions on these categories.
We prove in this article that the question has an affirmative answer for abelian varieties of arbitrary dimensions, generalizing the result of Fan-Filip. One of the key steps is to establish an alternative definition of the shifting numbers via bounded t-structures on triangulated categories. In particular, the full package of a Bridgeland stability condition (a bounded t-structure, and a central charge on a charge lattice) is not necessary for the purpose of computing the shifting numbers.
△ Less
Submitted 20 July, 2023; v1 submitted 23 June, 2023;
originally announced June 2023.
-
Formation Control with Unknown Directions and General Coupling Coefficients
Authors:
Zhen Li,
Yang Tang,
Yongqing Fan,
Tingwen Huang
Abstract:
Generally, the normal displacement-based formation control has a sensing mode that requires the agent not only to have certain knowledge of its direction, but also to gather its local information characterized by nonnegative coupling coefficients. However, the direction may be unknown in the sensing processes, and the coupling coefficients may also involve negative ones due to some circumstances.…
▽ More
Generally, the normal displacement-based formation control has a sensing mode that requires the agent not only to have certain knowledge of its direction, but also to gather its local information characterized by nonnegative coupling coefficients. However, the direction may be unknown in the sensing processes, and the coupling coefficients may also involve negative ones due to some circumstances. This paper introduces these phenomena into a class of displacement-based formation control problem. Then, a geometric approach have been employed to overcome the difficulty of analysis on the introduced phenomena. The purpose of this approach is to construct some convex polytopes for containing the effects caused by the unknown direction, and to analyze the non-convexity by admitting the negative coupling coefficients in a certain range. Under the actions of these phenomena, the constructed polytopes are shown to be invariant in view of the contractive set method. It means that the convergence of formation shape can be guaranteed. Subsequently, an example is given to examine the applicability of derived result.
△ Less
Submitted 3 June, 2023;
originally announced June 2023.
-
A Shape-Newton Method for Free-boundary Problems Subject to The Bernoulli Boundary Condition
Authors:
Yiyun Fan,
John Billingham,
Kristoffer van der Zee
Abstract:
We develop a shape-Newton method for solving generic free-boundary problems where one of the free-boundary conditions is governed by the Bernoulli equation. The Newton-like scheme is developed by employing shape derivatives in the weak forms, which allows us to update the position of the free surface and the potential on the free boundary by solving a boundary-value problem at each iteration. To v…
▽ More
We develop a shape-Newton method for solving generic free-boundary problems where one of the free-boundary conditions is governed by the Bernoulli equation. The Newton-like scheme is developed by employing shape derivatives in the weak forms, which allows us to update the position of the free surface and the potential on the free boundary by solving a boundary-value problem at each iteration. To validate the effectiveness of the approach, we apply the scheme to solve a problem involving the flow over a submerged triangular obstacle.
△ Less
Submitted 23 May, 2023;
originally announced May 2023.
-
Multiplicities of Representations in Algebraic Families
Authors:
Li Cai,
Yangyu Fan
Abstract:
In this short notes, we consider multiplicities of representations in general algebraic families, especially the upper semi-continuity of homological multiplicities and the locally constancy of Euler-Poincare numbers. This generalizes the main result of Aizenbud-Sayag for unramified twisting families.
In this short notes, we consider multiplicities of representations in general algebraic families, especially the upper semi-continuity of homological multiplicities and the locally constancy of Euler-Poincare numbers. This generalizes the main result of Aizenbud-Sayag for unramified twisting families.
△ Less
Submitted 20 February, 2024; v1 submitted 19 May, 2023;
originally announced May 2023.
-
Spherical Characters in Families: the unitary Gan-Gross-Prasad case
Authors:
Li Cai,
Yangyu Fan
Abstract:
We consider the variation of spherical characters in families. We formulate conjectures for the rationality and meromorphic property of spherical characters. As an example, we establish these conjectures in the unitary Gan-Gross-Prasad case.
We consider the variation of spherical characters in families. We formulate conjectures for the rationality and meromorphic property of spherical characters. As an example, we establish these conjectures in the unitary Gan-Gross-Prasad case.
△ Less
Submitted 19 May, 2023;
originally announced May 2023.
-
Intermediate Service Facility Planning in a Stochastic and Competitive Market: Incorporating Agent-infrastructure Interactions over Networks
Authors:
Sina Baghali,
Julio Deride,
Yueyue Fan,
Zhaomiao Guo
Abstract:
This paper presents a network-based multi-agent optimization model for the strategic planning of service facilities in a stochastic and competitive market. We focus on the type of service facilities that are of intermediate nature, i.e., users may need to deviate from the shortest path to receive/provide services in between the users' planned origins and destinations. This problem has many applica…
▽ More
This paper presents a network-based multi-agent optimization model for the strategic planning of service facilities in a stochastic and competitive market. We focus on the type of service facilities that are of intermediate nature, i.e., users may need to deviate from the shortest path to receive/provide services in between the users' planned origins and destinations. This problem has many applications in emerging transportation mobility, including dynamic ride-sharing hub design and competitive facility location and allocation problems for alternative fuel vehicle refueling stations. The main contribution of this paper is establishing a new multi-agent optimization framework considering decentralized decision makings of facility investors and users over a transportation network and providing rigorous analyses of its mathematical properties, such as uniqueness and existence of system equilibrium. In addition, we develop an exact convex reformulation of the original multi-agent optimization problems to overcome computational challenges brought by non-convexity. Extensive analysis on case studies showed how the proposed model can capture the complex interaction between different stakeholders in an uncertain environment. Additionally, our model allowed quantifying the value of stochastic modeling and information availability by exploring stochastic metrics, including value of stochastic solution (VSS) and expected value of perfect information (EVPI), in a multi-agent framework.
△ Less
Submitted 2 April, 2023;
originally announced April 2023.
-
Higher Ext-groups in the triple product case
Authors:
Li Cai,
Yangyu Fan
Abstract:
In this short note, we compute higher extension groups for all irreducible representations and deduce the multiplicity formula for finite length representations in triple product case.
In this short note, we compute higher extension groups for all irreducible representations and deduce the multiplicity formula for finite length representations in triple product case.
△ Less
Submitted 10 March, 2023;
originally announced March 2023.
-
Nielsen realization problem for derived automorphisms of generic K3 surfaces
Authors:
Yu-Wei Fan,
Kuan-Wen Lai
Abstract:
We prove that all nontrivial finite subgroups of derived automorphisms of K3 surfaces of Picard number one have order two and give formulas for the numbers of their conjugacy classes. We also obtain a similar result for the subgroups which are finite up to shifts. This in turn shows that such a K3 surface admits an associated cubic fourfold if and only if it has a derived automorphism of order thr…
▽ More
We prove that all nontrivial finite subgroups of derived automorphisms of K3 surfaces of Picard number one have order two and give formulas for the numbers of their conjugacy classes. We also obtain a similar result for the subgroups which are finite up to shifts. This in turn shows that such a K3 surface admits an associated cubic fourfold if and only if it has a derived automorphism of order three up to shifts. These results are achieved by proving that such a subgroup fixes a Bridgeland stability condition up to $\mathbb{C}$-actions. We also establish similar existence results for curves, twisted abelian surfaces, generic twisted K3 surfaces, and standard autoequivalences on surfaces.
△ Less
Submitted 16 May, 2023; v1 submitted 24 February, 2023;
originally announced February 2023.
-
Testing mean and variance by e-processes
Authors:
Yixuan Fan,
Zhanyi Jiao,
Ruodu Wang
Abstract:
We address the problem of testing conditional mean and conditional variance for non-stationary data. We build e-values and p-values for four types of non-parametric composite hypotheses with specified mean and variance as well as other conditions on the shape of the data-generating distribution. These shape conditions include symmetry, unimodality, and their combination. Using the obtained e-value…
▽ More
We address the problem of testing conditional mean and conditional variance for non-stationary data. We build e-values and p-values for four types of non-parametric composite hypotheses with specified mean and variance as well as other conditions on the shape of the data-generating distribution. These shape conditions include symmetry, unimodality, and their combination. Using the obtained e-values and p-values, we construct tests via e-processes, also known as testing by betting, as well as some tests based on combining p-values for comparison. Although we mainly focus on one-sided tests, the two-sided test for the mean is also studied. Simulation and empirical studies are conducted under a few settings, and they illustrate features of the methods based on e-processes.
△ Less
Submitted 25 June, 2024; v1 submitted 29 January, 2023;
originally announced January 2023.
-
Virtual representations for closed graph manifolds and Seifert geometry
Authors:
Yao Fan
Abstract:
In this paper, we mainly discuss the representations of closed graph manifolds to the Seifert motion group. Then we prove that there exist graph manifolds virtually having no faithful representations to the Seifert motion group.
In this paper, we mainly discuss the representations of closed graph manifolds to the Seifert motion group. Then we prove that there exist graph manifolds virtually having no faithful representations to the Seifert motion group.
△ Less
Submitted 5 January, 2023;
originally announced January 2023.
-
Score-based Generative Modeling Secretly Minimizes the Wasserstein Distance
Authors:
Dohyun Kwon,
Ying Fan,
Kangwook Lee
Abstract:
Score-based generative models are shown to achieve remarkable empirical performances in various applications such as image generation and audio synthesis. However, a theoretical understanding of score-based diffusion models is still incomplete. Recently, Song et al. showed that the training objective of score-based generative models is equivalent to minimizing the Kullback-Leibler divergence of th…
▽ More
Score-based generative models are shown to achieve remarkable empirical performances in various applications such as image generation and audio synthesis. However, a theoretical understanding of score-based diffusion models is still incomplete. Recently, Song et al. showed that the training objective of score-based generative models is equivalent to minimizing the Kullback-Leibler divergence of the generated distribution from the data distribution. In this work, we show that score-based models also minimize the Wasserstein distance between them under suitable assumptions on the model. Specifically, we prove that the Wasserstein distance is upper bounded by the square root of the objective function up to multiplicative constants and a fixed constant offset. Our proof is based on a novel application of the theory of optimal transport, which can be of independent interest to the society. Our numerical experiments support our findings. By analyzing our upper bounds, we provide a few techniques to obtain tighter upper bounds.
△ Less
Submitted 12 December, 2022;
originally announced December 2022.
-
OBMeshfree: An optimization-based meshfree solver for nonlocal diffusion and peridynamics models
Authors:
Yiming Fan,
Huaiqian You,
Yue Yu
Abstract:
We present OBMeshfree, an Optimization-Based Meshfree solver for compactly supported nonlocal integro-differential equations (IDEs) that can describe material heterogeneity and brittle fractures. OBMeshfree is developed based on a quadrature rule calculated via an equality constrained least square problem to reproduce exact integrals for polynomials. As such, a meshfree discretization method is ob…
▽ More
We present OBMeshfree, an Optimization-Based Meshfree solver for compactly supported nonlocal integro-differential equations (IDEs) that can describe material heterogeneity and brittle fractures. OBMeshfree is developed based on a quadrature rule calculated via an equality constrained least square problem to reproduce exact integrals for polynomials. As such, a meshfree discretization method is obtained, whose solution possesses the asymptotically compatible convergence to the corresponding local solution. Moreover, when fracture occurs, this meshfree formulation automatically provides a sharp representation of the fracture surface by breaking bonds, avoiding the loss of mass. As numerical examples, we consider the problem of modeling both homogeneous and heterogeneous materials with nonlocal diffusion and peridynamics models. Convergences to the analytical nonlocal solution and to the local theory are demonstrated. Finally, we verify the applicability of the approach to realistic problems by reproducing high-velocity impact results from the Kalthoff-Winkler experiments. Discussions on possible immediate extensions of the code to other nonlocal diffusion and peridynamics problems are provided. OBMeshfree is freely available on GitHub.
△ Less
Submitted 27 November, 2022;
originally announced November 2022.
-
Linear spectral Turan problems for expansions of graphs with given chromatic number
Authors:
Chuan-Ming She,
Yi-Zheng Fan,
Liying Kang,
Yaoping Hou
Abstract:
An $r$-uniform hypergraph is linear if every two edges intersect in at most one vertex. The $r$-expansion $F^{r}$ of a graph $F$ is the $r$-uniform hypergraph obtained from $F$ by enlarging each edge of $F$ with a vertex subset of size $r-2$ disjoint from the vertex set of $F$ such that distinct edges are enlarged by disjoint subsets. Let $ex_{r}^{lin}(n,F^{r})$ and $spex_{r}^{lin}(n,F^{r})$ be th…
▽ More
An $r$-uniform hypergraph is linear if every two edges intersect in at most one vertex. The $r$-expansion $F^{r}$ of a graph $F$ is the $r$-uniform hypergraph obtained from $F$ by enlarging each edge of $F$ with a vertex subset of size $r-2$ disjoint from the vertex set of $F$ such that distinct edges are enlarged by disjoint subsets. Let $ex_{r}^{lin}(n,F^{r})$ and $spex_{r}^{lin}(n,F^{r})$ be the maximum number of edges and the maximum spectral radius of all $F^{r}$-free linear $r$-uniform hypergraphs with $n$ vertices, respectively. In this paper, we present the sharp (or asymptotic) bounds of $ex_{r}^{lin}( n,F^{r})$ and $spex_{r}^{lin}(n,F^{r})$ by establishing the connection between the spectral radii of linear hypergraphs and those of their shadow graphs, where $F$ is a $(k+1)$-color critical graph or a graph with chromatic number $k$.
△ Less
Submitted 4 June, 2023; v1 submitted 24 November, 2022;
originally announced November 2022.
-
Sharpened Uncertainty Principle
Authors:
Yun Fan
Abstract:
For any finite group $G$, any finite $G$-set $X$ and any field $F$, we consider the vector space $F^X$ of all functions from $X$ to $F$. When the group algebra $FG$ is semisimple and splitting, we find a specific basis $\widehat X$ of $F^X$, construct the Fourier transform: $F^X\to F^{\widehat X}$, $f\mapsto\widehat f$, and define the rank support $\mbox{rk-supp}(\widehat f)$; we prove that…
▽ More
For any finite group $G$, any finite $G$-set $X$ and any field $F$, we consider the vector space $F^X$ of all functions from $X$ to $F$. When the group algebra $FG$ is semisimple and splitting, we find a specific basis $\widehat X$ of $F^X$, construct the Fourier transform: $F^X\to F^{\widehat X}$, $f\mapsto\widehat f$, and define the rank support $\mbox{rk-supp}(\widehat f)$; we prove that $\mbox{rk-supp}(\widehat f)=\dim FGf$, where $FGf$ is the submodule of the permutation module $FX$ generated by the element $f=\sum_{x\in X}f(x)x$. Next, we extend a sharpened uncertainty principle for abelian finite groups by Feng, Hollmann, and Xiang [9] to the following extensive framework: for any field $F$, any transitive $G$-set $X$ and $0\neq f\in F^X$ we prove that $$
|{\rm supp}(f)|\cdot \dim FGf \geq |X|+|{\rm supp}(f)|-|X_{{\rm supp}(f)}|, $$ where ${\rm supp}(f)$ is the support of $f$, and $X_{{\rm supp}(f)}$ is a block of $X$ associated with the subset ${\rm supp}(f)$ such that ${\rm supp}(f)$ is a disjoint union of some translations of the block. Then many (sharpened or classical) versions of finite-dimensional uncertainty principle are derived as corollaries.
△ Less
Submitted 3 January, 2023; v1 submitted 21 November, 2022;
originally announced November 2022.
-
Pieri and Murnaghan--Nakayama type Rules for Chern classes of Schubert Cells
Authors:
Neil J. Y. Fan,
Peter L. Guo,
Rui Xiong
Abstract:
We develop Pieri type as well as Murnaghan--Nakayama type formulas for equivariant Chern--Schwartz--MacPherson classes of Schubert cells in the classical flag variety. These formulas include as special cases many previously known multiplication formulas for Chern--Schwartz--MacPherson classes or Schubert classes. We apply the equivariant Murnaghan--Nakayama formula to the enumeration of rim hook t…
▽ More
We develop Pieri type as well as Murnaghan--Nakayama type formulas for equivariant Chern--Schwartz--MacPherson classes of Schubert cells in the classical flag variety. These formulas include as special cases many previously known multiplication formulas for Chern--Schwartz--MacPherson classes or Schubert classes. We apply the equivariant Murnaghan--Nakayama formula to the enumeration of rim hook tableaux.
△ Less
Submitted 12 November, 2022;
originally announced November 2022.
-
SIMPLE-RC: Group Network Inference with Non-Sharp Nulls and Weak Signals
Authors:
Jianqing Fan,
Yingying Fan,
Jinchi Lv,
Fan Yang
Abstract:
Large-scale network inference with uncertainty quantification has important applications in natural, social, and medical sciences. The recent work of Fan, Fan, Han and Lv (2022) introduced a general framework of statistical inference on membership profiles in large networks (SIMPLE) for testing the sharp null hypothesis that a pair of given nodes share the same membership profiles. In real applica…
▽ More
Large-scale network inference with uncertainty quantification has important applications in natural, social, and medical sciences. The recent work of Fan, Fan, Han and Lv (2022) introduced a general framework of statistical inference on membership profiles in large networks (SIMPLE) for testing the sharp null hypothesis that a pair of given nodes share the same membership profiles. In real applications, there are often groups of nodes under investigation that may share similar membership profiles at the presence of relatively weaker signals than the setting considered in SIMPLE. To address these practical challenges, in this paper we propose a SIMPLE method with random coupling (SIMPLE-RC) for testing the non-sharp null hypothesis that a group of given nodes share similar (not necessarily identical) membership profiles under weaker signals. Utilizing the idea of random coupling, we construct our test as the maximum of the SIMPLE tests for subsampled node pairs from the group. Such technique reduces significantly the correlation among individual SIMPLE tests while largely maintaining the power, enabling delicate analysis on the asymptotic distributions of the SIMPLE-RC test. Our method and theory cover both the cases with and without node degree heterogeneity. These new theoretical developments are empowered by a second-order expansion of spiked eigenvectors under the $\ell_\infty$-norm, built upon our work for random matrices with weak spikes. Our theoretical results and the practical advantages of the newly suggested method are demonstrated through several simulation and real data examples.
△ Less
Submitted 31 October, 2022;
originally announced November 2022.
-
The trace of uniform hypergraphs with application to Estrada index
Authors:
Yi-Zheng Fan,
Jian Zheng,
Ya Yang
Abstract:
In this paper we investigate the traces of the adjacency tensor of hypergraphs (simply called the traces of hypergraphs). We give new expressions for the traces of hypertrees and linear unicyclic hypergraphs by the weight function assigned to their connected sub-hypergraphs, and provide some perturbation results for the traces of a hypergraph with cut vertices. As applications we determine the uni…
▽ More
In this paper we investigate the traces of the adjacency tensor of hypergraphs (simply called the traces of hypergraphs). We give new expressions for the traces of hypertrees and linear unicyclic hypergraphs by the weight function assigned to their connected sub-hypergraphs, and provide some perturbation results for the traces of a hypergraph with cut vertices. As applications we determine the unique hypertree with maximum Estrada index among all hypertrees with fixed number of edges and perfect matchings, and the unique unicyclic hypergraph with maximum Estrada index among all unicyclic hypergraph with fixed number of edges and girth $3$.
△ Less
Submitted 6 October, 2022;
originally announced October 2022.
-
High-ordered spectral characterization of unicyclic graphs
Authors:
Yi-Zheng Fan,
Hong-Xia Yang,
Jian Zheng
Abstract:
In this paper we will apply the tensor and its traces to investigate the spectral characterization of unicyclic graphs. Let $G$ be a graph and $G^m$ be the $m$-th power (hypergraph) of $G$. The spectrum of $G$ is referring to its adjacency matrix, and the spectrum of $G^m$ is referring to its adjacency tensor. The graph $G$ is called determined by high-ordered spectra (DHS for short) if, whenever…
▽ More
In this paper we will apply the tensor and its traces to investigate the spectral characterization of unicyclic graphs. Let $G$ be a graph and $G^m$ be the $m$-th power (hypergraph) of $G$. The spectrum of $G$ is referring to its adjacency matrix, and the spectrum of $G^m$ is referring to its adjacency tensor. The graph $G$ is called determined by high-ordered spectra (DHS for short) if, whenever $H$ is a graph such that $H^m$ is cospectral with $G^m$ for all $m$, then $H$ is isomorphic to $G$. In this paper we first give formulas for the traces of the power of unicyclic graphs, and then provide some high-ordered cospectral invariants of unicyclic graphs. We prove that a class of unicyclic graphs with cospectral mates is DHS, and give two examples of infinitely many pairs of cospectral unicyclic graphs but with different high-ordered spectra.
△ Less
Submitted 28 August, 2022;
originally announced August 2022.
-
Center of the Yangian double in type A
Authors:
Yang Fan,
Naihuan Jing
Abstract:
We prove the R-matrix and Drinfeld presentations of the Yangian double in type A are isomorphic. The central elements of the completed Yangian double in type A at the critical level are constructed. The images of these elements under a Harish-Chandra-type homomorphism are calculated by applying a version of the Poincaré-Birkhoff-Witt theorem for the R-matrix presentation. These images coincide wit…
▽ More
We prove the R-matrix and Drinfeld presentations of the Yangian double in type A are isomorphic. The central elements of the completed Yangian double in type A at the critical level are constructed. The images of these elements under a Harish-Chandra-type homomorphism are calculated by applying a version of the Poincaré-Birkhoff-Witt theorem for the R-matrix presentation. These images coincide with the eigenvalues of the central elements in the Wakimoto modules.
△ Less
Submitted 20 April, 2023; v1 submitted 4 July, 2022;
originally announced July 2022.
-
FACT: High-Dimensional Random Forests Inference
Authors:
Chien-Ming Chi,
Yingying Fan,
Jinchi Lv
Abstract:
Quantifying the usefulness of individual features in random forests learning can greatly enhance its interpretability. Existing studies have shown that some popularly used feature importance measures for random forests suffer from the bias issue. In addition, there lack comprehensive size and power analyses for most of these existing methods. In this paper, we approach the problem via hypothesis t…
▽ More
Quantifying the usefulness of individual features in random forests learning can greatly enhance its interpretability. Existing studies have shown that some popularly used feature importance measures for random forests suffer from the bias issue. In addition, there lack comprehensive size and power analyses for most of these existing methods. In this paper, we approach the problem via hypothesis testing, and suggest a framework of the self-normalized feature-residual correlation test (FACT) for evaluating the significance of a given feature in the random forests model with bias-resistance property, where our null hypothesis concerns whether the feature is conditionally independent of the response given all other features. Such an endeavor on random forests inference is empowered by some recent developments on high-dimensional random forests consistency. Under a fairly general high-dimensional nonparametric model setting with dependent features, we formally establish that FACT can provide theoretically justified feature importance test with controlled type I error and enjoy appealing power property. The theoretical results and finite-sample advantages of the newly suggested method are illustrated with several simulation examples and an economic forecasting application.
△ Less
Submitted 12 November, 2023; v1 submitted 4 July, 2022;
originally announced July 2022.
-
Distribution of zeros of matching polynomials of hypergraphs
Authors:
Jiang-Chao Wan,
Yi Wang,
Yi-zheng Fan
Abstract:
Let $\h$ be a connected $k$-graph with maximum degree $Δ\geq 2$ and let $μ(\h, x)$ be the matching polynomial of $\h$. In this paper, we focus on studying the distribution of zeros of the matching polynomials of $k$-graphs. We prove that the zeros (with multiplicities) of $μ(\h, x)$ are invariant under a rotation of an angle $2π/{\ell}$ in the complex plane for some positive integer $\ell$ and…
▽ More
Let $\h$ be a connected $k$-graph with maximum degree $Δ\geq 2$ and let $μ(\h, x)$ be the matching polynomial of $\h$. In this paper, we focus on studying the distribution of zeros of the matching polynomials of $k$-graphs. We prove that the zeros (with multiplicities) of $μ(\h, x)$ are invariant under a rotation of an angle $2π/{\ell}$ in the complex plane for some positive integer $\ell$ and $k$ is the maximum integer with this property. Let $λ(\h)$ denote the maximum modulus of all zeros of $μ(\h, x)$. We show that $λ(\h)$ is a simple root of $μ(\h, x)$ and $$Δ^{1\over k} \leq λ(\h)< \frac{k}{k-1}\big((k-1)(Δ-1)\big)^{1\over k}.$$ To achieve these, we introduce the path tree $\T(\h,u)$ of $\h$ with respect to a vertex $u$ of $\h$, which is a $k$-tree, and prove that $$\frac{μ(\h-u,x)}{μ(\h, x)} = \frac{μ(\T(\h,u)-u,x) }{μ(\T(\h,u),x)},$$ which generalizes the celebrated Godsil's identity on the matching polynomial of graphs.
△ Less
Submitted 30 July, 2023; v1 submitted 19 June, 2022;
originally announced June 2022.
-
The trace and Estrada index of uniform hypergraphs with cut vertices
Authors:
Yi-Zheng Fan,
Ya Yang,
Chuan-Ming She,
Jian Zheng,
Yi-Min Song,
Hong-Xia Yang
Abstract:
Let $\mathcal{H}$ be an $m$-uniform hypergraph, and let $\mathcal{A}(\mathcal{H})$ be the adjacency tensor of $\mathcal{H}$ which can be viewed as a system of homogeneous polynomials of degree $m-1$. Morozov and Shakirov generalized the traces of linear systems to nonlinear homogeneous polynomial systems and obtained explicit formulas for multidimensional resultants. Sun, Zhou and Bu introduced th…
▽ More
Let $\mathcal{H}$ be an $m$-uniform hypergraph, and let $\mathcal{A}(\mathcal{H})$ be the adjacency tensor of $\mathcal{H}$ which can be viewed as a system of homogeneous polynomials of degree $m-1$. Morozov and Shakirov generalized the traces of linear systems to nonlinear homogeneous polynomial systems and obtained explicit formulas for multidimensional resultants. Sun, Zhou and Bu introduced the Estrada index of uniform hypergraphs which is closely related to the traces of their adjacency tensors. In this paper we give formulas for the traces of $\mathcal{A}(\mathcal{H})$ when $\mathcal{H}$ contains cut vertices, and obtain results on the traces and Estrada index when $\mathcal{H}$ is perturbed under local changes. We prove that among all hypertrees with fixed number of edges, the hyperpath is the unique one with minimum Estrada index and the hyperstar is the unique one with maximum Estrada index.
△ Less
Submitted 30 May, 2022;
originally announced May 2022.
-
Towards Programming Adaptive Linear Neural Networks Through Chemical Reaction Networks
Authors:
Yuzhen Fan,
Xiaoyu Zhang,
Chuanhou Gao
Abstract:
This paper is concerned with programming adaptive linear neural networks (ALNNs) using chemical reaction networks (CRNs) equipped with mass-action kinetics. Through individually programming the forward propagation and the backpropagation of ALNNs, and also utilizing the permeation walls technique, we construct a powerful CRN possessing the function of ALNNs, especially having the function of autom…
▽ More
This paper is concerned with programming adaptive linear neural networks (ALNNs) using chemical reaction networks (CRNs) equipped with mass-action kinetics. Through individually programming the forward propagation and the backpropagation of ALNNs, and also utilizing the permeation walls technique, we construct a powerful CRN possessing the function of ALNNs, especially having the function of automatic computation. We also provide theoretical analysis and a case study to support our construction. The results will have potential implications for the developments of synthetic biology, molecular computer and artificial intelligence.
△ Less
Submitted 13 April, 2022; v1 submitted 6 April, 2022;
originally announced April 2022.
-
A Meshfree Peridynamic Model for Brittle Fracture in Randomly Heterogeneous Materials
Authors:
Yiming Fan,
Huaiqian You,
Xiaochuan Tian,
Xiu Yang,
Xingjie Li,
Naveen Prakash,
Yue Yu
Abstract:
In this work we aim to develop a unified mathematical framework and a reliable computational approach to model the brittle fracture in heterogeneous materials with variability in material microstructures, and to provide statistic metrics for quantities of interest, such as the fracture toughness. To depict the material responses and naturally describe the nucleation and growth of fractures, we con…
▽ More
In this work we aim to develop a unified mathematical framework and a reliable computational approach to model the brittle fracture in heterogeneous materials with variability in material microstructures, and to provide statistic metrics for quantities of interest, such as the fracture toughness. To depict the material responses and naturally describe the nucleation and growth of fractures, we consider the peridynamics model. In particular, a stochastic state-based peridynamic model is developed, where the micromechanical parameters are modeled by a finite-dimensional random vector, or a combination of random variables truncating the Karhunen-Loève decomposition or the principle component analysis (PCA). To solve this stochastic peridynamic problem, probabilistic collocation method (PCM) is employed to sample the random field representing the micromechanical parameters. For each sample, the deterministic peridynamic problem is discretized with an optimization-based meshfree quadrature rule. We present rigorous analysis for the proposed scheme and demonstrate its convergence for a number of benchmark problems, showing that it sustains the asymptotic compatibility spatially and achieves an algebraic or sub-exponential convergence rate in the random space as the number of collocation points grows. Finally, to validate the applicability of this approach on real-world fracture problems, we consider the problem of crystallization toughening in glass-ceramic materials, in which the material at the microstructural scale contains both amorphous glass and crystalline phases. The proposed stochastic peridynamic solver is employed to capture the crack initiation and growth for glass-ceramics with different crystal volume fractions, and the averaged fracture toughness are calculated. The numerical estimates of fracture toughness show good consistency with experimental measurements.
△ Less
Submitted 14 February, 2022;
originally announced February 2022.
-
High-Dimensional Knockoffs Inference for Time Series Data
Authors:
Chien-Ming Chi,
Yingying Fan,
Ching-Kang Ing,
Jinchi Lv
Abstract:
The model-X knockoffs framework provides a flexible tool for achieving finite-sample false discovery rate (FDR) control in variable selection in arbitrary dimensions without assuming any dependence structure of the response on covariates. It also completely bypasses the use of conventional p-values, making it especially appealing in high-dimensional nonlinear models. Existing works have focused on…
▽ More
The model-X knockoffs framework provides a flexible tool for achieving finite-sample false discovery rate (FDR) control in variable selection in arbitrary dimensions without assuming any dependence structure of the response on covariates. It also completely bypasses the use of conventional p-values, making it especially appealing in high-dimensional nonlinear models. Existing works have focused on the setting of independent and identically distributed observations. Yet time series data is prevalent in practical applications in various fields such as economics and social sciences. This motivates the study of model-X knockoffs inference for time series data. In this paper, we make some initial attempt to establish the theoretical and methodological foundation for the model-X knockoffs inference for time series data. We suggest the method of time series knockoffs inference (TSKI) by exploiting the ideas of subsampling and e-values to address the difficulty caused by the serial dependence. We also generalize the robust knockoffs inference to the time series setting and relax the assumption of known covariate distribution required by model-X knockoffs, because such an assumption is overly stringent for time series data. We establish sufficient conditions under which TSKI achieves the asymptotic FDR control. Our technical analysis reveals the effects of serial dependence and unknown covariate distribution on the FDR control. We conduct power analysis of TSKI using the Lasso coefficient difference knockoff statistic under linear time series models. The finite-sample performance of TSKI is illustrated with several simulation examples and an economic inflation study.
△ Less
Submitted 19 May, 2023; v1 submitted 18 December, 2021;
originally announced December 2021.