Biomedical applications for a generalized linear functional Poisson model

Many medical and biological studies involve response in the form of Poisson counts which can bemodelled using explanatory variables which also arise from count data. If the explanatory variables are observable without error (also as Poisson counts) we have a generalized linear model with a logarithmic link function and Poisson error structure. If,however, some of the explanatory variables are not directly observable, but arise with superimposed errors (again of Poisson form), the model is of a new type:a generalised linear functional Poisson model. In this paper,maximum likelihood estimates of the parameters of this model are determined along with the information matrix which (on noting its particular patterned form) is amenable to inversion in explicit form. Methods are proposed of an iterative type for computing estimates of the parameters and of their variational properties (e.g. standard errors) for this model, which also has application in other fields such as road traffic studies.     

Modelling the travel time of transit vehicles in real‐time through a GTFS‐based road network using GPS vehicle locations

Predicting the arrival time of a transit vehicle involves not only knowledge of its current position and schedule adherence, but also traffic conditions along the remainder of the route. Road networks are dynamic and can quickly change from free‐flowing to highly congested, which impacts the arrival time of transit vehicles, particularly buses which often share the road with other vehicles, so reliable predictions need to account for real‐time and future traffic conditions. The first step in this process is to construct a framework with which road state (traffic conditions) can be estimated using real‐time transit vehicle position data. Our proposed framework implements a vehicle model using a particle filter to estimate road travel times, which are used in a second model to estimate real‐time traffic conditions. Although development and testing took place in Auckland, New Zealand, we generalised each component to make the framework compatible with other public transport systems around the world. We demonstrate the real‐time feasibility and performance of our approach in real‐time, where a combination of R and C++ was used to obtain the necessary performance results. Future work will use these estimated traffic conditions in combination with historical data to obtain reliable arrival time predictions of transit vehicles.     

Data gold mining

The Motorway Incident Detection and Automatic Signalling (MIDAS) system data are summaries of traffic volumes, speeds and vehicle types, derived from loop detectors that sense traffic passing over them. These sensors are being rolled out across the UK's strategic highway network (motorways and major A-roads) and data from portions of one of the busiest motorways, the M25, now go back 10 years. Richard Gibbens and Wiebke Werft describe a journey time predictor based on MIDAS and explain how the data can shed light on topical issues such as road pricing.     

基于MixedLogit模型的道路交通统计生命价值评价研究

道路交通统计生命价值是道路交通安全项目经济评价的重要指标。运用意愿选择法和正交试验法设计了出行路径选择的调查问卷。假定时间、死亡风险为常数,费用服从对数正态分布,建立了基于Mixed Logit模型的统计生命价值评价模型。以大连市私家车出行者为调查对象获得调查数据,利用Monte Carlo方法并借助GAUSS软件对模型进行了150次仿真实验。研究表明：模型各参数的估计值具有较强集中性,t检验值具有较强显著性,模型优度比较高。统计生命价值服从参数为0．922和0．814的对数正态分布,其数学期望是35万元,可以作为道路交通安全项目经济评价的参考数据。     

Mixed inverse problems arising in the estimation of calibration factors in positron emission tomography

A penalized likelihood approach to the estimation of calibration factors in positron emission tomography (PET) is considered, in particular the problem of estimating the efficiency of PET detectors. Varying efficiencies among the detectors create a non-uniform performance and failure to account for the non-uniformities would lead to streaks in the image, so efficient estimation of the non-uniformities is desirable to reduce the propagation of noise to the final image. The relevant data set is provided by a blank scan, where a model may be derived that depends only on the sources affecting non-uniformities: inherent variation among the detector crystals and geometric effects. Physical considerations suggest a novel mixed inverse model with random crystal effects and smooth geometric effects. Using appropriate penalty terms, the penalized maximum likelihood estimates are derived and an efficient computational algorithm utilizing the fast Fourier transform is developed. Data-driven shrinkage and smoothing parameters are chosen to minimize an estimate of the predictive loss function. Various examples indicate that the approach proposed works well computationally and compares well with the standard method.     

