Diagnostics for generalized Poisson regression models with errors in variables

In this paper, we develop diagnostic methods for generalized Poisson regression (GPR) models with errors in variables based on the corrected likelihood. The one-step approximations of the estimates in the case-deletion model are given and case-deletion and local influence measures are presented. Meanwhile, based on a corrected score function, the testing statistics for the significance of dispersion parameters in GPR models with measurement errors are investigated. Finally, illustration of our methodology is given through numerical examples.     

Remedying the Neyman–Scott phenomenon in model discrimination

The objective of this paper is to investigate through simulation the possible presence of the incidental parameters problem when performing frequentist model discrimination with stratified data. In this context, model discrimination amounts to considering a structural parameter taking values in a finite space, with k points, k≥2. This setting seems to have not yet been considered in the literature about the Neyman–Scott phenomenon. Here we provide Monte Carlo evidence of the severity of the incidental parameters problem also in the model discrimination setting and propose a remedy for a special class of models. In particular, we focus on models that are scale families in each stratum. We consider traditional model selection procedures, such as the Akaike and Takeuchi information criteria, together with the best frequentist selection procedure based on maximization of the marginal likelihood induced by the maximal invariant, or of its Laplace approximation. Results of two Monte Carlo experiments indicate that when the sample size in each stratum is fixed and the number of strata increases, correct selection probabilities for traditional model selection criteria may approach zero, unlike what happens for model discrimination based on exact or approximate marginal likelihoods. Finally, two examples with real data sets are given.     

Computations of Bayesian t-values for multiple comparisons

Results from a simulation study of the power of eight statistics for testing that a sample is form a uniform distribution on the unit interval are reported. Power is given for each statistic against four classes if alternatives. The statistics studied include the discrete Pearson chi-square with ten and twenty cells, X² ₁₀ and X² ₂₀; Kolmogorov-smirov, D; Cramer-Von Mises, W²; Watson, U²; Anderson-Darling, A; Greenwood. G;and a new statistic called O A modified form of each of these statistic is also studied by first transforming the sample using a transformation given by Durbin. On the basis of the results observed in this study, the Watson U² statistic is recommended as a general test for uniformity.     

Inference for Box–Cox Transformed Threshold GARCH Models with Nuisance Parameters

Abstract. Generalized autoregressive conditional heteroscedastic (GARCH) models have been widely used for analyzing financial time series with time‐varying volatilities. To overcome the defect of the Gaussian quasi‐maximum likelihood estimator (QMLE) when the innovations follow either heavy‐tailed or skewed distributions, Berkes & Horváth (Ann. Statist., 32, 633, 2004) and Lee & Lee (Scand. J. Statist. 36, 157, 2009) considered likelihood methods that use two‐sided exponential, Cauchy and normal mixture distributions. In this paper, we extend their methods for Box–Cox transformed threshold GARCH model by allowing distributions used in the construction of likelihood functions to include parameters and employing the estimated quasi‐likelihood estimators (QELE) to handle those parameters. We also demonstrate that the proposed QMLE and QELE are consistent and asymptotically normal under regularity conditions. Simulation results are provided for illustration.     

Small-sample comparison of the exact and asymptotic upper tail probabilities of chi-squared goodness-of-fit statistics: the binomila and the mixture binomial

The exact and asymptotic upper tail probabilities (α = .10, .05, .01, .001) of the three chi-squared goodness-of-fit statistics Pearson's X ², likelihood ratioG ², and powerdivergence statisticD ²(λ), with λ= 2/3 are compared by complete enumeration for the binomial and the mixture binomial. For the two-component mixture binomial, three cases have been distinguished. 1. Both success probabilities and the mixing weights are unknwon. 2. One of the two success probabilities is known. And 3., the mixing weights are known. The binomial was investigated for the number of cellsk, being between 3 and 6 with sample sizes between 5 and 100, for k = 7 with sample sizes between 5 and 45, and for k = 10 with sample sizes ranging from 5 to 20. For the mixture binomial, solely k = 5 cells were considered with sample sizes from 5 to 100 and k = 8 cells with sample sizes between 4 and 20. Rating the relative accuracy of the chi-squared approximation in terms of ±10% and ±20% intervals around α led to the following conclusions for the binomial: 1. Using G2 is not recommendable. 2. At the significance levels α=.10 and α=.05X ² should be preferred over D ²; D ² is the best choice at α = .01. 3. Cochran's (1954; Biometrics, 10, 417-451) rule for the minimum expectation when using X ² seems to generalize to the binomial for G ² and D ² ; as a compromise, it gives a rather strong lower limit for the expected cell frequencies in some circumstances, but a rather liberal in others. To draw similar conclusions concerning the mixture binomial was not possible, because in that case, the accuracy of the chi-squared approximation is not only a function of the chosen test statistic and of the significance level, but also heavily depends on the numerical value of theinvolved unknown parameters and on the hypothesis to be tested. Thereto, the present study may give rise only to warnings against the application of mixture models to small samples.     

