Optimal designs for a mixed interference model

This paper generalizes Kunert and Martin’s (Ann Stat 28:1728–1742, 2000) method for finding optimal designs under a fixed interference model, to find optimal designs under a mixed interference model. The results are based on the properties of information matrices in fixed and mixed models given in Markiewicz (J Stat Plan Inference 59:127–137, 1997). The method is applied to find a design which is optimal for any given variances of random neighbor effects. Research partially supported by the KBN Grant Number 5 P03A 041 21.     

Estimating sensitive proportions from Warner’s randomized responses in alternative ways restricting to only distinct units sampled

Warner’s (J Am Stat Assoc 60:63–69, 1965) pioneering ‘randomized response’ (RR) technique (RRT), useful in unbiasedly estimating the proportion of people bearing a sensitive characteristic, is based exclusively on simple random sampling (SRS) with replacement (WR). Numerous developments that follow from it also use SRSWR and employ in estimation the mean of the RR’s yielded by each person everytime he/she is drawn in the sample. We examine how the accuracy in estimation alters when instead we use the mean of the RR’s from each distinct person sampled and also alternatively if the Horvitz and Thompson’s (HT, J Am Stat Assoc 47:663–685, 1952) method is employed for an SRSWR eliciting only, a single RR from each. Arnab’s (1999) approach of using repeated RR’s from each distinct person sampled is not taken up here to avoid algebraic complexity. Pathak’s (Sánkhya 24:287–302, 1962) results for direct response (DR) from SRSWR are modified with the use of such RR’s. Through simulations we present relative levels in accuracy in estimation by the above three alternative methods.     

On the detectability of different forms of interaction in regression models

We consider the ability to detect interaction structure from data in a regression context. We derive an asymptotic power function for a likelihood-based test for interaction in a regression model, with possibly misspecified alternative distribution. This allows a general investigation of different types of interactions which are poorly or well detected via data. Principally we contrast pairwise-interaction models with ‘diffuse interaction models’ as introduced in Gustafson et al. (Stat Med 24:2089–2104, 2005).     

On the calibration of design weights using a displacement function

In the present investigation, we propose a new method to calibrate the estimator of the general parameter of interest in survey sampling. We demonstrate that the linear regression estimator due to Hansen et al. (Sample Survey Method and Theory. Wiley, NY, 1953) is a special case of this. We reconfirm that the sum of calibrated weights has to be set equal to sum of the design weights within a given sample as shown in Singh (Advanced sampling theory with applications: How Michael ‘selected’ Amy, Vol. 1 and 2. Kluwer, The Netherlands, pp 1–1247, 2003; Proceedings of the American Statistical Association, Survey Method Section [CD-ROM], Toronto, Canada: American Statistical Association, pp 4382–4389, 2004; Metrika:1–18, 2006a; Presented at INTERFACE 2006, Pasadena, CA, USA, 2006b) and Stearns and Singh (Presented at Joint Statistical Meeting, MN, USA (Available on the CD), 2005; Comput Stat Data Anal 52:4253–4271, 2008). Thus, it shows that the Sir. R.A. Fisher’s brilliant idea of keeping sum of observed frequencies equal to that of expected frequencies leads to a “Honest-Balance” while weighing design weights in survey sampling. The major benefit of the proposed new estimator is that it always works unlike the pseudo empirical likelihood estimators listed in Owen (Empirical Likelihood. Chapman & Hall, London, 2001), Chen and Sitter (Stat Sin 9:385–406, 1999) and Wu (Sur Methodol 31(2):239–243, 2005). The main endeavor of this paper is to bring a change in the existing calibration technology, which is based on only positive distance functions, with a displacement function that has the flexibility of taking positive, negative, or zero value. At the end, the proposed technology has been compared with its competitors under several kinds of linear and non-linear non-parametric models using an extensive simulation study. A couple of open questions are raised.     

Robust estimation for the order of finite mixture models

In this paper, we study a robust and efficient estimation procedure for the order of finite mixture models based on the minimizing a penalized density power divergence estimator. For this task, we use the locally conic parametrization approach developed by Dacunha-Castelle and Gassiate (ESAIM Probab Stat 285–317, 1997a; Ann Stat 27:1178–1209, 1999), and verify that the minimizing a penalized density power divergence estimator is consistent. Simulation results are provided for illustration.     

