Bayesian estimation of renewal function for inverse Gaussian renewal process

Two approximation methods are used to obtain the Bayes estimate for the renewal function of inverse Gaussian renewal process. Both approximations use a gamma-type conditional prior for the location parameter, a non-informative marginal prior for the shape parameter, and a squared error loss function. Simulations compare the accuracy of the estimators and indicate that the Tieney and Kadane (T–K)-based estimator out performs Maximum Likelihood (ML)- and Lindley (L)-based estimator. Computations for the T–K-based Bayes estimate employ the generalized Newton's method as well as a recent modified Newton's method with cubic convergence to maximize modified likelihood functions. The program is available from the author.     

A Simulation Based Study for Comparing Tests Associated With Receiver Operating Characteristic (ROC) Curves

Receiver Operating Characteristic curves and the Area Under Curve (AUC) are widely used to evaluate the predictive accuracy of diagnostic tests. The parametric methods of estimating AUCs are well established while nonparametric methods, such as Wilcoxon's method, lack proper research. This study considered three standard error techniques, namely, Hanley and McNeil, Hanley and Tilaki, and DeLong methods. Several parameters were considered, while measuring the predictor on a binary scale. The normality and type I error rate was violated for Hanley and McNeil's method while asymptotically DeLong's method performed better. Hanley and Tilaki's Jackknife method and DeLong's method performed equally well.     

Inferences on regression coefficients in a regression model under heteroscedasticity and robustness with respect to departure from normality

In this paper we consider a simple linear regression model under heteroscedasticity and nonnormality. A statistical test for testing the regression coefficient is then derived by assuming normality for the random disturbances and by applying Welch's method. Some Monte Carlo studies are generated for assessing robustness of this test. By combining Tiku's robust procedure with the new test, a robust but more powerful test is developed.     

Parameter Estimation for First-Order Random Coefficient Autoregressive (RCA) Models Based on Kalman Filter

In this article, we propose a new estimate algorithm for the parameters of a first-order Random Coefficient Autoregressive (RCA) Model. This algorithm turns out to be very reliable in estimating the true parameter values of a given model. It combines quasi-maximum likelihood method, the Kalman filter algorithm, and the Powell's method. Simulation results demonstrate that the algorithm is viable and promising.     

An efficient algorithm for computing the noncentrality parameters of chi-squared tests

Use of Newton's method for computing the noncentrality parameter based on the specified power in sample size problems of chi-squared tests requires that we evaluate both the noncentral chi-squareddistribution function and its derivative with respect to the noncentrality parameter. A close relationship between computing formulas for them is revealed, by which their evaluations can be performed jointly. This property greatly reduces the amount of computation involved. The corresponding algorithm is provided in a step-by-step form.     

On Confidence Intervals for Process Capability Indices in a One-Way Random Model

In this article, we investigated the bootstrap calibrated generalized confidence limits for process capability indices C _pk for the one-way random effect model. Also, we derived Bissell's approximation formula for the lower confidence limit using Satterthwaite's method and calculated its coverage probabilities and expected values. Then we compared it with standard bootstrap (SB) method and generalized confidence interval method. The simulation results indicate that the confidence limit obtained offers satisfactory coverage probabilities. The proposed method is illustrated with the help of simulation studies and data sets.     

A method for testing interaction in unreplicated two-way tables: using all pairwise interaction contrasts

Several methods exist for testing interaction in unreplicated two-way layouts. Some are based on specifying a functional form for the interaction term and perform well provided that the functional form is appropriate. Other methods do not require such a functional form to be specified but only test for the presence of non-additivity and do not provide a suitable estimate of error variance for a non-additive model. This paper presents a method for testing for interaction in unreplicated two-way tables that is based on testing all pairwise interaction contrasts. This method (i) is easy to implement, (ii) does not assume a functional form for the interaction term, (iii) can find a sub-table of data which may be free from interaction and to base the estimate of unknown error variance, and (iv) can be used for incomplete two-way layouts. The proposed method is illustrated using examples and its power is investigated via simulation studies. Simulation results show that the proposed method is competitive with existing methods for testing for interaction in unreplicated two-way layouts.     

On estimating the dimensionality in discriminant analysis

In discriminant analysis, the dimension of the hyperplane which population mean vectors span is called the dimensionality. The procedures commonly used to estimate this dimension involve testing a sequence of dimensionality hypotheses as well as model fitting approaches based on (consistent) Akaike's method, (modified) Mallows' method and Schwarz's method. The marginal log-likelihood (MLL) method is developed and the asymptotic distribution of the dimensionality estimated by this method for normal populations is derived. Furthermore a modified marginal log-likelihood (MMLL) method is also considered. The MLL method is not consistent for large samples and two modified criteria are proposed which attain asymptotic consistency. Some comments are made with regard to the robustness of this method to departures from normality. The operating characteristics of the various methods proposed are examined and compared.     

On the ratio of two independent skewnesses

In many fields, the researchers are interested in making inferences about the ratio of skewnesses in two independent populations. In the present paper, the asymptotic distribution for the ratio of the sample skewnesses in two independent populations is established. Then the asymptotic distribution is used to derive the asymptotic confidence interval and to test the hypothesis for the ratio of population's skewnesses. Finally, the applicability of the proposed method is investigated through Monte Carlo simulations.     

Test and prediction in factorial models with independent variance estimates

