Estimating selective advantage of two alleles in discrete time

A class of estimators for the selective advantage, s, in a Wright-Fisher model with two alleles, variable population size, and genic selection is derived via martingale theory. Explicit expressions are given for these estimators which only involve simple computation. The optimal estimate among this class of estimators is obtained. Asymptotic results are readily established by an application of a martingale central limit theorem. The performance of this optimal estimator is compared to known estimators by means of a simulation study.     

On optimal Bayesian classification and risk estimation under multiple classes

A recently proposed optimal Bayesian classification paradigm addresses optimal error rate analysis for small-sample discrimination, including optimal classifiers, optimal error estimators, and error estimation analysis tools with respect to the probability of misclassification under binary classes. Here, we address multi-class problems and optimal expected risk with respect to a given risk function, which are common settings in bioinformatics. We present Bayesian risk estimators (BRE) under arbitrary classifiers, the mean-square error (MSE) of arbitrary risk estimators under arbitrary classifiers, and optimal Bayesian risk classifiers (OBRC). We provide analytic expressions for these tools under several discrete and Gaussian models and present a new methodology to approximate the BRE and MSE when analytic expressions are not available. Of particular note, we present analytic forms for the MSE under Gaussian models with homoscedastic covariances, which are new even in binary classification.     

Ratio estimation with measurement error in the auxiliary variate

Summary .  With auxiliary information that is well correlated with the primary variable of interest, ratio estimation of the finite population total may be much more efficient than alternative estimators that do not make use of the auxiliary variate. The well-known properties of ratio estimators are perturbed when the auxiliary variate is measured with error. In this contribution we examine the effect of measurement error in the auxiliary variate on the design-based statistical properties of three common ratio estimators. We examine the case of systematic measurement error as well as measurement error that varies according to a fixed distribution. Aside from presenting expressions for the bias and variance of these estimators when they are contaminated with measurement error we provide numerical results based on a specific population. Under systematic measurement error, the biasing effect is asymmetric around zero, and precision may be improved or degraded depending on the magnitude of the error. Under variable measurement error, bias of the conventional ratio-of-means estimator increased slightly with increasing error dispersion, but far less than the increased bias of the conventional mean-of-ratios estimator. In similar fashion, the variance of the mean-of-ratios estimator incurs a greater loss of precision with increasing error dispersion compared with the other estimators we examine. Overall, the ratio-of-means estimator appears to be remarkably resistant to the effects of measurement error in the auxiliary variate.     

Statistical Analysis of Nonrectangular Family Data

MINQUE (Minimum Norm Quadratic Unbiased Estimators) theory is applied to the problem of estimation of variance components in family data (siblings) with variable family size. Using this approach, the traditional iterative maximum likelihood estimators are shown to be asymptotically normal, even though the data come from non-identical parent distributions. Asymptotic expressions are also obtained for the variance of the MINQUE estimators which hold even if the data are decidedly non-normal (e.g. a mixture of normals). In the case of normal data, exact small-sample variance estimates are derived. Simulations demonstrate the fast rate of convergence to asymptotic properties as the number of families increases. These desirable qualities suggest that the easy to compute MINQUE class of estimators may provide a useful alternative method for modelling familial aggregation.     

A Combined Product Estimator in Sample Survey

For the estimation of the population mean in stratified random sampling a ‘Combined Product Estimator’ is proposed which is more efficient than the ‘Combined Ratio’ and ‘Separate Ratio’ estimators. Also, the proposed estimator have exact expressions for bias and mean square error. An empirical illustration is given to compare the efficiencies of different estimators.     

Bayesian inference for disease prevalence using negative binomial group testing

Group testing, also known as pooled testing, and inverse sampling are both widely used methods of data collection when the goal is to estimate a small proportion. Taking a Bayesian approach, we consider the new problem of estimating disease prevalence from group testing when inverse (negative binomial) sampling is used. Using different distributions to incorporate prior knowledge of disease incidence and different loss functions, we derive closed form expressions for posterior distributions and resulting point and credible interval estimators. We then evaluate our new estimators, on Bayesian and classical grounds, and apply our methods to a West Nile Virus data set.     

Non-Iterative Least Squares Estimation of Missing Values in Graeco-Latin Square Designs

