The double chain markov model

In the class of discrete time Markovian processes, two models are widely used, the Markov chain and the hidden Markov model. A major difference between these two models lies in the relation between successive outputs of the observed variable. In a visible Markov chain, these are directly correlated while in hidden models they are not. However, in some situations it is possible to observe both a hidden Markov chain and a direct relation between successive observed outputs. Unfortunately, the use of either a visible or a hidden model implies the suppression of one of these hypothesis. This paper prsents a Markovian model under random environment called the Double Chain Markov Model which takes into account the maijn features of both visible and hidden models. Its main purpose is the modeling of non-homogeneous time-series. It is very flexible and can be estimated with traditional methods. The model is applied on a sequence of wind speeds and it appears to model data more successfully than both the usual Markov chains and hidden Markov models.     

A bayesian comparison of randomized and voluntary response sampling models

This paper is concerned with making a Bayesian comparison between the two types of survey sampling procedures where the individuals are asked sensitive questions. The randomized response procedure is designed to obtain information in such a way that the respondent is protected against revealing his/her status about the sensitive issue. In voluntary response model each individual chooses either to respond or not to respond. It is assumed that the respondents do not falsify their answers. The problem of interest is to estimate the proportion of the sampled population who belong to the sensitive group. A beta prior distribution is used for the proportion of interest. The loss function is taken to be the quadratic loss function.     

Dynamic ordering design for dose finding in drug-combination trials

Drug-combination studies have become increasingly popular in oncology. One of the critical concerns in phase I drug-combination trials is the uncertainty in toxicity evaluation. Most of the existing phase I designs aim to identify the maximum tolerated dose (MTD) by reducing the two-dimensional searching space to one dimension via a prespecified model or splitting the two-dimensional space into multiple one-dimensional subspaces based on the partially known toxicity order. Nevertheless, both strategies often lead to complicated trials which may either be sensitive to model assumptions or induce longer trial durations due to subtrial split. We develop two versions of dynamic ordering design (DOD) for dose finding in drug-combination trials, where the dose-finding problem is cast in the Bayesian model selection framework. The toxicity order of dose combinations is continuously updated via a two-dimensional pool-adjacent-violators algorithm, and then the dose assignment for each incoming cohort is selected based on the optimal model under the dynamic toxicity order. We conduct extensive simulation studies to evaluate the performance of DOD in comparison with four other commonly used designs under various scenarios. Simulation results show that the two versions of DOD possess competitive performances in terms of correct MTD selection as well as safety, and we apply both versions of DOD to two real oncology trials for illustration.     

On goodness of fit for time series regression models

We Propose a Bayesian approach to chech the goodness of fit for time series regression models. The test statistics is proposed by Smith (1985) based on a sequence of random variables which are independently distributed standard normal if the model is correct. We estimate this sequence of random variables using several methods. The tests of goodness of fit are performed when either the error terms violate the Gaussian assumption, or the order is incorrect, or the model is misspecified. The methodology is illustrated using both a simulation study and three real date sets.     

Alternative models for the heterogeneity of mortality risks among the aged

The authors examine how sensitive the estimates of heterogeneity in the mortality risks in a population are to the choices of two types of function, "one describing the age-specific rate of increase of mortality risks for individuals and the other describing the distribution of mortality risks across individuals." U.S. data from published Medicare mortality rates for the period 1968-1978 are used to analyze total mortality among the aged. "In addition, national vital statistics data for the period 1950-1977 were used to analyze adult lung cancer mortality. For these data, the estimates of structural parameters were less sensitive to reasonable choices of the heterogeneity distribution (gamma vs. inverse Gaussian) than to reasonable choices of the hazard rate function (Gompertz vs. Weibull)."     

Criteria for Linear Model Selection Based on Kullback's Symmetric Divergence

Model selection criteria are frequently developed by constructing estimators of discrepancy measures that assess the disparity between the 'true' model and a fitted approximating model. The Akaike information criterion (AIC) and its variants result from utilizing Kullback's directed divergence as the targeted discrepancy. The directed divergence is an asymmetric measure of separation between two statistical models, meaning that an alternative directed divergence can be obtained by reversing the roles of the two models in the definition of the measure. The sum of the two directed divergences is Kullback's symmetric divergence. In the framework of linear models, a comparison of the two directed divergences reveals an important distinction between the measures. When used to evaluate fitted approximating models that are improperly specified, the directed divergence which serves as the basis for AIC is more sensitive towards detecting overfitted models, whereas its counterpart is more sensitive towards detecting underfitted models. Since the symmetric divergence combines the information in both measures, it functions as a gauge of model disparity which is arguably more balanced than either of its individual components. With this motivation, the paper proposes a new class of criteria for linear model selection based on targeting the symmetric divergence. The criteria can be regarded as analogues of AIC and two of its variants: 'corrected' AIC or AICc and 'modified' AIC or MAIC. The paper examines the selection tendencies of the new criteria in a simulation study and the results indicate that they perform favourably when compared to their AIC analogues.     

