Fréchet differentiability in statistical inference for time series

It is shown how the method of Fréchet differentiability can simplify the asymptotic derivations in an important range of robust inferential problems for stationary and related time series models. The uniform root-n consistency of the empirical distribution function for the Cramer von Mises norm under a weak mixing condition is indicated. Various regularity conditions naturally implemented and leading to the differentiability are discussed. A simulation study supplementing the theoretical discussion is included.     

Pseudo-likelihood for case–cohort studies under length-biased sampling

A case–cohort design was proposed by Prentice (1986) Prentice, R.L. (1986). A case-cohort design for epidemiologic cohort studies and disease prevention trials. Biometrika 73:1–11.[Crossref], [Web of Science ®] , [Google Scholar] in order to reduce costs. It involves the collection of covariate data from all subjects who experience the event of interest, and from the members of a random subcohort. This case–cohort design has been extensively studied, but is exclusively considered for right-censored data. In this article, we propose case–cohort designs adapted to length-biased data under the proportional hazards assumption. A pseudo-likelihood procedure is described for estimating parameters and the corresponding cumulative hazard function. The large sample properties, such as consistency and weak convergence, for such pseudo-likelihood estimators are presented. We also conduct simulation studies to show that the proposed estimators are appropriate for practical use. A real Oscar Awards data is provided.     

Hidden Markov model for discrete circular–linear wind data time series

ABSTRACT

In this work, we deal with a bivariate time series of wind speed and direction. Our observed data have peculiar features, such as informative missing values, non-reliable measures under a specific condition and interval-censored data, that we take into account in the model specification. We analyse the time series with a non-parametric Bayesian hidden Markov model, introducing a new emission distribution, suitable to model our data, based on the invariant wrapped Poisson, the Poisson and the hurdle density. The model is estimated on simulated datasets and on the real data example that motivated this work.     

On testes for threshold–type nonlinearity in irregulaly spaced time series

Formulating the model first in continuous time, we have developed a state space approach to the problem of testing for threshold-type nonlinearity when the data are irregularly spaced.     

Fitting semiparametric accelerated failure time models for nested case–control data

A nested case–control (NCC) study is an efficient cohort-sampling design in which a subset of controls are sampled from the risk set at each event time. Since covariate measurements are taken only for the sampled subjects, time and efforts of conducting a full scale cohort study can be saved. In this paper, we consider fitting a semiparametric accelerated failure time model to failure time data from a NCC study. We propose to employ an efficient induced smoothing procedure for rank-based estimating method for regression parameters estimation. For variance estimation, we propose to use an efficient resampling method that utilizes the robust sandwich form. We extend our proposed methods to a generalized NCC study that allows a sampling of cases. Finite sample properties of the proposed estimators are investigated via an extensive stimulation study. An application to a tumor study illustrates the utility of the proposed method in routine data analysis.     

Constant-partially accelerated life tests for Marshall–Olkin exponential series system with dependent masked data

This paper considers the constant-partially accelerated life tests for series system products, where dependent M-O bivariate exponential distribution is assumed for the components.

Based on progressive type-II censored and masked data, the maximum likelihood estimates for the parameters and acceleration factors are obtained by using the decomposition approach. In addition, this method can also be applied to the Bayes estimates, which are too complex to obtain as usual way. Finally, a Monte Carlo simulation study is carried out to verify the accuracy of the methods under different masking probabilities and censoring schemes.     

Inference using attributable risk for a 2×2 case–control study

Asymptotic variance plays an important role in the inference using interval estimate of attributable risk. This paper compares asymptotic variances of attributable risk estimate using the delta method and the Fisher information matrix for a 2×2 case–control study due to the practicality of applications. The expressions of these two asymptotic variance estimates are shown to be equivalent. Because asymptotic variance usually underestimates the standard error, the bootstrap standard error has also been utilized in constructing the interval estimates of attributable risk and compared with those using asymptotic estimates. A simulation study shows that the bootstrap interval estimate performs well in terms of coverage probability and confidence length. An exact test procedure for testing independence between the risk factor and the disease outcome using attributable risk is proposed and is justified for the use with real-life examples for a small-sample situation where inference using asymptotic variance may not be valid.     

