Empirical likelihood for heteroscedastic partially linear models

We make empirical-likelihood-based inference for the parameters in heteroscedastic partially linear models. Unlike the existing empirical likelihood procedures for heteroscedastic partially linear models, the proposed empirical likelihood is constructed using components of a semiparametric efficient score. We show that it retains the double robustness feature of the semiparametric efficient estimator for the parameters and shares the desirable properties of the empirical likelihood for linear models. Compared with the normal approximation method and the existing empirical likelihood methods, the empirical likelihood method based on the semiparametric efficient score is more attractive not only theoretically but empirically. Simulation studies demonstrate that the proposed empirical likelihood provides smaller confidence regions than that based on semiparametric inefficient estimating equations subject to the same coverage probabilities. Hence, the proposed empirical likelihood is preferred to the normal approximation method as well as the empirical likelihood method based on semiparametric inefficient estimating equations, and it should be useful in practice.     

RECOVERY BASED FINITE ELEMENT METHOD FOR BIHARMONIC EQUATION IN 2D

We design and numerically validate a recovery based linear finite element method for solving the biharmonic equation.The main idea is to replace the gradient operator▽on linear finite element space by G(▽)in the weak formulation of the biharmonic equation,where G is the recovery operator which recovers the piecewise constant function into the linear finite element space.By operator G,Laplace operator△is replaced by▽·G(▽).Furthermore,the boundary condition on normal derivative▽u-n is treated by the boundary penalty method.The explicit matrix expression of the proposed method is also introduced.Numerical examples on the uniform and adaptive meshes are presented to illustrate the correctness and effectiveness of the proposed method.     

Homotopy analysis method with a non-homogeneous term in the auxiliary linear operator

We demonstrate the efficiency of a modification of the normal homotopy analysis method (HAM) proposed by Liao [2] by including a non-homogeneous term in the auxiliary linear operator (this can be considered as a special case of “further generalization” of HAM given by Liao in [2]). We then apply the modified method to a few examples. It is observed that including a non-homogeneous term gives faster convergence in comparison to normal HAM. We also prove a convergence theorem, which shows that our technique yields the convergent solution.     

Spline-based sieve maximum likelihood estimation in the partly linear model under monotonicity constraints

We study a spline-based likelihood method for the partly linear model with monotonicity constraints. We use monotone B-splines to approximate the monotone nonparametric function and apply the generalized Rosen algorithm to compute the estimators jointly. We show that the spline estimator of the nonparametric component achieves the possible optimal rate of convergence under the smooth assumption and that the estimator of the regression parameter is asymptotically normal and efficient. Moreover, a spline-based semiparametric likelihood ratio test is established to make inference of the regression parameter. Also an observed profile information method to consistently estimate the standard error of the spline estimator of the regression parameter is proposed. A simulation study is conducted to evaluate the finite sample performance of the proposed method. The method is illustrated by an air pollution study.     

An equivalent linearization technique for nonlinear piezoelectric energy harvesters under Gaussian white noise

An equivalent linearization technique is proposed to determine approximately the output voltage a nonlinear piezoelectric energy harvester excited by Gaussian white noise excitations. Equivalent linear system is derived from minimizing the mean-squared of the error. The linear equivalent coefficients are presented by the method of normal truncation. The exact solution of equivalent linear system is derived obtained. The effectiveness of the method is demonstrated by numerical simulations.     

On possibility of obtaining linear accuracy evaluation of approximate solutions to inverse problems

We consider a concept of linear a priori estimate of the accuracy for approximate solutions to inverse problems with perturbed data. We establish that if the linear estimate is valid for a method of solving the inverse problem, then the inverse problem is well-posed according to Tikhonov. We also find conditions, which ensure the converse for the method of solving the inverse problem independent on the error levels of data. This method is well-known method of quasi-solutions by V. K. Ivanov. It provides for well-posed (according to Tikhonov) inverse problems the existence of linear estimates. If the error levels of data are known, a method of solving well-posed according to Tikhonov inverse problems is proposed. This method called the residual method on the correctness set (RMCS) ensures linear estimates for approximate solutions. We give an algorithm for finding linear estimates in the RMCS.     

An efficient sampling scheme for dynamic generalized models