Identification of the important parameters in thermal building simulation models

Two parameter screening techniques, a sequential bifurcation technique and a factorial sampling method, have been applied to a building thermal model, used to predict thermal comfort performance of a building in its design stage. Combined application of both screening methods revealed a set of 12 important model parameters out of a total of 81, explaining 94% of the variability in the model output. These important parameters were identified by the factorial sampling method on the basis of 246 model evaluations, while sequential bifurcation only needed 52 evaluations. However, the factorial sampling scheme was effective in identifying of not only the important parameters, but also the directions of parameter main effects and the severity of interaction effects. This additional information showed that isolated application of the sequential bifurcation method would have been unreliable, as satisfaction of the inherent assumptions could not be guaranteed. Only on the basis of proper knowledge of the sign of the parameter main effects, adequate clustering of important parameters and transformation of the model output, all obtained from the results of the factorial sampling scheme, reliable and economic application of sequential bifurcation was possible.     

The Computer Analysis of Factorial Experiments: In Memoriam—Frank Yates

The responses to a recent paper by Dallal in this journal are evaluated by reference to the ideas of Frank Yates. It is concluded that much unnecessary complication has been introduced into the computer analysis of linear models by (1) the imposition of constraints on parameters, (2) neglect of marginality relations in forming hypotheses, and (3) confusion between the form of noncentrality parameters and hypotheses.     

Loycell纤维织物设计探讨

通过对Loycell纤维性能的分析,阐述其织物的性能,提出了这类织物设计的主要结构参数.     

A Statistical Test of Change‐Point in Mean that Almost Surely Has Zero Error Probabilities

In this paper we develop a non‐conventional statistical test for the change‐point in a mean model by making use of an almost‐sure (a.s.) convergence (or strong convergence) result that we obtain, in respect of the difference between the sums of squared residuals under the null and alternative hypotheses. We prove that both types of error probabilities of the new test converge to zero almost surely when the sample size goes to infinity. This result does not hold for any conventional statistical test where the type I error probability, i.e. the significance level or the size, is prescribed at a low but non‐zero level (e.g. 0.05). The test developed is easy to use in practice, and is ready to be generalised to other change‐point models provided that the relevant almost‐sure convergence results are available. We also provide a simulation study in the paper to compare the new and conventional tests under different data scenarios. The results obtained are consistent with our asymptotic study. In addition we provide least squares estimators of those parameters used in the change‐point test together with their almost‐sure convergence properties.     

Efficient Regularization Parameter Selection Via Information Criteria

In a nonlinear regression model based on a regularization method, selection of appropriate regularization parameters is crucial. Information criteria such as generalized information criterion (GIC) and generalized Bayesian information criterion (GBIC) are useful for selecting the optimal regularization parameters. However, the optimal parameter is often determined by calculating information criterion for all candidate regularization parameters, and so the computational cost is high. One simple method by which to accomplish this is to regard GIC or GBIC as a function of the regularization parameters and to find a value minimizing GIC or GBIC. However, it is unclear how to solve the optimization problem. In the present article, we propose an efficient Newton–Raphson type iterative method for selecting optimal regularization parameters with respect to GIC or GBIC in a nonlinear regression model based on basis expansions. This method reduces the computational time remarkably compared to the grid search and can select more suitable regularization parameters. The effectiveness of the method is illustrated through real data examples.