Bayesian shrinkage estimates for regression coefficients in m populations

For a linear regression model over m populations with separate regression coefficients but a common error variance, a Bayesian model is employed to obtain regression coefficient estimates which are shrunk toward an overall value. The formulation uses Normal priors on the coefficients and diffuse priors on the grand mean vectors, the error variance, and the between-to-error variance ratios. The posterior density of the parameters which were given diffuse priors is obtained. From this the posterior means and variances of regression coefficients and the predictive mean and variance of a future observation are obtained directly by numerical integration in the balanced case, and with the aid of series expansions in the approximately balanced case. An example is presented and worked out for the case of one predictor variable. The method is an extension of Box & Tiao's Bayesian estimation of means in the balanced one-way random effects model.     

Estimation of gene/haplotype frequencies in genetic marker systems based on phenotype data

Two equivalent methods (gene counting and maximum likelihood) for estimating gene frequencies in a general genetic marker system based on observed phenotype data are derived. Under the maximum likelihood approach, an expression is given for the estimated covariance matrix from which estimated standard errors of the estimators can be found. In addition, consideration is given to the problem of estimating gene frequencies when there are available several independent population data sets.     

Analysing exponential family based messy data

Recently exponential family based random effects models have received considerable attention. These models usually arise from an unobservable random process added to the independent exponential family models. An unobservable correlated process, however, would cause correlations among the exponential family based data. This paper, first, develops an asymptotically optimal test for testing the appropriateness of a fixed effects model for the exponential family based independent data versus a random effects model for the exponential family based independent or correlated data. The paper, then, provides a general framework on regression analysis for the exponential family based data generated under the random effects models.     

On the efficiency of using the sample kurtosis in selecting optimal lpestimators

This paper examines the efficiency of thesample kurtosisin obtaining L_P estimates as an estimates of central tendency for symmetric distributions. Moreover, guidelines are established for determining an optimal value of P based on the kurtosis of the error distribution.     

ASYMPTOTIC AND SMALL SAMPLE STATISTICAL PROPERTIES OF RANDOM FRAILTY VARIANCE ESTIMATES FOR SHARED GAMMA FRAILTY MODELS

This paper concerns maximum likelihood estimation for the semiparametric shared gamma frailty model; that is the Cox proportional hazards model with the hazard function multiplied by a gamma random variable with mean 1 and variance θ. A hybrid ML-EM algorithm is applied to 26 400 simulated samples of 400 to 8000 observations with Weibull hazards. The hybrid algorithm is much faster than the standard EM algorithm, faster than standard direct maximum likelihood (ML, Newton Raphson) for large samples, and gives almost identical results to the penalised likelihood method in S-PLUS 2000. When the true value θ₀ of θ is zero, the estimates of θ are asymptotically distributed as a 50–50 mixture between a point mass at zero and a normal random variable on the positive axis. When θ₀ > 0, the asymptotic distribution is normal. However, for small samples, simulations suggest that the estimates of θ are approximately distributed as an x ? (100 ? x)% mixture, 0 ≤ x ≤ 50, between a point mass at zero and a normal random variable on the positive axis even for θ₀ > 0. In light of this, p-values and confidence intervals need to be adjusted accordingly. We indicate an approximate method for carrying out the adjustment.     

Bayesian Estimation of Change Point in Inverse Weibull Sequence

A sequence of independent lifetimes X ₁,…, X _m , X _m+1,…, X _n were observed from inverse Weibull distribution with mean stress θ₁ and reliability R ₁(t ₀) at time t ₀ but later it was found that there was a change in the system at some point of time m and it is reflected in the sequence after X _m by change in mean stress θ₁ and in reliability R ₂(t ₀) at time t ₀. The Bayes estimators of m, R ₁(t ₀) and R ₂(t ₀) are derived when a poor and a more detailed prior information is introduced into the inferential procedure. The effects of correct and wrong prior information on the Bayes estimators are studied.     

On sufficient conditions for the strong consistency of least-squares estimates

The strong consistency of the least-squares estimates in regression models is obtained when the errors are i.i.d. with absolute moment of order r, 0<r? 2. The assumptions presented for the random error sequence will permit us to obtain improvements of the conditions on the regressors in order to obtain the strong consistency of the least-squares estimates in linear and nonlinear regression models.