ESTIMATING POPULATION SIZE FROM CAPTURE-RECAPTURE EXPERIMENTS IN CONTINUOUS TIME

Estimates for the size of a closed population are given for multiple recapture studies in continuous time. The estimates are derived by a method of moments for martingales. An estimate and associated standard error of the population size are derived for a homogeneous population when the capture rates are permitted to depend on time in an unspecified manner. Corresponding results are obtained when the capture rates vary among individuals as well. Explicit expressions are given for these estimates and standard errors which involve only simple computation.     

a bayesian analysis of the multivariate mixed-linear model

The Bayesian analysis of the multivariate mixed linear model is considered. The exact posterior distribution for the fixed effects matrix and the error covariance matrix are obtained. The exact posterior means and variances of the Bayesian estimators for the covariance matrices of random effects are also derived. These posterior moments are computed without constrained optimization and numerical integration. The calculations are feasible for arbitrary models. Reasonable approximations for the posterior distributions for the covariance matrices associated with the random effects are obtained also. Results are illustrated with a numerical example.     

Measuring non-linear dependence for two random variables distributed along a curve

We propose new dependence measures for two real random variables not necessarily linearly related. Covariance and linear correlation are expressed in terms of principal components and are generalized for variables distributed along a curve. Properties of these measures are discussed. The new measures are estimated using principal curves and are computed for simulated and real data sets. Finally, we present several statistical applications for the new dependence measures.     

Bayesian credibility limits for small finite population proportions:stratified and cluster sampling

Bayes credibility limits for small proportions from stratified and fixed size cluster samples are discussed. Ericson’s (JRSS B (1969)) Beta Binomial and Dirichlet-Multinomial priors are used. Approximate limits that are appropriate for large samples and small proportions are derived in both cases. These allow asymptotic comparisons of the efficacy of stratified and cluster sampling relative to simple random sampling for estimating small proportions. Procedures for the selection of hyper parameters are also presented.     

Circular Binary Block Second- and Higher-Order Neighbor Designs

Some new neighbor designs are presented here. Second-order neighbor designs for different configurations are generated in circular binary blocks. Third-order and fourth-order neighbor designs for some cases are also constructed. In all cases, circular blocks are well separated and these designs are obtained through initial block/s. At the end of the study, some models for analysis of these designs are also presented.     

A new way of quantifying the symmetry of a random variable: Estimation and hypothesis testing

New measures of skewness for real-valued random variables are proposed. The measures are based on a functional representation of real-valued random variables. Specifically, the expected value of the transformed random variable can be used to characterize the distribution of the original variable. Firstly, estimators of the proposed skewness measures are analyzed. Secondly, asymptotic tests for symmetry are developed. The tests are consistent for both discrete and continuous distributions. Bootstrap versions improving the empirical results for moderated and small samples are provided. Some simulations illustrate the performance of the tests in comparison to other methods. The results show that our procedures are competitive and have some practical advantages.     

Goodness—of—fit for the two—parameter weibull distribution with estimated parameters

Goodness—of—fit statistics based on the empirical distribution function (EDF) are not distribution—free when parameters for the hypothesized distribution are estimated. Tables are percentile values of several EDF statistics are available for the two—parameter Weibull distribution when parameters are estimated by maximum likelihood. To determine how these tabled values change when simpler estimators are employed, percentile scores for EDF goodness—of—fit tests were obtained by Monte—Carlo simulation for maximum likelihood estimators (MLEs), good linear unbiased estimators (GLUEs), and modified Cramer—von Mises, Anderson—Darling, and Watson statistics are presented for GLUEs for both complete and censored samples. Critical values for Kolmogorov—Smirnov statistics were less affected by the method of estimation than were closer for MLEs and MGLUEs than for MGLUEs and GLUEs. On the other hand, MGLUE and GLUE results were much more similar to each other than to the MLE results when censoring was light and sample sizes were large.     

Asymptotic efficiencies of estiamtors of location a comutational study

