Search

A kappa-like classification statistic is proposed for assessing the fit of GEE regression models with a categorical response. The proposed statistic is a summary measure depicting how well categorical responses are predicted from the fitted GEE model. The statistic takes on a value of 1 if prediction is perfect and a value of 0 if the fitted model fares no better than random chance, i.e., fitting the repeated categorical responses with an intercept-only model. To demonstrate the usefulness of the classification statistic, we present simulation results as well as two examples from biomedical studies.

We develop a likelihood ratio test statistic, based on the beta- binomial distribution, for comparing a single treated group with dichoto- mous data to dual control groups. This statistic is useful in cases where there is overdispersion or extra-binomial variation. We apply the statis- tic to data from a two year rodent carcinogenicity study with dual control groups. The test statistic we developed is similar to others that have been developed for incorporation of historical control groups with rodent carcino- genicity experiments. However, for the small sample case we considered, large sample theory used by the other test statistics did not apply. We de- termined the critical values of this statistic by enumerating its distribution. A small Monte Carlo study shows the new test statistic controls the signifi- cance level much better than Fisher's exact test when there is overdispersion and that it has adequate power.

For a two-group comparative study, a stratified inference procedure is routinely used to estimate an overall group contrast to increase the precision of the simple two-sample estimator. Unfortunately, most commonly used methods including the Cochran-Mantel-Haenszel statistic for a binary outcome and the stratified Cox procedure for the event time endpoint do not serve this purpose well. In fact, these procedures may be worse than their two-sample counterparts even when the observed treatment allocations are imbalanced across strata. Various procedures beyond the conventional stratified methods have been proposed to increase the precision of estimation when the naive estimator is consistent. In this paper, we are interested in the case when the treatment allocation proportions vary markedly across strata. We study the stochastic properties of the two-sample naive estimator conditional on the ancillary statistics, the observed treatment allocation proportions and/or the stratum sizes, and present a biased-adjusted estimator. This adjusted estimator is asymptotically equivalent to the augmentation estimators proposed under the unconditional setting. Moreover, this consistent estimation procedure is also equivalent to a rather simple procedure, which estimates the mean response of each treatment group first via a stratum-size weighted average and then constructs the group contrast estimate. This simple procedure is flexible and readily applicable to any target patient population by choosing appropriate stratum weights. All the proposals are illustrated with the data from a cardiovascular clinical trial, whose treatment allocations are imbalanced.

Detection of multiple outliers or subset of influential points has been rarely considered in the linear measurement error models. In this paper a new influence statistic for one or a set of observations is generalized and characterized based on the corrected likelihood in the linear measurement error models. This influence statistic can be expressed in terms of the residuals and the leverages of linear measurement error regression. Unlike Cook’s statistic, this new measure of influence has asymptotically normal distribution and is able to detect a subset of high leverage outliers which is not identified by Cook’s statistic. As an illustrative example, simulation studies and a real data set are analysed.

New statistics are proposed to estimate and test the structural change when the data dimension is comparable to or larger than the sample size. Consistency of the new statistic in estimating the change point position is established under the alternative hypothesis. The asymptotic distribution of the new statistic in testing the existence of a change point is obtained under the null hypothesis. Some simulation results are presented which show that the numerical performance of our method is satisfactory. The method is illustrated via the analysis of the house price index of US.

In this study, a simple diagnostic statistic for determining the size of random forest is proposed. This method is based on MV (margin of victory), a scaled difference in the votes at the infinite forest between the first and second most popular categories of the current random forest. We can note that if MV is negative then there is discrepancy between the current and infinite forests. More precisely, our method is based on the proportion of cases that -MV is greater than a fixed small positive number (say, 0.03). We derive an appropriate diagnostic statistic for our method and then calculate the distribution of the statistic. A simulation study is performed to compare our method with a recently proposed diagnostic statistic.

A test for ordered categorical variables is of considerable importance, because they are frequently encountered in biomedical studies. This paper introduces a simple ordering test approach for the two-way r × c contingency tables with incomplete counts by developing six test statistics, i.e., the likelihood ratio test statistic, score test statistic, global score test statistic, Hausman-Wald test statistic, Wald test statistic and distance-based test statistic. Bootstrap resampling methods are also presented. The performance of the proposed tests is evaluated with respect to their empirical type I error rates and empirical powers. The results show that the likelihood ratio test statistic based on the bootstrap resampling methods perform satisfactorily for small to large sample sizes. A real example from a wheeze study in six cities is used to illustrate the proposed methodologies.