基于车牌识别数据的城市交通流时空特性分析

了解城市道路交通流的时空特性有助于提高交通流的预测精度。基于实测的车牌识别数据,运用相似系数、快速傅氏变换和混沌理论分析交通流相似性、周期性和混沌性的时间特性;运用相关系数及互相关函数作为度量标准对交通流空间特性进行描述,论证交通流检测断面之间存在空间相互作用并具有时滞性;根据交通流的空间相关性,使用多维标度法对实际路网中检测断面进行聚类和分组,为交通流预测中的多断面分组以及构建交通流时空预测模型提供理论基础。     

Testing Equality of Mean Vectors in Two Sample Problem with Missing Data

Finding the influence of traffic accident on the road is helpful to analyze the characteristics of traffic flow, and take reasonable and effective control measures. Here, the detrended fluctuation analysis method is applied to investigate the complexity of time series in mixed traffic flow with a blockage induced by an accident. As a parameter to depict the long-term evolutionary behavior of the time series in traffic flow, the scaling exponent is analyzed. According to the scaling exponent, it is shown that the traffic flow time series can display long-range correlation characteristics, short-range correlation characteristics, and non-power-law relation in the long-range correlation characteristics, which is strongly dependent on the entering probability of vehicle, the ratio of slow vehicle and the blockage duration time.     

ELICITING A DIRECTED ACYCLIC GRAPH FOR A MULTIVARIATE TIME SERIES OF VEHICLE COUNTS IN A TRAFFIC NETWORK

The problem of modelling multivariate time series of vehicle counts in traffic networks is considered. It is proposed to use a model called the linear multiregression dynamic model (LMDM). The LMDM is a multivariate Bayesian dynamic model which uses any conditional independence and causal structure across the time series to break down the complex multivariate model into simpler univariate dynamic linear models. The conditional independence and causal structure in the time series can be represented by a directed acyclic graph (DAG). The DAG not only gives a useful pictorial representation of the multivariate structure, but it is also used to build the LMDM. Therefore, eliciting a DAG which gives a realistic representation of the series is a crucial part of the modelling process. A DAG is elicited for the multivariate time series of hourly vehicle counts at the junction of three major roads in the UK. A flow diagram is introduced to give a pictorial representation of the possible vehicle routes through the network. It is shown how this flow diagram, together with a map of the network, can suggest a DAG for the time series suitable for use with an LMDM.     

On MLEs of the parameters of a modified Weibull distribution for progressively type-2 censored samples

Lifetimes of modern mechanic or electronic units usually exhibit bathtub-shaped failure rates. An appropriate probability distribution to model such data is the modified Weibull distribution proposed by Lai et al. [15]. This distribution has both the two-parameter Weibull and type-1 extreme value distribution as special cases. It is able to model lifetime data with monotonic and bathtub-shaped failure rates, and thus attracts some interest among researchers because of this property. In this paper, the procedure of obtaining the maximum likelihood estimates (MLEs) of the parameters for progressively type-2 censored and complete samples are studied. Existence and uniqueness of the MLEs are proved.     

A test for lack-of-fit of zero-inflated negative binomial models

When a count data set has excessive zero counts, nonzero counts are overdispersed, and the effect of a continuous covariate might be nonlinear, for analysis a semiparametric zero-inflated negative binomial (ZINB) regression model is proposed. The unspecified smooth functional form for the continuous covariate effect is approximated by a cubic spline. The semiparametric ZINB regression model is fitted by maximizing the likelihood function. The likelihood ratio procedure is used to evaluate the adequacy of a postulated parametric functional form for the continuous covariate effect. An extensive simulation study is conducted to assess the finite-sample performance of the proposed test. The practicality of the proposed methodology is demonstrated with data of a motorcycle survey of traffic regulations conducted in 2007 in Taiwan by the Ministry of Transportation and Communication.     

Nonparametric modeling of the gap time in recurrent event data

