CENSOR procedure

Pre-processes censored data before analysis by ANOVA (P.W. Lane).


Options

PRINT = string
Whether to monitor convergence (monitor); default * implies no monitoring

TERM = formula
Formula for lowest stratum residual term; no default - this option must be set

DESIGN = pointer
Identifier specifying design information for ANOVA, or to save design information; default *

MAXCYCLE = scalar
Maximum number of iterations; default 20


Parameters

Y = variates
Observed variate with censored values represented by values greater than or equal to the bound; no default - this parameter must be set

BOUND = scalars or variates
Upper bound for censoring for each unit; no default - this parameter must be set

DF = scalars
Estimated residual d.f. for lowest stratum, adjusting for censoring; default *


Description

An observation is said to be censored if it is known only that it is less than (or greater than) a particular value. Such observations can occur in designed experiments when the observed variable is the time until some event takes place for each experimental unit. For example, if the observed measurement is the lifetime of electric light-bulbs, it may happen that some bulbs are still alight when the experiment has to be concluded.

   The response variate should be specified using the Y parameter, representing the censored values as values that are greater than or equal to the censoring bound. The bound is specified using the BOUND parameter, either as a scalar - if the bound is constant over the experiment - or as a variate of the same length as the response variate. Missing values in Y will be treated as usual, not as censored values. The procedure deals with the case of censoring with an upper bound. For a problem involving a lower bound, the structures in Y and BOUND should be multiplied by -1 before using the procedure (and the analysis interpreted accordingly).

   The results for any experiment analysable by the ANOVA directive may be processed by the procedure. You must give BLOCKSTRUCTURE, TREATMENTSTRUCTURE and COVARIATE statements, as relevant, before using the procedure. If the analysis of the experiment requires a setting of the WEIGHTS, FACTORIAL, CONTRASTS or DEVIATIONS options of the ANOVA directive, you should give an ANOVA statement with these settings before using the procedure, setting the DESIGN option and then using the same identifier in the DESIGN option of CENSOR. The lowest stratum of the experiment must be identified explicitly in the BLOCKSTRUCTURE statement, rather than being implicitly taken as the *units* stratum by ANOVA; the model term representing this stratum must be specified using the TERM option of CENSOR. For example, a split-plot experiment with blocks might be specified by

BLOCKS block/plot/subplot

CENSOR [TERM=block.plot.subplot; ...

   If you set the option PRINT=monitor, the procedure will print the values of the standard error of the lowest stratum at each cycle of the iterative estimation process. The maximum number of oterations is specified by the MAXCYCLE option, with a default of 20. After using the procedure, the censored values in the variate in the Y parameter will have been replaced by estimated values. The analysis of this variate can be displayed with ADISPLAY, or results saved with AKEEP. The analysis will not be exact, due to the estimation of the censored values. However, the residual degrees of freedom in the lowest stratum may be corrected to the value output by the DF parameter (this is the absolute correction referred to by Taylor 1973; see the Method section).


Options: PRINT, TERM, DESIGN, MAXCYCLE. Parameters: Y, BOUND, DF.


Method

The censored units in the observed variate are replaced by estimated values, using the method outlined by Taylor (1973). This method estimates the expected value of each censored observation iteratively conditional on the fact that the value must be greater than the fixed bound, and using the relevant information from the other observations in the experiment. The iterative process is deemed to have converged when the relative change in each fitted value, and in the stratum standard error, is less than 0.1%. If convergence is not reached within the number of iteractions specified by the MAXCYCLE option, a message will be printed and the process will terminate. This should occur only if there is a high proportion of censored values, or if the design affords no information about one or more of the censored values.


Action with RESTRICT

The Y variate can be restricted, but the BOUNDS variate must not be.


Reference

Taylor, J. (1973). The analysis of designed experiments with censored observations. Biometrics, 29, 35-43.