Detections of changes in return by a wavelet smoother with conditional heteroscedastic volatility

In this paper, we propose two estimators, an integral estimator and a discretized estimator, for the wavelet coefficient of regression functions in nonparametric regression models with heteroscedastic variance. These estimators can be used to test the jumps of the regression function. The model allows for lagged-dependent variables and other mixing regressors. The asymptotic distributions of the statistics are established, and the asymptotic critical values are analytically obtained from the asymptotic distribution. We also use the test to determine consistent estimators for the locations of change points. The jump sizes and locations of change points can be consistently estimated using wavelet coefficients, and the convergency rates of these estimators are derived. We perform some Monte Carlo simulations to check the powers and sizes of the test statistics. Finally, we give practical examples in finance and economics to detect changes in stock returns and short-term interest rates using the empirical wavelet method.     

Efficient rank tests for semiparametric competing risk models

We consider a semiparametric competing risk model given by k independent survival times. The paper offers an asymptotic treatment of tests for the semiparametric null hypothesis of equality of the underlying risks. It turns out that modified rank tests are asymptotically efficient for certain semiparametric submodels, where the baseline hazard is a nuisance parameter. In addition, the asymptotic relative efficiency of the present tests is derived. A comparison of asymptotic power functions can then be used to classify various tests proposed earlier in the literature. For instance a chi-square type test is efficient for proportional hazards. Data driven tests of likelihood ratio type are proposed for cones of alternatives. We will consider certain stochastically increasing alternatives as a special example. The paper shows how the concept of local asymptotic normality of Le Cam works for hazard oriented models.     

Estimation and model adequacy checking for multivariate seasonal autoregressive time series models with periodically varying parameters

We introduce a class of multivariate seasonal time series models with periodically varying parameters, abbreviated by the acronym SPVAR. The model is suitable for multivariate data, and combines a periodic autoregressive structure and a multiplicative seasonal time series model. The stationarity conditions (in the periodic sense) and the theoretical autocovariance functions of SPVAR stochastic processes are derived. Estimation and checking stages are considered. The asymptotic normal distribution of the least squares estimators of the model parameters is established, and the asymptotic distributions of the residual autocovariance and autocorrelation matrices in the class of SPVAR time series models are obtained. In order to check model adequacy, portmanteau test statistics are considered and their asymptotic distributions are studied. A simulation study is briefly discussed to investigate the finite-sample properties of the proposed test statistics. The methodology is illustrated with a bivariate quarterly data set on travelers entering in to Canada.     

A product Pareto distribution

The Pareto distributions are becoming increasing prominent in several applied areas. In this note, a new Pareto distribution is introduced. It takes the form of the product of two Pareto probability density functions. Various structural properties of this distribution are derived, including its cumulative distribution function, moments, mean deviation about the mean, mean deviation about the median, entropy, asymptotic distribution of the extreme order statistics, method of moments estimates, maximum likelihood estimates and the Fisher information matrix. The calculations involve the use of several special functions.     

Frequency Domain Tests for Residual Serial Correlation in Cointegration Regressions

This paper introduces tests for residual serial correlation in cointegrating regressions. The tests are devised in the frequency domain by using the spectral measure estimates. The asymptotic distributions of the tests are derived and test consistency is established. The asymptotic distributions are obtained by using the assumptions and methods that are different from those used in Grenander and Rosenblatt (1957) and Durlauf (1991). Small-scale simulation results are reported to illustrate the finite sample performance of the tests under various distributional assumptions on the data generating process. The distributions considered are normal and t-distributions. The tests are shown to have stable size at sample sizes as large as 50 or 100. Additionally, it is shown that the tests are reasonably powerful against the ARMA residuals. An empirical application of the tests to investigate the ‘weak-form’ efficiency in the foreign exchange market is also reported.     

Marginal likelihood and unit roots

We develop new tests for the hypothesis of unit roots that are based on the marginal likelihood of the general linear model. The marginal likelihood allows the incorporation of invariance arguments in the likelihood function. It turns out that marginal likelihood tests for unit roots appear to be more powerful than other unit root tests. For some basic models power functions almost coincide with the power envelopes, even in small samples. General correlation structures can be incorporated, either by standard likelihood procedures or by adjustments of the test statistics on the basis of asymptotic distributions.     

Maximum entropy autoregressive conditional heteroskedasticity model