In practice, it is often not possible to find an appropriate family of distributions which can be used for fitting the sample distribution with high precision. In these cases, it seems to be opportune to search for the best approximation by a family of distributions instead of an exact fit. In this paper, we consider the Anderson–Darling statistic with plugged-in minimum distance estimator for the parameter vector. We prove asymptotic normality of the Anderson–Darling statistic which is used for a test of goodness of approximation. Moreover, we introduce a measure of discrepancy between the sample distribution and the model class.

The statistician often needs to test whether or not a time series has a seasonal first moment. The problem often arises in environmental series, where most time-ordered data display some type of periodic structure. This paper reviews the problem, proposing new statistics in both the time and frequency domains. Our new time domain statistic has an analysis of variance form that is based on the one-step-ahead prediction errors of the series. This statistic inherits the classic traits of the F-distribution arising in one-way analysis of variance tests, is easy to use, and is asymptotically equivalent to the likelihood ratio test. The statistics asymptotic distribution is quantified when time series parameters are estimated. In the frequency domain, a statistic modifying Fisher's classical test for a sinusoidal mean superimposed on independent and identically distributed Gaussian noise is devised. The performance and comparison of these statistics are studied via simulation. Implementation of the methods merely requires sample means, autocovariances, and periodograms of the series. Application to a data set of monthly temperatures from Tuscaloosa, Alabama, is given. Copyright © 2016 John Wiley & Sons, Ltd.

A new method is proposed for identifying clusters in continuous data indexed by time or by space. The scan statistic we introduce is derived from the well-known Mann–Whitney statistic. It is completely non parametric as it relies only on the ranks of the marks. This scan test seems to be very powerful against any clustering alternative. These results have applications in various fields, such as the study of climate data or socioeconomic data.

This article reexamines ANOVA problems for autocorrelated data. Using linear prediction techniques for stationary time series, a new test statistic that assesses a null hypothesis of equal means is proposed and investigated. Our test statistic mimics the classical F-type ratio form used with independent data, but substitutes estimated prediction residuals in for the errors. This simple tactic departs from past studies that adjust the quadratic forms in the numerator and denominator in the F ratio for autocorrelation. One of the advantages is that our statistic retains the classical null hypothesis F distribution (now as a limit) with the customary degrees of freedom. The statistic is shown to perform well in simulations. Asymptotic proofs are given in the case of autoregressive random errors; a sports application is supplied.[Received December 2014. Revised August 2015.]

In this paper, we consider the consistency of Cp-type statistics for selecting variables in a normality-assumed linear regression with multiple responses when the dimension of the vector of the response variables may be large. We propose a new consistent Cp-type statistic for which consistency can be achieved whenever the dimension of the response variables vector is fixed or goes to infinity. A high probability of selecting the true subset of explanatory variables can be expected under a moderate sample size when the proposed Cp-type statistic is used to select variables, even when there is a high-dimensional response variables vector.

The C statistic is a commonly reported measure of screening test performance. Optimistic estimation of the C statistic is a frequent problem because of overfitting of statistical models in small data sets, and methods exist to correct for this issue. However, many studies do not use such methods, and those that do correct for optimism use diverse methods, some of which are known to be biased. We used clinical data sets (United Kingdom Down syndrome screening data from Glasgow (1991-2003), Edinburgh (1999-2003), and Cambridge (1990-2006), as well as Scottish national pregnancy discharge data (2004-2007)) to evaluate different approaches to adjustment for optimism. We found that sample splitting, cross-validation without replication, and leave-1-out cross-validation produced optimism-adjusted estimates of the C statistic that were biased and/or associated with greater absolute error than other available methods. Cross-validation with replication, bootstrapping, and a new method (leave-pair-out cross-validation) all generated unbiased optimism-adjusted estimates of the C statistic and had similar absolute errors in the clinical data set. Larger simulation studies confirmed that all 3 methods performed similarly with 10 or more events per variable, or when the C statistic was 0.9 or greater. However, with lower events per variable or lower C statistics, bootstrapping tended to be optimistic but with lower absolute and mean squared errors than both methods of cross-validation.