The present paper deals with the estimation of several missing values in Graeco-Latin Square-designs. When the observations are missing in a particular pattern, explicit computable expressions are presented for the estimators of the missing values. This procedure is illustrated with the help of a numerical example.     

The similarity index and DNA fingerprinting

DNA-fingerprint similarity is being used increasingly to make inferences about levels of genetic variation within and between natural populations. It is shown that the similarity index--the average fraction of shared restriction fragments--provides upwardly biased estimates of population homozygosity but nearly unbiased estimates of the average identity-in-state for random pairs of individuals. A method is suggested for partitioning the DNA-fingerprint dissimilarity into within- and between-population components. Some simple expressions are given for the sampling variances of these estimators.     

Non-Iterative Least Squares Estimation of Missing Values in Cross-Over Designs Without Residual Effect

In the present study an attempt has been made to esttimate several missing values in cross-over designs. When the observations are missing in a particular pattern, explicit expressions are given for the estimators of the missing values. This procedure is illustrated with the help of a numerical example.     

Estimation of Several Missing Values in F-Square Designs

The present study deals with the estimation of several missing values in F-square designs. The estimating equations for the non-iterative least squares estimation of Missing Values and explicit expressions for the estimators of the particular patterns of Missing Values are presented. This procedure is illustrated with the help of a numerical example.     

An Estimate of the Variance of Estimators for Lead Time and Screening Benefit Time in Randomized Cancer Screening Trials

Variance estimators are derived for estimators of the average lead time and average benefit time due to screening in a randomized screening trial via influence functions. The influence functions demonstrate that these estimators are asymptotically equivalent to the mean difference, between the study and control case groups, in the appropriate survival times. For estimating benefit time, the survival time is measured since start of study; for estimating lead time, the survival time is measured since time of diagnosis. Asymptotic variances of these estimators can be calculated in a straightforward manner from the influence functions, and these variances can be estimated from actual trial data. The performance of the variance estimators is assessed via a simulated screening trial. The situation involving censored data is also discussed.     

Mixed Moment Estimation in Hermite Distribution and its Comparison with Other Estimation Procedures from Sample Size Considerations

In this paper (1) expressions (correct to n^?2 terms) for biases, variances, and covariances of the estimators a and b of Hermite distribution with probability generating function Exp[a(t–1) + b(t–1)] are obtained for two mixed moment estimates; (2) for the biases and variance-covariances, approximate regions of the parameter space (a>0, b>0) have been outlined where a sample of size 100 can be considered as “safe” in the sense that contribution of second order terms in them is 5% of that from the first order term; (3) comparison of the biases and variance-covariances of these two sets of estimators are made with those for the moment estimators, maximum likelihood estimates and the even point estimators for a sample of size 100 using the terms up to order n^?2; (4) the comparisons based on n^?2 terms in (3) have not only provided information on the estimation procedures included in the Hermite distribution, but also demonstrated the importance of higher order terms in the sampling properties of the various alternative techniques for the Hermite distribution.     

Shrinkage Pre-Test Estimator of the Intercept Parameter for a Regression Model with Multivariate Student-t Errors

In the presence of an uncertain prior information about the value of the slope parameter, the estimation of the intercept parameter of a simple regression model with a multivariate Student-t error distribution is investigated. The unrestricted, restricted and shrinkage preliminary test maximum likelihood estimators are defined. The expressions for the bias and the mean square error of the three estimators are provided and the relative efficiences are analyzed. A maximin criterion is established, and graphs are constructed for an arbitrary number of degrees of freedom (D.F.) as well as sample sizes. A criterion to select optimal significance level is also discussed.     

Discrete-time nonlinear observer-based estimators for the on-line estimation of the kinetic rates in bioreactors

Simple discrete-time estimators which allow the on-line estimation of the kinetic rates from the measurements of components' concentrations inside a bioreactor are proposed. In fact, the proposed estimators are obtained by a direct forward Euler discretization of continuous-time estimators. The design of the estimators in the continuous as well as in the discrete-time does not require or assume any model for the kinetic rates. One of the main characteristics of these estimators lies in the easiness of their calibration. We here emphasize on the performances of the discrete version of these estimators, whose stability and convergence are proved under the same conditions as in the continuous case with an additional mild assumption on the sampling time. Simulation and real-life experiments results corresponding to the discrete estimation are given. The accuracy of the obtained estimates as well as the easiness of the estimators' implementation do constitute reliable and powerful arguments for their use, in particular in adaptive control schemes.     