In many applications, it has been found that the autoregressive conditional heteroskedasticity (ARCH) model under the conditional normal or Student’s t distributions are not general enough to account for the excess kurtosis in the data. Moreover, asymmetry in the financial data is rarely modeled in a systematic way. In this paper, we suggest a general density function based on the maximum entropy (ME) approach that takes account of asymmetry, excess kurtosis and also of high peakedness. The ME principle is based on the efficient use of available information, and as is well known, many of the standard family of distributions can be derived from the ME approach. We demonstrate how we can extract information functional from the data in the form of moment functions. We also propose a test procedure for selecting appropriate moment functions. Our procedure is illustrated with an application to the NYSE stock returns. The empirical results reveal that the ME approach with a fewer moment functions leads to a model that captures the stylized facts quite effectively.     

Relative Efficiency of the Fuzzy p-Value Approach to Hypothesis Testing

In missing data problems, it is often the case that there is a natural test statistic for testing a statistical hypothesis had all the data been observed. A fuzzy  p -value approach to hypothesis testing has recently been proposed which is implemented by imputing the missing values in the "complete data" test statistic by values simulated from the conditional null distribution given the observed data. We argue that imputing data in this way will inevitably lead to loss in power. For the case of scalar parameter, we show that the asymptotic efficiency of the score test based on the imputed "complete data" relative to the score test based on the observed data is given by the ratio of the observed data information to the complete data information. Three examples involving probit regression, normal random effects model, and unidentified paired data are used for illustration. For testing linkage disequilibrium based on pooled genotype data, simulation results show that the imputed Neyman Pearson and Fisher exact tests are less powerful than a Wald-type test based on the observed data maximum likelihood estimator. In conclusion, we caution against the routine use of the fuzzy  p -value approach in latent variable or missing data problems and suggest some viable alternatives.     

Shrinkage estimation of the proportion in randomized response

This article considers the asymptotic estimation theory for the proportion in randomized response survey usinguncertain prior information (UPI) about the true proportion parameter which is assumed to be available on the basis of some sort of realistic conjecture. Three estimators, namely, the unrestricted estimator, the shrinkage restricted estimator and an estimator based on a preliminary test, are proposed. Their asymptotic mean squared errors are derived and compared. The relative dominance picture of the estimators is presented.     

Likelihood Ratio Test for and Against Nonlinear Inequality Constraints

In applied statistics a finite dimensional parameter involved in the distribution function of the observed random variable is very often constrained by a number of nonlinear inequalities. This paper is devoted to studying the likelihood ratio test for and against the hypothesis that the parameter is restricted by some nonlinear inequalities. The asymptotic null distributions of the likelihood ratio statistics are derived by using the limits of the related optimization problems. The author also shows how to compute critical values for the tests.     

Goodness-of-fit tests based on series estimators in nonparametric instrumental regression

This paper proposes several tests of restricted specification in nonparametric instrumental regression. Based on series estimators, test statistics are established that allow for tests of the general model against a parametric or nonparametric specification as well as a test of exogeneity of the vector of regressors. The tests’ asymptotic distributions under correct specification are derived and their consistency against any alternative model is shown. Under a sequence of local alternative hypotheses, the asymptotic distributions of the tests are derived. Moreover, uniform consistency is established over a class of alternatives whose distance to the null hypothesis shrinks appropriately as the sample size increases. A Monte Carlo study examines finite sample performance of the test statistics.     

An exact test for a column of the covariance matrix based on a single observation

In this paper, we derive an exact test for a column of the covariance matrix. The test statistic is calculated by using a single observation. The exact distributions of the test statistic are derived under both the null and alternative hypotheses. We also obtain an analytical expression of the power function of the test for the equality of a column of the covariance matrix to a given vector. It is shown that the information contained in a single vector is large enough to ensure a good performance of the test. Moreover, the suggested test can be applied for time-dependent multivariate Gaussian processes.     

Asymptotic efficiency in estimation with conditional moment restrictions

In this paper, bounds on asymptotic efficiency are derived for a class of non-parametric models. The data are independent and identically distributed according to some unknown distribution F. There is a given function of the data and a parameter. The restrictions are that a conditional expectation of this function is zero at some point in the parameter space; this point is to be estimated. If F is assumed to be a multinomial distribution with known (finite) support, then the problem becomes parametric and the bound can be obtained from the information matrix. This bound turns out to depend only upon certain conditional moments, and not upon the support of the distribution. Since a general F can be approximated by a multinomial distribution, the multinomial bound applies to the general case.     

A parametric test of the negativity of the substitution matrix

The negativity of the substitution matrix implies that its latent roots are non-positive. When inequality restrictions are tested, standard test statistics such as a likelihood ratio or a Wald test are not X²-distributed in large samples. We propose a Wald test for testing the negativity of the substitution matrix. The asymptotic distribution of the statistic is a mixture of X²-distributions. The Wald test provides an exact critical value for a given significance level. The problems involved in computing the exact critical value can be avoided by using the upper and lower bound critical values derived by Kodde and Palm (1986). Finally the methods are applied to the empirical results obtained by Barten and Geyskens (1975).     