Some New Results on the Multinomial Randomized Response Model

ABSTRACT

The randomized response technique is an effective survey method designed to elicit sensitive information while ensuring the privacy of the respondents. In this article, we present some new results on the randomization response model in situations wherein one or two response variables are assumed to follow a multinomial distribution. For a single sensitive question, we use the well-known Hopkins randomization device to derive estimates, both under the assumption of truthful and untruthful responses, and present a technique for making pairwise comparisons. When there are two sensitive questions of interest, we derive a Pearson product moment correlation estimator based on the multinomial model assumption. This estimator may be used to quantify the linear relationship between two variables when multinomial response data are observed according to a randomized-response protocol.     

Inference about parameters in Binomial-Poisson distribution with additional incomplete data

ABSTRACT

This article considers the distribution of Binomial-Poisson random vector which has two components and includes two parameters: one is the rate of a Poisson distribution, the other is the proportion in a Binomial distribution. The inference about the two parameters is usually made based on only paired observations. However, the number of paired observations is, in general, not large enough because of either technical difficulty or budget limitation, and so one can not make efficient inferences with only paired data. Instead, it is often much easier and not too costly to have incomplete observation on only one component independently. In this article we will combine both the paired complete data and unpaired incomplete data for estimating the two parameters. The performances of various estimators are compared both analytically and numerically. It is observed that fully using the unpaired incomplete data can always improve the inference, and the improvement is very significant in the case when there are only a few paired complete observations.     

An EWMA chart for monitoring process mean

In the statistical process control literature, there exists several improved quality control charts based on cost-effective sampling schemes, including the ranked set sampling (RSS) and median RSS (MRSS). A generalized cost-effective RSS scheme has been recently introduced for efficiently estimating the population mean, namely varied L RSS (VLRSS). In this article, we propose a new exponentially weighted moving average (EWMA) control chart for monitoring the process mean using VLRSS, named the EWMA-VLRSS chart, under both perfect and imperfect rankings. The EWMA-VLRSS chart encompasses the existing EWMA charts based on RSS and MRSS (named the EWMA-RSS and EWMA-MRSS charts). We use extensive Monte Carlo simulations to compute the run length characteristics of the EWMA-VLRSS chart. The proposed chart is then compared with the existing EWMA charts. It is found that, with either perfect or imperfect rankings, the EWMA-VLRSS chart is more sensitive than the EWMA-RSS and EWMA-MRSS charts in detecting small to large shifts in the process mean. A real dataset is also used to explain the working of the EWMA-VLRSS chart.     

Effects of non-normality on the power function in a one-way random model

The effects of non-normality on type-I and type-II errors in a one-way random model are investigated for moderate departures from normality. It is found that the probabilities of both errors are more sensitive to the kurtosis of between group effects than that of within group effects.     

Improved estimation methods for unrelated question randomized response techniques

The randomized response technique (RRT) is an important tool, commonly used to avoid biased answers in survey on sensitive issues by preserving the respondents’ privacy. In this paper, we introduce a data collection method for survey on sensitive issues combining both the unrelated-question RRT and the direct question design. The direct questioning method is utilized to obtain responses to a non sensitive question that is related to the innocuous question from the unrelated-question RRT. These responses serve as additional information that can be used to improve the estimation of the prevalence of the sensitive behavior. Furthermore, we propose two new methods for the estimation of the proportion of respondents possessing the sensitive attribute under a missing data setup. More specifically, we develop the weighted estimator and the weighted conditional likelihood estimator. The performances of our estimators are studied numerically and compared with that of an existing one. Both proposed estimators are more efficient than the Greenberg's estimator. We illustrate our methods using real data from a survey study on illegal use of cable TV service in Taiwan.     

Forecasting future values of changing sequences

Sequences of independent random variables are observed and on the basis of these observations future values of the process are forecast. The Bayesian predictive density of k future observations for normal, exponential, and binomial sequences which change exactly once are analyzed for several cases. It is seen that the Bayesian predictive densities are mixtures of standard probability distributions. For example, with normal sequences the Bayesian predictive density is a mixture of either normal or t-distributions, depending on whether or not the common variance is known. The mixing probabilities are the same as those occurring in the corresponding posterior distribution of the mean(s) of the sequence. The predictive mass function of the number of future successes that will occur in a changing Bernoulli sequence is computed and point and interval predictors are illustrated.     

Detection and Interpretation of a Multivariate Signal Using Combined Charts

A control procedure is presented in this article that is based on jointly using two separate control statistics in the detection and interpretation of signals in a multivariate normal process. The procedure detects the following three situations: (i) a mean vector shift without a shift in the covariance matrix; (ii) a shift in process variation (covariance matrix) without a mean vector shift; and (iii) both a simultaneous shift in the mean vector and covariance matrix as the result of a change in the parameters of some key process variables. It is shown that, following the occurrence of a signal on either of the separate control charts, the values from both of the corresponding signaling statistics can be decomposed into interpretable elements. Viewing the two decompositions together helps one to specifically identify the individual components and associated variables that are being affected. These components may include individual means or variances of the process variables as well as the correlations between or among variables. An industrial data set is used to illustrate the procedure.     