The multiple inference character of several tests in the same application is usually taken into consideration by requiring that the tests have a multiple level of significance. Also, a prediction problem in an application with several possible predictor variables requires that the multiple inference character of the problem be considered. This is not being done in the methods commonly used to choose predictor variables. Here, we discuss both the test and prediction methods in two-level factorial designs and suggest a principle for choosing variables which is based on multiple inference thinking. By an example use demonstrated that the principle proposed leads to the use of fewer prediction variables than does the Akaike method.     

Non Uniform Bounds on Geometric Approximation Via Stein's Method and w-Functions

In this article, we use Stein's method and w-functions to give uniform and non uniform bounds in the geometric approximation of a non negative integer-valued random variable. We give some applications of the results of this approximation concerning the beta-geometric, Pólya, and Poisson distributions.     

Laplace approximations to means and variances with asymptotic modes

The moment-generating function method, which is proposed by Tierney et al. [1989a. Fully exponential Laplace approximations to expectations and variances of nonpositive functions. J. Amer. Statist. Assoc. 84, 710–716], is an asymptotic technique of approximating a posterior mean of a general function by approximating the moment-generating function (MGF), and then differentiating it. In this article, we give approximations to the posterior means and variances by combining the MGF method and the Laplace approximations with asymptotic modes. We prove that asymptotic errors of the approximate means and variances are of order n^-2$n^{-2}$ and of order n^-3$n^{-3}$, respectively. Our approximation is closely related to a standard-form approximation, and is given without evaluating the exact posterior mode and third derivatives of the log-likelihood function. The MGF method also improves numerical instability of the fully exponential Laplace approximation for a predictive mean in logistic regression.     

Simultaneous testing of parameter subsets in less than full rank models

If the experimental design (or lack of design) results in a nodel which is not of f u l l rank, the problem of variable selection becomes rather complex. Tkis is due to the fact that, in less than f u l l rank models, not every linear combination of regression parameters is estimable. In this paper, we present a procedure fortesting all “testable” subsets of a complete set of regression parameters, using a technique based on Scheffe's method (1959). A class of “adequate” subsets of regression parameters is obtained in a manner similar to that of Aitkin (1974). The proposed procedure is illustraced with an example.     

Asymptotic hypothesis test to simultaneously compare the weighted kappa coefficients of multiple binary diagnostic tests in the presence of ignorable missing data

The weighted kappa coefficient of a binary diagnostic test is a measure of the beyond-chance agreement between the diagnostic test and the gold standard, and is a measure that allows us to assess and compare the performance of binary diagnostic tests. In the presence of partial disease verification, the comparison of the weighted kappa coefficients of two or more binary diagnostic tests cannot be carried out ignoring the individuals with an unknown disease status, since the estimators obtained would be affected by verification bias. In this article, we propose a global hypothesis test based on the chi-square distribution to simultaneously compare the weighted kappa coefficients when in the presence of partial disease verification the missing data mechanism is ignorable. Simulation experiments have been carried out to study the type I error and the power of the global hypothesis test. The results have been applied to the diagnosis of coronary disease.     

Testing the difference between two independent regression models

ABSTRACT

In some situations, for example, in biology or psychology studies, we wish to determine whether the linear relationship between response variable and predictor variables differs in two populations. The analysis of the covariance (ANCOVA) or, equivalently, the partial F-test approaches are the commonly used methods. In this study, the asymptotic distribution for the difference between two independent regression coefficients was established. The proposed method was used to derive the asymptotic confidence set for the difference between coefficients and hypothesis testing for the equality of the two regression models. Then a simulation study was conducted to compare the proposed method with the partial F method. The performance of the new method was comparable with that of the partial F method.     

Robust estimation for the coefficient of a first order autoregressive process

Nonparametric methods, Theil's method and Hussain's method have been applied to simple linear regression problems for estimating the slope of the regression line.We extend these methods and propose a robust estimator to estimate the coefficient of a first order autoregressive process under various distribution shapes, A simulation study to compare Theil's estimator, Hus-sain's estimator, the least squares estimator, and the proposed estimator is also presented.     

A family of IDMRL tests with unknown turning point

In this paper we propose a family of tests for exponentiality against the IDMRL alternative. Here we assume that the turning point or the proportion before the turning point is unknown. We derive the asymptotic null distributions of the test statistics and obtain their asymptotic critical values based on Durbin's approximation method. A simulation study is conducted to evaluate the proposed tests.     

ALGORITHMS FOR RETURN PROBABILITIES FOR STOCHASTIC FLUID FLOWS

Abstract

We consider several known algorithms and introduce some new algorithms that can be used to calculate the probability of return to the initial level in the Markov stochastic fluid flow model. We give the physical interpretations of these algorithms within the fluid flow environment. The rates of convergence are explained in terms of the physical properties of the fluid flow processes. We compare these algorithms with respect to the number of iterations required and their complexity. The performance of the algorithms depends on the nature of the process considered in the analysis. We illustrate this with examples and give appropriate recommendations.     

A note on scheffé's method of deriving confidence sets for directions and ratios of normals

Scheffé (1970) introduced a method for deriving confidence sets for directions and ratios of normals. The procedure requires use of an approximation and Scheffé provided evidence that the method performs well for cases in which the variances of the random deviates are known. This paper extends Scheffé's numerical integrations to the case of unknown variances. Our results indicate that Scheffé's method works well when variances are unknown