Recurrent event data arise in many biomedical and engineering studies when failure events can occur repeatedly over time for each study subject. In this article, we are interested in nonparametric estimation of the hazard function for gap time. A penalized likelihood model is proposed to estimate the hazard as a function of both gap time and covariate. Method for smoothing parameter selection is developed from subject-wise cross-validation. Confidence intervals for the hazard function are derived using the Bayes model of the penalized likelihood. An eigenvalue analysis establishes the asymptotic convergence rates of the relevant estimates. Empirical studies are performed to evaluate various aspects of the method. The proposed technique is demonstrated through an application to the well-known bladder tumor cancer data.     

Penalized likelihood regression: General formulation and efficient approximation

The authors consider a formulation of penalized likelihood regression that is sufficiently general to cover canonical and noncanonical links for exponential families as well as accelerated life models with censored survival data. They present an asymptotic analysis of convergence rates to justify a simple approach to the lower‐dimensional approximation of the estimates. Such an approximation allows for much faster numerical calculation, paving the way to the development of algorithms that scale well with large data sets.     

Estimation of exponential regression parameters using binary data

Exponential regression model is important in analyzing data from heterogeneous populations. In this paper we propose a simple method to estimate the regression parameters using binary data. Under certain design distributions, including ellipticaily symmetric distributions, for the explanatory variables, the estimators are shown to be consistent and asymptotically normal when sample size is large. For finite samples, the new estimates were shown to behave reasonably well. They are competitive with the maximum likelihood estimates and more importantly, according to our simulation results, the cost of CPU time for computing new estimates is only 1/7 of that required for computing the usual maximum likelihood estimates. We expect the savings in CPU time would be more dramatic with larger dimension of the regression parameter space.     

Missing data techniques for multilevel data: implications of model misspecification

When modeling multilevel data, it is important to accurately represent the interdependence of observations within clusters. Ignoring data clustering may result in parameter misestimation. However, it is not well established to what degree parameter estimates are affected by model misspecification when applying missing data techniques (MDTs) to incomplete multilevel data. We compare the performance of three MDTs with incomplete hierarchical data. We consider the impact of imputation model misspecification on the quality of parameter estimates by employing multiple imputation under assumptions of a normal model (MI/NM) with two-level cross-sectional data when values are missing at random on the dependent variable at rates of 10%, 30%, and 50%. Five criteria are used to compare estimates from MI/NM to estimates from MI assuming a linear mixed model (MI/LMM) and maximum likelihood estimation to the same incomplete data sets. With 10% missing data (MD), techniques performed similarly for fixed-effects estimates, but variance components were biased with MI/NM. Effects of model misspecification worsened at higher rates of MD, with the hierarchical structure of the data markedly underrepresented by biased variance component estimates. MI/LMM and maximum likelihood provided generally accurate and unbiased parameter estimates but performance was negatively affected by increased rates of MD.     

Point and Interval Estimations of Marginal Risk Difference by Logistic Model

We use logistic model to get point and interval estimates of the marginal risk difference in observational studies and randomized trials with dichotomous outcome. We prove that the maximum likelihood estimate of the marginal risk difference is unbiased for finite sample and highly robust to the effects of dispersing covariates. We use approximate normal distribution of the maximum likelihood estimates of the logistic model parameters to get approximate distribution of the maximum likelihood estimate of the marginal risk difference and then the interval estimate of the marginal risk difference. We illustrate application of the method by a real medical example.     

Hierarchical likelihood approach for the Weibull frailty model

In this article, the proportional hazard model with Weibull frailty, which is outside the range of the exponential family, is used for analysing the right-censored longitudinal survival data. Complex multidimensional integrals are avoided by using hierarchical likelihood to estimate the regression parameters and to predict the realizations of random effects. The adjusted profile hierarchical likelihood is adopted to estimate the parameters in frailty distribution, during which the first- and second-order methods are used. The simulation studies indicate that the regression-parameter estimates in the Weibull frailty model are accurate, which is similar to the gamma frailty and lognormal frailty models. Two published data sets are used for illustration.     

