Estimating the fundamental frequency of a periodic function

In this paper we consider the problem of estimation of the fundamental frequency of a periodic function, which has several applications in Speech Signal Processing. The problem was originally proposed by Hannan (1974) and later on Quinn and Thomson (1991) provided an estimation procedure of the unknown parameters. It is observed that the estimation procedure of Quinn and Thomson (1991) is quite involved numerically. In this paper we propose to use two simple estimators and it is observed that their performance are quite satisfactory. Asymptotic properties of the proposed estimators are obtained. The large sample properties of the estimators are compared theoretically. We present some simulation results to compare their small sample performance. One speech data is analyzed using this particular model.     

On Markov-switching periodic ARMA models

In this work, we propose a generalization of the classical Markov-switching ARMA models to the periodic time-varying case. Specifically, we propose a Markov-switching periodic ARMA (MS-PARMA) model. In addition of capturing regime switching often encountered during the study of many economic time series, this new model also captures the periodicity feature in the autocorrelation structure. We first provide some probabilistic properties of this class of models, namely the strict periodic stationarity and the existence of higher-order moments. We thus propose a procedure for computing the autocovariance function where we show that the autocovariances of the MS-PARMA model satisfy a system of equations similar to the PARMA Yule–Walker equations. We propose also an easily implemented algorithm which can be used to obtain parameter estimates for the MS-PARMA model. Finally, a simulation study of the performance of the proposed estimation method is provided.     

Periodic integration: further results on model selection and forecasting

This paper considers model selection and forecasting issues in two closely related models for nonstationary periodic autoregressive time series [PAR]. Periodically integrated seasonal time series [PIAR] need a periodic differencing filter to remove the stochastic trend. On the other hand, when the nonperiodic first order differencing filter can be applied, one can have a periodic model with a nonseasonal unit root [PARI]. In this paper, we discuss and evaluate two testing strategies to select between these two models. Furthermore, we compare the relative forecasting performance of each model using Monte Carlo simulations and some U.K. macroeconomic seasonal time series. One result is that forecasting with PARI models while the data generating process is a PIAR process seems to be worse thanvice versa.     

An alternative circular smoothing method to nonparametric estimation of periodic functions

This article provides alternative circular smoothing methods in nonparametric estimation of periodic functions. By treating the data as ‘circular’, we solve the “boundary issue” in the nonparametric estimation treating the data as ‘linear’. By redefining the distance metric and signed distance, we modify many estimators used in the situations involving periodic patterns. In the perspective of ‘nonparametric estimation of periodic functions’, we present the examples in nonparametric estimation of (1) a periodic function, (2) multiple periodic functions, (3) an evolving function, (4) a periodically varying-coefficient model and (5) a generalized linear model with periodically varying coefficient. In the perspective of ‘circular statistics’, we provide alternative approaches to calculate the weighted average and evaluate the ‘linear/circular–linear/circular’ association and regression. Simulation studies and an empirical study of electricity price index have been conducted to illustrate and compare our methods with other methods in the literature.     

Estimating nonlinear additive models with nonstationarities and correlated errors

In this paper, we study a nonparametric additive regression model suitable for a wide range of time series applications. Our model includes a periodic component, a deterministic time trend, various component functions of stochastic explanatory variables, and an AR(p) error process that accounts for serial correlation in the regression error. We propose an estimation procedure for the nonparametric component functions and the parameters of the error process based on smooth backfitting and quasimaximum likelihood methods. Our theory establishes convergence rates and the asymptotic normality of our estimators. Moreover, we are able to derive an oracle‐type result for the estimators of the AR parameters: Under fairly mild conditions, the limiting distribution of our parameter estimators is the same as when the nonparametric component functions are known. Finally, we illustrate our estimation procedure by applying it to a sample of climate and ozone data collected on the Antarctic Peninsula.     

Asymptotic Inefficiency of Mean-Correction on Parameter Estimation for a Periodic First-Order Autoregressive Model

A common practice in time series analysis is to fit a centered model to the mean-corrected data set. For stationary autoregressive moving-average (ARMA) processes, as far as the parameter estimation is concerned, fitting an ARMA model without intercepts to the mean-corrected series is asymptotically equivalent to fitting an ARMA model with intercepts to the observed series. We show that, related to the parameter least squares estimation of periodic ARMA models, the second approach can be arbitrarily more efficient than the mean-corrected counterpart. This property is illustrated by means of a periodic first-order autoregressive model. The asymptotic variance of the estimators for both approaches is derived. Moreover, empirical experiments based on simulations investigate the finite sample properties of the estimators.     

SEMIPARAMETRIC COMPARISON OF REGRESSION CURVES VIA NORMAL LIKELIHOODS

