Paper deep dive
Rethinking Irregular Time Series Forecasting from the Perspective of Basis Functions
Rongwen Li, Changjian Chen
Intelligence
Status: succeeded | Model: Gemma-4-26B-A4B | Prompt: intel-v1 | Confidence: 90%
Last extracted: 8/19/2026, 4:40:40 AM
Summary
The paper introduces DNBNet, a Debiased Neural Basis-Function Network for irregular time series forecasting. It addresses two main limitations of existing methods: asymptotic bias caused by ignoring sampling density and limited adaptability of predefined basis functions. DNBNet employs a debiased neural basis-function response mechanism using importance sampling and Kernel Density Estimation (KDE) to correct bias, while parameterizing basis functions with neural networks for adaptability. It also utilizes a multi-scale decomposition module with mass-aware fusion and a dual-branch decoder for improved prediction.
Entities (11)
Relation Signals (11)
DNBNet → addresses → Limited Adaptability of Predefined Basis Functions
confidence 95% · parameterizing basis functions with neural networks to adapt to diverse temporal patterns
DNBNet → addresses → Asymptotic Bias
confidence 95% · DNBNet... corrects asymptotic bias through importance sampling
Existing Methods → suffersfrom → Limited Adaptability of Predefined Basis Functions
confidence 92% · the limited adaptability of predefined basis functions to diverse temporal patterns
Existing Methods → suffersfrom → Asymptotic Bias
confidence 92% · existing methods... suffers from... a potential non-vanishing asymptotic bias caused by ignoring the sampling density
DNBNet → uses → Importance Sampling
confidence 90% · corrects asymptotic bias through importance sampling
DNBNet → uses → Kernel Density Estimation
confidence 90% · estimates the timestamp density via Kernel Density Estimation (KDE)
DNBNet → contains → Mass-aware Fusion
confidence 88% · together with a mass-aware fusion mechanism
DNBNet → contains → Multi-scale Decomposition
confidence 88% · design a novel multi-scale decomposition module based on average pooling
Cypher Suggestions (0)
No Cypher suggestions yet.
Abstract
Abstract:Irregular time series forecasting is crucial in many domains, such as healthcare and meteorological observation. However, due to the inherent characteristics of irregular time series, including sparse observations and non-uniform sampling, accurately predicting future dynamics remains challenging. In light of these two characteristics, many existing methods aggregate irregular observations into fixed-dimensional estimated response coefficients through predefined basis functions and use these coefficients as sequence representations. Nevertheless, this modeling paradigm still suffers from two key limitations: (i) a potential non-vanishing asymptotic bias caused by ignoring the sampling density of timestamps; and (ii) the limited adaptability of predefined basis functions to diverse temporal patterns. In this study, we propose a Debiased Neural Basis-Function Network (DNBNet) to address these challenges. Its core is a debiased neural basis-function response mechanism, which corrects asymptotic bias through importance sampling while parameterizing basis functions with neural networks to adapt to diverse temporal patterns. In addition, considering the sparsity of irregular data, we design a novel multi-scale decomposition module based on average pooling, together with a mass-aware fusion mechanism, to obtain richer representations. Finally, a dual-branch decoder is employed for forecasting. Extensive experiments on multiple real-world datasets demonstrate the effectiveness of DNBNet and its strong generalizability across diverse irregular time series scenarios. Our code can be obtained at this https URL.
Tags
Links
- Source: https://arxiv.org/abs/2608.17284v1
- Canonical: https://arxiv.org/abs/2608.17284v1
Trouble viewing inline? Open PDF directly →
Full Text
70,147 characters extracted from source content.
Expand or collapse full text
Rethinking Irregular Time Series Forecasting from the Perspective of Basis Functions Rongwen Li Affiliation: College of Computer Science and Electronic Engineering, Hunan University, Changsha, Hunan, China Changjian Chen Affiliation: College of Computer Science and Electronic Engineering, Hunan University, Changsha, Hunan, China Correspondence to: changjianchen@hnu.edu.cn Abstract Irregular time series forecasting is crucial in many domains, such as healthcare and meteorological observation. However, due to the inherent characteristics of irregular time series, including sparse observations and non-uniform sampling, accurately predicting future dynamics remains challenging. In light of these two characteristics, many existing methods aggregate irregular observations into fixed-dimensional estimated response coefficients through predefined basis functions and use these coefficients as sequence representations. Nevertheless, this modeling paradigm still suffers from two key limitations: (i) a potential non-vanishing asymptotic bias caused by ignoring the sampling density of timestamps; and (i) the limited adaptability of predefined basis functions to diverse temporal patterns. In this study, we propose a Debiased Neural Basis-Function Network (DNBNet) to address these challenges. Its core is a debiased neural basis-function response mechanism, which corrects asymptotic bias through importance sampling while parameterizing basis functions with neural networks to adapt to diverse temporal patterns. In addition, considering the sparsity of irregular data, we design a novel multi-scale decomposition module based on average pooling, together with a mass-aware fusion mechanism, to obtain richer representations. Finally, a dual-branch decoder is employed for forecasting. Extensive experiments on multiple real-world datasets demonstrate the effectiveness of DNBNet and its strong generalizability across diverse irregular time series scenarios. Our code can be obtained at https://github.com/hnu-vis/DNBNet. Keywords: Machine Learning, ICML 1 Introduction Irregular Multivariate Time Series (IMTS) consist of observations from multiple variables collected asynchronously at non-uniform, variable-specific timestamps, often with unequal sampling intervals and missing values. Such data are widely encountered in various domains, including medical monitoring, sensor networks, and meteorological observations (6; 4; 1). They are typically generated by underlying continuous-time dynamic processes, but can only be observed through sparse and irregular measurements in practice. As a result, learning reliable representations and making accurate predictions remain challenging. Table 1: Comparison of the types of predefined basis functions adopted by different methods, including their specific mathematical forms and function curves, reflecting their distinct modeling assumptions and interpretations of the temporal dynamics in the data. Basis Type Formulation Shape Related Work RBF basis function exp(−(t−ck)22σk2) \! (- (t-c_k)^22 _k^2 ) (29) Fourier basis function cos(2πωkt),sin(2πωkt) aligned & (2π _kt),\\ & (2π _kt) aligned (5; 19) Attention basis function exp(ckQ⊤Kt)d (c_kQ Kt) d (22; 9) Soft-window basis function σ(tkright−tsoftplus(δ))⋅σ(t−tkleftsoftplus(δ)) σ\! ( t_k^right-tsoftplus(δ) )·σ\! ( t-t_k^leftsoftplus(δ) ) (11) To address this issue, many recent methods (22; 29; 11; 9; 19) attempt to map irregular observations into regular fixed-dimensional representations for subsequent prediction. Although these methods differ substantially in model design, we observe that most of them can be interpreted from a unified basis function perspective: the model computes the responses of an irregular time series to a set of predefined temporal basis functions, and uses the resulting coefficients as the sequence representation. These coefficients summarize how strongly the observed irregular time series aligns with each temporal basis, thereby encoding potential temporal patterns in a compact form. From this perspective, the key distinction among different methods lies in the form of the adopted basis functions. As shown in Table 1, existing methods use different predefined temporal basis functions and their corresponding response curves to encode temporal information. For example, radial basis functions (RBFs) typically emphasize local temporal responses and are therefore well suited for capturing information around specific time regions (29). In contrast, Fourier basis functions provide global temporal support and are more naturally aligned with modeling long-range structures such as periodicity and seasonality (19). Therefore, the choice of basis functions is not merely a technical detail; it implicitly determines which temporal patterns are easier for the model to represent, thereby introducing different structural priors and inductive biases. This unified perspective also reveals two key limitations of existing methods. 1) Potential asymptotic bias under irregular sampling: existing discrete approximations for response coefficients are generally reasonable under regular or nearly uniform sampling, but they ignore the distribution of observation timestamps in irregular sampling scenarios, which may lead to a non-vanishing asymptotic bias. 2) Limited adaptability of predefined bases: although predefined basis functions introduce useful inductive biases, their fixed forms restrict the model’s ability to adapt to different datasets and sampling patterns for irregular time series. To overcome the above limitations, we propose a Debiased Neural Basis-Function Network (DNBNet). Specifically, we design a Debiased Neural Basis-function Response mechanism, which estimates the timestamp density via Kernel Density Estimation (KDE) to correct the asymptotic bias caused by non-uniform sampling. Moreover, it parameterizes temporal basis functions with learnable neural networks, thereby improving the model’s adaptability to complex temporal patterns. In addition, considering the sparsity of irregular time series data, we further develop a novel multi-scale decomposition based on average pooling and mass-aware fusion mechanism to obtain richer feature representations. Finally, the fused representations are fed into a dual-branch decoder for more accurate prediction. Extensive experiments on multiple real-world datasets demonstrate the effectiveness of the proposed framework. • A unified basis-function perspective of irregular time series forecasting, which reveals two key limitations of existing methods: potential asymptotic estimation bias caused by irregular sampling, and the limited ability of predefined basis functions to adaptively model diverse temporal patterns. • A DNBNet with a neural basis-response mechanism that corrects the asymptotic bias and produces more expressive response coefficients. Building on this mechanism, we also incorporate multi-scale decomposition and mass-aware fusion to learn richer feature representations. • Extensive experiments on multiple irregular time series forecasting benchmarks show that DNBNet consistently achieves superior predictive performance compared with existing competitive methods. 2 Related Work 2.1 Irregular Multivariate Time Series Forecasting Recent IMTS forecasting methods can be roughly divided into three categories. (i) Some methods directly model the underlying continuous time based ODE or CDE (14; 16). For instance, Latent ODE (20) leverages ODEs to construct continuous-time sequences in the latent space, thereby enabling interpolation and extrapolation. To avoid inefficient numerical integration for ODEs, Neural Flow (2) turn to directly model the solution of ODEs. (i) Methods based on specific data structures model message passing among observation points (24; 27). For example, SeFT (7) utilizes set, whereas GraFITi (24) and HyperIMTS (10) use bipartite graphs and hypergraphs, respectively. (i) Some methods learn compact representations through temporal alignment and aggregation (26; 12). tPatchGNN (26) integrates region-level information by a fixed-span patching strategy, and APN (11) designs an adaptive patch-span learning mechanism to avoid the meaningless and sparse feature representation from rigid patch size. Despite these advancements, many representation learning methods ignore the influence of timestamp density, thereby resulting in asymptotic biases. 2.2 Basis Functions for Irregular Multivariate Time Series Basis functions provide a classical perspective for continuous signals, i.e., decomposing and characterizing complex signals through a set of response coefficients over temporal patterns. This idea has been widely used in regularly sampled time series analysis (17). For example, some methods (28; 18) capture periodic components of sequences via the Fast Fourier Transform (FFT). Similarly, many IMTS methods also utilize basis functions to extract temporal information. For instance, mTAN (22) proposes a time-aware attention mechanism to aggregate irregular observations using attention-based basis functions with predefined reference points. KAFNet (29) utilizes Gaussian temporal kernel aggregation constructed from radial basis functions (RBFs) to efficiently model irregular time series. LSCD (5) performs imputation with better consideration of the global structure of irregular time series through Fourier bases. However, previous methods usually ignore the inherent non-uniform sampling of irregular time series when estimating basis responses, thereby introducing asymptotic bias. 3 Preliminaries and Motivation In this section, we first formalize the irregular time series forecasting problem and then present the discrete computation form of basis-function response coefficients adopted by existing methods. Based on this formulation, we further analyze the bias introduced by discrete response estimation under irregular sampling, which motivates us to design a debiased response estimation mechanism. 3.1 Problem Definition Given an irregular multivariate time series =nn=1NO=\O^n\_n=1^N with N variates, where n=(tin,xin)i=1LnO^n=\(t_i^n,x_i^n)\_i=1^L_n denotes the observations of the n-th variate at non-uniform timestamps, and (tin,xin)(t_i^n,x_i^n) represents the value xinx_i^n recorded at time tint_i^n. Irregular time series forecasting aims to predict future values at a set of query timestamps =nn=1NQ=\Q^n\_n=1^N, with n=qjnj=1QnQ^n=\q_j^n\_j=1^Q_n. Formally, the goal is to learn a model fθf_θ that maps the historical observations and query timestamps to the corresponding predictions: fθ(,)→^=y^jnj=1Qnn=1N, f_θ(O,Q)→ Y= \\ y_j^n\_j=1^Q_n \_n=1^N, (1) where y^jn y_j^n denotes the predicted value of the n-th variate at timestamp qjnq_j^n. 3.2 Basis-Function Response Representation Given a continuous-time signal x(t)x(t) and a set of temporal basis functions ϕk(t)k=1K\ _k(t)\_k=1^K, the ideal response coefficient with respect to the k-th basis is defined by the following integral: ck=⟨x,ψk⟩=∫x(t)ψk(t)t, c_k= x, _k = _Tx(t) _k(t)\,dt, (2) where ψk(t)=ϕk(t)∫ϕk(t)t _k(t)\!\!=\!\! _k(t) _T _k(t)\,dt or ψk(t)=ϕk(t)∫ϕk2(t)t _k(t)\!\!=\!\! _k(t) _T _k^2(t)\,dt correspond to the weighted average response and projection response for nonnegative and orthogonal basis functions, respectively. For example, a commonly used weighted average response coefficient is given by ck=∫x(t)ϕk(t)t∫ϕk(t)tc_k\!\!=\!\! _Tx(t) _k(t)\,dt _T _k(t)\,dt. In practice, only discrete observations (ti,xi)i=1L\(t_i,x_i)\_i=1^L are available. Existing methods therefore usually approximate the above continuous integral by summing over the observed data, e.g., c^k=∑i=1Lxiϕk(ti)∑i=1Lϕk(ti) c_k\!=\! _i=1^Lx_i _k(t_i) _i=1^L _k(t_i) (22; 29; 11). In the following, we focus on the weighted average response formulation for clarity, while the analysis for the projection response is similar. 3.3 Asymptotic Bias Analysis of Discrete Response Estimation Existing methods approximate the continuous-time response ckc_k using the discrete computation form c^k c_k, which is reasonable under regular sampling. However, when the observations are irregularly sampled, this discrete approximation may lead to an asymptotic bias. To facilitate the subsequent analysis, we consider a normalized time domain =[0,1]T=[0,1] and introduce the following mild assumption: Assumption 3.1. The continuous signal x(t)x(t) and the nonnegative basis functions ϕk(t) _k(t) are bounded, i.e., |x(t)|≤1|x(t)|≤ 1 and |ϕk(t)|≤m| _k(t)|≤ m for any t and k. Moreover, each basis function has a nontrivial mass under both the uniform time density and the timestamp sampling density: min∫ϕk(t)t,∫ϕk(t)p(t)t≥a. \ _T _k(t)dt, _T _k(t)p(t)dt \≥ a. (3) Assumption 3.1 is reasonable in practice. Data normalization and the construction of basis functions typically ensure the existence of upper bounds for x(t)x(t) and ϕk(t) _k(t), while the lower bound on the basis mass guarantees that each basis function carries sufficient information. Under Assumption 3.1, we give the following theorem, which provides an estimation error bound between the discrete coefficient c^k c_k and the ideal response coefficient ckc_k: Theorem 3.2. Given observations (ti,xi)i=1L\(t_i,x_i)\_i=1^L, where ti∼i.i.d.p(t)t_i .i.d. p(t), for any δ∈(0,1)δ∈(0,1), when L is sufficiently large, with probability at least 1−δ1-δ, the following inequality holds: maxk∈[K]|c^k−ck|≤4ma2log(4K/δ)L+2ma‖p(t)−1‖L1, _k∈[K]| c_k-c_k|≤ 4ma 2 (4K/δ)L+ 2ma\|p(t)-1\|_L^1, where ‖p(t)−1‖L1=∫|p(t)−1|t\|p(t)-1\|_L^1= _T|p(t)-1|dt. Proof. See Appendix A. ∎ Figure 1: (a) The overall pipeline of DNBNet. (b) The debiased neural basis-function response computes debiased coefficients and basis masses using KDE and learnable basis functions. (c) The basis mass-aware fusion module aggregates coefficients from different scales. (d) The final prediction is obtained by combining the feature decoder with basis-function reconstruction. The first term on the right-hand side of the above bound decreases as the number of observations L increases, reflecting the standard finite-sample estimation error. In contrast, the second term is independent of L and does not vanish with more observations unless the timestamp density p(t)p(t) approaches the uniform distribution (p(t)=1p(t)=1), i.e., the observations are obtained through regular sampling. Therefore, this term characterizes the potential asymptotic bias induced by irregular sampling, which does not vanish as L increases. This observation motivates the need for a debiased response estimation mechanism that explicitly considers the sampling density of observed timestamps. 4 Debiased Neural Basis-Function Network In this section, we introduce the proposed DNBNet in detail. As shown in Fig. 1, we first generate observation subsequences at different scales through time-aware average pooling. We then design a debiased neural basis response mechanism to compute debiased and more expressive response coefficients. Next, we introduce a mass-aware fusion mechanism to integrate multi-scale response information. Finally, a dual decoder is employed to produce the final predictions. 4.1 Time-aware Average Pool The underlying dynamics of irregular time series often exhibit heterogeneous temporal patterns across different time scales, while the inherent sparsity of the original sequence limits its ability to provide sufficient information for feature learning. This motivates us to construct subsequences at different scales to capture temporal patterns of varying granularities and obtain richer representations through fusion. Therefore, we first utilize time-aware average pooling to obtain the subsequence of the original input sequence =(ti,xi)i=1LO=\(t_i,x_i)\_i=1^L at the s-th scale: (s) ^(s) =AvgPool(,ws,rs), =AvgPool (O;w_s,r_s ), (4) (s) ^(s) =(ti(s),xi(s))i=1Ls, =\(t_i^(s),x_i^(s))\_i=1^L_s, (5) where wsw_s and rsr_s respectively denote the window size and stride defined by absolute time spans at scale s, and we let (1)O^(1) denotes the sequence of original scale. Note that our model learns representation mainly in a channel-independent manner. For convenience, we temporarily omit the variate index n in this section. 4.2 Debiased Neural Basis-function Response According to the previous analysis, existing discrete computation methods suffer from asymptotic bias under irregular sampling. To correct this bias, we employ importance sampling to construct a density-corrected estimator of the original integral: ∫x(t)ϕk(t)p(t)p(t)t∫ϕk(t)p(t)p(t)t=t[x(t)ϕk(t)p(t)]t[ϕk(t)p(t)]≈∑i=1Lxiϕk(ti)p(ti)∑i=1Lϕk(ti)p(ti)=c~k. aligned _T x(t) _k(t)p(t)p(t)dt _T _k(t)p(t)p(t)dt\!=\! E_t [ x(t) _k(t)p(t) ]E_t [ _k(t)p(t) ]\!≈\! _i=1^L x_i _k(t_i)p(t_i) _i=1^L _k(t_i)p(t_i)\!=\! c_k. aligned (6) The proposed estimator asymptotically corrects the bias. A detailed explanation of the debiasing property of our estimates is provided in the Appendix A.2. Note that c^k c_k of the existing methods is actually a special case of our formulation c~k c_k under the uniform timestamp distribution, i.e., p(t)=1p(t)=1. In irregular sampling scenarios, 1/p(t)1/p(t) can be interpreted as an importance weight that explicitly corrects the sampling bias: it suppresses the overemphasis on densely sampled regions while compensating for contributions from sparsely sampled regions. In our implementation, we use Kernel Density Estimation (KDE) to estimate the timestamp density: p(ti)=1Lh∑j=1Lκ(ti−tjh), p(t_i)= 1Lh _j=1^Lκ ( t_i-t_jh ), (7) where κ(⋅)κ(·) denotes the kernel function, for which we adopt the Gaussian kernel. We further parameterize the bandwidth as a learnable variable, i.e., h=softplus(ρ)h=softplus(ρ), to adapt to different temporal scales. In addition, as shown in Table 1, many existing methods rely on predefined basis functions. The choice of basis functions usually reflects prior assumptions about the characteristics of the data. However, in the real scenarios, a fixed set of basis functions may be difficult to adapt to diverse datasets due to the changes of data distribution. Meanwhile, selecting appropriate basis functions for different datasets is highly dependent on domain expert knowledge, further limiting the generalization and applicability of the method. Therefore, considering the strong approximation ability of neural networks, we directly learn basis functions themselves from raw data by fully connected neural networks. Specifically, let ϕ(t)=[ϕ1(t),…,ϕK(t)]∈ℝK φ(t)= [ _1(t),…, _K(t) ] ^K denote the vector formed by concatenating K basis functions. We parameterize it with a two-layer MLP: ϕ(t)=softmax(MLPbasis(t)), φ(t)=softmax(MLP_basis(t)), (8) where softmax is to bound basis functions, consistent with Assumption 3.1. It also normalizes the responses across basis functions, thereby improving the stability of computation. Next, we incorporate the learnable basis-function ϕ(t)φ(t) into the debiased estimation in Eq. (37). For each scale-specific sequence (s)O^(s), we compute the corresponding basis response coefficient with respect to the k-th basis function: c~k(s)=∑i=1Lsxi(s)ϕk(ti(s))/p(ti(s))∑i=1Lsϕk(ti(s))/p(ti(s)), c_k^(s)= _i=1^L_sx_i^(s) _k(t_i^(s))/p(t_i^(s)) _i=1^L_s _k(t_i^(s))/p(t_i^(s)), (9) In this way, each scale produces a set of basis response coefficients c~k(s)k=1K\ c_k^(s)\_k=1^K. 4.3 Basis Mass-aware Fusion To obtain richer representations, we further fuse the basis response coefficients across different scales. Note that the denominator in Eq. (9) can be viewed as the basis mass of the k-th basis function at scale s: Mk(s)=∑i=1Lsϕk(ti(s))/p(ti(s)). M_k^(s)= _i=1^L_s _k(t_i^(s))/p(t_i^(s)). (10) A larger Mk(s)M_k^(s) indicates that the corresponding basis function receives stronger support at scale s, and thus its coefficient may contain more informative signals. Based on this observation, we design a mass-aware fusion mechanism to adaptively aggregate coefficients from different scales: rs,k r_s,k =τs,klog(1+Mk(s))+βs,k, = _s,k (1+M_k^(s))+ _s,k, (11) αs,k _s,k =softmax(rs,ks=1S), =softmax (\r_s,k\_s=1^S ), (12) c~k c_k =∑s=1Sαs,kc~k(s), = _s=1^S _s,k c_k^(s), (13) where τs,k _s,k and βs,k _s,k are learnable scale and bias parameters. The logarithmic transformation prevents excessively large basis masses from dominating the fusion weights. Finally, we concatenate the fused coefficients of all basis functions as =[c~1,…,c~K]∈ℝKc=[ c_1,…, c_K] ^K, and further project it into a latent space to enhance its representation capacity: =Linear(), =Linear(c), (14) =LayerNorm(+MLP()). =LayerNorm (h+MLP(h) ). (15) The representation ∈ℝdz ^d encodes debiased multi-scale basis responses for the subsequent forecasting decoder. 4.4 Dual-Branch Forecasting Decoder To exploit both the latent representation and basis coefficients, we design a dual-branch forecasting decoder consisting of a feature decoder and a basis reconstruction branch. Given a future query timestamp qlnq_l^n from the n-th variate at the l-th prediction step, following existing methods (26), we construct its time embedding as follows: TE(qln)[d]=ω0⋅qln+b0,if d=0,sin(ωd⋅qln+bd),if 0<d≤Dt. (q_l^n)[d]= cases _0· q_l^n+b_0,&if d=0,\\ ( _d· q_l^n+b_d),&if 0<d≤ D_t. cases (16) where ω0,b0\ _0,b_0\ and ωd,bdd=1Dt\ _d,b_d\_d=1^D_t are learnable parameters for the linear and periodic components, respectively. To distinguish different variates, we add a variate embedding ne_n to obtain n=+nz^n=z+e_n. The feature-branch prediction is then computed by y^ln,fea=MLPfea([n,TE(qln)]). y_l^n,fea=MLP_fea ( [z^n,TE(q_l^n) ] ). (17) The basis functions provide a natural way to reconstruct the signal from the learned coefficients. This allows us to explicitly extrapolate future values by evaluating the learned basis functions at query timestamps and combining them with the basis coefficients. Therefore, in the basis branch, the prediction for the n-th variate at query timestamp qlnq_l^n is computed as y^ln,basis=∑k=1Kc~kϕk(qln), y_l^n,basis= _k=1^K c_k _k(q_l^n), (18) where c~k c_k denotes the fused coefficient in Eq. (13) and ϕ(t) φ(t) is the same learnable basis-function in Eq. (8). Finally, we fuse the two branches via a learnable gate: y^ln=λy^ln,fea+(1−λ)y^ln,basis, y_l^n=λ\, y_l^n,fea+(1-λ)\, y_l^n,basis, (19) where λ=sigmoid(γ)λ=sigmoid(γ) is a learnable parameter that balances the contribution of the feature and basis branch. 5 Experiments 5.1 Experimental Setup 5.1.1 Datasets and Baselines To evaluate our method on irregular multivariate time series forecasting, we conduct experiments on five widely used benchmark datasets across multiple real-world domains. Specifically, we consider two healthcare datasets, PhysioNet and MIMIC, containing ICU clinical records; two biomechanics-related datasets, Human Activity and Student Life, consisting of sensor data from subjects performing various activities; and one climate dataset, USHCN, containing historical meteorological observations from U.S. weather stations. Following common practice, each dataset is split into training, validation, and test sets with proportions of 80%80\%, 10%10\%, and 10%10\%, respectively. We compare our method with 12 baseline models covering IMTS forecasting, classification, and interpolation, including PrimeNet (3), SeFT (7), mTAN (22), CRU (21), GNeuralFlow (14), Raindrop (27), tPatchGNN (26), GraFITi (24), Warpformer (25), Hi-Patch (12), KAFNet (29), and APN (11). Additional details on datasets and baselines are provided in Appendix B.1. Table 2: MSE results of DNBNet and baselines on five irregular time-series datasets. Best results are in bold and second-best results are underlined. Method USHCN Human Activity Student Life PhysioNet MIMIC PrimeNet 0.7328±0.00000.7328± 0.0000 4.2564±0.00074.2564± 0.0007 0.8799±0.00010.8799± 0.0001 0.7952±0.00000.7952± 0.0000 0.9728±0.00010.9728± 0.0001 SeFT 0.6647±0.00080.6647± 0.0008 1.3750±0.00361.3750± 0.0036 0.8747±0.00260.8747± 0.0026 0.7801±0.00260.7801± 0.0026 0.9905±0.00030.9905± 0.0003 mTAN 0.4197±0.01610.4197± 0.0161 0.1017±0.00490.1017± 0.0049 0.6555±0.00640.6555± 0.0064 0.3787±0.00570.3787± 0.0057 0.9719±0.01690.9719± 0.0169 CRU 0.5179±0.01130.5179± 0.0113 0.1452±0.00090.1452± 0.0009 0.7487±0.00550.7487± 0.0055 0.6169±0.00360.6169± 0.0036 0.7370±0.00590.7370± 0.0059 GNeuralFlow 0.4971±0.00930.4971± 0.0093 0.2550±0.03730.2550± 0.0373 0.7168±0.04220.7168± 0.0422 0.3516±0.00280.3516± 0.0028 0.6631±0.01120.6631± 0.0112 Raindrop 0.4841±0.00640.4841± 0.0064 0.0939±0.00170.0939± 0.0017 0.6672±0.00540.6672± 0.0054 0.3491±0.00150.3491± 0.0015 0.6219±0.00400.6219± 0.0040 tPatchGNN 0.5499±0.10680.5499± 0.1068 0.0578±0.00070.0578± 0.0007 0.6332±0.00230.6332± 0.0023 0.3097±0.00140.3097± 0.0014 0.4673±0.00160.4673± 0.0016 GraFITi 0.4270±0.00930.4270± 0.0093 0.0606±0.00290.0606± 0.0029 0.6366±0.00010.6366± 0.0001 0.3021±0.00110.3021± 0.0011 0.4349±0.0070¯ 0.4349± 0.0070 Warpformer 0.4263±0.00730.4263± 0.0073 0.0576±0.00090.0576± 0.0009 0.6267±0.0066¯ 0.6267± 0.0066 0.3153±0.00110.3153± 0.0011 0.4358±0.00180.4358± 0.0018 Hi-Patch 0.4108±0.0049¯ 0.4108± 0.0049 0.0592±0.00030.0592± 0.0003 0.6277±0.00100.6277± 0.0010 0.3226±0.00150.3226± 0.0015 0.4724±0.00670.4724± 0.0067 KAFNet 0.4132±0.00600.4132± 0.0060 0.0572±0.00020.0572± 0.0002 0.6312±0.00090.6312± 0.0009 0.3328±0.00150.3328± 0.0015 0.4612±0.00730.4612± 0.0073 APN 0.4262±0.00300.4262± 0.0030 0.0562±0.0004¯ 0.0562± 0.0004 0.6512±0.00150.6512± 0.0015 0.3133±0.00050.3133± 0.0005 0.4527±0.00930.4527± 0.0093 DNBNet (Ours) 0.4013±0.01080.4013± 0.0108 0.0551±0.00010.0551± 0.0001 0.6136±0.00020.6136± 0.0002 0.3089±0.0005¯ 0.3089± 0.0005 0.4289±0.00120.4289± 0.0012 Table 3: Ablation study of DNBNet on three irregular time-series datasets, with MSE reported. Variant USHCN PhysioNet MIMIC w/o p(t)p(t) 0.4449±0.02910.4449± 0.0291 0.3161±0.00030.3161± 0.0003 0.4485±0.00120.4485± 0.0012 Rep. RBF 0.4144±0.00910.4144± 0.0091 0.3231±0.00070.3231± 0.0007 0.4507±0.00090.4507± 0.0009 Rep. Fourier 0.4101±0.00880.4101± 0.0088 0.3267±0.00020.3267± 0.0002 0.4869±0.00480.4869± 0.0048 w/o AvgPool 0.4195±0.01660.4195± 0.0166 0.3119±0.00060.3119± 0.0006 0.4494±0.00120.4494± 0.0012 w/o BasDec 0.4552±0.04240.4552± 0.0424 0.3092±0.00020.3092± 0.0002 0.4613±0.00230.4613± 0.0023 DNBNet (Full) 0.4013±0.01080.4013± 0.0108 0.3089±0.00050.3089± 0.0005 0.4289±0.00120.4289± 0.0012 (a) MSE on three datasets w.r.t. different K. (b) MSE on three datasets w.r.t. different S. Figure 2: Parameter sensitivity study of number of basis functions K and Scales S on three datasets. 5.1.2 Implementation Details Following the prior work (2; 12), for USHCN, we use 3 years of historical observations to predict the next 1 year; for Human Activity, the look-back window and prediction horizon are set to 3000 ms and 1000 ms, respectively; for Student Life, we use observations from the past 30 days to predict the next 30 days; and for PhysioNet and MIMIC, we use observations from the past 36 hours to predict the next 3 observation points. All experiments are implemented in PyTorch 2.5.1+cu121 and conducted on a single NVIDIA 4090 GPU. During training, all models are optimized using the AdamW optimizer, with mean squared error (MSE) adopted as the training objective. We set the maximum number of training epochs to 200 and use an early stopping patience of 10 to prevent overfitting. Mean squared error and mean absolute error (MAE) are used as evaluation metrics. Each experiment is repeated with five random seeds, and we report the mean and standard deviation of the results. 5.2 Main results Due to space limitations, we report the MSE results of our method and all baselines across all datasets in Table 2, while the MAE results are provided in Appendix B.2. As shown in the table, DNBNet achieves the best performance on four of the five datasets and obtains the second-best result on PhysioNet. In particular, compared with Warpformer, the average second-best baseline, DNBNet achieves an average relative MSE reduction of approximately 3.18%3.18\%. Furthermore, compared with methods that also utilize basis-function based modeling, such as mTAN, KAFNet, and APN, DNBNet consistently achieves leading results on datasets from different domains, including healthcare and climate field. This indicates that, compared with directly relying on predefined basis functions, the proposed learnable basis functions can more effectively adapt to differences in sampling frequency or temporal irregularity across datasets, which is consistent with our motivation. In addition, we provide experimental results and analyses under different look-back windows and prediction horizons in the Appendix B.3, further validating the stable advantages of DNBNet. 5.3 Ablation Study We conduct ablation studies on three irregular time-series datasets to verify the effectiveness of our model components. The main results are reported in Table 3, with complete results provided in the Appendix B.4. Specifically, (1) w/o p(t)p(t) removes the density factor in Eq. (9), leading to consistent performance degradation and confirming the effectiveness of our debiased response mechanism. (2) Rep. RBF and Rep. Fourier replace our learnable basis functions ϕ(t)φ(t) with RBF and Fourier bases, respectively, both yielding inferior performance on all datasets. This suggests that our learnable basis-function is more robust and generalizable than predefined bases across different datasets. (3) w/o AvgPool and w/o BasDec remove the average pooling module and the basis-function prediction branch, respectively. The consistent performance degradation validates the effectiveness of leveraging multi-scale features and further demonstrates the potential of basis-function coefficients in representing future continuous dynamics. (4) To further validate the effectiveness of the debiased mechanism, we introduce the density debiasing factor p(t)p(t) into two existing basis-function-based methods, KAFNet and APN. Specifically, we replace their biased weighted average responses with the estimator in Eq. (37). Table 4 reports the original results of KAFNet and APN, as well as their performance after incorporating p(t)p(t) on three datasets. Compared with the original biased weighted average aggregation, simply dividing by p(t)p(t) yields consistent and stable improvements with almost no additional complexity. This verifies the effectiveness of our debiasing analysis. 5.4 Parameter Sensitivity In this subsection, we study the impact of the number of basis functions and the multi-scale aggregation setting on model performance. As shown in Figure 2, we report the results with different numbers of basis functions K and scales S on three datasets. We observe that using more basis functions does not necessarily lead to better performance; relatively small values, such as K=16K=16 or K=32K=32, are often sufficient to achieve the best or competitive results. Similarly, aggregating features from 22 or 33 scales already yields strong predictive performance, while further increasing the number of scales does not bring consistent gains. These results indicate that DNBNet is not sensitive to key hyperparameters and can maintain stable performance over a wide range of settings, providing practical guidance for transferring it to other datasets or application scenarios. 5.5 Visualization and Efficiency Analysis To better understand the basis functions learned from data, we visualize the top four basis functions ϕ(t)φ(t) with the largest mass on PhysioNet and MIMIC in Figure 3. For PhysioNet, the model learns localized basis functions similar to RBFs. Notably, each learned basis function has a distinct center, bandwidth, and peak location, suggesting that the model can adaptively adjust its local response range, which is difficult to achieve with predefined RBF bases. For the more challenging MIMIC dataset, DNBNet learns more complex and diverse basis-function shapes, including the soft window function ϕ12 _12 and sigmoid-like functions ϕ0 _0, ϕ8 _8, and ϕ13 _13. These basis functions capture response patterns over different temporal regions, further showing that DNBNet can learn temporal basis representations suited to the data distribution and dynamics, rather than being constrained by fixed predefined forms. Figure 4 further compares the computational efficiency of DNBNet with six competitive baselines on MIMIC, the largest dataset, under the same batch size. DNBNet achieves the best forecasting performance while introducing nearly the lowest computational cost, except for APN. This efficiency mainly comes from avoiding quadratic-complexity operations and instead extracting temporal dynamics through learnable basis functions, resulting in a lightweight and efficient model design. Table 4: Effectiveness of the density debiasing factor 1/p(t)1/p(t) on existing basis-function-based methods. Method USHCN PhysioNet MIMIC KAFNet w/o p(t)p(t) 0.4132±0.00600.4132± 0.0060 0.3328±0.00150.3328± 0.0015 0.4612±0.00730.4612± 0.0073 KAFNet w/ p(t)p(t) 0.4050±0.01730.4050± 0.0173 0.3117±0.00070.3117± 0.0007 0.4488±0.00240.4488± 0.0024 APN w/o p(t)p(t) 0.4262±0.00300.4262± 0.0030 0.3133±0.00050.3133± 0.0005 0.4527±0.00930.4527± 0.0093 APN w/ p(t)p(t) 0.4193±0.01560.4193± 0.0156 0.3094±0.00200.3094± 0.0020 0.4409±0.00190.4409± 0.0019 (a) The learned basis functions with the four largest basis function mass on PhysioNet dataset. (b) The learned basis functions with the four largest basis function mass on MIMIC dataset. Figure 3: Visualization of the four learned basis functions with the largest basis function mass on two datasets. Figure 4: Model efficiency comparison with six competitive baselines on MIMIC dataset. 6 Conclusion We show that existing basis-based irregular time series forecasting methods can exhibit a non-vanishing asymptotic estimation bias when timestamp density is ignored, while their reliance on predefined basis functions may limit adaptability. To address these limitations, we propose the Debiased Neural Basis Function Network (DNBNet), a flexible framework that integrates density correction, learnable neural basis functions, multi-scale response extraction, and mass-aware fusion to enable robust and adaptive representation learning. DNBNet further employs a dual-branch decoder that combines implicit latent prediction with explicit basis-function reconstruction, thereby improving both predictive flexibility and reconstruction interpretability. Experiments on five real-world benchmarks demonstrate that DNBNet consistently achieves superior performance against strong baselines. References Afrifa-Yamoah et al. (2020) E. Afrifa-Yamoah, U. A. Mueller, S. M. Taylor, and A. J. Fisher Missing data imputation of high-resolution temporal climate time series data. Meteorological Applications 27 (1), p. e1873. Cited by: §1. Bilos et al. (2021) M. Bilos, J. Sommer, S. S. Rangapuram, T. Januschowski, and S. Günnemann Neural flows: efficient alternative to neural odes. In Advances in Neural Information Processing Systems 34: Annual Conference on Neural Information Processing Systems 2021, NeurIPS 2021, December 6-14, 2021, virtual, p. 21325–21337. Cited by: §2.1, §5.1.2. Chowdhury et al. (2023) R. R. Chowdhury, J. Li, X. Zhang, D. Hong, R. K. Gupta, and J. Shang Primenet: pre-training for irregular multivariate time series. In Proceedings of the AAAI Conference on Artificial Intelligence, p. 7184–7192. Cited by: §C.2, §5.1.1. Decorte et al. (2024) T. Decorte, S. Mortier, J. J. Lembrechts, F. J. R. Meysman, S. Latré, E. Mannens, and T. Verdonck Missing value imputation of wireless sensor data for environmental monitoring. Sensors 24 (8), p. 2416. Cited by: §1. Fons et al. (2025) E. Fons, A. Sztrajman, Y. El-Laham, L. Ferrer, S. Vyetrenko, and M. Veloso LSCD: lomb-scargle conditioned diffusion for time series imputation. In Forty-second International Conference on Machine Learning, ICML 2025, Vancouver, BC, Canada, July 13-19, 2025, Vol. 267. Cited by: Table 1, §2.2. Ghassemi et al. (2015) M. Ghassemi, M. A. F. Pimentel, T. Naumann, T. Brennan, D. A. Clifton, P. Szolovits, and M. Feng A multivariate timeseries modeling approach to severity of illness assessment and forecasting in ICU with sparse, heterogeneous clinical data. In Proceedings of the Twenty-Ninth AAAI Conference on Artificial Intelligence, January 25-30, 2015, Austin, Texas, USA, p. 446–453. Cited by: §1. Horn et al. (2020) M. Horn, M. Moor, C. Bock, B. Rieck, and K. Borgwardt Set functions for time series. In International Conference on Machine Learning, p. 4353–4363. Cited by: §C.2, §2.1, §5.1.1. Johnson et al. (2016) A. E. Johnson, T. J. Pollard, L. Shen, L. H. Lehman, M. Feng, M. Ghassemi, B. Moody, P. Szolovits, L. Anthony Celi, and R. G. Mark MIMIC-i, a freely accessible critical care database. Scientific data 3 (1), p. 1–9. Cited by: §C.1. Lee et al. (2026) S. Lee, K. Min, Y. Son, and H. Do Adaptive time encoding for irregular multivariate time-series classification. Advances in neural information processing systems 38, p. 101795–101824. Cited by: Table 1, §1. Li et al. (2025) B. Li, Y. Luo, Z. Liu, J. Zheng, J. Lv, and Q. Ma HyperIMTS: hypergraph neural network for irregular multivariate time series forecasting. In Forty-second International Conference on Machine Learning, ICML 2025, Vancouver, BC, Canada, July 13-19, 2025, Cited by: §2.1. Liu et al. (2026) X. Liu, X. Qiu, X. Wu, Z. Li, C. Guo, J. Hu, and B. Yang Rethinking irregular time series forecasting: a simple yet effective baseline. In Proceedings of the AAAI Conference on Artificial Intelligence, p. 23873–23881. Cited by: §C.2, Table 1, §1, §2.1, §3.2, §5.1.1. Luo et al. (2025) Y. Luo, B. Zhang, Z. Liu, and Q. Ma Hi-patch: hierarchical patch GNN for irregular multivariate time series. In Forty-second International Conference on Machine Learning, ICML 2025, Vancouver, BC, Canada, July 13-19, 2025, Cited by: §C.2, §2.1, §5.1.1, §5.1.2. Menne et al. (2015) M. J. Menne, C. N. Williams Jr, and R. S. Vose Long-term daily and monthly climate records from stations across the contiguous united states (us historical climatology network). Technical report Environmental System Science Data Infrastructure for a Virtual Ecosystem; CDIAC. Cited by: §C.1. Mercatali et al. (2024) G. Mercatali, A. Freitas, and J. Chen Graph neural flows for unveiling systemic interactions among irregularly sampled time series. In Advances in Neural Information Processing Systems 37: Annual Conference on Neural Information Processing Systems 2024, NeurIPS 2024, Vancouver, BC, Canada, December 10 - 15, 2024, Cited by: §2.1, §5.1.1. Nepal et al. (2024) S. Nepal, W. Liu, A. Pillai, W. Wang, V. Vojdanovski, J. F. Huckins, C. Rogers, M. L. Meyer, and A. T. Campbell Capturing the college experience: a four-year mobile sensing study of mental health, resilience and behavior of college students during the pandemic. Proceedings of the ACM on interactive, mobile, wearable and ubiquitous technologies 8 (1), p. 1–37. Cited by: §C.1. Oh et al. (2026) Y. Oh, D. Lim, and S. Kim FlowPath: learning data-driven manifolds with invertible flows for robust irregularly-sampled time series classification. In Proceedings of the AAAI Conference on Artificial Intelligence, p. 24594–24603. Cited by: §2.1. Oreshkin et al. (2020) B. N. Oreshkin, D. Carpov, N. Chapados, and Y. Bengio N-BEATS: neural basis expansion analysis for interpretable time series forecasting. In 8th International Conference on Learning Representations, ICLR 2020, Addis Ababa, Ethiopia, April 26-30, 2020, Cited by: §2.2. Qiu et al. (2025) X. Qiu, X. Wu, Y. Lin, C. Guo, J. Hu, and B. Yang DUET: dual clustering enhanced multivariate time series forecasting. In Proceedings of the 31st ACM SIGKDD Conference on Knowledge Discovery and Data Mining, V.1, KDD 2025, Toronto, ON, Canada, August 3-7, 2025, p. 1185–1196. Cited by: §2.2. Qiu et al. (2026) X. Qiu, K. Yan, X. Liu, X. Wu, and J. Hu Bridging time and frequency: a joint modeling framework for irregular multivariate time series forecasting. arXiv preprint arXiv:2602.00582. Cited by: Table 1, §1, §1. Rubanova et al. (2019) Y. Rubanova, R. T. Chen, and D. K. Duvenaud Latent ordinary differential equations for irregularly-sampled time series. Advances in neural information processing systems 32. Cited by: §2.1. Schirmer et al. (2022) M. Schirmer, M. Eltayeb, S. Lessmann, and M. Rudolph Modeling irregular time series with continuous recurrent units. In International conference on machine learning, p. 19388–19405. Cited by: §C.2, §5.1.1. Shukla and Marlin (2021) S. N. Shukla and B. M. Marlin Multi-time attention networks for irregularly sampled time series. In 9th International Conference on Learning Representations, ICLR 2021, Virtual Event, Austria, May 3-7, 2021, Cited by: §C.2, Table 1, §1, §2.2, §3.2, §5.1.1. Silva et al. (2012) I. Silva, G. Moody, D. J. Scott, L. A. Celi, and R. G. Mark Predicting in-hospital mortality of icu patients: the physionet/computing in cardiology challenge 2012. In 2012 computing in cardiology, p. 245–248. Cited by: §C.1. Yalavarthi et al. (2024) V. K. Yalavarthi, K. Madhusudhanan, R. Scholz, N. Ahmed, J. Burchert, S. Jawed, S. Born, and L. Schmidt-Thieme Grafiti: graphs for forecasting irregularly sampled time series. In Proceedings of the AAAI Conference on Artificial Intelligence, p. 16255–16263. Cited by: §C.2, §2.1, §5.1.1. Zhang et al. (2023) J. Zhang, S. Zheng, W. Cao, J. Bian, and J. Li Warpformer: a multi-scale modeling approach for irregular clinical time series. In Proceedings of the 29th ACM SIGKDD Conference on Knowledge Discovery and Data Mining, p. 3273–3285. Cited by: §C.2, §5.1.1. Zhang et al. (2024) W. Zhang, C. Yin, H. Liu, X. Zhou, and H. Xiong Irregular multivariate time series forecasting: A transformable patching graph neural networks approach. In Forty-first International Conference on Machine Learning, ICML 2024, Vienna, Austria, July 21-27, 2024, R. Salakhutdinov, Z. Kolter, K. A. Heller, A. Weller, N. Oliver, J. Scarlett, and F. Berkenkamp (Eds.), p. 60179–60196. Cited by: §C.2, §2.1, §4.4, §5.1.1. Zhang et al. (2022) X. Zhang, M. Zeman, T. Tsiligkaridis, and M. Zitnik Graph-guided network for irregularly sampled multivariate time series. In The Tenth International Conference on Learning Representations, ICLR 2022, Virtual Event, April 25-29, 2022, Cited by: §C.2, §2.1, §5.1.1. Zhou et al. (2022) T. Zhou, Z. Ma, Q. Wen, X. Wang, L. Sun, and R. Jin FEDformer: frequency enhanced decomposed transformer for long-term series forecasting. In International Conference on Machine Learning, ICML 2022, 17-23 July 2022, Baltimore, Maryland, USA, p. 27268–27286. Cited by: §2.2. Zhou et al. (2026) Z. Zhou, Y. Huang, Y. Wang, Y. Wu, J. Kwok, and Y. Liang Revitalizing canonical pre-alignment for irregular multivariate time series forecasting. In Proceedings of the AAAI Conference on Artificial Intelligence, p. 29115–29123. Cited by: §C.2, Table 1, §1, §1, §2.2, §3.2, §5.1.1. Appendix A Proof of Theorem 1 For each basis function ϕk _k, the ideal continuous-time response coefficient is defined as ck=∫x(t)ϕk(t)t∫ϕk(t)t. c_k= _Tx(t) _k(t)dt _T _k(t)dt. (20) The uncorrected discrete estimator is defined as c^k=∑i=1Lxiϕk(ti)∑i=1Lϕk(ti). c_k= _i=1^Lx_i _k(t_i) _i=1^L _k(t_i). (21) Since the timestamps are sampled from density p(t)p(t), the discrete estimator does not directly approximate ckc_k. Instead, it approximates the response coefficient under the timestamp sampling density: ckp=∫x(t)ϕk(t)p(t)t∫ϕk(t)p(t)t. c_k^p= _Tx(t) _k(t)p(t)dt _T _k(t)p(t)dt. (22) Therefore, we decompose the error as |c^k−ck|≤|c^k−ckp|+|ckp−ck|. | c_k-c_k|≤| c_k-c_k^p|+|c_k^p-c_k|. (23) The first term corresponds to the finite-sample approximation error, while the second term corresponds to the sampling bias induced by the non-uniform timestamp density. We first bound the finite-sample error. Define Yik=x(ti)ϕk(ti),Zik=ϕk(ti). Y_i^k=x(t_i) _k(t_i), Z_i^k= _k(t_i). (24) Then c^k=Y¯kZ¯k,ckp=p[Yik]p[Zik], c_k= Y_k Z_k, c_k^p= E_p[Y_i^k]E_p[Z_i^k], (25) where Y¯k=1L∑i=1LYik,Z¯k=1L∑i=1LZik. Y_k= 1L _i=1^LY_i^k, Z_k= 1L _i=1^LZ_i^k. (26) Since |x(t)|≤1|x(t)|≤ 1 and 0≤ϕk(t)≤m0≤ _k(t)≤ m, we have |Yik|≤m,0≤Zik≤m. |Y_i^k|≤ m, 0≤ Z_i^k≤ m. (27) By Hoeffding’s inequality and a union bound over all k∈[K]k∈[K] for both Y¯k Y_k and Z¯k Z_k, with probability at least 1−δ1-δ, the following inequalities hold simultaneously: maxk∈[K]|Y¯k−p[Yik]|≤ϵL,maxk∈[K]|Z¯k−p[Zik]|≤ϵL, _k∈[K]| Y_k-E_p[Y_i^k]|≤ _L, _k∈[K]| Z_k-E_p[Z_i^k]|≤ _L, (28) where ϵL=m2log(4K/δ)L. _L=m 2 (4K/δ)L. (29) When L is sufficiently large such that ϵL≤a/2 _L≤ a/2, i.e., L≥8m2a2log4KδL≥ 8m^2a^2 4Kδ, we have Z¯k≥p[Zik]−ϵL≥a2. Z_k _p[Z_i^k]- _L≥ a2. (30) Now we compare the two ratios: |c^k−ckp| | c_k-c_k^p| =|Y¯kZ¯k−p[Yik]p[Zik]| = | Y_k Z_k- E_p[Y_i^k]E_p[Z_i^k] | ≤|Y¯k−p[Yik]|Z¯k+|p[Yik]||1Z¯k−1p[Zik]|. ≤ | Y_k-E_p[Y_i^k]| Z_k+|E_p[Y_i^k]| | 1 Z_k- 1E_p[Z_i^k] |. ≤|Y¯k−p[Yik]|Z¯k+p[Zik]|Z¯k−p[Zik]|Z¯kp[Zik] ≤ | Y_k-E_p[Y_i^k]| Z_k+E_p[Z_i^k] | Z_k-E_p[Z_i^k]| Z_kE_p[Z_i^k] ≤2ϵLa+2ϵLa ≤ 2 _La+ 2 _La =4ϵLa. = 4 _La. (31) The second inequality holds because |p[Yik]|=|∫x(t)ϕk(t)p(t)t|≤∫ϕk(t)p(t)t=p[Zik]|E_p[Y_i^k]|= | _Tx(t) _k(t)p(t)dt |≤ _T _k(t)p(t)dt=E_p[Z_i^k]. Taking the maximum over all k∈[K]k∈[K], we obtain maxk∈[K]|c^k−ckp|≤4ma2log(4K/δ)L. _k∈[K]| c_k-c_k^p|≤ 4ma 2 (4K/δ)L. (32) Next, we bound the sampling bias term |ckp−ck||c_k^p-c_k|. Define Ak=∫ϕk(t)t,Bk=∫x(t)ϕk(t)t, A_k= _T _k(t)dt, B_k= _Tx(t) _k(t)dt, (33) and Akp=∫ϕk(t)p(t)t,Bkp=∫x(t)ϕk(t)p(t)t. A_k^p= _T _k(t)p(t)dt, B_k^p= _Tx(t) _k(t)p(t)dt. (34) Then ck=BkAk,ckp=BkpAkp. c_k= B_kA_k, c_k^p= B_k^pA_k^p. (35) Thus, |ckp−ck| |c_k^p-c_k| =|BkpAkp−BkAk| = | B_k^pA_k^p- B_kA_k | ≤|Bkp−Bk|Akp+|Bk||1Akp−1Ak| ≤ |B_k^p-B_k|A_k^p+|B_k| | 1A_k^p- 1A_k | ≤|Bkp−Bk|Akp+Ak|Akp−Ak|AkpAk ≤ |B_k^p-B_k|A_k^p+A_k |A_k^p-A_k|A_k^pA_k ≤|Bkp−Bk|+|Akp−Ak|a ≤ |B_k^p-B_k|+|A_k^p-A_k|a =|∫x(t)ϕk(t)(p(t)−1)t|+|∫ϕk(t)(p(t)−1)t|a = | _Tx(t) _k(t)(p(t)-1)dt |+ | _T _k(t)(p(t)-1)dt |a ≤∫|x(t)|ϕk(t)|p(t)−1|t+∫ϕk(t)|p(t)−1|ta ≤ _T|x(t)| _k(t)|p(t)-1|dt+ _T _k(t)|p(t)-1|dta ≤2m∫|p(t)−1|ta ≤ 2m _T|p(t)-1|dta =2m‖p(t)−1‖L1a = 2m\|p(t)-1\|_L^1a (36) The second inequality holds due to |Bk|=|∫x(t)ϕk(t)t|≤∫ϕk(t)t=Ak|B_k|= | _Tx(t) _k(t)dt |≤ _T _k(t)dt=A_k. Taking the maximum over k∈[K]k∈[K], we obtain maxk∈[K]|ckp−ck|≤2ma‖p(t)−1‖L1. _k∈[K]|c_k^p-c_k|≤ 2ma\|p(t)-1\|_L^1. (37) Finally, combining the finite-sample approximation error in Eq. (32) and the sampling bias term in Eq. (37) yields maxk∈[K]|c^k−ck| _k∈[K]| c_k-c_k| ≤maxk∈[K]|c^k−ckp|+maxk∈[K]|ckp−ck| ≤ _k∈[K]| c_k-c_k^p|+ _k∈[K]|c_k^p-c_k| ≤4ma2log(4K/δ)L+2ma‖p(t)−1‖L1. ≤ 4ma 2 (4K/δ)L+ 2ma\|p(t)-1\|_L^1. (38) This completes the proof. Appendix B Debiasing Property of Density-Corrected Estimation In this section, we further explain why the density-corrected estimator can eliminate the sampling bias. Recall when the observed timestamps are sampled from a non-uniform density p(t)p(t), the uncorrected discrete estimator converges to the response coefficient under the timestamp sampling density, which leads to the non-vanishing bias term in Theorem 1. To remove this bias, we apply importance weighting to compensate for the timestamp sampling density. The density-corrected estimator is defined as c~k=∑i=1Lxiϕk(ti)p(ti)∑i=1Lϕk(ti)p(ti). c_k= _i=1^L x_i _k(t_i)p(t_i) _i=1^L _k(t_i)p(t_i). (39) The key difference from the uncorrected estimator is that the numerator and denominator now directly estimate the corresponding continuous-time integrals under the uniform time domain. Specifically, since tit_i is sampled from density p(t)p(t), we have p[x(t)ϕk(t)p(t)] _p [ x(t) _k(t)p(t) ] =∫x(t)ϕk(t)p(t)p(t)t=∫x(t)ϕk(t)t, = _T x(t) _k(t)p(t)p(t)dt= _Tx(t) _k(t)dt, (40) p[ϕk(t)p(t)] _p [ _k(t)p(t) ] =∫ϕk(t)p(t)p(t)t=∫ϕk(t)t. = _T _k(t)p(t)p(t)dt= _T _k(t)dt. (41) Therefore, the density-corrected estimator directly approximates the desired coefficient ckc_k, rather than the sampling-density-weighted coefficient ckpc_k^p. To obtain a finite-sample bound, we additionally assume that the sampling density is lower bounded, i.e., there exists a constant ρ>0ρ>0 such that p(t)≥ρp(t)≥ρ, for ∀t∈∀ t . This condition ensures that the importance weights 1/p(t)1/p(t) are bounded. Under this assumption, together with |x(t)|≤1|x(t)|≤ 1 and 0≤ϕk(t)≤m0≤ _k(t)≤ m, we have |x(ti)ϕk(ti)p(ti)|≤mρ,0≤ϕk(ti)p(ti)≤mρ. | x(t_i) _k(t_i)p(t_i) |≤ mρ, 0≤ _k(t_i)p(t_i)≤ mρ. (42) Following the same concentration and ratio-estimation arguments as in Appendix A, when L is sufficiently large, with probability at least 1−δ1-δ, the following bound holds: maxk∈[K]|c~k−ck|≤4maρ2log(4K/δ)L. _k∈[K]| c_k-c_k|≤ 4maρ 2 (4K/δ)L. (43) Compared with the bound for the uncorrected estimator in Theorem 1, the above bound does not contain the non-vanishing bias term related to ‖p(t)−1‖L1\|p(t)-1\|_L^1. Therefore, the density-corrected estimator removes the systematic bias induced by irregular timestamp sampling. Appendix C Datasets and Baseline Model Details C.1 Datasets In this section, we provide a detailed description of the five public datasets used in our experiments, including their sources and characteristics. The basic statistics are summarized in Table 5. Dataset # Variables # samples Avg # Obs. Max Length PhysioNet 36 11,981 308.6 47 MIMIC 96 21,250 144.6 96 HumanActivity 12 1,359 362.2 131 StudentLife 9 20 7,680.5 1,370 USHCN 5 1,114 313.5 337 Table 5: Dataset statistics. Table 6: MAE results of DNBNet and baselines on five irregular time-series datasets. Best results are in bold and second-best results are underlined. Method USHCN Human Activity Student Life PhysioNet MIMIC PrimeNet 0.5094±0.00020.5094± 0.0002 1.7054±0.00051.7054± 0.0005 0.7231±0.00010.7231± 0.0001 0.6859±0.00000.6859± 0.0000 0.6632±0.00010.6632± 0.0001 SeFT 0.4981±0.00210.4981± 0.0021 0.9733±0.00060.9733± 0.0006 0.7201±0.00220.7201± 0.0022 0.6766±0.00050.6766± 0.0005 0.6689±0.00080.6689± 0.0008 mTAN 0.3412±0.00760.3412± 0.0076 0.2302±0.00710.2302± 0.0071 0.5742±0.00330.5742± 0.0033 0.4271±0.00410.4271± 0.0041 0.6749±0.00230.6749± 0.0023 CRU 0.3948±0.00690.3948± 0.0069 0.2673±0.00260.2673± 0.0026 0.6403±0.00430.6403± 0.0043 0.5782±0.00400.5782± 0.0040 0.5828±0.00460.5828± 0.0046 GNeuralFlow 0.3842±0.00240.3842± 0.0024 0.3802±0.02170.3802± 0.0217 0.6085±0.01630.6085± 0.0163 0.4012±0.00090.4012± 0.0009 0.5336±0.00580.5336± 0.0058 Raindrop 0.3515±0.00700.3515± 0.0070 0.2159±0.00270.2159± 0.0027 0.5816±0.00720.5816± 0.0072 0.4057±0.00060.4057± 0.0006 0.5047±0.00350.5047± 0.0035 tPatchGNN 0.4105±0.07250.4105± 0.0725 0.1444±0.00200.1444± 0.0020 0.5603±0.00090.5603± 0.0009 0.3665±0.00180.3665± 0.0018 0.4069±0.00150.4069± 0.0015 GraFITi 0.3437±0.01050.3437± 0.0105 0.1559±0.00560.1559± 0.0056 0.5644±0.00050.5644± 0.0005 0.3538±0.00150.3538± 0.0015 0.4016±0.0040¯ 0.4016± 0.0040 Warpformer 0.3215±0.0017¯ 0.3215± 0.0017 0.1427±0.00050.1427± 0.0005 0.5510±0.00160.5510± 0.0016 0.3718±0.00130.3718± 0.0013 0.4069±0.00140.4069± 0.0014 Hi-Patch 0.3364±0.01140.3364± 0.0114 0.1450±0.00080.1450± 0.0008 0.5580±0.00110.5580± 0.0011 0.3715±0.00360.3715± 0.0036 0.4112±0.00220.4112± 0.0022 KAFNet 0.3219±0.00280.3219± 0.0028 0.1501±0.00130.1501± 0.0013 0.5674±0.00110.5674± 0.0011 0.3809±0.00080.3809± 0.0008 0.4098±0.00220.4098± 0.0022 APN 0.3368±0.01470.3368± 0.0147 0.1401±0.0007¯ 0.1401± 0.0007 0.5682±0.00050.5682± 0.0005 0.3669±0.00020.3669± 0.0002 0.4145±0.00560.4145± 0.0056 DNBNet (Ours) 0.3053±0.00460.3053± 0.0046 0.1376±0.00030.1376± 0.0003 0.5515±0.0003¯ 0.5515± 0.0003 0.3644±0.0011¯ 0.3644± 0.0011 0.3981±0.00090.3981± 0.0009 Table 7: Complete ablation study of DNBNet on five irregular time-series datasets. Best results are highlighted in bold. MSE Variant USHCN Human Activity Student Life PhysioNet MIMIC w/o p(t)p(t) 0.4449±0.02910.4449± 0.0291 0.0564±0.00000.0564± 0.0000 0.6291±0.00040.6291± 0.0004 0.3161±0.00030.3161± 0.0003 0.4485±0.00120.4485± 0.0012 Rep. RBF 0.4144±0.00910.4144± 0.0091 0.0568±0.00020.0568± 0.0002 0.6242±0.00070.6242± 0.0007 0.3231±0.00070.3231± 0.0007 0.4507±0.00090.4507± 0.0009 Rep. Fourier 0.4101±0.00880.4101± 0.0088 0.0562±0.00010.0562± 0.0001 0.6257±0.00090.6257± 0.0009 0.3267±0.00020.3267± 0.0002 0.4869±0.00480.4869± 0.0048 w/o AvgPool 0.4195±0.01660.4195± 0.0166 0.0569±0.00010.0569± 0.0001 0.6269±0.00040.6269± 0.0004 0.3119±0.00060.3119± 0.0006 0.4494±0.00120.4494± 0.0012 w/o BasDec 0.4552±0.04240.4552± 0.0424 0.0559±0.00020.0559± 0.0002 0.6281±0.00110.6281± 0.0011 0.3092±0.00020.3092± 0.0002 0.4613±0.00230.4613± 0.0023 DNBNet (Full) 0.4013±0.01080.4013± 0.0108 0.0551±0.00010.0551± 0.0001 0.6136±0.00020.6136± 0.0002 0.3089±0.00050.3089± 0.0005 0.4289±0.00120.4289± 0.0012 MAE Variant USHCN Human Activity Student Life PhysioNet MIMIC w/o p(t)p(t) 0.3276±0.02000.3276± 0.0200 0.1386±0.00060.1386± 0.0006 0.5575±0.00050.5575± 0.0005 0.3694±0.00020.3694± 0.0002 0.4078±0.00140.4078± 0.0014 Rep. RBF 0.3079±0.00640.3079± 0.0064 0.1382±0.00030.1382± 0.0003 0.5549±0.00140.5549± 0.0014 0.3798±0.00120.3798± 0.0012 0.4116±0.00090.4116± 0.0009 Rep. Fourier 0.3139±0.00520.3139± 0.0052 0.1399±0.00050.1399± 0.0005 0.5571±0.00110.5571± 0.0011 0.3794±0.00010.3794± 0.0001 0.4239±0.00330.4239± 0.0033 w/o AvgPool 0.3147±0.01100.3147± 0.0110 0.1388±0.00030.1388± 0.0003 0.5569±0.00050.5569± 0.0005 0.3686±0.00100.3686± 0.0010 0.4080±0.00090.4080± 0.0009 w/o BasDec 0.3395±0.02780.3395± 0.0278 0.1387±0.00060.1387± 0.0006 0.5565±0.00110.5565± 0.0011 0.3652±0.00130.3652± 0.0013 0.4169±0.00220.4169± 0.0022 DNBNet (Full) 0.3053±0.00460.3053± 0.0046 0.1376±0.00030.1376± 0.0003 0.5515±0.00030.5515± 0.0003 0.3644±0.00110.3644± 0.0011 0.3981±0.00090.3981± 0.0009 MIMIC. MIMIC is a large, freely accessible critical care database (8). It contains de-identified health records from patients who stayed in intensive care units (ICUs). The dataset is highly detailed, including vital signs, medications, and laboratory measurements. In our experiments, we use the clinical time series collected from the first 48 hours of each patient’s ICU stay. The resulting benchmark contains 21,250 samples with 96 variables. PhysioNet. PhysioNet is another widely used clinical benchmark for irregular time series analysis (23). It was originally released for predicting in-hospital mortality of ICU patients. Each record consists of multivariate measurements collected during the first 48 hours after admission. The processed benchmark used in our experiments contains 11,981 samples and 36 variables, such as serum glucose and heart rate. Human Activity. Human Activity is a biomechanics-related dataset from the UCI Machine Learning Repository. It contains sensor measurements collected from five subjects performing various activities. The benchmark includes 1,359 samples and 12 variables representing 3D positions. The data are naturally irregular because the sensors record information at non-uniform time intervals. StudentLife. StudentLife is a human behavior sensing dataset collected from college students using smartphones and wearable devices (15). In our setting, we use its numerical time series modality, which includes activity-related features such as biking, walking, and sleeping. The irregularity mainly arises from human scheduling, device usage patterns, and asynchronous sensor recording, making it a realistic benchmark for irregular multivariate time series forecasting. USHCN. USHCN (United States Historical Climatology Network) is a climate dataset containing long-term meteorological observations from weather stations across the United States (13). The processed benchmark contains 1,114 samples and 5 variables, such as daily maximum temperature, daily minimum temperature, and precipitation. Although the data are recorded daily, missing observations are common in practice, which makes this dataset suitable for irregular time series forecasting. Following previous work, we use a subset covering the period from 1996 to 2000. C.2 Baseline Model Details We compare our method with 12 representative baseline models, covering irregular time series forecasting, classification, and interpolation methods. These baselines include PrimeNet, SeFT, mTAN, CRU, GNeuralFlow, Raindrop, tPatchGNN, GraFITi, Warpformer, Hi-Patch, KAFNet, and APN. PrimeNet (3) is a pretraining framework for irregular multivariate time series, which learns generalizable temporal representations from irregular observations and can be adapted to downstream forecasting tasks. SeFT (7) is a set-function-based model that treats irregular observations as an unordered set and uses permutation-invariant aggregation to learn representations without requiring time alignment. mTAN (22) employs multi-time attention to map irregular observations onto a set of reference time points, thereby obtaining a regularized latent representation for downstream prediction. CRU (21) is a continuous-time recurrent model that combines recurrent updates with latent uncertainty modeling, making it suitable for irregularly sampled observations with varying time gaps. GNeuralFlow (mercatali2024graph) extends neural flow modeling with graph structures to capture both continuous-time dynamics and inter-variable dependencies in irregular multivariate time series. Raindrop (27) is a graph attention model for IMTS, which propagates information across variables and timestamps to handle missingness and asynchronous observations. tPatchGNN (26) is a patch-based graph neural network that partitions irregular sequences into temporal patches and models temporal and variable dependencies jointly. GraFITi (24) represents irregular time series as a bipartite graph over observations, timestamps, and variables, allowing the model to operate directly on sparse observations without dense alignment. Warpformer (25) is a Transformer-based model that introduces a temporal warping mechanism for multi-scale modeling, aiming to better capture irregular and non-stationary temporal patterns. Hi-Patch (12) adopts a hierarchical patching strategy to model irregular time series at multiple temporal granularities, which helps capture both local dynamics and long-range dependencies. KAFNet (29) is a recent IMTS forecasting model that uses kernel-based aggregation to encode irregular observations into compact latent representations for future prediction. APN (11) is a recent basis-function-based forecasting method that models irregular observations through adaptive temporal aggregation, providing a strong and efficient baseline for IMTS forecasting. Appendix D Additional Experiments Results D.1 Complete MAE Results Due to space limitations, we report only the MSE results in the main paper and provide the complete MAE results in Table 6. DNBNet achieves the best MAE on USHCN, Human Activity, and MIMIC, and obtains the second-best results on Student Life and PhysioNet, which is generally consistent with the MSE results. This demonstrates that DNBNet remains effective under different evaluation metrics and across diverse application domains. D.2 Varying Lookback Lengths and Forecast Horizons To further evaluate the robustness of DNBNet under different temporal settings, we conduct additional experiments by varying the lookback lengths and forecast horizons. As shown in Fig. 5, DNBNet consistently achieves stable and competitive performance across different lookback windows. This indicates that the proposed model does not rely on a specific input length, and can effectively extract useful temporal patterns from both shorter and longer historical observations. We also report the results under different forecast horizons in Table 8. The lookback lengths follow the settings in Table 5 of main paper, while the prediction horizons are set to the remaining length of each sequence, including 12 hours for MIMIC-I and PhysioNet, 300 milliseconds for Human Activity, the next 3 observations for USHCN, and 10 days for Student Life. The results show that DNBNet maintains strong performance under these horizon settings, further demonstrating the effectiveness of the debiased neural basis-function response and the dual-branch forecasting decoder for irregular time-series forecasting. Table 8: Experimental results on five irregular multivariate time-series datasets evaluated by MSE and MAE. The look-back lengths follow Table 2, and the prediction horizons are set to the remaining length of each sequence: 12 hours for MIMIC-I and PhysioNet, 300 milliseconds for Human Activity, the next 3 observations for USHCN, and 10 days for Student Life. MSE Method USHCN Human Activity Student Life PhysioNet MIMIC GraFTi 0.2089±0.02180.2089±0.0218 0.0431±0.00060.0431±0.0006 0.5789±0.00180.5789±0.0018 0.3553±0.00160.3553 0.0016 0.4957±0.00320.4957±0.0032 Warpformer 0.1738±0.00690.1738±0.0069 0.0446±0.00130.0446±0.0013 0.5861±0.00320.5861±0.0032 0.3558±0.00160.3558±0.0016 0.5295±0.00350.5295±0.0035 Hi-Patch 0.2416±0.02050.2416±0.0205 0.0472±0.00010.0472±0.0001 0.5902±0.00200.5902±0.0020 0.3666±0.00010.3666±0.0001 0.5352±0.00620.5352±0.0062 KAFNet 0.2045±0.00910.2045±0.0091 0.0677±0.00060.0677±0.0006 0.5980±0.00240.5980±0.0024 0.3647±0.00100.3647±0.0010 0.5145±0.00760.5145±0.0076 APN 0.1816±0.01330.1816±0.0133 0.0436±0.00040.0436±0.0004 0.6057±0.01400.6057±0.0140 0.3689±0.00260.3689±0.0026 0.5396±0.00530.5396±0.0053 NBFN (Ours) 0.1730±0.00970.1730 0.0097 0.0429±0.00010.0429 0.0001 0.5763±0.00370.5763 0.0037 0.3584±0.00060.3584±0.0006 0.4902±0.00850.4902 0.0085 MAE Method USHCN Human Activity Student Life PhysioNet MIMIC GraFTi 0.2813±0.00310.2813±0.0031 0.1203±0.00300.1203±0.0030 0.5308±0.00120.5308±0.0012 0.4051±0.00030.4051 0.0003 0.4369±0.00230.4369±0.0023 Warpformer 0.2682±0.00260.2682±0.0026 0.1221±0.00220.1221±0.0022 0.5317±0.00150.5317±0.0015 0.4063±0.00190.4063±0.0019 0.4412±0.00230.4412±0.0023 Hi-Patch 0.2923±0.01050.2923±0.0105 0.1250±0.00010.1250±0.0001 0.5318±0.00110.5318±0.0011 0.4150±0.00110.4150±0.0011 0.4544±0.00620.4544±0.0062 KAFNet 0.2714±0.01350.2714±0.0135 0.1582±0.00150.1582±0.0015 0.5386±0.00160.5386±0.0016 0.4138±0.00140.4138±0.0014 0.4419±0.00460.4419±0.0046 APN 0.2818±0.00980.2818±0.0098 0.1211±0.00110.1211±0.0011 0.5480±0.01560.5480±0.0156 0.4153±0.00270.4153±0.0027 0.4564±0.00440.4564±0.0044 NBFN (Ours) 0.2518±0.00730.2518 0.0073 0.1172±0.00060.1172 0.0006 0.5291±0.00220.5291 0.0022 0.4084±0.00120.4084±0.0012 0.4338±0.00670.4338 0.0067 Figure 5: Forecasting performance with varying lookback lengths and fixed forecast horizons. D.3 Complete Ablation Results We further report the complete ablation results on all five datasets in Tables 7. The full DNBNet consistently outperforms its variants, verifying the contribution of each component. In particular, removing the density correction term p(t)p(t), replacing learnable basis functions with predefined RBF or Fourier bases, or removing the average pooling module or basis-function decoder all lead to performance degradation. These results confirm the effectiveness of debiased response estimation, learnable basis functions, multi-scale representation, and basis-function-based prediction. langley00