On fitting the conditional difference asymmetry models to square contingency tables with nominal categories

In this paper, we implement the conditional difference asymmetry model (CDAS) for square tables with nominal categories proposed by Tomizawa et al. (J. Appl. Stat. 31(3): 271–277, 2004) with the use of the non-standard log-linear model formulation approach. The implementation is carried out by refitting the model in the 3 ×  3 table in (Tomizawa et al. J. Appl. Stat. 31(3): 271–277, 2004). We extend this approach to a larger 4 ×  4 table of religious affiliation. We further calculated the measure of asymmetry along with its asymptotic standard error and confidence bounds. The procedure is implemted with SAS PROC GENMOD but can also be implemented in SPSS by following the discussion in (Lawal, J. Appl. Stat. 31(3): 279–303, 2004; Lawal, Qual. Quant. 38(3): 259–289, 2004).     

Asymptotics of SIMEX-based variance estimation

In this paper, we study the asymptotic properties of simulation extrapolation (SIMEX) based variance estimation that was proposed by Wang et al. (J R Stat Soc Series B 71:425–445, 2009). We first investigate the asymptotic normality of the parameter estimator in general parametric variance function and the local linear estimator for nonparametric variance function when permutation SIMEX (PSIMEX) is used. The asymptotic optimal bandwidth selection with respect to approximate mean integrated squared error (AMISE) for nonparametric estimator is also studied. We finally discuss constructing confidence intervals/bands of the parameter/function of interest. Other than applying the asymptotic results so that normal approximation can be used, we recommend a nonparametric Monte Carlo algorithm to avoid estimating the asymptotic variance of estimator. Simulation studies are carried out for illustration.     

Order Effects in Batteries of Questions

Batteries of questions with identical response items are commonly used in survey research. This paper suggests that question order has the potential to cause systematic positive or negative bias on responses to all questions in a battery. Whilst question order effects have been studied for many decades, almost no attention has been given to this topic. The primary aim is to draw attention to this effect, to demonstrate its possible magnitude, and to discuss a range of mechanisms through which it might occur. These include satisficing, anchoring and cooperativeness. The effect seems apparent in the results of a recent survey. This was a survey of Emergency Department patients presenting to Wollongong Hospital (Australia) with apparently less urgent conditions in 2004. Two samples were taken. Question order was fixed in the first sample (n = 104; response rate RR2 = 94%), but randomised in the second sample (n = 46; response rate RR2 = 96%). Respondents were asked to indicate whether each of 18 reasons for presenting to the ED was a ‘very important reason’ a ‘moderately important reason’ or ‘not a reason’ The mean number of very important reasons selected was 56% higher in the first sample as compared to the second sample.     

Multidimensional strong large deviation results

We establish strong large deviation results for an arbitrary sequence of random vectors under some assumptions on the normalized cumulant generating function. In other words, we give asymptotic approximations for a multivariate tail probability of the same kind as the one obtained by Bahadur and Rao (Ann Math Stat 31:1015–1027, 1960) for the sample mean (in the one-dimensional case). The proof of our results follows the same lines as in Chaganty and Sethuraman (J Stat Plan Inference, 55:265–280, 1996). We also present three statistical applications to illustrate our results, the first one dealing with a vector of independent sample variances, the second one with a Gaussian multiple linear regression model and the third one with the multivariate Nadaraya–Watson estimator. Some numerical results are also presented for the first two applications.     

Response of double-auction markets to instantaneous Selling�CBuying signals with stochastic Bid�CAsk spread

Statistical properties of order-driven double-auction markets with Bid–Ask spread are investigated through the dynamical quantities such as response function. We first attempt to utilize the so-called Madhavan–Richardson–Roomans model (MRR for short) to simulate the stochastic process of the price-change in empirical data sets (say, EUR/JPY or USD/JPY exchange rates) in which the Bid–Ask spread fluctuates in time. We find that the MRR theory apparently fails to simulate so much as the qualitative behaviour (‘non-monotonic’ behaviour) of the response function R(l) (l denotes the difference of times at which the response function is evaluated) calculated from the data. Especially, we confirm that the stochastic nature of the Bid–Ask spread causes apparent deviations from a linear relationship between the R(l) and the auto-correlation function C(l), namely, R(l) μ -C(l){R(l) \propto -C(l)}. To make the microscopic model of double-auction markets having stochastic Bid–Ask spread, we use the minority game with a finite market history length and find numerically that appropriate extension of the game shows quite similar behaviour of the response function to the empirical evidence. We also reveal that the minority game modeling with the adaptive (‘annealed’) look-up table reproduces the non-linear relationship R(l) μ -f(C(l)){R(l) \propto -f(C(l))} (f(x) stands for a non-linear function leading to ‘λ-shapes’) more effectively than the fixed (‘quenched’) look-up table does.     

