Numerical algorithms for solving nonlinear L р-norm estimation problems: part II - a mixture method for large residual and illo-conditioned problems

The nonlinear least squares algorithm of Gill and Murray (1978) is extended and modified to solve nonlinear L _р-norm estimation problems efficiently. The new algorithm uses a mixture of 1st-order derivative (Guass-Newton) and 2nd-order derivative (Newton) search directions. A new rule for selecting the “grade” r of the p-jacobiab matrix J_p was also incorporated. This brought about rapid convergence of the algorithm on previously reported test examples.     

Inference procedures for the variance gamma model and applications

Goodness-of-fit tests for the family of the four-parameter normal–variance gamma distribution are constructed. The tests are based on a weighted integral incorporating the empirical characteristic function of suitably standardized data. Non-standard algorithms are employed for the computation of the maximum-likelihood estimators of the parameters involved in the test statistic, while Monte Carlo results are used in order to compare the new test with some classical goodness-of-fit methods. A real-data application is also included.     

A control in stratified sampling

A new expression of the variance of the units from a stratified finite population wit L strata, is obtained in funcions of the means L–1 strata, and in conseqeuence we obtained formulae which relates the mean of a stratum in a statistical study with the population mean and the ones of the remaining strate. As an application, we obtained a useful checking of the estimation consistency by stratified sampling in any precies survey     

Markov Random Fields with Higher-order Interactions

Discrete-state Markov random fields on regular arrays have played a significant role in spatial statistics and image analysis. For example, they are used to represent objects against background in computer vision and pixel-based classification of a region into different crop types in remote sensing. Convenience has generally favoured formulations that involve only pairwise interactions. Such models are in themselves unrealistic and, although they often perform surprisingly well in tasks such as the restoration of degraded images, they are unsatisfactory for many other purposes. In this paper, we consider particular forms of Markov random fields that involve higher-order interactions and therefore are better able to represent the large-scale properties of typical spatial scenes. Interpretations of the parameters are given and realizations from a variety of models are produced via Markov chain Monte Carlo. Potential applications are illustrated in two examples. The first concerns Bayesian image analysis and confirms that pairwise-interaction priors may perform very poorly for image functionals such as number of objects, even when restoration apparently works well. The second example describes a model for a geological dataset and obtains maximum-likelihood parameter estimates using Markov chain Monte Carlo. Despite the complexity of the formulation, realizations of the estimated model suggest that the representation is quite realistic.     

相干分布式目标一维波达方向搜索迭代估计方法

在信号源角分布函数具体的数学形式未知的情况下,基于角分布函数共轭对称性约束条件,提出了一种相干分布式目标一维波达方向搜索迭代估计方法,该方法不用多维参数搜索,并适用于角信号分布函数形式不同的分布式目标同时存在的情况。     

A semiparametric estimator of the distribution function of a variable measured with error

The estimation of the distribution functon of a random variable X measured with error is studied. Let the i-th observation on X be denoted by Y_iX_i+ε_i where ε_i is the measuremen error. Let {Y_i} (i=1,2,…,n) be a sample of independent observations. It is assumed that {X_i} and {∈_i} are mutually independent and each is identically distributed. As is standard in the literature for this problem, the distribution of e is assumed known in the development of the methodology. In practice, the measurement error distribution is estimated from replicate observations.

The proposed semiparametric estimator is derived by estimating the quantises of X on a set of n transformed V-values and smoothing the estimated quantiles using a spline function. The number of parameters of the spline function is determined by the data with a simple criterion, such as AIC. In a simulation study, the semiparametric estimator dominates an optimal kernel estimator and a normal mixture estimator for a wide class of densities.

The proposed estimator is applied to estimate the distribution function of the mean pH value in a field plot. The density function of the measurement error is estimated from repeated measurements of the pH values in a plot, and is treated as known for the estimation of the distribution function of the mean pH value.     

An algorithm for computing nonparametric estimators of change‐points

An efficient FORTRAN subroutine for computing three nonparametric point estimators of change-points is offered. Estimates computed from the subroutine are obtained for two classical data sets     

Pareto Poisson–Lindley distribution with applications

A new lifetime distribution is introduced based on compounding Pareto and Poisson–Lindley distributions. Several statistical properties of the distribution are established, including behavior of the probability density function and the failure rate function, heavy- and long-right tailedness, moments, the Laplace transform, quantiles, order statistics, moments of residual lifetime, conditional moments, conditional moment generating function, stress–strength parameter, Rényi entropy and Song's measure. We get maximum-likelihood estimators of the distribution parameters and investigate the asymptotic distribution of the estimators via Fisher's information matrix. Applications of the distribution using three real data sets are presented and it is shown that the distribution fits better than other related distributions in practical uses.     

End-Point Estimation for Decreasing Densities: Asymptotic Behaviour of the Penalized Likelihood Ratio

Abstract.  We consider the problem of estimating the modal value of a decreasing density on the positive real line. This has application in several interesting phenomena arising, for example, in renewal theory, and in biased and distance samplings. We use a penalized likelihood ratio-based approach for inference and derive the scale-free universal large sample null distribution of the log-likelihood ratio, using a suitably chosen penalty parameter. We present simulation results and a real data analysis to corroborate our findings, and compare the performance of the confidence sets with the existing results.     

Testing in semiparametric models with interaction, with applications to gene–environment interactions

Summary.  Motivated from the problem of testing for genetic effects on complex traits in the presence of gene–environment interaction, we develop score tests in general semiparametric regression problems that involves Tukey style 1 degree-of-freedom form of interaction between parametrically and non-parametrically modelled covariates. We find that the score test in this type of model, as recently developed by Chatterjee and co-workers in the fully parametric setting, is biased and requires undersmoothing to be valid in the presence of non-parametric components. Moreover, in the presence of repeated outcomes, the asymptotic distribution of the score test depends on the estimation of functions which are defined as solutions of integral equations, making implementation difficult and computationally taxing. We develop profiled score statistics which are unbiased and asymptotically efficient and can be performed by using standard bandwidth selection methods. In addition, to overcome the difficulty of solving functional equations, we give easy interpretations of the target functions, which in turn allow us to develop estimation procedures that can be easily implemented by using standard computational methods. We present simulation studies to evaluate type I error and power of the method proposed compared with a naive test that does not consider interaction. Finally, we illustrate our methodology by analysing data from a case–control study of colorectal adenoma that was designed to investigate the association between colorectal adenoma and the candidate gene NAT2 in relation to smoking history.