A multimove sampling scheme for the state parameters of non-Gaussian and nonlinear dynamic models for univariate time series is proposed. This procedure follows the Bayesian framework, within a Gibbs sampling algorithm with steps of the Metropolis–Hastings algorithm. This sampling scheme combines the conjugate updating approach for generalized dynamic linear models, with the backward sampling of the state parameters used in normal dynamic linear models. A quite extensive Monte Carlo study is conducted in order to compare the results obtained using our proposed method, conjugate updating backward sampling (CUBS), with those obtained using some algorithms previously proposed in the Bayesian literature. We compare the performance of CUBS with other sampling schemes using two real datasets. Then we apply our algorithm in a stochastic volatility model. CUBS significantly reduces the computing time needed to attain convergence of the chains, and is relatively simple to implement.     

Local linear regression for data with AR errors

In many statistical applications, data are collected over time, and they are likely correlated. In this paper, we investigate how to incorporate the correlation information into the local linear regression. Under the assumption that the error process is an auto-regressive process, a new estimation procedure is proposed for the nonparametric regression by using local linear regression method and the profile least squares techniques. We further propose the SCAD penalized profile least squares method to determine the order of auto-regressive process. Extensive Monte Carlo simulation studies are conducted to examine the finite sample performance of the proposed procedure, and to compare the performance of the proposed procedures with the existing one. From our empirical studies, the newly proposed procedures can dramatically improve the accuracy of naive local linear regression with working-independent error structure. We illustrate the proposed methodology by an analysis of real data set.     

Instrumental variable approach to covariate measurement error in generalized linear models

The paper presents the method of moments estimation for generalized linear measurement error models using the instrumental variable approach. The measurement error has a parametric distribution that is not necessarily normal, while the distributions of the unobserved covariates are nonparametric. We also propose simulation-based estimators for the situation where the closed forms of the moments are not available. The proposed estimators are strongly consistent and asymptotically normally distributed under some regularity conditions. Finite sample performances of the estimators are investigated through simulation studies.     

Solving boundary value problems,integral, and integro-differential equations using Gegenbauer integration matrices

We introduce a hybrid Gegenbauer (ultraspherical) integration method (HGIM) for solving boundary value problems (BVPs), integral and integro-differential equations. The proposed approach recasts the original problems into their integral formulations, which are then discretized into linear systems of algebraic equations using Gegenbauer integration matrices (GIMs). The resulting linear systems are well-conditioned and can be easily solved using standard linear system solvers. A study on the error bounds of the proposed method is presented, and the spectral convergence is proven for two-point BVPs (TPBVPs). Comparisons with other competitive methods in the recent literature are included. The proposed method results in an efficient algorithm, and spectral accuracy is verified using eight test examples addressing the aforementioned classes of problems. The proposed method can be applied on a broad range of mathematical problems while producing highly accurate results. The developed numerical scheme provides a viable alternative to other solution methods when high-order approximations are required using only a relatively small number of solution nodes.     

A new construction for skew multivariate distributions

This paper considers a new approach to develop a very general class of skew multivariate distributions. The approach is based on a linear combination of an elliptically distributed random variable with a linear constraint. Using this approach two different classes of multivariate distributions are constructed based on original distribution. These new classes include different types of skew normal (type A and type B) and other skew elliptical distributions, exist in the literature. We also derive the moment generating function, marginal and conditional density of our proposed classes of distributions. Straightforward explanations are applied to demonstrate the relationships among previous approaches by others with our proposed class of skew distributions.     

Numerical discretization-based kernel type estimation methods for ordinary differential equation models

We consider the problem of parameter estimation in both linear and nonlinear ordinary differential equation(ODE) models. Nonlinear ODE models are widely used in applications. But their analytic solutions are usually not available. Thus regular methods usually depend on repetitive use of numerical solutions which bring huge computational cost. We proposed a new two-stage approach which includes a smoothing method(kernel smoothing or local polynomial fitting) in the first stage, and a numerical discretization method(Eulers discretization method, the trapezoidal discretization method,or the Runge–Kutta discretization method) in the second stage. Through numerical simulations, we find the proposed method gains a proper balance between estimation accuracy and computational cost.Asymptotic properties are also presented, which show the consistency and asymptotic normality of estimators under some mild conditions. The proposed method is compared to existing methods in term of accuracy and computational cost. The simulation results show that the estimators with local linear smoothing in the first stage and trapezoidal discretization in the second stage have the lowest average relative errors. We apply the proposed method to HIV dynamics data to illustrate the practicability of the estimator.     

