Cross-validation Bandwidth Matrices for Multivariate Kernel Density Estimation

Abstract.  The performance of multivariate kernel density estimates depends crucially on the choice of bandwidth matrix, but progress towards developing good bandwidth matrix selectors has been relatively slow. In particular, previous studies of cross-validation (CV) methods have been restricted to biased and unbiased CV selection of diagonal bandwidth matrices. However, for certain types of target density the use of full (i.e. unconstrained) bandwidth matrices offers the potential for significantly improved density estimation. In this paper, we generalize earlier work from diagonal to full bandwidth matrices, and develop a smooth cross-validation (SCV) methodology for multivariate data. We consider optimization of the SCV technique with respect to a pilot bandwidth matrix. All the CV methods are studied using asymptotic analysis, simulation experiments and real data analysis. The results suggest that SCV for full bandwidth matrices is the most reliable of the CV methods. We also observe that experience from the univariate setting can sometimes be a misleading guide for understanding bandwidth selection in the multivariate case.     

Goodness-of-fit Tests Based on the Kernel Density Estimator

Abstract.  Given an i.i.d. sample drawn from a density f on the real line, the problem of testing whether f is in a given class of densities is considered. Testing procedures constructed on the basis of minimizing the L ₁-distance between a kernel density estimate and any density in the hypothesized class are investigated. General non-asymptotic bounds are derived for the power of the test. It is shown that the concentration of the data-dependent smoothing factor and the 'size' of the hypothesized class of densities play a key role in the performance of the test. Consistency and non-asymptotic performance bounds are established in several special cases, including testing simple hypotheses, translation/scale classes and symmetry. Simulations are also carried out to compare the behaviour of the method with the Kolmogorov-Smirnov test and an L ₂ density-based approach due to Fan [ Econ. Theory 10 (1994) 316].     

Kernel spline regression

The authors propose «kernel spline regression,» a method of combining spline regression and kernel smoothing by replacing the polynomial approximation for local polynomial kernel regression with the spline basis. The new approach retains the local weighting scheme and the use of a bandwidth to control the size of local neighborhood. The authors compute the bias and variance of the kernel linear spline estimator, which they compare with local linear regression. They show that kernel spline estimators can succeed in capturing the main features of the underlying curve more effectively than local polynomial regression when the curvature changes rapidly. They also show through simulation that kernel spline regression often performs better than ordinary spline regression and local polynomial regression.     

DDS阵列频率源技术研究

杂散抑制差和输出频率低是DDS的两大缺点,并且随着输出带宽的增加杂散性能更加恶化。传统的解决方法是采用更好的算法和改进DDS芯片本身的结构。该文在分析了DDS杂散性能和输出带宽相互制约的DDS阵列方法基础上,尝试了一种新的解决办法棗DDS阵列方法。同时还分析了这种方法的性能,并通过实验验证了其可行性。     

Data Sharpening for Hazard Rate Estimation

Data sharpening is a general tool for enhancing the performance of statistical estimators, by altering the data before substituting them into conventional methods. In one of the simplest forms of data sharpening, available for curve estimation, an explicit empirical transformation is used to alter the data. The attraction of this approach is diminished, however, if the formula has to be altered for each different application. For example, one could expect the formula for use in hazard rate estimation to differ from that for straight density estimation, since a hazard rate is a ratio–type functional of a density. This paper shows that, in fact, identical data transformations can be used in each case, regardless of whether the data involve censoring. This dramatically simplifies the application of data sharpening to problems involving hazard rate estimation, and makes data sharpening attractive.     

IPv6可用带宽测量方法的设计与实现

提出一种与现有带宽测量方法(包对法)不同的探测理论——基于链路状态的带宽测量方法。源端向目的端发送一系列小的探测报文,通过对探测报文的时延参数进行数理统计得到链路忙闲状态,进而获得链路可用带宽。采用此方法也可用于计算各种链路容量和链路空闲率。在实际测量中为解决探测报文与背景业务之间的相互影响,利用IPv6基本报头中流标签和业务类型字段设计专门的测试流和测试级。仿真证明方法正确。     

树的圈上带宽和下界

给出了一些以顶点数、直径或独立数表示的树的圈上带宽和的下界，并以此计算了K1，n(星）和Wn 1（轮）的圈上带宽和。     

Variable kernel density estimation

This paper considers the problem of selecting optimal bandwidths for variable (sample‐point adaptive) kernel density estimation. A data‐driven variable bandwidth selector is proposed, based on the idea of approximating the log‐bandwidth function by a cubic spline. This cubic spline is optimized with respect to a cross‐validation criterion. The proposed method can be interpreted as a selector for either integrated squared error (ISE) or mean integrated squared error (MISE) optimal bandwidths. This leads to reflection upon some of the differences between ISE and MISE as error criteria for variable kernel estimation. Results from simulation studies indicate that the proposed method outperforms a fixed kernel estimator (in terms of ISE) when the target density has a combination of sharp modes and regions of smooth undulation. Moreover, some detailed data analyses suggest that the gains in ISE may understate the improvements in visual appeal obtained using the proposed variable kernel estimator. These numerical studies also show that the proposed estimator outperforms existing variable kernel density estimators implemented using piecewise constant bandwidth functions.     

Large sample properties and confidence bands for component-wise varying-coefficient regression with longitudinal dependent variable

Longitudinal studies with repeatedly measured dependent variable (out-come) and time-invariant covariates are common in biomedical and epidemi-ological studies. A useful statistical tool to evaluate the effects of covariates on the outcome variable over time is the varying-coefficient regression, which considers a linear relationship between the covariates and the outcome at a specific time point but assumes the linear coefficients to be smooth curves over time. In order to provide adequate smoothing for each coefficient curve, Wu and Chiang ( 1999 ) proposed a class of component-wise kernel estimators and determined the large sample convergence rates and some of the constant terms of the mean squared errors of their estimators. In this paper we calcu¬late the explicit large sample mean squared errors, including the convergence rates and ail the constant terms, and the asymptotic distributions of the kernel estimators of Wu and Chiang ( 1999 ). These asymptotic distributions are used to construct point-wise confidence intervals and Bonferroni-type confidence bands for the coefficient curves. Through a Monte Carlo simulation, wre show that our confidence regions have adequate coverage probabilities. Applying our procedures to a NIH fetal growth study, we show that our procedures are useful to determine the effects of maternal height, cigarette smoking and al¬cohol consumption on the growth of fetal abdominal circumference over time during pregnancy.     

Spectral analysis with replicated time series

A doubly stochastic process {x(b,t);b?B,t?Z} is considered, with (B,β,P_β) being a probability space so that for each b, {X(b,t);t ? Z} is a stationary process with an absolutely continuous spectral distribution. The population spectrum is defined as f(ω) = E_B[Q(b,ω)] with Q(b,ω) being the spectral density function of X(b,t). The aim of this paper is to estimate f(ω) by means of a random sample b₁,…,b_r from (B,β,P_β). For each b₁? B, the processes X(b₁,t) are observed at the same times t=1,…,N. Thus, r time series (x(b₁,t)} are available in order to estimate f(ω). A model for each individual periodogram, which involves f(ω), is formulated. It has been proven that a certain family of linear stationary processes follows the above model In this context, a kernel estimator is proposed in order to estimate f(ω). The bias, variance and asymptotic distribution of this estimator are investigated under certain conditions.