Maximum likelihood estimation of stochastic frontier models by the Fourier transform

The paper is concerned with several kinds of stochastic frontier models whose likelihood function is not available in closed form. First, with output-oriented stochastic frontier models whose one-sided errors have a distribution other than the standard ones (exponential or half-normal). The gamma and beta distributions are leading examples. Second, with input-oriented stochastic frontier models which are common in theoretical discussions but not in econometric applications. Third, with two-tiered stochastic frontier models when the one-sided error components follow gamma distributions. Fourth, with latent class models with gamma distributed one-sided error terms. Fifth, with models whose two-sided error component is distributed as stable Paretian and the one-sided error is gamma. The principal aim is to propose approximations to the density of the composed error based on the inversion of the characteristic function (which turns out to be manageable) using the Fourier transform. Procedures that are based on the asymptotic normal form of the log-likelihood function and have arbitrary degrees of asymptotic efficiency are also proposed, implemented and evaluated in connection with output-oriented stochastic frontiers. The new methods are illustrated using data for US commercial banks, electric utilities, and a sample from the National Youth Longitudinal Survey.     

A more flexible joint latent model for longitudinal and survival time data

In this article, a more flexible and easily explained joint latent model with time-varying coefficients is used to characterize time-dependent responses and a failure time. The dependence within time-dependent responses and between time-dependent responses and a failure time, and the heterogeneity in both processes are established through partially non-parametric latent variables. Based on longitudinal and survival time data, an estimation procedure is proposed for the parameter functions of the joint latent model. In our estimation, the approximated likelihood is constructed via substituting the basis function expansions for parameter functions. The expectation and maximization (EM) algorithm is then implemented to obtain the maximizer of the approximated likelihood function, and, hence, the estimated parameter functions. The validity of the considered joint latent model enables us to derive the asymptotic properties of the estimated functions. Moreover, the corresponding finite sample properties and the usefulness of our methods are demonstrated through a Monte Carlo simulation and the AIDS Clinical Trials Group (ACTG) 175 data. A possible extension of our joint latent model and some additional topics of interest are also discussed herein.     

Locally Most Powerful and Other Rank Tests for Independence – with a Contaminated Weighted Alternative

Several authors in the literature have attempted the quantification of the concept of stochastic dependence for bivariate distribution. Two weighted rank tests for testing independence against a weighted contamination alternative is proposed and their distributional properties are studied. We also derived a locally most powerful rank test for the alternative setting. The rank tests proposed are shown to be asymptotic locally most powerful for specific distributions.     

Efficiency of iterated WLS in the linear model with completely unknown heteroskedasticity

Iterated weighted least squares (IWLS) is investigated for estimating the regression coefficients in a linear model with symmetrically distributed errors. The variances of the errors are not specified; it is not assumed that they are unknown functions of the explanatory variables nor that they are given in some parametric way.
IWLS is carried out in a random number of steps, of which the first one is OLS. In each step the error variance at time t is estimated with a weighted sum of m squared residuals in the neighbourhood of t and the coefficients are estimated using WLS. Furthermore an estimate of the co-variance matrix is obtained. If this estimate is minimal in some way the iteration process is stopped.
Asymptotic properties of IWLS are derived for increasing sample size n . Some particular cases show that the asymptotic efficiency can be increased by allowing more than two steps. Even asymptotic efficiency with respect to WLS with the true error variances can be obtained if m is not fixed but tends to infinity with n and if the heteroskedasticity is smooth.     

Vergleich zweier Testverfahren zur Symmetrieprüfung in 4-Felder-Tafeln

In order to jugde the success of an election campaign, two opinion surveys are to be carried out at two different times, each timen people being asked about their opinion with respect to the partyx. Now the question arises whether we should ask the same sample ofn people both times (test of McNemar for the comparison of dependent frequencies) or whether it would be more suitable to carry out the second survey independently of the first one (test for the comparison of independent frequencies). In the present paper we calculate the asymptotic power functions of these two test procedures and derive the asymptotic relative efficiency (ARE).     

Inverse-probability-weighted logrank test for stratified survival data with missing measurements

The stratified logrank test can be used to compare survival distributions of several groups of patients, while adjusting for the effect of some discrete variable that may be predictive of the survival outcome. In practice, it can happen that this discrete variable is missing for some patients. An inverse-probability-weighted version of the stratified logrank statistic is introduced to tackle this issue. Its asymptotic distribution is derived under the null hypothesis of equality of the survival distributions. A simulation study is conducted to assess behavior of the proposed test statistic in finite samples. An analysis of a medical dataset illustrates the methodology.