The function may be called by the names: g03dcc or nag_mv_discrim_group.
Discriminant analysis is concerned with the allocation of observations to groups using information from other observations whose group membership is known, ; these are called the training set. Consider variables observed on populations or groups. Let be the sample mean and the within-group variance-covariance matrix for the th group; these are calculated from a training set of observations with observations in the th group, and let be the th observation from the set of observations to be allocated to the groups. The observation can be allocated to a group according to a selected rule. The allocation rule or discriminant function will be based on the distance of the observation from an estimate of the location of the groups, usually the group means. A measure of the distance of the observation from the th group mean is given by the Mahalanobis distance, :
If the pooled estimate of the variance-covariance matrix is used rather than the within-group variance-covariance matrices, then the distance is:
Instead of using the variance-covariance matrices and , g03dcc uses the upper triangular matrices and supplied by g03dac such that and . can then be calculated as where or as appropriate.
In addition to the distances, a set of prior probabilities of group membership, , for , may be used, with . The prior probabilities reflect your view as to the likelihood of the observations coming from the different groups. Two common cases for prior probabilities are , that is, equal prior probabilities, and , for , that is, prior probabilities proportional to the number of observations in the groups in the training set.
g03dcc uses one of four allocation rules. In all four rules the variables are assumed to follow a multivariate Normal distribution with mean and variance-covariance matrix if the observation comes from the th group. The different rules depend on whether or not the within-group variance-covariance matrices are assumed equal, i.e., , and whether a predictive or estimative approach is used. If is the probability of observing the observation from group , then the posterior probability of belonging to group is:
In the estimative approach, the arguments and in (3) are replaced by their estimates calculated from . In the predictive approach, a non-informative prior distribution is used for the arguments and a posterior distribution for the arguments, , is found. A predictive distribution is then obtained by integrating over the argument space. This predictive distribution then replaces in (3). See Aitchison and Dunsmore (1975), Aitchison et al. (1977) and Moran and Murphy (1979) for further details.
The observation is allocated to the group with the highest posterior probability. Denoting the posterior probabilities, , by , the four allocation rules are:
(i)Estimative with equal variance-covariance matrices – Linear Discrimination.
(ii)Estimative with unequal variance-covariance matrices – Quadratic Discrimination.
(iii)Predictive with equal variance-covariance matrices.
(iv)Predictive with unequal variance-covariance matrices
In the above the appropriate value of from (1) or (2) is used. The values of the are standardized so that,
Moran and Murphy (1979) show the similarity between the predictive methods and methods based upon likelihood ratio tests.
In addition to allocating the observation to a group, g03dcc computes an atypicality index, . This represents the probability of obtaining an observation more typical of group than the observed (see Aitchison and Dunsmore (1975) and Aitchison et al. (1977)). The atypicality index is computed as:
where is the lower tail probability from a beta distribution where, for unequal within-group variance-covariance matrices,
and for equal within-group variance-covariance matrices,
If is close to 1 for all groups it indicates that the observation may come from a grouping not represented in the training set. Moran and Murphy (1979) provide a frequentist interpretation of .
Aitchison J and Dunsmore I R (1975) Statistical Prediction Analysis Cambridge
Aitchison J, Habbema J D F and Kay J W (1977) A critical comparison of two methods of statistical discrimination Appl. Statist.26 15–25
Kendall M G and Stuart A (1976) The Advanced Theory of Statistics (Volume 3) (3rd Edition) Griffin
Krzanowski W J (1990) Principles of Multivariate Analysis Oxford University Press
Moran M A and Murphy B J (1979) A closer look at two alternative methods of statistical discrimination Appl. Statist.28 223–232
Morrison D F (1967) Multivariate Statistical Methods McGraw–Hill
1: – Nag_DiscrimMethodInput
On entry: indicates whether the estimative or predictive approach is to be used.
The estimative approach is used.
The predictive approach is used.
2: – Nag_GroupCovarsInput
On entry: indicates whether or not the within-group variance-covariance matrices are assumed to be equal and the pooled variance-covariance matrix used.
The within-group variance-covariance matrices are assumed equal and the matrix stored in the first elements of gc is used.
The within-group variance-covariance matrices are assumed to be unequal and the matrices , for , stored in the remainder of gc are used.
3: – Nag_PriorProbabilityInput
On entry: indicates the form of the prior probabilities to be used.
Equal prior probabilities are used.
Prior probabilities proportional to the group sizes in the training set, , are used.
An internal error has occurred in this function. Check the function call
and any array sizes. If the call is correct then please contact NAG for
On entry, .
Constraint: must be within machine precision of 1 when .
On entry, .
Constraint: , when .
The number of variables, nvar in the analysis , while number of variables included in the analysis via array .
Constraint: these two numbers must be the same.
The accuracy of the returned posterior probabilities will depend on the accuracy of the input or matrices. The atypicality index should be accurate to four significant places.
8Parallelism and Performance
Background information to multithreading can be found in the Multithreading documentation.
g03dcc is not threaded in any implementation.
The distances can be computed using g03dbc if other forms of discrimination are required.
The data, taken from Aitchison and Dunsmore (1975), is concerned with the diagnosis of three ‘types’ of Cushing's syndrome. The variables are the logarithms of the urinary excretion rates (mg/24hr) of two steroid metabolites. Observations for a total of 21 patients are input and the group means and matrices are computed by g03dac. A further six observations of unknown type are input and allocations made using the predictive approach and under the assumption that the within-group covariance matrices are not equal. The posterior probabilities of group membership, , and the atypicality index are printed along with the allocated group. The atypicality index shows that observations 5 and 6 do not seem to be typical of the three types present in the initial 21 observations.