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

??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.     

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.     

基于径向基的自适应惩罚样条回归模型

传统惩罚样条回归模型中惩罚项的设置未考虑数据的空间异质性,因而对复杂数据的拟合缺乏自适应性.文章通过对径向基函数的几何意义分析,以节点两侧相邻区域内数据点的纵向极差为基础,构造局部惩罚权重向量并加入到约束回归模型的惩罚项中,构造了基于径向基的自适应惩罚样条回归模型.新模型在观测数据波动较大的区域,给予拟合曲线较小的惩罚,而在观测数据波动较小的区域,给予拟合曲线较大的惩罚,从而使拟合曲线能自适应地反映观测数据的局部变化特征.模拟和应用结果显示新模型的拟合效果显著优于传统的惩罚样条回归模型.     

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

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

Computing bivariate splines in scattered data fitting and the finite-element method

A number of useful bivariate spline methods are global in nature, i.e., all of the coefficients of an approximating spline must be computed at one time. Typically this involves solving a system of linear equations. Examples include several well-known methods for fitting scattered data, such as the minimal energy, least-squares, and penalized least-squares methods. Finite-element methods for solving boundary-value problems are also of this type. It is shown here that these types of globally-defined splines can be efficiently computed, provided we work with spline spaces with stable local minimal determining sets.     

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.     

Trivariate near-best blending spline quasi-interpolation operators

A method to define trivariate spline quasi-interpolation operators (QIOs) is developed by blending univariate and bivariate operators whose linear functionals allow oversampling. In this paper, we construct new operators based on univariate B-splines and bivariate box splines, exact on appropriate spaces of polynomials and having small infinity norms. An upper bound of the infinity norm for a general blending trivariate spline QIO is derived from the Bernstein-Bézier coefficients of the fundamental functions associated with the operators involved in the construction. The minimization of the resulting upper bound is then proposed and the existence of a solution is proved. The quadratic and quartic cases are completely worked out and their exact solutions are explicitly calculated.     

Robust penalized regression spline fitting with application to additive mixed modeling

An increasingly popular method for smoothing noisy data is penalized regression spline fitting. In this paper a new procedure is proposed for fitting robust penalized regression splines. This procedure is computationally fast, straightforward to implement, and can be paired with any smoothing parameter selection method. In addition, it can also be extended to other settings, such as additive mixed modeling. Both simulated and real data examples are used to illustrate the effectiveness of the procedure.     

Principal Component Analysis of Two-Dimensional Functional Data

This article presents and compares two approaches of principal component (PC) analysis for two-dimensional functional data on a possibly irregular domain. The first approach applies the singular value decomposition of the data matrix obtained from a fine discretization of the two-dimensional functions. When the functions are only observed at discrete points that are possibly sparse and may differ from function to function, this approach incorporates an initial smoothing step prior to the singular value decomposition. The second approach employs a mixed effects model that specifies the PC functions as bivariate splines on triangulations and the PC scores as random effects. We apply the thin-plate penalty for regularizing the function estimation and develop an effective expectation–maximization algorithm for calculating the penalized likelihood estimates of the parameters. The mixed effects model-based approach integrates scatterplot smoothing and functional PC analysis in a unified framework and is shown in a simulation study to be more efficient than the two-step approach that separately performs smoothing and PC analysis. The proposed methods are applied to analyze the temperature variation in Texas using 100 years of temperature data recorded by Texas weather stations. Supplementary materials for this article are available online.     

固定效应部分线性可加面板数据模型的惩罚二次推断估计

具有较强解释力和灵活性的部分线性可加面板数据模型在各学科领域应用广泛.针对个体内存在相关结构的固定效应部分线性可加面板数据模型,本文在结合幂样条函数和最小二乘虚拟变量(LSDV)法的基础上,利用惩罚二次推断函数(PQIF)法对模型进行估计,在一定的正则条件下,证明了参数估计的渐近正态性和非参数估计的收敛性,Monte Carlo数值模拟显示所述估计方法具有良好的有限样本表现,同时,我们还将估计技术应用于实际数据分析中.     

Bivariate Polynomial Natural Spline Interpolation to Scattered Data

By means of the theory of spline interpolation in Hilbert spaces, the bivariate polynomial natural spline interpolation to scattered data is constructed. The method can easily be carried out on a computer, and parallelly generalized to high dimensional cases as well. The results can be used for numerical integration in higher dimensions and numerical solution of partial differential equations, and so on.     