Asymptotic efficiencies of four classes of estimators of location are evaluated for a family of distributions consisting of t, lambda and contaminated normal densities. For the class of estimators derived from signed rank tests, maximin efficiencies between pairs of distributions in the family are computed using a formula of Gastwirth ( 1966 ). Asymptotic efficiencies are also evaluated for the scale dependent estimators of the form proposed by Hubcr ( 1964 ) and the efficiencies of procedures utilizing interquantiie ranges.are evaluated. Efficiencies of linear estimators such as trimmed means, BLUE's for the lambda family are computed for each density considered. Efficiencies of linear, polynomial and trigonometric approximations to BLUE weight functions are determined. Using the method of Birnbaum and Laska ( 1967 ) maximin efficiencies are computed using four linear or polynomial terms. On the basis of comparisons of these numerical results, suggestions for robust estimators are given     

Generalized integer-valued random coefficient for a first order structure autoregressive (RCINAR) process

A random coefficient autoregressive process for count data based on a generalized thinning operator is presented. Existence and weak stationarity conditions for these models are established. For the particular case of the (generalized) binomial thinning, it is proved that the necessary and sufficient conditions for weak stationarity are the same as those for continuous-valued AR(1) processes. These kinds of processes are appropriate for modelling non-linear integer-valued time series. They allow for over-dispersion and are appropriate when including covariates. Model parameters estimators are calculated and their properties studied analytically and/or through simulation.     

Kolmogorov-smirnov confidence sets: the case of more than one sample

For the two-sample problem, distribution-free confidence sets for the shift parameter when the scale parameters are equal and for both the shift and the ratio of scale parameters are derived. Multiple comparisons for the k sample location problem are constructed when all scale parameters are equal. Examples are given. Procedures may be completed with only pencil and paper.

     

Exact confidence regions for scale and location parameters based on sample quartiles

In this paper we present relatively simple (ruler, paper, and pencil) nonparametric procedures for constructing joint confidence regions for (i) the median and the inner quartile range for the symmetric one-sample problem and (ii) the shift and ratio of scale parameters for the two-sample case. Both procedures are functions of the sample quartiles and have exact confidence levels when the populations are continuous. The one-sample case requires symmetry of first and third quartiles about the median.

The confidence regions we propose are always convex, nested for decreasing confidence levels and are compact for reasonably large sample sizes. Both exact small sample and approximate large sample distributions are given.     

One-sided normal prediction intervals and sample sizes for a third-stage sample

One-sided two-stage prediction intervals for a normal population are extended to a third sampling stage. Procedures and tables are given for two situations. In the first situation, methods for obtaining such intervals are presented, and tables for calculating such prediction intervals are provided. In the second situation, a two-stage prediction interval has been applied, and a third stage is now required. Sample sizes are given for the third stage.     

A-Optimal Block Designs with Additional Singly Replicated Treatments

Block designs to which have been added a number of singly-replicated treatments, known as secondary treatments, are particularly useful for experiments where only small amounts of material are available for some treatments, for example new plant varieties. The designs are of particular use in the microarray situation. Such designs are known as 'augmented designs'. This paper obtains the properties of these designs and shows that, with an equal number of secondary treatments in each block, the A-optimal design is obtained by using the A-optimal design for the original block design. It develops formulae for the variance of treatment comparisons, for both the primary and the secondary treatments. A number of examples are used to illustrate the results.     

Asymptotic Inference on Cointegrating Rank in Partial Systems

The likelihood ratio test for cointegrating rank is analyzed for partial (or conditional) systems in the vector autoregressive error-correction model. Under the assumption of weak exogeneity for the cointegrating parameters, the asymptotic distributions are given and tables of critical values are provided. A discussion is given of some of the assumptions of the model, why they are needed, and how they are tested.     

Non-homogeneous Markov models in the analysis of survival after breast cancer

A cohort of 300 women with breast cancer who were submitted for surgery is analysed by using a non-homogeneous Markov process. Three states are onsidered: no relapse, relapse and death. As relapse times change over time, we have extended previous approaches for a time homogeneous model to a non omogeneous multistate process. The trends of the hazard rate functions of transitions between states increase and then decrease, showing that a changepoint can be considered. Piecewise Weibull distributions are introduced as transition intensity functions. Covariates corresponding to treatments are incorporated in the model multiplicatively via these functions. The likelihood function is built for a general model with k changepoints and applied to the data set, the parameters are estimated and life-table and transition probabilities for treatments in different periods of time are given. The survival probability functions for different treatments are plotted and compared with the corresponding function for the homogeneous model. The survival functions for the various cohorts submitted for treatment are fitted to the mpirical survival functions.     

