纵向数据边际模型非参数光滑方法的比较

对于纵向数据边际模型的均值函数, 有很多非参数估计方法, 其中回归样条, 光滑样条, 似乎不相关(SUR)核估计等方法在工作协方差阵正确指定时具有最小的渐近方差. 回归样条的渐近偏差与工作协方差阵无关, 而SUR核估计和光滑样条估计的渐近偏差却依赖于工作协方差阵. 本文主要研究了回归样条, 光滑样条和SUR核估计的效率问题. 通过模拟比较发现回归样条估计的表现比较稳定, 在大多数情况下比光滑样条估计和SUR核估计的效率高.     

Spline Multiscale Smoothing to Control FDR for Exploring Features of Regression Curves

SiZer (significant zero crossing of the derivatives) is a multiscale smoothing method for exploring trends, maxima, and minima in data. In this article, a regression spline version of SiZer is proposed in a nonparametric regression setting by the fiducial method. The number of knots for spline interpolation is used as the scale parameter of the new SiZer, which controls the smoothness of estimate. In the construction of the new SiZer, multiple testing adjustment is made to control the row-wise false discovery rate (FDR) of SiZer. This adjustment is appealing for exploratory data analysis and has potential to increase the power. A special map is also produced on a continuous scale using p-values to assess the significance of features. Simulations and a real data application are carried out to investigate the performance of the proposed SiZer, in which several comparisons with other existing SiZers are presented. Supplementary materials for this article are available online.     

Iterative weighted partial spline least squares estimation in semiparametric modeling of longitudinal data

In this paper we consider the estimating problem of a semiparametric regression modelling whenthe data are longitudinal.An iterative weighted partial spline least squares estimator(IWPSLSE)for the para-metric component is proposed which is more efficient than the weighted partial spline least squares estimator(WPSLSE)with weights constructed by using the within-group partial spline least squares residuals in the sense     

Zur numerischen Lösung gewöhnlicher Differential-gleichungen mit Splines in einem Sonderfall

In an earlier paper [1] a general procedure has been presented to obtain polynomial spline approximations for the solution of the initial value problem for ordinary differential equations. In this paper the general procedure is described by an equivalent one step method. Furthermore two convergence theorems are proved for a special case which is not included in the general convergence or divergence theory given in [1].     

On the Construction of Optimal Monotone Cubic Spline Interpolations

In this paper we derive necessary optimality conditions for an interpolating spline function which minimizes the Holladay approximation of the energy functional and which stays monotone if the given interpolation data are monotone. To this end optimal control theory for state-restricted optimal control problems is applied. The necessary conditions yield a complete characterization of the optimal spline. In the case of two or three interpolation knots, which we call thelocalcase, the optimality conditions are treated analytically. They reduce to polynomial equations which can very easily be solved numerically. These results are used for the construction of a numerical algorithm for the optimal monotone spline in the general (global) case via Newton's method. Here, the local optimal spline serves as a favourable initial estimation for the additional grid points of the optimal spline. Some numerical examples are presented which are constructed by FORTRAN and MATLAB programs.     

������ڵľֲ��ͷ������ع鷽��

??Inspired by intuitive meanings of truncated power basis's coefficients, the local penalization based on range's linear decreasing function is given in penalized spline regression model. This method gives less penalization to fitting curve where data is with more volatility, which makes fitted curve controls tradeoff between goodness-of-fit and smoothness better. Simulations show that regression models with local penalized spline obtain lower information rules' scores than global penalized spline when the data is with heteroskedasticity.     

Geometrically designed,variable knot regression splines

A new method of Geometrically Designed least squares (LS) splines with variable knots, named GeDS, is proposed. It is based on the property that the spline regression function, viewed as a parametric curve, has a control polygon and, due to the shape preserving and convex hull properties, it closely follows the shape of this control polygon. The latter has vertices whose x-coordinates are certain knot averages and whose y-coordinates are the regression coefficients. Thus, manipulation of the position of the control polygon may be interpreted as estimation of the spline curve knots and coefficients. These geometric ideas are implemented in the two stages of the GeDS estimation method. In stage A, a linear LS spline fit to the data is constructed, and viewed as the initial position of the control polygon of a higher order (\(n>2\)) smooth spline curve. In stage B, the optimal set of knots of this higher order spline curve is found, so that its control polygon is as close to the initial polygon of stage A as possible and finally, the LS estimates of the regression coefficients of this curve are found. The GeDS method produces simultaneously linear, quadratic, cubic (and possibly higher order) spline fits with one and the same number of B-spline coefficients. Numerical examples are provided and further supplemental materials are available online.     

