A class of semiparametric transormation frailty models with application to estimating treatment effects and heterogenity in aids clinical traials

Muitivariate failure time data are common in medical research; com¬monly used statistical models for such correlated failure-time data include frailty and marginal models. Both types of models most often assume pro¬portional hazards (Cox, 1972); but the Cox model may not fit the data well This article presents a class of linear transformation frailty models that in¬cludes, as a special case, the proportional hazards model with frailty. We then propose approximate procedures to derive the best linear unbiased es¬timates and predictors of the regression parameters and frailties. We apply the proposed methods to analyze results of a clinical trial of different dose levels of didansine (ddl) among HIV-infected patients who were intolerant of zidovudine (ZDV). These methods yield estimates of treatment effects and of frailties corresponding to patient groups defined by clinical history prior to entry into the trial.     

Randomization and blocking in the design of experiments

It is suggested that, whenever possible, an experiment be run in a completely randomized fashion. One reason for randomizing Is to protect against violations in the usual linear model assump¬tions. The protection has always been argued on qualitative grounds. This paper quantitatively demonstrates the protection by hypothesizing models in violation of the usual assumptions, mathe¬matically representing the physical act of randomization, and algebraically deriving expected mean squares, EMS, and F tests. It is shown that randomization offers considerable but not com¬plete protection against model violations.

The same methodology is also applied to blocked experiments, i.e. to experiments performed under a specific type of incomplete randomization commonly referred to as blocking. It is shown that blocking offers little protection against certain model viola¬tions. The common practice of representing blocks as a treatment factor applied to the experimental units approximates the form of the EMS derived under the violated assumptions model.     

Comparative performance of tests of normality in detecting mixtures of parallel regression lines

This paper addresses the problem of detecting a mixture of parallel regression lines when information about group member¬ship of individual cases is not given. The problem is approached as a missing variable problem, with the missing variables being the dummy variables that code for groups. If a mixture of par¬allel regression lines with normally distributed error terms is present, a simple regression model without dummy variables will produce residuals that follow approximately a mixed normal dis¬tribution. In a simulation studyr several goodness-of-fit tests of normality were used to test the residuals obtained from mis-specified models that excluded dummy variables, Factors varied in the simulation included the number and the separation of the parallel lines and the sample size, The goodness-of-fit test based on the sample kurtosis (82) was overall most powerful in detecting mixtures of parallel regression lines, Applications are discussed.     

Simple correlated arma processes

The generalization of the simple correlated autoregressive processes introduced by RISAGEE ( 1980, 1981 ) is presented in the paper. The correlation structure of the model is investigated for the purpose of its identification. Various methods of estimation are discussed. Verification of the model is based on the portmanteau statistics whose behaviour is derived using the results for the simple correlated autoregressive processes and the prin¬ciple 0f Box and PIEBCE concerning the relation of estimated residual correlations in auto-regressive and ARMA models.     

Generalized additive mixed models

Following the extension from linear mixed models to additive mixed models, extension from generalized linear mixed models to generalized additive mixed models is made, Algorithms are developed to compute the MLE's of the nonlinear effects and the covariance structures based on the penalized marginal likelihood. Convergence of the algorithms and selection of the smooth param¬eters are discussed.     

For censored and truncated data

In this paper we develop nonparametric methods for regression analysis when the response variable is subject to censoring and/or truncation. The development is based on a data completion princple that enables us to apply, via an iterative scheme, nonparametric regression techniques to iteratively com¬pleted data from a given sample with censored and/or truncated observations. In particular, locally weighted regression smoothers and additive regression models are extended to left-truncated and right-censored data Nonparamet¬ric regression analysis is applied to the Stanford heart transplant data, which have been analyzed by previous authors using semiparametric regression meth¬ods. and provides new insights into the relationship between expected survival time after a heart transplant and explanatory variables.     

Estimation of the parameters of a regression model with a multivariate t error variable

This paper deals with a regression model for several vari¬ables under the assumption that the errors have a multivariate t-distribution. The parameters of the model, the regression parameters, as well as the scale parameters and the degress of freedom of the error variable are estimated and the estimation procedure is illustrated by a numerical example, Also, the prop¬erties of the estimators and tests for the regression parameters are discussed.     

A generalization of the compound rayleigh distribution: using a bayesian method on cancer survival times

In this paper the generalized compound Rayleigh model, exhibiting flexible hazard rate, is high¬lighted. This makes it attractive for modelling survival times of patients showing characteristics of a random hazard rate. The Bayes estimators are derived for the parameters of this model and some survival time parameters from a right censored sample. This is done with respect to conjugate and discrete priors on the parameters of this model, under the squared error loss function, Varian's asymmetric linear-exponential (linex) loss function and a weighted linex loss function. The future survival time of a patient is estimated under these loss functions. A Monte Carlo simu¬lation procedure is used where closed form expressions of the estimators cannot be obtained. An example illustrates the proposed estimators for this model.     

Multiyear,through the season crop acreage estimation using estimated acreage in sample segments

Large scale crop surveys can be made frequently and inex¬pensively during a crop growing season using Landsat data. A crop's at-harvest acreage in a stratum can be estimated from the crop's estimated at-harvest acreage in a small sample of the stratum's segments. The stratum estimate can utilize Landsat imagery obtained during the current crop grow¬ing season and in previous years. A mixed effects analysis of variance model is used to generate a weighted least squares es¬timate of the stratum at-harvest acreage proportion for the cur¬rent year. Similar Landsat based stratum crop proportion esti¬mates can be combined with historical information on non-sampled (or unsuccessfully sampled) strata to provide crop acreage estimates for large regions. These regional estimates of the at-harvest acreage can be determined early in the crop growing season, at different intermediate points, and at har¬vest time     

