Inverse Sampling with Unequal Selection Probabilities

In some real situations the population of interest is divided into two groups, of which one contains only a few units. In other cases, the population may be considered as subdivided into two group', for example, if only a few units display a value of the variable of interest which is highly different from zero, while all the other units show a value equal to or near zero. In both cases, inverse sampling is more efficient than classical fixed sample-size designs to obtain the parameter estimators for the whole population as well as for its groups (e.g., Salehi and Seber, 2004 Salehi , M. M. , Seber , G. A. F. ( 2004 ). A general inverse sampling scheme and its application to adaptive cluster sampling . Austral. NZ J. Statist. 46 : 483 – 494 .[Crossref], [Web of Science ®] , [Google Scholar]). In fact, in this design the procedure selection continues until a prefixed number of units with the characteristic of interest is sampled. Since it is not known a priori to which group the population units belong, the sample size is a random variable. Christman and Lan (2001 Christman , M. C. , Lan , F. ( 2001 ). Inverse adaptive cluster sampling . Biometrics 57 : 1096 – 1105 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) and Salehi and Seber (2001 Salehi , M. M. , Seber , G. A. F. ( 2001 ). A new proof of Murthy's estimator which applies to sequential sampling . Austral. NZ J. Statist. 43 : 281 – 286 . [Google Scholar] 2004 Salehi , M. M. , Seber , G. A. F. ( 2004 ). A general inverse sampling scheme and its application to adaptive cluster sampling . Austral. NZ J. Statist. 46 : 483 – 494 .[Crossref], [Web of Science ®] , [Google Scholar]) considered inverse sampling designs when all the population units have equal selection probabilities. In this article, we consider the general case in which the units may have unequal probabilities of being included in the sample. In fact, in many real situations different units may have different selection probabilities because of some inherent features of the sampling procedure, or in order to obtain better estimates. We derive unbiased estimators of the totals of the two groups, their variance and the corresponding unbiased variance estimators in inverse sampling with replacement. Finally, we derive similar results for more complex designs, where the selection procedure stops before observing the prefixed number of units from the rare group.     

Statistical Modeling of Temporal Dependence in Financial Data via a Copula Function

In financial analysis it is useful to study the dependence between two or more time series as well as the temporal dependence in a univariate time series. This article is concerned with the statistical modeling of the dependence structure in a univariate financial time series using the concept of copula. We treat the series of financial returns as a first order Markov process. The Archimedean two-parameter BB7 copula is adopted to describe the underlying dependence structure between two consecutive returns, while the log-Dagum distribution is employed to model the margins marked by skewness and kurtosis. A simulation study is carried out to evaluate the performance of the maximum likelihood estimates. Furthermore, we apply the model to the daily returns of four stocks and, finally, we illustrate how its fitting to data can be improved when the dependence between consecutive returns is described through a copula function.     

A Permutation Approach to Split-Plot Experiments

Split-plot design may be refer to a common experimental setting where a particular type of restricted randomization has occurred during a planned experiment. The aim of this article is to suggest a new method to perform inference on split-plot experiments by combination-based permutation tests. This novel nonparametric approach has been studied and validated using a Monte Carlo simulation study where we compared it with the parametric and nonparametric procedures proposed in the literature. Results suggest that in each experimental situation where normality is hard to justify and especially when errors have heavy-tailed distribution, the proposed nonparametric procedure can be considered as a valid solution.     

On Bayesian learning via loss functions

We provide a decision theoretic approach to the construction of a learning process in the presence of independent and identically distributed observations. Starting with a probability measure representing beliefs about a key parameter, the approach allows the measure to be updated via the solution to a well defined decision problem. While the learning process encompasses the Bayesian approach, a necessary asymptotic consideration then actually implies the Bayesian learning process is best. This conclusion is due to the requirement of posterior consistency for all models and of having standardized losses between probability distributions. This is shown considering a specific continuous model and a very general class of discrete models.     

Multivariate nonparametric methods in two-way balanced designs: performances and limitations in small samples

Investigations of multivariate population are pretty common in applied researches, and the two-way crossed factorial design is a common design used at the exploratory phase in industrial applications. When assumptions such as multivariate normality and covariance homogeneity are violated, the conventional wisdom is to resort to nonparametric tests for hypotheses testing. In this paper we compare the performances, and in particular the power, of some nonparametric and semi-parametric methods that have been developed in recent years. Specifically, we examined resampling methods and robust versions of classical multivariate analysis of variance (MANOVA) tests. In a simulation study, we generate data sets with different configurations of factor''s effect, number of replicates, number of response variables under null hypothesis, and number of response variables under alternative hypothesis. The objective is to elicit practical advice and guides to practitioners regarding the sensitivity of the tests in the various configurations, the tradeoff between power and type I error, the strategic impact of increasing number of response variables, and the favourable performance of one test when the alternative is sparse. A real case study from an industrial engineering experiment in thermoformed packaging production is used to compare and illustrate the application of the various methods.     

Fractional Moments and Maximum Entropy: Geometric Meaning

We consider the problem of recovering a probability density on a bounded or unbounded subset D of [0, ∞), from the knowledge of its sequence of fractional moments within a maximum entropy (MaxEnt) setup. Based upon entropy convergence results previously formulated, the fractional moments are selected so that the entropy of the MaxEnt approximation be minimum. A geometric interpretation of the reconstruction procedure is formulated as follows: the two moment curves generated by the unknown density and its MaxEnt approximation are interpolating in Hermite-Birkoff sense; that is, they are both interpolating and tangent at the selected nodes.     

An R package for the simulation of correlated discrete variables

A package for the stochastic simulation of discrete variables with assigned marginal distributions and correlation matrix is presented and discussed. The simulating mechanism relies upon the Gaussian copula, linking the discrete distributions together, and an iterative scheme recovering the correlation matrix for the copula that ensures the desired correlations among the discrete variables. Examples of its use are provided as well as three possible applications (related to probability, sampling, and inference), which illustrate the utility of the package as an efficient and easy-to-use tool both in statistical research and for didactic purposes.