Härdle & Marron (1990) treated the problem of semiparametric comparison of nonparametric regression curves by proposing a kernel-based estimator derived by minimizing a version of weighted integrated squared error. The resulting estimators of unknown transformation parameters are n-consistent, which prompts a consideration of issues. of optimality. We show that when the unknown mean function is periodic, an optimal nonparametric estimator may be motivated by an elegantly simple argument based on maximum likelihood estimation in a parametric model with normal errors. Strikingly, the asymptotic variance of an optimal estimator of θ does not depend at all on the manner of estimating error variances, provided they are estimated n-consistently. The optimal kernel-based estimator derived via these considerations is asymptotically equivalent to a periodic version of that suggested by Härdle & Marron, and so the latter technique is in fact optimal in this sense. We discuss the implications of these conclusions for the aperiodic case.     

Student-t censored regression model: properties and inference

In statistical analysis, particularly in econometrics, it is usual to consider regression models where the dependent variable is censored (limited). In particular, a censoring scheme to the left of zero is considered here. In this article, an extension of the classical normal censored model is developed by considering independent disturbances with identical Student-t distribution. In the context of maximum likelihood estimation, an expression for the expected information matrix is provided, and an efficient EM-type algorithm for the estimation of the model parameters is developed. In order to know what type of variables affect the income of housewives, the results and methods are applied to a real data set. A brief review on the normal censored regression model or Tobit model is also presented.     

A note on Läuter-type robust design for polynomial regression models

In this paper, we consider the problem of model robust design for simultaneous parameter estimation among a class of polynomial regression models with degree up to k. A generalized D-optimality criterion, the Ψ_α‐optimality criterion, first introduced by Läuter (1974) is considered for this problem. By applying the theory of canonical moments and the technique of maximin principle, we derive a model robust optimal design in the sense of having highest minimum Ψ_α‐efficiency. Numerical comparison indicates that the proposed design has remarkable performance for parameter estimation in all of the considered rival models.     

On periodic autoregressive stochastic volatility models: structure and estimation

To capture both the volatility evolution and the periodicity feature in the autocorrelation structure exhibited by many nonlinear time series, a Periodic AutoRegressive Stochastic Volatility (PAR-SV ) model is proposed. Some probabilistic properties, namely the strict and second-order periodic stationarity, are provided. Furthermore, conditions for the existence of higher-order moments are established. The autocovariance structure of the squares and higher order powers of the PAR-SV process is studied. Its dynamic properties are shown to be consistent with financial time series empirical findings. Ways in which the model may be estimated are discussed. Finally, a simulation study of the performance of the proposed estimation methods is provided and the PAR-SV is applied to model the spot rates of the euro and US dollar both against the Algerian dinar. The empirical analysis shows that the proposed PAR-SV model can be considered as a viable alternative to the periodic generalized autoregressive conditionally heteroscedastic (PGARCH) model.     

A tutorial on rank-based coefficient estimation for censored data in small- and large-scale problems

The analysis of survival endpoints subject to right-censoring is an important research area in statistics, particularly among econometricians and biostatisticians. The two most popular semiparametric models are the proportional hazards model and the accelerated failure time (AFT) model. Rank-based estimation in the AFT model is computationally challenging due to optimization of a non-smooth loss function. Previous work has shown that rank-based estimators may be written as solutions to linear programming (LP) problems. However, the size of the LP problem is O(n ²+p) subject to n ² linear constraints, where n denotes sample size and p denotes the dimension of parameters. As n and/or p increases, the feasibility of such solution in practice becomes questionable. Among data mining and statistical learning enthusiasts, there is interest in extending ordinary regression coefficient estimators for low-dimensions into high-dimensional data mining tools through regularization. Applying this recipe to rank-based coefficient estimators leads to formidable optimization problems which may be avoided through smooth approximations to non-smooth functions. We review smooth approximations and quasi-Newton methods for rank-based estimation in AFT models. The computational cost of our method is substantially smaller than the corresponding LP problem and can be applied to small- or large-scale problems similarly. The algorithm described here allows one to couple rank-based estimation for censored data with virtually any regularization and is exemplified through four case studies.     

Evaluating hospital readmission rates in dialysis facilities; adjusting for hospital effects