Ridge estimator in singulah oesiun with application to age-period-cohort analysis of disease rates

Ridge estimator of a singular design is considered for linear and gener¬alized linear models. Ridge penalty helps determine a unique estimator in singmar uesign. me tuning parameter o± tue penalty is seiecteu via gener¬alized cross-validation (GCV) method. It is proven that the ridge estimator lies in a special sub-parameter space and converges to the intrinsic estimator, an estimable function in singular design, as the shrinkage penalty diminishes. The expansion of the ridge estimator and its variance are also obtained. Thismethod is demonstrated through an application to age-period-cohort (APC) analysis of the incidence rates of cervical cancer in Ontario women 1980-1994     

Estimation of parameters of the gped

Janardan (1973) introduced the generalized Polya-Eggenberger distribution as a limiting form of the generalized Markov-Polya distribution (GMPD), Ja¬nardan (1998) derived GPED formally by means of Lagrange's expansion and discussed its various properties systematically. Here, a new urn model is pro¬vided for the GPED. Moment estimators of the parameters are given in closed form. Maximum hkelihood estimators are also given. Some apphcations are provided.     

Some results on simultaneous confidence intervals for monotone contrasts in one-way anova model

The one-way ANOVA model with common variance is considered. Simultaneous confidence Intervals (SCI) for monotone contrasts in the means are derived and compared to alternative intervals gene¬rated by Williams (1977)     

The mean of the conditional,truncated multinormal distribution

An expression is derived for the mean of the conditional, truncated multinormal distribution in the general case, and for the situation where the conditioning variates have iden¬tical correlation. This equicorrelated case occurs in some models of duration and in panel data. An example is taken from economics, where a cross section of firms consider the adoption of several technological innovations. The results of the paper are used to estimate the firm-specific.     

Hypothesis tests for linear functional relationships

The analysis of linear functional relationships is considered. Expressions useful for the estimation of both unconstrained and con¬strained parameters are presented. The testing of hypotheses which consist of linear constraints on the parameters is discussed. Wald type and Wilks type test statistics and their asymptotic null dis¬tributions are derived. It appears that these test statistics are not asymptotically equivalent.     

An evaluation of the huber estimator

The asymptotic efficiency for the Huber estimator is evalu¬ated for distributions in the class of normal scale mixtures, gamma power, and Student's t-distributions. Rather high efficien¬cies are found for small deviations from normality, while intract¬able efficiencies occur for larger deviations.     

Parameter estimation based upon nonparametric function estimators

By considering the solution to a linear approximation of a nonlinear regression problem, a procedure for developing a para¬meter estimator, based upon a nonpammetric estimator of a para¬metric function, is given. The resulting estimators, which are determinable in closed form, are asymptotically normally distri¬buted and are optimal among the class of estimators based upon the function estimator. Further, in many cases, the estimator will have the same asymptotic distribution theory as the correspond¬ing maximum likelihood estimator. Estimators based upon the Kaplan-Meier quantile function are developed for randomly censored samples.     

Testing for a change lit the regression matrix of a multivariate linear model

A Bayesian test procedure Is developed to test; the null hypothesis of no change In the regression matrix of a multivariate lin¬ear model against the alternative hypothesis of exactly one change The resulting test is based on the marginal posterior distribution of the change point; To illustrate the test procedure a numerical example using a bivariate regression model is considered.     

Approximate conditional inference and the linear functional model

Approximate conditional inference is developed for the slope parameter of the linear functional model with two variables. It is shown that the model can be transformed so that the slope parameter becomes an angle and nuisance parameters are radial distances. If the nuisance parameters are known an exact confidence interval based on a location-type conditional distribution is available for the angle. More gen¬erally, confidence distributions are used to average the conditional distribution over the nuisance parameters yielding an approximate conditional confidence interval that reflects the precision indicated by the data. An example is analyzed.     

Consistent estimators of the variance-covariance matrix of the gmanova model with missing data

A common problem in multivariate general linear models is partially missing response data. The simplest method of analysis in the presence of missing data has been to delete all observations on any individual with any missing data(listwise deletion) and utilize a traditional complete data approach. However: this can result in a great loss of information: and perhaps inconsistencies in the estimation of the variance-covariance matrix. In the generalized multivariate analysis of variance(GMANOVA) model with missing data: Kleinbaum(1973) proposed an estimated generalized least squares approach. In order to apply this: however: a consistent estimate of the variance-covariance matrix is needed. Kleinbaum proposed an estimator which is unbiased and consistent: but it does not take advantage of the fact that the underlying model is GMANOVA and not MANOVA. Using the fact that the underlying model is GMANOVA we have constructed four other con¬sistent estimators. A Monte Carlo simulation experiment is conducted tto further examine how well these estimators compare to the estimator proposed by Kleinbaum.     

On the asymptotic behaviour of selector statistics

Adaptive procedures proposed by Hogg are based on selector statistics for the skewness and the tails. The asymptotic properties of several proposed selector statistics are investigated. Since ail these statistics have under some assumptions asymptoti¬cally a normal distribution, their properties depend on the asymptotic bias and variance. The reasonable concept to compare the different selector statistics is based on the selection probabil¬ities in discriminating the type of the underlying distribution. These values are numerically calculated and analyzed in detail for a number of underlying distributions.