RANDOMISED RESPONSE VS. DIRECT QUESTIONING: ESTIMATING THE PREVALENCE OF ALCOHOL RELATED PROBLEMS IN A FIELD SURVEY

Fifteen hundred adults selected for a population survey from the City of Edinburgh Electoral Registers by probability sampling were allocated to randomised response interview or direct interview at random. These alternative interviewing strategies were adopted for two questions relating to attributes which might be considered stigmatising. There were no significant differences in reporting between the two methods either at aggregate or at subgroup level. The apparent failure of the randomised response method to increase estimates of the proportions with the attributes hypothesised indicates either that traditional methods are satisfactory for the collection of this sensitive information, or that randomised response was not an acceptable alternative.     

Exact distributions of constrained (k, ℓ) strings of failures between subsequent successes

Consider a sequence of binary (success–failure) random variables (RVs) ordered on a line. The number of strings with a constrained number of consecutive failures between two subsequent successes is studied under an overlapping enumeration scheme. The respective waiting time is examined as well. The study is first developed on sequences of independent and identically distributed RVs. It is extended then on sequences of dependent, exchangeability and Markovian dependency is considered, and independent, not necessarily identically distributed, RVs. Exact probabilities and moments are obtained by means of combinatorial analysis and via recursive schemes. An explicit expression of the mean value of the number of strings for both independent and dependent sequences is derived. An application in system reliability is provided.     

Cost‐efficient higher‐order crossover designs for two‐treatment clinical trials

Higher‐order crossover designs have drawn considerable attention in clinical trials, because of their ability to test direct treatment effects in the presence of carry‐over effects. The important question, when applying higher‐order crossover designs in practice, is how to choose a design with both statistical and cost efficiencies from various alternatives. In this paper, we propose a general cost function and compare five statistically optimal or near‐optimal designs with this cost function for a two‐treatment study under different carry‐over models. Based on our study, to achieve both statistical and cost efficiencies, a four‐period, four‐sequence crossover design is generally recommended under the simple carry‐over or no carry‐over models, and a three‐period, two‐sequence crossover design is generally recommended under the steady‐state carry‐over models. Copyright © 2005 John Wiley & Sons, Ltd.     

On Kronecker block designs for factorial experiments

In this paper the use of Kronecker designs for factorial experiments is considered. The two-factor Kronecker design is considered in some detail and the efficiency factors of the main effects and interaction in such a design are derived. It is shown that the efficiency factor of the interaction is at least as large as the product of the efficiency factors of the two main effects and when both the component designs are totally balanced then its efficiency factor will be higher than the efficiency factor of either of the two main effects. If the component designs are nearly balanced then its efficiency factor will be approximately at least as large as the efficiency factor of either of the two main effects. It is argued that these designs are particularly useful for factorial experiments.Extensions to the multi-factor design are given and it is proved that the two-factor Kronecker design will be connected if the component designs are connected.     

Run statistics in a sequence of arbitrarily dependent binary trials

Let {Z _i } _i≥1 be an arbitrary sequence of trials with two possible outcomes either success (1) or failure (0). General expressions for the exact distributions of runs, both success and failure, in Z ₁, . . . , Z _n are presented. Our method is based on the use of joint distribution of success and failure run lengths and unifies the results on distribution of runs. As a special case of our results we obtain the distributions of runs for various binary sequences. As illustrated in the paper the results enable us to derive the distribution of runs for binary trials arising in urn models.     

Variable sampling intervals for multiparameter shewhart charts

This paper considers the problem of using control charts to simultaneously monitor more than one parameter with emphasis on simultaneously monitoring the mean and variance. Fixed sampling interval control charts are modified to use variable sampling intervals depending on what is being observed from the data. Two basic strategies are investigated. One strategy uses separate control charts for each parameter, A second strategy uses a proposed single combined statistic which is sensitive to shifts in both the mean and variance. Each procedure is compared to corresponding fixed interval procedures. It is seen that for both strategies the variable sampling interval approach is substantially more efficient than fixed interval procedures.     

New tests for stochastic order with application to case control studies

It is shown that if a binary regression function is increasing then retrospective sampling induces a stochastic ordering of the covariate distributions among the responders, which we call cases, and the non-responders, which we call controls. We also show that if the covariate distributions are stochastically ordered then the regression function must be increasing. This means that testing whether the regression function is monotone is equivalent to testing whether the covariate distributions are stochastically ordered. Capitalizing on these new probabilistic observations we proceed to develop two new non-parametric tests for stochastic order. The new tests are based on either the maximally selected, or integrated, chi-bar statistic of order one. The tests are easy to compute and interpret and their large sampling distributions are easily found. Numerical comparisons show that they compare favorably with existing methods in both small and large samples. We emphasize that the new tests are applicable to any testing problem involving two stochastically ordered distributions.