Functional principal component analysis for the explorative analysis of multisite–multivariate air pollution time series with long gaps

The knowledge of the urban air quality represents the first step to face air pollution issues. For the last decades many cities can rely on a network of monitoring stations recording concentration values for the main pollutants. This paper focuses on functional principal component analysis (FPCA) to investigate multiple pollutant datasets measured over time at multiple sites within a given urban area. Our purpose is to extend what has been proposed in the literature to data that are multisite and multivariate at the same time. The approach results to be effective to highlight some relevant statistical features of the time series, giving the opportunity to identify significant pollutants and to know the evolution of their variability along time. The paper also deals with missing value issue. As it is known, very long gap sequences can often occur in air quality datasets, due to long time failures not easily solvable or to data coming from a mobile monitoring station. In the considered dataset, large and continuous gaps are imputed by empirical orthogonal function procedure, after denoising raw data by functional data analysis and before performing FPCA, in order to further improve the reconstruction.     

What matters most? Different stakeholder perspectives on estimands for an invented case study in COPD

In drug development, we ask ourselves which population, endpoint and treatment comparison should be investigated. In this context, we also debate what matters most to the different stakeholders that are involved in clinical drug development, for example, patients, physicians, regulators and payers. With the publication of draft ICH E9 addendum on estimands in 2017, we now have a common framework and language to discuss such questions in an informed and transparent way. This has led to the estimand discussion being a key element in study development, including design, analysis and interpretation of a treatment effect. At an invited session at the 2018 PSI annual conference, PSI hosted a role‐play debate where the aim of the session was to mimic a regulatory and payer scientific advice discussion for a COPD drug. Including role‐play views from an industry sponsor, a patient, a regulator and a payer. This paper presents the invented COPD case‐study design and considerations relating to appropriate estimands are discussed by each of the stakeholders from their differing viewpoints with the additional inclusion of a technical (academic) perspective. The rationale for each perspective on approaches for handling intercurrent events is presented, with a key emphasis on the application of while‐on‐treatment and treatment policy estimands in this context. It is increasingly recognised that the treatment effect estimated by the treatment policy approach may not always be of primary clinical interest and may not appropriately communicate to patients the efficacy they can expect if they take the treatment as directed.     

Nonlinear filters based on taylor series expansions ∗

The nonlinear filters based on Taylor series approximation are broadly used for computational simplicity, even though their filtering estimates are clearly biased. In this paper, first, we analyze what is approximated when we apply the expanded nonlinear functions to the standard linear recursive Kalman filter algorithm. Next, since the state variable α_t and α_t-t are approximated as a conditional normal distribution given information up to time t - 1 (i.e., I_t-1) in approximation of the Taylor series expansion, it might be appropriate to evaluate each expectation by generating normal random numbers of α_t and α_t-1 given I_t-1 and those of the error terms θ and η_t. Thus, we propose the Monte-Carlo simulation filter using normal random draws. Finally we perform two Monte-Carlo experiments, where we obtain the result that the Monte-Carlo simulation filter has a superior performance over the nonlinear filters such as the extended Kalman filter and the second-order nonlinear filter.     

Moving average models—time series in m-dimensions

Stochastic models for discrete time series in the time domain are well known but such models lack consideration of spatial dependency I We expand on their work by constructing spatially dependent moving average models. Definitions of order, stationarity, invertibility, autocorrelation function, and spectrum are made as natural extensions of those in zero dimensions and are implemented in the one and two-space dimensional models.     

Comment on “recent developments in bootstrapping time series”