Nonnegative splines, in particular of degree five

Summary. We investigate splines from a variational point of view, which have the following properties: (a) they interpolate given data, (b) they stay nonnegative, when the data are positive, (c) for a given integer they minimize the functional for all nonnegative, interpolating . We extend known results for to larger , in particular to and we find general necessary conditions for solutions of this restricted minimization problem. These conditions imply that solutions are splines in an augmented grid. In addition, we find that the solutions are in and consist of piecewise polynomials in with respect to the augmented grid. We find that for general, odd there will be no boundary arcs which means (nontrivial) subintervals in which the spline is identically zero. We show also that the occurrence of a boundary arc in an interval between two neighboring knots prohibits the existence of any further knot in that interval. For we show that between given neighboring interpolation knots, the augmented grid has at most two additional grid points. In the case of two interpolation knots (the local problem) we develop polynomial equations for the additional grid points which can be used directly for numerical computation. For the general (global) problem we propose an algorithm which is based on a Newton iteration for the additional grid points and which uses the local spline data as an initial guess. There are extensions to other types of constraints such as two-sided restrictions, also ones which vary from interval to interval. As an illustration several numerical examples including graphs of splines manufactured by MATLAB- and FORTRAN-programs are given. Received November 16, 1995 / Revised version received February 24, 1997     

Parallel mesh methods for tension splines

This paper addresses the problem of shape preserving spline interpolation formulated as a differential multipoint boundary value problem (DMBVP for short). Its discretization by mesh method yields a five-diagonal linear system which can be ill-conditioned for unequally spaced data. Using the superposition principle we split this system in a set of tridiagonal linear systems with a diagonal dominance. The latter ones can be stably solved either by direct (Gaussian elimination) or iterative methods (SOR method and finite-difference schemes in fractional steps) and admit effective parallelization. Numerical examples illustrate the main features of this approach.     

Parallel mesh methods for tension splines

This paper addresses the problem of shape preserving spline interpolation formulated as a differential multipoint boundary value problem (DMBVP for short). Its discretization by mesh method yields a five-diagonal linear system which can be ill-conditioned for unequally spaced data. Using the superposition principle we split this system in a set of tridiagonal linear systems with a diagonal dominance. The latter ones can be stably solved either by direct (Gaussian elimination) or iterative methods (SOR method and finite-difference schemes in fractional steps) and admit effective parallelization. Numerical examples illustrate the main features of this approach.     

MARS: selecting basis functions and knots with an empirical Bayes method

An empirical Bayes method to select basis functions and knots in multivariate adaptive regression spline (MARS) is proposed, which takes both advantages of frequentist model selection approaches and Bayesian approaches. A penalized likelihood is maximized to estimate regression coefficients for selected basis functions, and an approximated marginal likelihood is maximized to select knots and variables involved in basis functions. Moreover, the Akaike Bayes information criterion (ABIC) is used to determine the number of basis functions. It is shown that the proposed method gives estimation of regression structure that is relatively parsimonious and more stable for some example data sets.     

Optimal design for smoothing splines

In the common nonparametric regression model we consider the problem of constructing optimal designs, if the unknown curve is estimated by a smoothing spline. A special basis for the space of natural splines is introduced and the local minimax property for these splines is used to derive two optimality criteria for the construction of optimal designs. The first criterion determines the design for a most precise estimation of the coefficients in the spline representation and corresponds to D-optimality, while the second criterion is the G-optimality criterion and corresponds to an accurate prediction of the curve. Several properties of the optimal designs are derived. In general, D- and G-optimal designs are not equivalent. Optimal designs are determined numerically and compared with the uniform design.     

广义线性模型中一类推广的拟似然估计与异方差性诊断检验

本文考虑了随机设计情形下一类普通的异方差回归模型，在这个模型中，假定回归函数与方差函数之间的关系服从推广的广义非线性模型，该模型在实际中很常见，广义线性模型便是其特例，首先，我们导出了均值函数的局部加权拟似然估计，然后，用它来得到方差函数的估计，并且证明了这些估计有较好的性质，最后，建立了异方差检验统计量，文中的方法很吸引人。     

Functional-coefficient partially linear regression model

