Unified Inference for Sparse and Dense Longitudinal Data in Time‐varying Coefficient Models

Time‐varying coefficient models are widely used in longitudinal data analysis. These models allow the effects of predictors on response to vary over time. In this article, we consider a mixed‐effects time‐varying coefficient model to account for the within subject correlation for longitudinal data. We show that when kernel smoothing is used to estimate the smooth functions in time‐varying coefficient models for sparse or dense longitudinal data, the asymptotic results of these two situations are essentially different. Therefore, a subjective choice between the sparse and dense cases might lead to erroneous conclusions for statistical inference. In order to solve this problem, we establish a unified self‐normalized central limit theorem, based on which a unified inference is proposed without deciding whether the data are sparse or dense. The effectiveness of the proposed unified inference is demonstrated through a simulation study and an analysis of Baltimore MACS data.     

Estimation and Inference Procedures for Semiparametric Distribution Models with Varying Linear‐Index

More flexible semiparametric linear‐index regression models are proposed to describe the conditional distribution. Such a model formulation captures varying effects of covariates over the support of a response variable in distribution, offers an alternative perspective on dimension reduction and covers a lot of widely used parametric and semiparameteric regression models. A feasible pseudo likelihood approach, accompanied with a simple and easily implemented algorithm, is further developed for the mixed case with both varying and invariant coefficients. By justifying some theoretical properties on Banach spaces, the uniform consistency and asymptotic Gaussian process of the proposed estimator are also established in this article. In addition, under the monotonicity of distribution in linear‐index, we develop an alternative approach based on maximizing a varying accuracy measure. By virtue of the asymptotic recursion relation for the estimators, some of the achievements in this direction include showing the convergence of the iterative computation procedure and establishing the large sample properties of the resulting estimator. It is noticeable that our theoretical framework is very helpful in constructing confidence bands for the parameters of interest and tests for the hypotheses of various qualitative structures in distribution. Generally, the developed estimation and inference procedures perform quite satisfactorily in the conducted simulations and are demonstrated to be useful in reanalysing data from the Boston house price study and the World Values Survey.     

Robust Poisson likelihood estimation for frailty Cox models: A simulation study

Frailty models can be fit as mixed-effects Poisson models after transforming time-to-event data to the Poisson model framework. We assess, through simulations, the robustness of Poisson likelihood estimation for Cox proportional hazards models with log-normal frailties under misspecified frailty distribution. The log-gamma and Laplace distributions were used as true distributions for frailties on a natural log scale. Factors such as the magnitude of heterogeneity, censoring rate, number and sizes of groups were explored. In the simulations, the Poisson modeling approach that assumes log-normally distributed frailties provided accurate estimates of within- and between-group fixed effects even under a misspecified frailty distribution. Non-robust estimation of variance components was observed in the situations of substantial heterogeneity, large event rates, or high data dimensions.     

On the complex dynamics of functional-coefficients nonlinear autoregressive time series models

In this article, a new chaotic functional-coefficient nonlinear autoregressive time series model is formulated. Asymptotic stability of the equilibria of the skeleton is studied. Complex dynamics of the proposed model are investigated by means of the bifurcations, time series diagrams, and phase portraits. The effects of noise intensity on its dynamics and the intermittency phenomenon are also discussed via simulation. Two chaotic indicators, namely, the fractal dimension and the Lyapunov exponent methods are investigated for the model.     

Model comparison of linear and nonlinear Bayesian structural equation models with dichotomous data

In this article, dichotomous variables are used to compare between linear and nonlinear Bayesian structural equation models. Gibbs sampling method is applied for estimation and model comparison. Statistical inferences that involve estimation of parameters and their standard deviations and residuals analysis for testing the selected model are discussed. Hidden continuous normal distribution (censored normal distribution) is used to solve the problem of dichotomous variables. The proposed procedure is illustrated by a simulation data obtained from R program. Analyses are done by using R2WinBUGS package in R-program.     

Variable selection in partially linear additive models for modal regression

Based on B-spline basis functions and smoothly clipped absolute deviation (SCAD) penalty, we present a new estimation and variable selection procedure based on modal regression for partially linear additive models. The outstanding merit of the new method is that it is robust against outliers or heavy-tail error distributions and performs no worse than the least-square-based estimation for normal error case. The main difference is that the standard quadratic loss is replaced by a kernel function depending on a bandwidth that can be automatically selected based on the observed data. With appropriate selection of the regularization parameters, the new method possesses the consistency in variable selection and oracle property in estimation. Finally, both simulation study and real data analysis are performed to examine the performance of our approach.     

Optimum designs for parameter estimation in a mixture experiment with two correlated responses

In this paper, we investigate a mixture problem with two responses, which are functions of the mixing proportions, and are correlated with known dispersion matrix. We obtain D- and A-optimal designs for estimating the parameters of the response functions, when none or some of the regression coefficients of the two functions are the same. It is shown that when no prior knowledge about the regression coefficients is available, the D-optimal design is independent of the dispersion matrix, while the A-optimal design depends on it, provided the response functions are of different degree. On the other hand, when some of the regression coefficients are known to be the same for both the functions, the D-optimal design depends on the dispersion matrix when the two response functions are not of the same degree.     

An evaluation of the reconstructed coefficient of determination and potential adjustments

Previously, a method was proposed for calculating a reconstructed coefficient of determination in the case of right-censored regression using the expectation–maximization (EM) algorithm. This measure is assessed via simulation study for the purpose of evaluating the utility of model fit. Further, several reconstructed adjusted coefficients of determination are proposed and compared via simulation study for the purpose of model selection. The application of these proposed measures is illustrated on a real dataset.     

Bootstrap forecast intervals for asymmetric volatilities via EGARCH model

Bootstrap forecast intervals are developed for volatilities having asymmetric features, which are accounted for by fitting EGARCH models. A Monte-Carlo simulation compares the proposed forecast intervals with those based on GARCH fittings which ignore asymmetry. The comparison reveals substantial advantage of addressing asymmetry through EGARCH fitting over ignoring it as the conventional GARCH forecast. The EGARCH forecast intervals have empirical coverage probabilities closer to the nominal level and/or have shorter average lengths than the GARCH forecast intervals. The finding is also supported by real dataset analysis of Dow–Jones index and financial times stock exchange (FTSE) 100 index.     

On the estimation of a sensitive quantitative mean using blank cards

In this article, we introduce for the first time, the blank card methods for estimation of finite population mean of a sensitive variable. Two generic randomization devices are suggested, and for each device we identify the choices of special models. We introduce additive, multiplicative, and combination of both additive and multiplicative scrambling models that require use of a non sensitive variable. We derive the basic statistical properties of each model. It is interesting to note that various existing estimators can be viewed as the special cases of those presented here. The statistical efficiency of proposed techniques is compared with Greenberg et al. (1971 Greenberg, B.G., Kuebler Jr, R.R., Abernathy, J.R., Horvitz, D.G. (1971). Application of the randomized response technique in obtaining quantitative data. J. Am. Stat. Assoc. 66(334):243–250.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) and modified Perri (2008 Perri, P.F. (2008). Modified randomized devices for Simmons’ model. Model Assisted Stat. Appl. 3(3):233–239. [Google Scholar]) model. The proposed devices can easily be adjusted to achieve the required efficiency level by making suitable choices of different design parameters.