时空模型的局部众数回归

时空数据经常含有奇异点或来自重尾分布,此时基于最小二乘的估计方法效果欠佳,需要更稳健的估计方法.本文提出时空模型的基于局部众数(local modal, LM)的局部线性估计方法.理论和数据分析结果都显示,若数据含有奇异点或来自重尾分布,基于局部众数的局部线性方法比基于最小二乘的局部线性方法有效;若数据无奇异点且来自正态分布,两种方法效率渐近一致.本文采用众数期望最大化(modal expectation-maximization, MEM)算法,并在数据相依情形下得出估计量的渐近正态性.     

Bernstein Polynomial Estimation in the Partially Linear Model under Monotonicity Constraints

In this paper, a Bernstein-polynomial-based likelihood method is proposed for the partially linear model under monotonicity constraints. Monotone Bernstein polynomials are employed to approximate the monotone nonparametric function in the model. The estimator of the regression parameter is shown to be asymptotically normal and efficient, and the rate of convergence of the estimator of the nonparametric component is established, which could be the optimal under the smooth assumptions. A simulation study and a real data analysis are conducted to evaluate the finite sample performance of the proposed method.     

Partially varying coefficient single-index additive hazard models

The partially linear additive hazards model has been proposed to study the interaction between some covariates and an exposure variable. In this paper, we extend it to the partially varying coefficient single-index additive hazard model where the high dimension covariates are collapsed to a single index, due to practical needs. Two sets of estimating equations were proposed to estimate the varying coefficient functions in the linear components: the link function for the single index and the single-index parameter vector separately. It was shown that the proposed local and global estimators are asymptotically normal. Simulation studies were conducted to examine the finite-sample performance of our method to compare the relative performance of our method with existing ones. A real data analysis was used to illustrate the proposed methods.     

响应变量存在缺失时部分线性模型的经验似然推断

考虑响应变量带有缺失的部分线性模型,采用借补的思想,研究了参数部分和非参数部分的经验似然推断,证明了所提出的经验对数似然比统计量依分布收敛到χ2分布,由此构造参数部分和函数部分的置信域和逐点置信区间.对参数部分,模拟比较了经验似然与正态逼近方法;对函数部分,模拟了函数的逐点置信区间.     

响应变量随机缺失下的半参变系数模型的均值估计

在响应变量随机缺失时,研究了半参数变系数模型响应变量均值的借补估计.首先利用完整个体估计模型中的参数与非参数部分,然后再用借补方法与加权借补方法估计响应变量的均值.最后求出了估计的渐近偏差与渐近方差,研究了所得到的估计的渐近性质,并进行模拟比较.     

Testing in linear composite quantile regression models

Composite quantile regression (CQR) can be more efficient and sometimes arbitrarily more efficient than least squares for non-normal random errors, and almost as efficient for normal random errors. Based on CQR, we propose a test method to deal with the testing problem of the parameter in the linear regression models. The critical values of the test statistic can be obtained by the random weighting method without estimating the nuisance parameters. A distinguished feature of the proposed method is that the approximation is valid even the null hypothesis is not true and power evaluation is possible under the local alternatives. Extensive simulations are reported, showing that the proposed method works well in practical settings. The proposed methods are also applied to a data set from a walking behavior survey.     

Modified integral current methods in electrodynamics of nonhomogeneous media

The main difficulty in numerical solution of integral equations of electrodynamics is associated with the need to solve a high-order system of linear equations with a dense matrix. It is therefore relevant to develop numerical methods that lead to linear equation systems of lower order at the cost of more complex evaluation of the coefficients. In this article we propose a method for solving linear equations of electrodynamics which is a modification of the integral current method. The main distinctive feature of the proposed method is double integration of the electric Green’s tensor in the process of algebraization of the original integral equation. The solutions of the system of linear equations are thus integral means of the electric field inside the anomaly constructed by the proposed transformation formula. We prove convergence and derive error bounds for both the solution of the integral equation and the electromagnetic field components evaluated from approximate transformation formulas.     

基于生存数据的半参数变换的删失回归估计

生存数据经过未知的单调变换后等于协变量的线性函数加上随机误差, 随机误差的分布函数已知或是带未知参数的已知函数\bd 本文先给出未知单调变换的一个相合估计, 再对删失数据做变换, 在此基础上给出了协变量系数的最小二乘估计, 并讨论它的大样本性质.