Estimation for u-shaped beta distributions: minimum hellinger distance and related methods

We compare minimum Hellinger distance and minimum Heiiinger disparity estimates for U-shaped beta distributions. Given suitable density estimates, both methods are known to be asymptotically efficient when the data come from the assumed model family, and robust to small perturbations from the model family. Most implementations use kernel density estimates, which may not be appropriate for U-shaped distributions. We compare fixed binwidth histograms, percentile mesh histograms, and averaged shifted histograms. Minimum disparity estimates are less sensitive to the choice of density estimate than are minimum distance estimates, and the percentile mesh histogram gives the best results for both minimum distance and minimum disparity estimates. Minimum distance estimates are biased and a bias-corrected method is proposed. Minimum disparity estimates and bias-corrected minimum distance estimates are comparable to maximum likelihood estimates when the model holds, and give better results than either method of moments or maximum likelihood when the data are discretized or contaminated, Although our re¬sults are for the beta density, the implementations are easily modified for other U-shaped distributions such as the Dirkhlet or normal generated distribution.     

A half circular distribution for modeling the posterior corneal curvature

The posterior corneal curvature and many other medical, environmental, and ecological variables are measured with angles where its range is less than π. Such data are so-called axial or half circular data. Half circular data modeling has not received much attention from researchers. This paper proposes a new half circular distribution model based on inverse stereographic projection technique of Burr-XII distribution. The maximum likelihood estimates of parameters are obtained and a simulation study to evaluate the performance of estimates was carried out. The application on the posterior corneal curvature of 23 patients shows that the proposed distribution fits the data well.     

Two Applications of Statistics to Traffic Models

Two statistical applications for estimation and prediction of flows in traffic networks are presented. In the first, the number of route users are assumed to be independent α-shifted gamma Γ(θ, λ₀) random variables denoted H(α, θ, λ₀), with common λ₀. As a consequence, the link, OD (origin-destination) and node flows are also H(α, θ, λ₀) variables. We assume that the main source of information is plate scanning, which permits us to identify, totally or partially, the vehicle route, OD and link flows by scanning their corresponding plate numbers at an adequately selected subset of links. A Bayesian approach using conjugate families is proposed that allows us to estimate different traffic flows. In the second application, a stochastic demand dynamic traffic model to predict some traffic variables and their time evolution in real networks is presented. The Bayesian network model considers that the variables are generalized Beta variables such that when marginally transformed to standard normal become multivariate normal. The model is able to provide a point estimate, a confidence interval or the density of the variable being predicted. Finally, the models are illustrated by their application to the Nguyen Dupuis network and the Vermont-State example. The resulting traffic predictions seem to be promising for real traffic networks and can be done in real time.     

Modelling a non-stationary BINAR(1) Poisson process

ABSTRACT

Non-stationarity in bivariate time series of counts may be induced by a number of time-varying covariates affecting the bivariate responses due to which the innovation terms of the individual series as well as the bivariate dependence structure becomes non-stationary. So far, in the existing models, the innovation terms of individual INAR(1) series and the dependence structure are assumed to be constant even though the individual time series are non-stationary. Under this assumption, the reliability of the regression and correlation estimates is questionable. Besides, the existing estimation methodologies such as the conditional maximum likelihood (CMLE) and the composite likelihood estimation are computationally intensive. To address these issues, this paper proposes a BINAR(1) model where the innovation series follow a bivariate Poisson distribution under some non-stationary distributional assumptions. The method of generalized quasi-likelihood (GQL) is used to estimate the regression effects while the serial and bivariate correlations are estimated using a robust moment estimation technique. The application of model and estimation method is made in the simulated data. The GQL method is also compared with the CMLE, generalized method of moments (GMM) and generalized estimating equation (GEE) approaches where through simulation studies, it is shown that GQL yields more efficient estimates than GMM and equally or slightly more efficient estimates than CMLE and GEE.