Exact power comparison of three criteria to assess the independence between two sets of variates

Powers of the three criteria are evaluated for testing the hypothesis of the independence between a -set and a q-set of variates in a (p + q) -variate normal population. They are: (1) the likelihood ratio type criterion, Wt (2) the largest root criterion, r1, and (3) criterion of the sum of roots, V. For p= 2, Pillai and Jayachandran, and others have studied for the restricted range of the alternative hypothesis. Recently the power of the largest root was investigated in detail by Sugiyama and %%. In this paper, their power functions are compared in a wide range of the alternative hypotheses. The powers of rl and V are locally optimum, but the W shows a large power in a wide range.     

Model Averaging Using Fractional Polynomials to Estimate a Safe Level of Exposure

Quantitative risk assessment involves the determination of a safe level of exposure. Recent techniques use the estimated dose-response curve to estimate such a safe dose level. Although such methods have attractive features, a low-dose extrapolation is highly dependent on the model choice. Fractional polynomials, basically being a set of (generalized) linear models, are a nice extension of classical polynomials, providing the necessary flexibility to estimate the dose-response curve. Typically, one selects the best-fitting model in this set of polynomials and proceeds as if no model selection were carried out. We show that model averaging using a set of fractional polynomials reduces bias and has better precision in estimating a safe level of exposure (say, the benchmark dose), as compared to an estimator from the selected best model. To estimate a lower limit of this benchmark dose, an approximation of the variance of the model-averaged estimator, as proposed by Burnham and Anderson, can be used. However, this is a conservative method, often resulting in unrealistically low safe doses. Therefore, a bootstrap-based method to more accurately estimate the variance of the model averaged parameter is proposed.     

A corrected likelihood method for the proportional hazards model with covariates subject to measurement error

There has been extensive interest in discussing inference methods for survival data when some covariates are subject to measurement error. It is known that standard inferential procedures produce biased estimation if measurement error is not taken into account. With the Cox proportional hazards model a number of methods have been proposed to correct bias induced by measurement error, where the attention centers on utilizing the partial likelihood function. It is also of interest to understand the impact on estimation of the baseline hazard function in settings with mismeasured covariates. In this paper we employ a weakly parametric form for the baseline hazard function and propose simple unbiased estimating functions for estimation of parameters. The proposed method is easy to implement and it reveals the connection between the naive method ignoring measurement error and the corrected method with measurement error accounted for. Simulation studies are carried out to evaluate the performance of the estimators as well as the impact of ignoring measurement error in covariates. As an illustration we apply the proposed methods to analyze a data set arising from the Busselton Health Study [Knuiman, M.W., Cullent, K.J., Bulsara, M.K., Welborn, T.A., Hobbs, M.S.T., 1994. Mortality trends, 1965 to 1989, in Busselton, the site of repeated health surveys and interventions. Austral. J. Public Health 18, 129–135].     

An Equivalence of Conditional and Unconditional Maximum Likelihood Estimators via Infinite Replication of Observations

Motivated by an application with complex survey data, we show that for logistic regression with a simple matched-pairs design, infinitely replicating observations and maximizing the conditional likelihood results in an estimator exactly identical to the unconditional maximum likelihood estimator based on the original sample, which is inconsistent. Therefore, applying conditional likelihood methods to a pseudosample with observations replicated a large number of times can lead to an inconsistent estimator; this casts doubt on one possible approach to conditional logistic regression with complex survey data. We speculate that for more general designs, an asymptotic equivalence holds.     

谱域法分析微带线电流分布特性

用改进的谱域法分析微带线的电流分布特性,并用切比雪夫多项式做基函数来表示切向电流分布,得到了较宽频率范围内的切向电流分布特性,再根据这些特性在不同频段上取适当数量的基函数对有效相对介电常数进行了计算,与其他数值计算方法相比,该方法不仅可以得到高精度的计算结果,还可以在相应频段减少计算时间,计算结果与HFSS仿真结果在h/λ0≤0.7范围内很接近。     

Extended stirling family of discrete probability distributions

In this paper. introduced is a family of discrete probability distributions. whose probability function includes explieitly the Striling-Carlitz polynomial of the first or the second kind. The new family extends the stirling family of distributions. Sibuya (1988) includes the conditional distributions of the orginal ones and enlarges the application area.     

Multilinear generating functions for Charlier polynomials

Charlier configurations provide a combinatorial model for Charlier polynomials. We use this model to give a combinatorial proof of a multilinear generating function for Charlier polynomials. As special cases of the multilinear generating function, we obtain the bilinear generating function for Charlier polynomials and formulas for derangements.     

Partitioning Pearson's Chi-squared Statistic for Singly Ordered Two-way Contingency Tables

This paper presents a partition of Pearson's chi-squared statistic for singly ordered two-way contingency tables. The partition involves using orthogonal polynomials for the ordinal variable while generalized basic vectors are used for the non-ordinal variable. The benefit of this partition is that important information about the structure of the ordered variable can be identified in terms of locations, dispersion and higher order components. For the non-ordinal variable, it is shown that the squared singular values from the singular value decomposition of the transformed dataset can be partitioned into location, dispersion and higher order components. The paper also uses the chi-squared partition to present an alternative to the maximum likelihood technique of parameter estimation for the log-linear analysis of the contingency table.     

Automatic Bayesian quantile regression curve fitting

Quantile regression, including median regression, as a more completed statistical model than mean regression, is now well known with its wide spread applications. Bayesian inference on quantile regression or Bayesian quantile regression has attracted much interest recently. Most of the existing researches in Bayesian quantile regression focus on parametric quantile regression, though there are discussions on different ways of modeling the model error by a parametric distribution named asymmetric Laplace distribution or by a nonparametric alternative named scale mixture asymmetric Laplace distribution. This paper discusses Bayesian inference for nonparametric quantile regression. This general approach fits quantile regression curves using piecewise polynomial functions with an unknown number of knots at unknown locations, all treated as parameters to be inferred through reversible jump Markov chain Monte Carlo (RJMCMC) of Green (Biometrika 82:711–732, 1995). Instead of drawing samples from the posterior, we use regression quantiles to create Markov chains for the estimation of the quantile curves. We also use approximate Bayesian factor in the inference. This method extends the work in automatic Bayesian mean curve fitting to quantile regression. Numerical results show that this Bayesian quantile smoothing technique is competitive with quantile regression/smoothing splines of He and Ng (Comput. Stat. 14:315–337, 1999) and P-splines (penalized splines) of Eilers and de Menezes (Bioinformatics 21(7):1146–1153, 2005).