The paper develops a general framework for identification, estimation, and hypothesis testing in cointegrated systems when the cointegrating coefficients are subject to (possibly) non-linear and cross-equation restrictions, obtained from economic theory or other relevant a priori information. It provides a proof of the consistency of the quasi maximum likelihood estimators (QMLE), establishes the relative rates of convergence of the QMLE of the short-run and the long-run parameters, and derives their asymptotic distributions; thus generalizing the results already available in the literature for the linear case. The paper also develops tests of the over-identifying (possibly) non-linear restrictions on the cointegrating vectors. The estimation and hypothesis testing procedures are applied to an Almost Ideal Demand System estimated on U.K. quarterly observations. Unlike many other studies of consumer demand this application does not treat relative prices and real per capita expenditures as exogenously given.     

Functional coefficient time series models with trending regressors

This paper studies a functional coe?cient time series model with trending regressors, where the coe?cients are unknown functions of time and random variables. We propose a local linear estimation method to estimate the unknown coe?cient functions, and establish the corresponding asymptotic theory under mild conditions. We also develop a test procedure to see if the functional coe?cients take particular parametric forms. For practical use, we further propose a Bayesian approach to select the bandwidths, and conduct several numerical experiments to examine the finite sample performance of our proposed local linear estimator and the test procedure. The results show that the local linear estimator works well and the proposed test has satisfactory size and power. In addition, our simulation studies show that the Bayesian bandwidth selection method performs better than the cross-validation method. Furthermore, we use the functional coe?cient model to study the relationship between consumption per capita and income per capita in United States, and it was shown that the functional coe?cient model with our proposed local linear estimator and Bayesian bandwidth selection method performs well in both in-sample fitting and out-of-sample forecasting.     

Approximate confidence intervals for the weighted kappa coefficient of a binary diagnostic test subject to a case–control design

Case–control design to assess the accuracy of a binary diagnostic test (BDT) is very frequent in clinical practice. This design consists of applying the diagnostic test to all of the individuals in a sample of those who have the disease and in another sample of those who do not have the disease. The sensitivity of the diagnostic test is estimated from the case sample and the specificity is estimated from the control sample. Another parameter which is used to assess the performance of a BDT is the weighted kappa coefficient. The weighted kappa coefficient depends on the sensitivity and specificity of the diagnostic test, on the disease prevalence and on the weighting index. In this article, confidence intervals are studied for the weighted kappa coefficient subject to a case–control design and a method is proposed to calculate the sample sizes to estimate this parameter. The results obtained were applied to a real example.     

Efficient von Mises–Fisher concentration parameter estimation using Taylor series

A simple method is outlined for constructing a Taylor series for the Maximum Likelihood Estimate of the von Mises–Fisher concentration parameter based around an initial heuristic estimate. While existing treatments require multiple computationally intensive calculations of a Bessel ratio, this method provides accurate results using only one such calculation. The accuracy of the method is tested extensively, and the reuse of the Taylor series for multiple calculations is explored.     

Testing homogeneity in semiparametric mixture case–control models

Parametric and semiparametric mixture models have been widely used in applications from many areas, and it is often of interest to test the homogeneity in these models. However, hypothesis testing is non standard due to the fact that several regularity conditions do not hold under the null hypothesis. We consider a semiparametric mixture case–control model, in the sense that the density ratio of two distributions is assumed to be of an exponential form, while the baseline density is unspecified. This model was first considered by Qin and Liang (2011 Qin, J., Liang, K.Y. (2011). Hypothesis testing in a mixture case–control model. Biometrics 67(1):182–198.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar], biometrics), and they proposed a modified score statistic for testing homogeneity. In this article, we consider alternative testing procedures based on supremum statistics, which could improve power against certain types of alternatives. We demonstrate the connection and comparison among the proposed and existing approaches. In addition, we provide a unified theoretical justification of the supremum test and other existing test statistics from an empirical likelihood perspective. The finite-sample performance of the supremum test statistics was evaluated in simulation studies.