Construction of Bivariate Compactly Supported Biorthogonal Box Spline Wavelets with Arbitrarily High Regularities

We give a simple formula for the duals of the filters associated with bivariate box spline functions. We show how to construct bivariate non-separable compactly supported biorthogonal wavelets associated with box spline functions which have arbitrarily high regularities.     

散乱数据带自然边界条件二元样条光顺及数值微分

考虑一种新的散乱数据带自然边界二元样条光顺问题.根据样条变分理论和Hilbert空间样条函数方法,构造出了显式的二元带自然边界光顺样条解,其表达式简单且系数可以由系数矩阵对称正定的线性方程组确定.证明了解的存在和唯一性,讨论了收敛性和误差估计.并由此得到一种新的基于散乱数据上的正则化二元数值微分的方法.最后,给出了一些数值例子对方法进行了验证.     

Mixed Model-Based Hazard Estimation

This article proposes a new method for estimation of the hazard function from a set of censored failure time data, with a view to extending the general approach to more complicated models. The approach is based on a mixed model representation of penalized spline hazard estimators. One payoff is the automation of the smoothing parameter choice through restricted maximum likelihood. Another is the option to use standard mixed model software for automatic hazard estimation.     

Estimation of Smooth Time-Varying Parameters in State Space Models

In this article, we propose a penalized likelihood method to estimate time-varying parameters in standard linear state space models. The time-varying parameter is modeled as a smoothing spline and then expressed as a state space model. The maximum likelihood method is used to estimate the smoothing parameter. The proposed method is assessed by a simulation study and applied to virological response data from an HIV-infected patient receiving antiretroviral treatment.     

Hazard Estimation With Bivariate Survival Data and Copula Density Estimation

Bivariate survival function can be expressed as the composition of marginal survival functions and a bivariate copula and, consequently, one may estimate bivariate hazard functions via marginal hazard estimation and copula density estimation. Leveraging on earlier developments on penalized likelihood density and hazard estimation, a nonparametric approach to bivariate hazard estimation is being explored in this article. The new ingredient here is the nonparametric estimation of copula density, a subject of interest by itself, and to accommodate survival data one needs to allow for censoring and truncation in the setting. A simple copularization process is implemented to convert density estimates into copula densities, and a cross-validation scheme is devised for density estimation under censoring and truncation. Empirical performances of the techniques are investigated through simulation studies, and potential applications are illustrated using real-data examples and open-source software.     

Bayesian binary regression involving two explanatory variables

Summary The purpose of the present paper is to propose a practical Bayesian procedure for the estimation of binary response probability where the explanatory variable is bivariate. The procedure is an extension of the procedure for univariate case which was proposed by the present authors [2] and is based on a model which approximates the logistic transformation of response probability by a quadratic orthogonal spline function on the two-dimensional space of explanatory variable. The flexibility of the model is guaranteed by assuming a spline function on sufficiently fine mesh. To obtain stable estimates we introduce a prior distribution of the parameters of the model. The prior distribution has several parameters (hyper-parameters) which are chosen to minimize an Bayesian information criterion ABIC. The procedure is applicable cable to cases where each explanatory variable takes continuous values provided that the probability of the occurrence changes smoothly. The practical utility of the procedure is demonstrated by examples of applications to five sets of data. The Institute of Statistical Mathematics     

Spline quasi-interpolation in the Bernstein basis and its application to digital elevation models

A nonstandard low-cost spline approximation method for approximating bivariate functions is constructed. It is applied for Digital Elevation approximation and then its accuracy in the downscaling process is studied.     

Extrema detection of bivariate spline functions

In this paper, we develop a systematic method for detecting the extrema of bivariate spline functions and of their derivatives. It is assumed that the splines are constituted by employing normalized, uniform B-splines as the basis functions, and the detection procedure fully utilizes the spline properties. All the extrema can be found except those with singular Hessian matrix. By numerical examples, we demonstrate the effectiveness of the method.     

Local Lagrange Interpolation with Bivariate Splines of Arbitrary Smoothness

We describe a method which can be used to interpolate function values at a set of scattered points in a planar domain using bivariate polynomial splines of any prescribed smoothness. The method starts with an arbitrary given triangulation of the data points, and involves refining some of the triangles with Clough-Tocher splits. The construction of the interpolating splines requires some additional function values at selected points in the domain, but no derivatives are needed at any point. Given n data points and a corresponding initial triangulation, the interpolating spline can be computed in just O(n) operations. The interpolation method is local and stable, and provides optimal order approximation of smooth functions.