Durbin’s random substitution and conditional Monte Carlo

Durbin (Biometrika 48:41–55, 1961) proposed a method called random substitution, by which a composite problem of goodness-of-fit can be reduced to a simple one. In this paper we provide a method of finding the p-value of any test statistic, for a composite goodness-of-fit problem, based on the simulation of a large number of conditional samples, using an analog of Durbin’s proposal in a reverse-type application. We analyze a Bayesian chi-square test proposed in Johnson (Ann Stat 32:2361–2384, 2004) which relies on a single randomization and relate it with Durbin’s original method. We also review a related proposal for conditional Monte-Carlo simulation in Lindqvist and Taraldsen (Biometrika 92:451–464, 2005) and compare it with our procedure. We show our method in a non-group example introduced in Lindqvist and Taraldsen (Biometrika 90:489–490, 2003).     

A note on multinomial maximum likelihood estimation under ordered restrictions and the EM algorithm

A frequently occurring problem is to find the maximum likelihood estimation (MLE) of p subject to p∈C (C⊂ P the probability vectors in R ^k ). The problem has been discussed by many authors and they mainly focused when p is restricted by linear constraints or log-linear constraints. In this paper, we construct the relationship between the the maximum likelihood estimation of p restricted by p∈C and EM algorithm and demonstrate that the maximum likelihood estimator can be computed through the EM algorithm (Dempster et al. in J R Stat Soc Ser B 39:1–38, 1997). Several examples are analyzed by the proposed method.     

Identification-commitment inventory (ICI model): confirmatory factor analysis and construct validity

The aim of this study is to confirm the factorial structure of the Identification-Commitment Inventory (ICI) developed within the frame of the Human System Audit (HSA) (Quijano et al. in Revist Psicol Soc Apl 10(2):27–61, 2000; Pap Psicól Revist Col Of Psicó 29:92–106, 2008). Commitment and identification are understood by the HSA at an individual level as part of the quality of human processes and resources in an organization; and therefore as antecedents of important organizational outcomes, such as personnel turnover intentions, organizational citizenship behavior, etc. (Meyer et al. in J Org Behav 27:665–683, 2006). The theoretical integrative model which underlies ICI Quijano et al. (2000) was tested in a sample (N = 625) of workers in a Spanish public hospital. Confirmatory factor analysis through structural equation modeling was performed. Elliptical least square solution was chosen as estimator procedure on account of non-normal distribution of the variables. The results confirm the goodness of fit of an integrative model, which underlies the relation between Commitment and Identification, although each one is operatively different.     

A new condition for identifiability of finite mixture distributions

In this paper a sufficient condition for the identifiability of finite mixtures is given. This condition is less restrictive than Teicher’s condition Teicher H, Ann Math Stat 34:1265–1269 (1963) and therefore it can be applied to a wider range of families of mixtures. In particular, it applies to the classes of all finite mixtures of Log-gamma and of reversed Log-gamma distributions. These families have been already studied by Henna J Jpn Stat Soc 24:193–200 (1994) using another condition, different from Teicher’s, but more difficult to check in many cases. Furthermore, the result given in this paper is very appropiated for the case of mixtures of the union of different distribution families. To illustrate this an application to the class of all finite mixtures generated by the union of Lognormal, Gamma and Weibull distributions is given, where Teicher’s and Henna’s conditions are not applicable     

A unified approach to bivariate discrete distributions