Approximate one-sided tolerance limits in random effects model and in some mixed models and comparisons

An approximate closed-form one-sided tolerance limit (TL) in a general mixed model is proposed. One-sided TLs for the distribution of observable random variable and for the distribution of unobservable random variable in one-way random model are obtained as special cases from the one for the general mixed model. Applications to a two-way nested random model are also given. The merits of the TLs are evaluated using Monte Carlo simulation and compared with the existing ones. Our comparison studies indicate that the approximate TLs are quite satisfactory for all parameter and sample size configurations, and better than the existing ones in some cases. Approximate confidence intervals for exceedance probabilities in one-way random effects model are also proposed. The procedures are illustrated using three examples.     

Weighing the Evidence for Hypotheses with Small Samples of Right-censored Exponential Data

The p-value evidence for an alternative to a null hypothesis regarding the mean lifetime can be unreliable if based on asymptotic approximations when there is only a small sample of right-censored exponential data. However, a guarded weight of evidence for the alternative can always be obtained without approximation, no matter how small the sample, and has some other advantages over p-values. Weights of evidence are defined as estimators of 0 when the null hypothesis is true and 1 when the alternative is true, and they are judged on the basis of the ensuing risks, where risk is mean squared error of estimation. The evidence is guarded in that a preassigned bound is placed on the risk under the hypothesis. Practical suggestions are given for choosing the bound and for interpreting the magnitude of the weight of evidence. Acceptability profiles are obtained by inversion of a family of guarded weights of evidence for two-sided alternatives to point hypotheses, just as confidence intervals are obtained from tests; these profiles are arguably more informative than confidence intervals, and are easily determined for any level and any sample size, however small. They can help understand the effects of different amounts of censoring. They are found for several small size data sets, including a sample of size 12 for post-operative cancer patients. Both singly Type I and Type II censored examples are included. An examination of the risk functions of these guarded weights of evidence suggests that if the censoring time is of the same magnitude as the mean lifetime, or larger, then the risks in using a guarded weight of evidence based on a likelihood ratio are not much larger than they would be if the parameter were known.     

Sparse group lasso for multiclass functional logistic regression models

Sparsity-inducing penalties are useful tools for variable selection and are also effective for regression problems where the data are functions. We consider the problem of selecting not only variables but also decision boundaries in multiclass logistic regression models for functional data, using sparse regularization. The parameters of the functional logistic regression model are estimated in the framework of the penalized likelihood method with the sparse group lasso-type penalty, and then tuning parameters for the model are selected using the model selection criterion. The effectiveness of the proposed method is investigated through simulation studies and the analysis of a gene expression data set.     

Tolerance intervals for symmetric location-scale families based on uncensored or censored samples

Exact methods for constructing two-sided tolerance intervals (TIs) and tolerance intervals that control percentages in both tails for a location-scale family of distributions are proposed. The proposed methods are illustrated by constructing TIs for a normal, logistic, and Laplace (double exponential) distributions based on type II singly censored samples. Factors for constructing one-sided and two-sided TIs for a logistic distribution are tabulated for the case of uncensored samples. Factors for constructing TIs based on censored samples for all three distributions are also tabulated. The factors for all cases are estimated by Monte Carlo simulation. An adjustment to the tolerance factors based on type II censored samples is proposed so that they can be used to find approximate TIs based on type I censored samples. Coverage studies of the approximate TIs based on type I censored samples indicate that the approximation is satisfactory as long as the proportion of censored observations is no more than 0.70. The methods are illustrated using some practical examples.     

Modified mantel-haenszel procedures for matched pairs

In the prospective study of a finely stratified population, one individual from each stratum is chosen at random for the “treatment” group and one for the “non-treatment” group. For each individual the probability of failure is a logistic function of parameters designating the stratum, the treatment and a covariate. Uniformly most powerful unbiased tests for the treatment effect are given. These tests are generally cumbersome but, if the covariate is dichotomous, the tests and confidence intervals are simple. Readily usable (but non-optimal) tests are also proposed for poly-tomous covariates and factorial designs. These are then adapted to retrospective studies (in which one “success” and one “failure” per stratum are sampled). Tests for retrospective studies with a continuous “treatment” score are also proposed.