Non-Iterative Least Squares Estimation of Missing Values in Hyper-Graeco-Latin Square Designs

In this paper an attempt has been made to obtain explicit expressions for the estimators of the several missing values in hyper-graeco-latin square designs. Further it has been shown that the estimates of the missing values in latin square designs and graeco-latin square designs are obtained as a particular case of the estimates of the missing values in hyper-graeco-latin square designs.     

Analysis of population genetic structure with RAPD markers

Recent advances in the application of the polymerase chain reaction make it possible to score individuals at a large number of loci. The RAPD (random amplified polymorphic DNA) method is one such technique that has attracted widespread interest. The analysis of population structure with RAPD data is hampered by the lack of complete genotypic information resulting from dominance, since this enhances the sampling variance associated with single loci as well as induces bias in parameter estimation. We present estimators for several population-genetic parameters (gene and genotype frequencies, within- and between-population heterozygosities, degree of inbreeding and population subdivision, and degree of individual relatedness) along with expressions for their sampling variances. Although completely unbiased estimators do not appear to be possible with RAPDs, several steps are suggested that will insure that the bias in parameter estimates is negligible. To achieve the same degree of statistical power, on the order of 2 to 10 times more individuals need to be sampled per locus when dominant markers are relied upon, as compared to codominant (RFLP, isozyme) markers. Moreover, to avoid bias in parameter estimation, the marker alleles for most of these loci should be in relatively low frequency. Due to the need for pruning loci with low-frequency null alleles, more loci also need to be sampled with RAPDs than with more conventional markers, and some problems of bias cannot be completely eliminated.     

Determination of Number and Size of Spherical Particles from Extraction Replicas

For spherical particles randomly dispersed in the space of a specimen the estimators of the parameters of the space structure from the measurements obtained from extraction replicas are given. First an arbitrary form of the probability density function f(x) of the diameter X and then the generalized RAYLEIGH and lognormal distributions of X are considered. Unbiased estimators of the space parameters and of parameters of these distributions are found. The variances of these estimators are given and unbiased estimators of these variances are determined.     

Order-restricted inference for clustered ROC data with application to fingerprint matching accuracy

Receiver operating characteristic (ROC) curve is commonly used to evaluate and compare the accuracy of classification methods or markers. Estimating ROC curves has been an important problem in various fields including biometric recognition and diagnostic medicine. In real applications, classification markers are often developed under two or more ordered conditions, such that a natural stochastic ordering exists among the observations. Incorporating such a stochastic ordering into estimation can improve statistical efficiency (Davidov and Herman, 2012). In addition, clustered and correlated data arise when multiple measurements are gleaned from the same subject, making estimation of ROC curves complicated due to within-cluster correlations. In this article, we propose to model the ROC curve using a weighted empirical process to jointly account for the order constraint and within-cluster correlation structure. The algebraic properties of resulting summary statistics of the ROC curve such as its area and partial area are also studied. The algebraic expressions reduce to the ones by Davidov and Herman (2012) for independent observations. We derive asymptotic properties of the proposed order-restricted estimators and show that they have smaller mean-squared errors than the existing estimators. Simulation studies also demonstrate better performance of the newly proposed estimators over existing methods for finite samples. The proposed method is further exemplified with the fingerprint matching data from the National Institute of Standards and Technology Special Database 4.     

A nonlinear approach for the on-line estimation of the kinetic rates in bioreactors

The objective of this paper is to propose simple nonlinear observer-based estimators which allow the on-line estimation of kinetic rates inside a bioreactor from the measurements of components' concentrations. The design of these estimators does not require any model for the kinetic rates and the corresponding gains do not necessitate the resolution of any dynamical system. One of the main characteristics of these estimators lies in the easiness of their implementation and in particular their calibration. Their performances are firstly illustrated in simulation through a typical example. Then, estimation results corresponding to a real-life experiment are given in order to confirm the theoretical results and to demonstrate the usefulness of the proposed estimators in practical situations.