A latent variable model for analyzing mixed longitudinal (k,l)-inflated count and ordinal responses

A random effects model for analyzing mixed longitudinal count and ordinal data is presented where the count response is inflated in two points (k and l) and an (k,l)-Inflated Power series distribution is used as its distribution. A full likelihood-based approach is used to obtain maximum likelihood estimates of parameters of the model. For data with non-ignorable missing values models with probit model for missing mechanism are used.The dependence between longitudinal sequences of responses and inflation parameters are investigated using a random effects approach. Also, to investigate the correlation between mixed ordinal and count responses of each individuals at each time, a shared random effect is used. In order to assess the performance of the model, a simulation study is performed for a case that the count response has (k,l)-Inflated Binomial distribution. Performance comparisons of count-ordinal random effect model, Zero-Inflated ordinal random effects model and (k,l)-Inflated ordinal random effects model are also given. The model is applied to a real social data set from the first two waves of the national longitudinal study of adolescent to adult health (Add Health study). In this data set, the joint responses are the number of days in a month that each individual smoked as the count response and the general health condition of each individual as the ordinal response. For the count response there is incidence of excess values of 0 and 30.     

Whither practice-near research in the modernization programme? Policy blunders in children's services

In this article, we lament the effects of practice-distant research and associated policy initiatives on contemporary children's services in England. In the last decade, as a result of high profile inquiries into non-accidental child deaths, statutory children's social care services in the UK have been subject to a wide-reaching ‘modernization’ programme. We studied decision-making in the high blame environment of local authority children's services. Our research sought to examine the relationship between performance management and the impact of anticipated blame within the decision-making practices of those providing, supervising and managing these services. We show that systems and technologies can be developed which both assist the users in their daily work and achieve desired organizational goals, but without an ethnographically informed, practice-near approach, unsafe work regimes and practices can ensue.     

To be or not to be valid in testing the significance of the slope in simple quantitative linear models with autocorrelated errors

In this article, the validity of procedures for testing the significance of the slope in quantitative linear models with one explanatory variable and first-order autoregressive [AR(1)] errors is analyzed in a Monte Carlo study conducted in the time domain. Two cases are considered for the regressor: fixed and trended versus random and AR(1). In addition to the classical t -test using the Ordinary Least Squares (OLS) estimator of the slope and its standard error, we consider seven t -tests with n-2\,\hbox{df} built on the Generalized Least Squares (GLS) estimator or an estimated GLS estimator, three variants of the classical t -test with different variances of the OLS estimator, two asymptotic tests built on the Maximum Likelihood (ML) estimator, the F -test for fixed effects based on the Restricted Maximum Likelihood (REML) estimator in the mixed-model approach, two t -tests with n - 2 df based on first differences (FD) and first-difference ratios (FDR), and four modified t -tests using various corrections of the number of degrees of freedom. The FDR t -test, the REML F -test and the modified t -test using Dutilleul's effective sample size are the most valid among the testing procedures that do not assume the complete knowledge of the covariance matrix of the errors. However, modified t -tests are not applicable and the FDR t -test suffers from a lack of power when the regressor is fixed and trended ( i.e. , FDR is the same as FD in this case when observations are equally spaced), whereas the REML algorithm fails to converge at small sample sizes. The classical t -test is valid when the regressor is fixed and trended and autocorrelation among errors is predominantly negative, and when the regressor is random and AR(1), like the errors, and autocorrelation is moderately negative or positive. We discuss the results graphically, in terms of the circularity condition defined in repeated measures ANOVA and of the effective sample size used in correlation analysis with autocorrelated sample data. An example with environmental data is presented.     

Finite Sample Performance of tests for Symmetry of the Errors in a Linear Model

The finite sample performance of a number of tests for symmetry of the distribution of the errors of a linear model is considered. The first family of tests is based on the discrepancy between two regression fits. The first fit is appropriate under symmetric errors while the second is appropriate for skewed as well as symmetric error distributions. The second family of procedures consists of tests for the univariate symmetry problem. Thus, in the linear model setting these tests are based on residuals. An extensive empirical study of the finite sample, null behavior of the tests is presented. The results of a power comparison among the tests is also discussed.     