Here we introduce a bivariate generalized hypergeometric factorial moment distribution (BGHFMD) through its probability generating function (p.g.f.) whose marginal distributions are the generalized hypergeometric factorial moment distributions introduced by Kemp and Kemp (Bull Int Stat Inst 43:336–338,1969). Well-known bivariate versions of distributions such as binomial, negative binomial and Poisson are special cases of this distribution. A genesis of the distribution and explicit closed form expressions for the probability mass function of the BGHFMD, its factorial moments and the p.g.f.’s of its conditional distributions are derived here. Certain recurrence relations for probabilities, moments and factorial moments of the bivariate distribution are also established.     

Estimating regression coefficients from survey data by asymptotic design-cum-model based approach

Postulating a linear regression of a variable of interest on an auxiliary variable with values of the latter known for all units of a survey population, we consider appropriate ways of choosing a sample and estimating the regression parameters. Recalling Thomsen’s (1978) results on non-existence of ‘design-cum-model’ based minimum variance unbiased estimators of regression coefficients we apply Brewer’s (1979) ‘asymptotic’ analysis to derive ‘asymptotic-design-cummodel’ based optimal estimators assuming large population and sample sizes. A variance estimation procedure is also proposed.     

Estimating a restricted normal mean

Let X ₁, X ₂, ..., X _n be a random sample from a normal distribution with unknown mean μ and known variance σ ². In many practical situations, μ is known a priori to be restricted to a bounded interval, say [−m, m] for some m > 0. The sample mean , then, becomes an inadmissible estimator for μ. It is also not minimax with respect to the squared error loss function. Minimax and other estimators for this problem have been studied by Casella and Strawderman (Ann Stat 9:870–878, 1981), Bickel (Ann Stat 9:1301–1309, 1981) and Gatsonis et al. (Stat Prob Lett 6:21–30, 1987) etc. In this paper, we obtain some new estimators for μ. The case when the variance σ ² is unknown is also studied and various estimators for μ are proposed. Risk performance of all estimators is numerically compared for both the cases when σ ² may be known and unknown.     

Of songs and men: a model for multiple choice with herding

We propose a generic model for multiple choice situations in the presence of herding and compare it with recent empirical results from a Web-based music market experiment. The model predicts a phase transition between a weak imitation phase and a strong imitation, ‘fashion’ phase, where choices are driven by peer pressure and the ranking of individual preferences is strongly distorted at the aggregate level. The model can be calibrated to reproduce the main experimental results of Salganik et al. (Science, 311, 854–856 (2006)); we show in particular that the value of the social influence parameter can be estimated from the data. In one of the experimental situation, this value is found to be close to the critical value of the model.     

Randomized response: Estimating mean square errors of linear estimators and finding optimal unbiased strategies

Summary General procedures are described to generate quantitative randomized response (RR) required to estimate the finite population total of a sensitive variable. Permitting sample selection with arbitrary probabilities a formula for the mean square error (MSE) of a linear estimator of total based on RR is noted indicating the simple modification over one that might be based on direct response (DR) if the latter were available. A general formula for an unbiased estimator of the MSE is presented. A simple approximation is proposed in case the RR ratio estimator is employed based on a simple random sample (SRS) taken without replacement (WOR). Among sampling strategies employing unbiased but not necessarily linear estimators based on RR, certain optimal ones are identified under two alternative models analogously to well-known counterparts based on DR, if available. Unlike Warner’s (1965) treatment of categorical RR we consider quantitative RR here.     

Factor analysis of the Big Five Questionnaire using polychoric correlations in children

Although the Big Five Questionnaire for children (BFQ-C) (C. Barbaranelli et al., Manuale del BFQ-C. Big Five Questionnaire Children, O.S. Organizazioni, Firenze, 1998) is an ordinal scale, its dimensionality has often been studied using factor analysis with Pearson correlations. In contrast, the present study takes this ordinal metric into account and examines the dimensionality of the scale using factor analysis based on a matrix of polychoric correlations. The sample comprised 852 subjects (40.90% boys and 59.10% girls). As in previous studies the results obtained through exploratory factor analysis revealed a five-factor structure (extraversion, agreeableness, conscientiousness, emotional instability and openness). However, the results of the confirmatory factor analysis were consistent with both a four and five-factor structure, the former showing a slightly better fit and adequate theoretical interpretation. These data confirm the need for further research as to whether the factor ‘Openness’ should be maintained as an independent factor (five-factor model), or whether it would be best to omit it and distribute its items among the factors ‘Extraversion’ and ‘Conscientiousness’ (four-factor model).