Motivated by the national evaluation of readmission rates among kidney dialysis facilities in the United States, we evaluate the impact of including discharging hospitals on the estimation of facility-level standardized readmission ratios (SRRs). The estimation of SRRs consists of two steps. First, we model the dependence of readmission events on facilities and patient-level characteristics, with or without an adjustment for discharging hospitals. Second, using results from the models, standardization is achieved by computing the ratio of the number of observed events to the number of expected events assuming a population norm and given the case-mix in that facility. A challenging aspect of our motivating example is that the number of parameters is very large and estimation of high-dimensional parameters is troublesome. To solve this problem, we propose a structured Newton-Raphson algorithm for a logistic fixed effects model and an approximate EM algorithm for the logistic mixed effects model. We consider a re-sampling and simulation technique to obtain p-values for the proposed measures. Finally, our method of identifying outlier facilities involves converting the observed p-values to Z-statistics and using the empirical null distribution, which accounts for overdispersion in the data. The finite-sample properties of proposed measures are examined through simulation studies. The methods developed are applied to national dialysis data. It is our great pleasure to present this paper in honor of Ross Prentice, who has been instrumental in the development of modern methods of modeling and analyzing life history and failure time data, and in the inventive applications of these methods to important national data problem.     

ESTIMATING PREDATOR-PREY SYSTEMS VIA ORDINARY DIFFERENTIAL EQUATIONS WITH CLOSED ORBITS

This paper considers periodic regression functions, which are solutions to a planar system of differential equations. In particular, it introduces a simple stochastic model which describes the interaction between predator and prey populations. The regression functions are solutions to the classical Lotka‐Volterra system of equations, which admits closed orbits. The proposed method of estimation can be applied whenever pairs of predator‐prey data are available, and the prey is the main source of food of the predator. Canadian mink‐muskrat data are analysed from this new viewpoint. The estimation method is based on the existence of closed trajectories that describe the relationship between the two population sizes, and the paper shows how it can be extended to other systems of differential equations which admit closed orbits (e.g. Hamiltonian systems).     

A note on estimating parameters of partial prior information

In an empirical Bayes model, examples for estimators of parameters of partial prior information are given. A typical application is the estimation of a probability ?(p≤a),p the “fraction defective”, which is used as prior information in Quality Control.     

Occupancy numbers for dynamic heterogeneous populations: Estimate of particles lifetimes

A dynamic model of a heterogeneous population is studied. Particles belonging to a population are divided, at every time t, into a finite number of classes according to their types and the partition changes over time. The role of the occupancy numbers, namely the cardinality of each class, is highlighted. The relationship between the stochastic process of occupancy numbers and the process of particle types is analyzed. The main goal of this paper is the estimation of the lifetime of each particle at a given time t, when the observed data are the history of the process of the number of dead particles up to t. Furthermore, a discrete time approximation of the filter is given.     

On the discrete analogues of continuous distributions

In this paper, a new method is proposed for generating discrete distributions. A special class of the distributions, namely, the T-geometric family contains the discrete analogues of continuous distributions. Some general properties of the T-geometric family of distributions are obtained. A member of the T-geometric family, namely, the exponentiated-exponential–geometric distribution is defined and studied. Various properties of the exponentiated-exponential–geometric distribution such as the unimodality, the moments and the probability generating function are discussed. The method of maximum likelihood estimation is proposed for estimating the model parameters. Three real data sets are used to illustrate the applications of the exponentiated-exponential–geometric distribution.     

A new alternative to the standard F test for clustered data

The data collection process and the inherent population structure are the main causes for clustered data. The observations in a given cluster are correlated, and the magnitude of such correlation is often measured by the intra-cluster correlation coefficient. The intra-cluster correlation can lead to an inflated size of the standard F test in a linear model. In this paper, we propose a solution to this problem. Unlike previous adjustments, our method does not require estimation of the intra-class correlation, which is problematic especially when the number of clusters is small. Our simulation results show that the new method outperforms the existing methods.     

On the discrete analogues of continuous distributions

In this paper, a new method is proposed for generating discrete distributions. A special class of the distributions, namely, the T-geometric family contains the discrete analogues of continuous distributions. Some general properties of the T-geometric family of distributions are obtained. A member of the T-geometric family, namely, the exponentiated-exponential–geometric distribution is defined and studied. Various properties of the exponentiated-exponential–geometric distribution such as the unimodality, the moments and the probability generating function are discussed. The method of maximum likelihood estimation is proposed for estimating the model parameters. Three real data sets are used to illustrate the applications of the exponentiated-exponential–geometric distribution.     

Rank-order choice-based conjoint experiments: Efficiency and design

In a rank-order choice-based conjoint experiment, the respondent is asked to rank a number of alternatives of a number of choice sets. In this paper, we study the efficiency of those experiments and propose a D-optimality criterion for rank-order experiments to find designs yielding the most precise parameter estimators. For that purpose, an expression of the Fisher information matrix for the rank-ordered conditional logit model is derived which clearly shows how much additional information is provided by each extra ranking step. A simulation study shows that, besides the Bayesian D-optimal ranking design, the Bayesian D-optimal choice design is also an appropriate design for this type of experiments. Finally, it is shown that considerable improvements in estimation and prediction accuracy are obtained by including extra ranking steps in an experiment.