Barndorff-Nielsen’s celebrated p*-formula and variations thereof have amongst their various attractions the ability to approximate bimodal distributions. In this paper we show that in general this requires a crucial adjustment to the basic formula. The adjustment is based on a simple idea and straightforward to implement, yet delivers important improvements. It is based on recognizing that certain outcomes are theoretically impossible and the density of the MLE should then equal zero, rather than the positive density that a straight application of p* would suggest. This has implications for inference and we show how to use the new p**-formula to construct improved confidence regions. These can be disjoint as a consequence of the bimodality. The degree of bimodality depends heavily on the value of an approximate ancillary statistic and conditioning on the observed value of this statistic is therefore desirable. The p**-formula naturally delivers the relevant conditional distribution. We illustrate these results in small and large samples using a simple nonlinear regression model and errors in variables model where the measurement errors in dependent and explanatory variables are correlated and allow for weak proxies.

Tests of goodness of fit of sparse multinomial models with non-canonical links is proposed by using approximations to the first three moments of the conditional distribution of a modified Pearson Chi-square statistic. The modified Pearson statistic is obtained using a supplementary estimating equation approach. Approximations to the first three conditional moments of the modified Pearson statistic are derived. A simulation study is conducted to compare, in terms of empirical size and power, the usual Pearson Chi-square statistic, the standardized modified Pearson Chi-square statistic using the first two conditional moments, a method using Edgeworth approximation of the p-values based on the first three conditional moments and a score test statistic. There does not seems to be any qualitative difference in size of the four methods. However, the standardized modified Pearson Chi-square statistic and the Edgeworth approximation method of obtaining p-values using the first three conditional moments show power advantages compared to the usual Pearson Chi-square statistic, and the score test statistic. In some situations, for example, for small nominal level, the standardized modified Pearson Chi-square statistic shows some power advantage over the method using Edgeworth approximation of the p-values using the first three conditional moments. Also, the former is easier to use and so is preferable. Two data sets are analyzed and a discussion is given.

In this article, several scan statistics are discussed for detecting a local change in variance for one dimensional normal data. When the length of the scanning window is known, a fixed window scan statistic based on moving sum of squares is proposed. Two approximations for the distribution of this scan statistic are investigated. When the length of the scanning window is unknown, a variable window scan statistic based on a generalized likelihood ratio test and a multiple window minimum P-value scan statistic are proposed for detecting the local change in variance. For a moderate or large shift in variance, numerical results indicate that both the variable and multiple window scan statistics perform well. For large data sets, considering the detection power and computing efficiency, the multiple window scan statistic is recommended.

Pearson's Chi-square statistic for frequency tables depends on what is hypothesized as the expected frequencies. Its partitions also depend on the hypothesis. Lancaster (J R Stat Soc B 13:242---249, 1951) proposed ANOVA-like partitions of Pearson's statistic under several representative hypotheses about the expected frequencies. His expositions were, however, not entirely clear. In this paper, we clarify his method of derivations, and extend it to more general situations. A comparison is made with analogous decompositions of the log likelihood ratio statistic associated with log-linear analysis of contingency tables.

This paper proposes a new simple panel unit-root test by extending the cross-sectionally augmented panel unit-root test (CIPS) developed by Pesaran et al. (2013) to allow for smoothing structural changes in deterministic terms, approximated by a Fourier series. The proposed statistic is the simple average of the individual statistics constructed from the breaks and cross-sectional dependence augmented Dickey-Fuller (BCADF) regression and is called the BCIPS statistic. We initially develop the tests by assuming that the number of factors in the model is known and show that the limiting distribution of the BCADF statistic is free of nuisance parameters. The nonstandard limiting distribution of the (truncated) BCIPS statistic is also shown to exist and its critical values are tabulated. Monte-Carlo experiments point out that the sizes and powers of the BCIPS statistic are generally satisfactory as long as T is greater than or equal to fifty and a hundred, respectively. By using two different methods to determine the number of factors, both the BCIPS and CIPS tests are applied to examine the validity of long-run purchasing power parity. The proposed test complements the panel unit-root tests with breaks using dummy variables.