In this paper, the functional-coefficient partially linear regression (FCPLR) model is proposed by combining nonparametric and functional-coefficient regression (FCR) model. It includes the FCR model and the nonparametric regression (NPR) model as its special cases. It is also a generalization of the partially linear regression (PLR) model obtained by replacing the parameters in the PLR model with some functions of the covariates. The local linear technique and the integrated method are employed to give initial estimators of all functions in the FCPLR model. These initial estimators are asymptotically normal. The initial estimator of the constant part function shares the same bias as the local linear estimator of this function in the univariate nonparametric model, but the variance of the former is bigger than that of the latter. Similarly, initial estimators of every coefficient function share the same bias as the local linear estimates in the univariate FCR model, but the variance of the former is bigger than that of the latter. To decrease the variance of the initial estimates, a one-step back-fitting technique is used to obtain the improved estimators of all functions. The improved estimator of the constant part function has the same asymptotic normality property as the local linear nonparametric regression for univariate data. The improved estimators of the coefficient functions have the same asymptotic normality properties as the local linear estimates in FCR model. The bandwidths and the smoothing variables are selected by a data-driven method. Both simulated and real data examples related to nonlinear time series modeling are used to illustrate the applications of the FCPLR model.     

Quadrature formulas descending from BS Hermite spline quasi-interpolation

Two new classes of quadrature formulas associated to the BS Boundary Value Methods are discussed. The first is of Lagrange type and is obtained by directly applying the BS methods to the integration problem formulated as a (special) Cauchy problem. The second descends from the related BS Hermite quasi-interpolation approach which produces a spline approximant from Hermite data assigned on meshes with general distributions. The second class formulas is also combined with suitable finite difference approximations of the necessary derivative values in order to define corresponding Lagrange type formulas with the same accuracy.     

Pseudoaffinity, de Boor algorithm, and blossoms

In order to ensure existence of a de Boor algorithm (hence of a B-spline basis) in a given spline space with (n+1)-dimensional sections, it is important to be able to generate each spline by restriction to the diagonal of a symmetric function of n variables supposed to be pseudoaffine w.r. to each variable. We proved that a way to obtain these three properties (symmetry, n-pseudoaffinity, diagonal property) is to suppose the existence of blossoms on the set of admissible n-tuples, given that blossoms are defined in a geometric way by means of intersections of osculating flats. In the present paper, we examine the converse: do symmetry, n-pseudoaffinity, and diagonal property imply existence of blossoms?     

加密网格点二元局部基插值样条函数

1.简介 由于在理论以及应用两方面的重要性,多元样条引起了许多人的注意([6],[7]),紧支撑光滑分片多项式函数对于曲面的逼近是一个十分有效的工具。由于它们的局部支撑性,它们很容易求值;由于它们的光滑性,它们能被应用到要满足一定光滑条件的情况下;由于它们是紧支撑的,它们的线性包有很大的逼近灵活性,而且用它们构造逼近方法来解决的系统是     

Bivariate Polynomial Natural Spline Interpolation Algorithms with Local Basis for Scattered Data

Because of its importance in both theory and applications, multivariate splines have attracted special attention in many fields. Based on the theory of spline functions in Hilbert spaces, bivariate polynomial natural splines for interpolating, smoothing or generalized interpolating of scattered data over an arbitrary domain are constructed with one-sided functions. However, this method is not well suited for large scale numerical applications. In this paper, a new locally supported basis for the bivariate polynomial natural spline space is constructed. Some properties of this basis are also discussed. Methods to order scattered data are shown and algorithms for bivariate polynomial natural spline interpolating are constructed. The interpolating coefficient matrix is sparse, and thus, the algorithms can be easily implemented in a computer.     

An EM Algorithm for Estimating Equations

This article presents an algorithm for accommodating missing data in situations where a natural set of estimating equations exists for the complete data setting. The complete data estimating equations can correspond to the score functions from a standard, partial, or quasi-likelihood, or they can be generalized estimating equations (GEEs). In analogy to the EM, which is a special case, the method is called the ES algorithm, because it iterates between an E-Step wherein functions of the complete data are replaced by their expected values, and an S-Step where these expected values are substituted into the complete-data estimating equation, which is then solved. Convergence properties of the algorithm are established by appealing to general theory for iterative solutions to nonlinear equations. In particular, the ES algorithm (and indeed the EM) are shown to correspond to examples of nonlinear Gauss-Seidel algorithms. An added advantage of the approach is that it yields a computationally simple method for estimating the variance of the resulting parameter estimates.     

一种Logistic模型参数估计的新方法及应用

介绍了logistic曲线参数估计的一种新方法,它是利用三次样条插值函数求导代替logistic曲线在这一点的导数值,进而利用最小二乘法得出参数的估计值,通过实例分析表明本文提出的方法比一般的三点法估计的参数值k再用线性化方法估计的参数值b,c,拟合精度更高.