美学:爱美之学——从辞源学的角度为美学正名

美学与审美学是两个存在差异的学术范畴,由于西语只有aesthetics一个术语,致使长期以来二者混淆在一起。为了厘清二者的边界,为美学的复兴开辟道路,为美学正名成为当务之急。学界通常所谓的aesthetics("美学")其实是"审美学","美学"应该另有其名。两百年来,美学家为"美学"创造了多个名字,其中以callology最为恰当。但是,真正最为恰切精当的"美学"一词是源自古希腊的philocaly(爱美之学)。     

Rank invariant tests with left truncated and interval censored data

Peto and Peto (1972) have studied rank invariant tests to compare two survival curves for right censored data. We apply their tests, including the logrank test and the generalized Wilcoxon test, to left truncated and interval censored data. The significance levels of the tests are approximated by Monte Carlo permutation tests. Simulation studies are conducted to show their size and power under different distributional differences. In particular, the logrank test works well under the Cox proportional hazards alternatives, as for the usual right censored data. The methods are illustrated by the analysis of the Massachusetts Health Care Panel Study dataset.     

A bootstrap control chart for inverse Gaussian percentiles

The conventional Shewhart-type control chart is developed essentially on the central limit theorem. Thus, the Shewhart-type control chart performs particularly well when the observed process data come from a near-normal distribution. On the other hand, when the underlying distribution is unknown or non-normal, the sampling distribution of a parameter estimator may not be available theoretically. In this case, the Shewhart-type charts are not available. Thus, in this paper, we propose a parametric bootstrap control chart for monitoring percentiles when process measurements have an inverse Gaussian distribution. Through extensive Monte Carlo simulations, we investigate the behaviour and performance of the proposed bootstrap percentile charts. The average run lengths of the proposed percentage charts are investigated.     

Optimal hedging through limit orders

This article applies the methods of stochastic dynamic programming to a risk management problem, where an agent hedges her derivative position by submitting limit orders. Therefore, this model is the first, in the literature on optimal trading with limit orders, to handle a problem of hedging options or other derivatives. A hedging strategy is developed where both the size and the limit price of each order is optimally set.     

Best-worst multi-criteria decision-making method: Some properties and a linear model

The Best Worst Method (BWM) is a multi-criteria decision-making method that uses two vectors of pairwise comparisons to determine the weights of criteria. First, the best (e.g. most desirable, most important), and the worst (e.g. least desirable, least important) criteria are identified by the decision-maker, after which the best criterion is compared to the other criteria, and the other criteria to the worst criterion. A non-linear minmax model is then used to identify the weights such that the maximum absolute difference between the weight ratios and their corresponding comparisons is minimized. The minmax model may result in multiple optimal solutions. Although, in some cases, decision-makers prefer to have multiple optimal solutions, in other cases they prefer to have a unique solution. The aim of this paper is twofold: firstly, we propose using interval analysis for the case of multiple optimal solutions, in which we show how the criteria can be weighed and ranked. Secondly, we propose a linear model for BWM, which is based on the same philosophy, but yields a unique solution.     

A Diversity Index for Model Space Selection in the Estimation of Benchmark and Infectious Doses via Model Averaging

In chemical and microbial risk assessments, risk assessors fit dose‐response models to high‐dose data and extrapolate downward to risk levels in the range of 1–10%. Although multiple dose‐response models may be able to fit the data adequately in the experimental range, the estimated effective dose (ED) corresponding to an extremely small risk can be substantially different from model to model. In this respect, model averaging (MA) provides more robustness than a single dose‐response model in the point and interval estimation of an ED. In MA, accounting for both data uncertainty and model uncertainty is crucial, but addressing model uncertainty is not achieved simply by increasing the number of models in a model space. A plausible set of models for MA can be characterized by goodness of fit and diversity surrounding the truth. We propose a diversity index (DI) to balance between these two characteristics in model space selection. It addresses a collective property of a model space rather than individual performance of each model. Tuning parameters in the DI control the size of the model space for MA.