A new characterization of the Pareto distribution is proposed, and new goodness-of-fit tests based on it are constructed. Test statistics are functionals of U-empirical processes. The first of these statistics is of integral type, it is similar to the classical statistics $$\omega _n^1$$ . The second one is a Kolmogorov type statistic. We show that the kernels of our statistics are non-degenerate. The limiting distribution and large deviations asymptotics of the new statistics under null hypothesis are described. Their local Bahadur efficiency for parametric alternatives is calculated. This type of efficiency is mostly appropriate for the solution of our problem since the Kolmogorov type statistic is not asymptotically normal, and the Pitman approach is not applicable to this statistic. For the second statistic we evaluate the critical values by using Monte-Carlo methods. Also conditions of local optimality of new statistics in the sense of Bahadur are discussed and examples of such special alternatives are given. For small sample size we compare the power of those tests with some common goodness-of-fit tests.

As an important extension of the regular kappa statistic, the weighted kappa statistic has been widely used to assess the agreement between two procedures for independent matched-pair ordinal data. For clustered matched-pair ordinal data, based on the delta method and sampling techniques, a non-parametric variance estimator for the weighted kappa statistic is proposed without within-cluster correlation structure or distributional assumptions. The results of an extensive Monte Carlo simulation study demonstrate that the proposed weighted kappa statistic provides consistent estimation, and the proposed variance estimator behaves reasonably well for at least a moderately large number of clusters (e.g., K ≥ 50 ). Compared with the variance estimator ignoring dependence within a cluster, the proposed variance estimator performs better in maintaining the nominal coverage probability when the intra-cluster correlation is fair ( ρ ≥ 0.3 ), with more pronounced improvement when ρ is further increased. Moreover, under the general analysis of variance setting with systematic variability between procedures and clusters being included as a component of total variation, the equivalence between weighted kappa statistic and intraclass correlation coefficient is established. To illustrate the practical application of the proposed estimator, two real medical research data examples of clustered matched-pair ordinal data are analyzed, including an agreement study to compare two methods for assessing cervical ectopy, and a physician–patients data example from the Enhancing Communication and HIV Outcomes study.

In normal theory exploratory factor analysis, likelihood ratio (LR) statistic plays an important role in evaluating the goodness-of-fit of the model. In this paper, we derive an approximation of the LR statistic. The approximation is then used to show explicitly that the expectation of the LR statistic agrees with the degrees of freedom of the asymptotic chi-square distribution.

A U-statistic for the tail index of a multivariate stable random vector is given as an extension of the univariate case introduced by Fan (2006). Asymptotic normality and consistency of the proposed U-statistic for the tail index are proved theoretically. The proposed estimator is used to estimate the spectral measure. The performance of both introduced tail index and spectral measure estimators is compared with the known estimators by comprehensive simulations and real datasets.

Comparing two suppliers for linear profiles is a very important task for supplier management. The difference test statistic based on the process-yield index is proposed to tackle the better process selection for linear profiles. A simple form of the sampling distribution of the process-yield index is derived by a simulation study. Therefore, the asymptotic normal distribution of the difference test statistic is established. The results provide useful information to practitioners. An example from the leather industry is presented to illustrate the applicability of the proposed method. Copyright © 2014 John Wiley & Sons, Ltd.

A 2-stage robust procedure as well as an R package, rsem, were recently developed for structural equation modeling with nonnormal missing data by Yuan and Zhang (2012). Several test statistics that have been used for complete data analysis are employed to evaluate model fit in the 2-stage robust method. However, properties of these statistics under robust procedures for incomplete nonnormal data analysis have never been studied. This study aims to systematically evaluate and compare 5 test statistics, including a test statistic derived from normal-distribution-based maximum likelihood, a rescaled chi-square statistic, an adjusted chi-square statistic, a corrected residual-based asymptotical distribution-free chi-square statistic, and a residual-based F statistic. These statistics are evaluated under a linear growth curve model by varying 8 factors: population distribution, missing data mechanism, missing data rate, sample size, number of measurement occasions, covariance between the latent intercept and s...

Single moments of order statistics from the modified Makeham distribution (MMD) are derived, an identity about the single moments of order statistics is given, and the specific expected value and variance of the single moments of order statistics from the MMD are calculated. In this study, the order statistic from the MMD was applied to the rank sum test in a two-sample problem. The exact critical values of the designated statistics were evaluated. Simulations were used to investigate the power of these statistics for the two-sided alternative with several population distributions. The powers of the statistics were compared with the Wilcoxon rank sum statistic, the Lepage statistic, the modified Baumgartner statistic, the Savage test and the normal score test. The Edgeworth expansion was used to evaluate the upper tail probability for the preferred statistic, given finite sample sizes.

This paper proposes an overlapping-based test statistic for testing the equality of two exponential distributions with different scale and location parameters. The test statistic is defined as the maximum likelihood estimate of the Weitzman's overlapping coefficient, which estimates the agreement of two densities. The proposed test statistic is derived in closed form. Simulated critical points are generated for the proposed test statistic for various sample sizes and significance levels via Monte Carlo Simulations. Statistical powers of the proposed test are computed via simulation studies and compared to those of the existing Log likelihood ratio test.

Hotelling's T2 statistic is the default control statistic for continuous multivariate data, but there are dangers in applying this statistic without the appropriate level of checks and balances. This paper discusses the potential issues with using the Hotelling's T2 statistic when the quality variable measures are highly correlated and provides some solutions that will help mitigate the risks with applying the Hotelling's T2 control charts in such practical examples. Copyright © 2014 John Wiley & Sons, Ltd.

Hotelling's T2 statistic has many applications in multivariate analysis. In particular, it can be used to measure the influence that a particular observation vector has on parameter estimation. For example, in the bivariate case, there exists a direct relationship between the ellipse generated using a T2 statistic for individual observations and the hyperbolae generated using Hampel's influence function for the corresponding correlation coefficient. In this paper, we jointly use the components of an orthogonal decomposition of the T2 statistic and some influence functions to identify outliers or influential observations. Since the conditional components in the T2 statistic are related to the possible changes in the correlation between a variable and a group of other variables, we consider the theoretical influence functions of the correlations and multiple correlation coefficients. Finite-sample versions of these influence functions are used to find the estimated influence function values.

The scan statistic is widely used in spatial cluster detection applications of inhomogeneous Poisson processes. However, real data may present substantial departure from the underlying Poisson process. One of the possible departures has to do with zero excess. Some studies point out that when applied to data with excess zeros, the spatial scan statistic may produce biased inferences. In this work, we develop a closed-form scan statistic for cluster detection of spatial zero-inflated count data. We apply our methodology to simulated and real data. Our simulations revealed that the Scan-Poisson statistic steadily deteriorates as the number of zeros increases, producing biased inferences. On the other hand, our proposed Scan-ZIP and Scan-ZIP+EM statistics are, most of the time, either superior or comparable to the Scan-Poisson statistic.

The discrete scan statistic is used in many areas of applied probability and statistics to study local clumping of patterns. Testing based on the statistic requires tail probabilities. Whereas the distribution has been studied extensively, most of the results are approximations, due to the difficulties associated with the computation. Results for exact p-values for the statistic have been given for a binary sequence that is independent or first-order Markovian. We give an algorithm to obtain probabilities for the statistic over multi-state trials that are Markovian of a general order of dependence, and explore the algorithm's usefulness.

The PARSCALE G2 statistic is arguably the most popular item fit statistic in operational testing. For long tests, the Type I error rates of the statistic have often been found to be satisfactory. However, the Type I error rates of the statistic have only been studied for sample sizes of up to several thousands. The authors examined the Type I error rates of the PARSCALE G2 statistic in a simulation study using sample sizes much larger than those considered in the literature. For any fixed test length, the Type I error rate of the PARSCALE G2 statistic is found to increase to 1 as the sample size increases. The findings contradict the claim in the PARSCALE software manual that the PARSCALE G2 statistic leads to a large-sample test and also contradict the common belief that the statistic has reasonable Type I error rates for long tests. Thus, this simulation study conveys the important practical message that the use of the PARSCALE G2 statistic cannot always be recommended even for long tests. The Type I error rates of the item fit statistics of Orlando and Thissen were found to be close to the nominal level for all simulation conditions considered here.

A wide literature is available on the asymptotic behavior of the Durbin-Watson statistic for autoregressive models. However, it is impossible to find results on the Durbin-Watson statistic for autoregressive models with adaptive control. Our purpose is to fill the gap by establishing the asymptotic behavior of the Durbin Watson statistic for ARX models in adaptive tracking. On the one hand, we show the almost sure convergence as well as the asymptotic normality of the least squares estimators of the unknown parameters of the ARX models. On the other hand, we establish the almost sure convergence of the Durbin-Watson statistic and its asymptotic normality. Finally, we propose a bilateral statistical test for residual autocorrelation in adaptive tracking.

This study deals with the small sample likelihood based inference for the ratio of twonormal variances. The small sample likelihood inference is an approximation method.The signed log-likelihood ratio statistic and the modi ed signed log-likelihood ratiostatistic, which converge to standard normal distribution, are proposed for the normalvariance ratio. Through the simulation study, the coverage probabilities of con denceinterval and power of the exact, the signed log-likelihood and the modi ed signed log-likelihood ratio statistic will be compared. A real data example will be provided.Keywords: Likelihood based inference, modi ed signed log-likelihood ratio statistic,normal variance ratio, signed log-likelihood ratio statistic. 1. Introduction The normal distribution plays an important role in statistical inference, and there are somany studies about this distribution. The statistical inference for the ratio of two normalvariances arises in the areas for comparing the precision of two independent normal popula-tions. About the ratio of two normal variances, there also exists an exact statistical inferencewith Fstatistic. This problem is simple, obvious and important.When comparing the dispersion of two normal population, one can use the distributiontable of F statistic. Since F distribution depends on degrees of freedom, the distributiontable of Fdistribution is more than several pages. For the degrees of freedom exceed 30, thetable does not have the value. In this case, one must use an interpolation to obtaining thequantile.An approximation of a statistic to standard normal distribution has been developed inmany statistical models. Since the percentile of standard normal distribution is well known,a statistic, which distributes as standard normal distribution asymptotically, may be veryuseful. But these statistics have the error rate depending on sample size. When the samplesize is small, these approximation to standard normal distribution is quite inaccurate. Forexample, a signed log-likelihood ratio statistic converges to standard normal distributionwith an error of O(n

In this paper, multivariate two-sample testing problems were examined based on the Jureckova–Kalina's ranks of distances. The multivariate two-sample rank test based on the modified Baumgartner statistic for the two-sided alternative was proposed. The proposed statistic was a randomized statistic. Simulations were used to investigate the power of the suggested statistic for various population distributions.

A great deal of inference in statistics is based on making the approximation that a statistic is normally distributed. The error in doing so is generally O(n−1/2), where n is the sample size and can be considered when the distribution of the statistic is heavily biased or skewed. This note shows how one may reduce the error to O(n−(j+1)/2), where j is a given integer. The case considered is when the statistic is the mean of the sample values of a continuous distribution with a scale or location change after the sample has undergone an initial transformation, which may depend on an unknown parameter. The transformation corresponding to Fisher's score function yields an asymptotically efficient procedure.

This paper focuses on inferential tools in the logistic regression model fitted by the Firth penalized likelihood. In this context, the Likelihood Ratio statistic is often reported to be the preferred choice as compared to the ‘traditional’ Wald statistic. In this work, we consider and discuss a wider range of test statistics, including the robust Wald, the Score, and the recently proposed Gradient statistic. We compare all these asymptotically equivalent statistics in terms of interval estimation and hypothesis testing via simulation experiments and analyses of two real datasets. We find out that the Likelihood Ratio statistic does not appear the best inferential device in the Firth penalized logistic regression.

An explicit statement of a statistic which is a nonparametric analogue to one-way MANOVA is presented. The statistic is a multivariate extension of the nonparametric Kruskal-Wallis test (1952). The large sample reference distribution of the test statistic is derived together with a set of computational formulas for the test statistic. In addition two post hoc procedures are developed and compared. The statistic and its post hoc procedures are illustrated with a data example from the behavioral sciences.

In this chapter, we introduce the new likelihood-based test statistic named as the gradient statistic . The test that uses the gradient statistic is a newly large-sample test which was introduced in the statistic literature by George R. Terrell. The gradient statistic is very simple to be computed and it can be an interesting alternative to the classical statistics, namely the likelihood ratio, Wald, and Rao score statistics. We provide some properties of the gradient statistic and present several examples.

Analysis of variance (ANOVA) is a common use parametric method to test the differences in means for more than two groups when the populations are normally distributed. ANOVA is highly inefficient under the influence of non- normal and heteroscedastic settings. When the assumptions are violated, researchers are looking for alternative such as Kruskal-Wallis under nonparametric or robust method. This study focused on flexible method, S1 statistic for comparing groups using median as the location estimator. S1 statistic was modified by substituting the median with Hodges-Lehmann and the default scale estimator with the variance of Hodges-Lehmann and MADn to produce two different test statistics for comparing groups. Bootstrap method was used for testing the hypotheses since the sampling distributions of these modified S1 statistics are unknown. The performance of the proposed statistic in terms of Type I error was measured and compared against the original S1 statistic, ANOVA and Kruskal-Wallis. The propose ...

In Chapter 6, the statistical methods of residual diagnostics are described first. A variety of residual types are specified, both generally and with specific regard to longitudinal data analysis. I also delineate semi-variogram, a popular residual diagnostic technique applied in longitudinal data analysis. This unique method is used to check whether serial correlation among repeated measurements for the same subject is present given the specified fixed effects and the covariance parameters. The method can be specified in both the random intercept and the random coefficient perspectives. Next, the statistical methods on influence diagnostics are displayed. The variety of the basic diagnostics to identify influential observations include the Cook’s distance statistic, leverage, the DFFITS score, the MDFFITS statistic, COVTRACE, COVRATIO, the likelihood displacement statistic (the LD statistic), and the LMAX standardized statistic. An empirical illustration is provided for displaying how to check whether there are any influential observations in fitting a linear mixed model.

We study the exact distribution of the likelihood-ratio statistic used for testing a normal sample for three upper (lower) outliers. We obtain recursive correlations for the integral distribution function of this statistic. We apply the obtained correlations for calculating critical values of the likelihood-ratio statistic which appear to be close to critical values of this statistic simulated by the Monte Carlo method. We give an example of the joint use of the likelihood-ratio statistic for testing a sample for more than one outlier.

The purpose of this paper is to assess goodness of fit properties of the conditional distribution of the Pearson statistic, for noncanonical generalized linear models with data that are extensive but sparse, by Edgeworth approximation of p-values using higher order moment corrections. In this paper we obtain approximations to the fourth moments of the unconditional and conditional distribution of the modified Pearson statistic with noncanonical links. This extends previous results where approximations to the first three moments are available and completes all usual higher order moment calculations of the modified Pearson statistic. We consider the asymptotic limit in which the data are extensive but sparse and a supplementary estimating equation for the dispersion parameter. Specific results for binomial and Poisson data are obtained separately. The methods for assessing goodness of fit using higher order moments are discussed. For testing goodness of fit of generalized linear models to sparse data some simulations are conducted to compare, in terms of empirical size and power, the performance of the classical Pearson statistic (X 2), a standardized modified Pearson statistic ( $X_{\ast}^2$ ), a standardized modified deviance statistic (D ∗ ), a modified Pearson statistic based on Edgeworth approximation with the first three conditional moments (Z 1), and a modified Pearson statistic based on Edgeworth approximation with the first four conditional moments (Z 2). The statistic Z 2 holds level most effectively in most situations and has some power advantage.

This paper studies the applications of extreme value theory on analysis for closing price data of the Dow-Jones industrial index and Danish fire insurance claims data. The generalized extreme value (GEV) distribution is considered in analyzing the real data, and the hypothesis testing problem for the shape parameter of GEV distribution is investigated based on a new test statistic—the $$L_q$$ -likelihood ratio ( $$L_q$$ R) statistic. The $$L_q$$ R statistic can be treated as a generalized form of the classical likelihood ratio (LR) statistic. We show that the asymptotic behavior of proposed statistic is characterized by the degree of distortion $$q$$ . For small and modest sample sizes, the $$L_q$$ R statistic is still available when $$q$$ is properly chosen. By simulation studies, the proposed statistic not only performs the asymptotic properties, but also outperforms the classical LR statistic as the sample sizes are modest or even small. Meanwhile, the test power based on the new statistic is also simulated by Monte Carlo methods. At last, the models are diagnosed by graphical methods.

Under model correctness, highly accurate inference on a scalar interest parameter in the presence of a nuisance parameter can be achieved by several routes, among them considering the bootstrap distribution of the signed root likelihood ratio statistic. The context of model mis-specification is considered and inference based on a robust form of the signed root statistic is discussed in detail. Stability of the distribution of the statistic allows accurate inference, outperforming that based on first-order asymptotic approximation, by considering the bootstrap distribution of the statistic under the incorrectly assumed distribution. Comparisons of this simple approach with alternative analytic and non-parametric inference schemes are discussed.