LifeTables (version 1.0)

hmd.DA: Discriminant analysis to determine life table family

Description

Returns the results of a discriminant analysis to find the appropriate life table family based on either a single indicator of child mortality (1m0, 5m0, 1q0, 5q0) or one of the four child indicators and 45q15

Usage

hmd.DA(x, child.mort=4, sex="female", adult.mort = NULL)

Arguments

x

A numeric value for the child mortality indicator to be classified expressed as a proportion (e.g. x=0.05)

sex

"male" or "female" indicates the sex for the indicators

child.mort

An integer (1-4) to indicate which child mortality indicator is being supplied (1 - 1m0; 2 - 5m0; 3 - 1q0; 4 - 5q0)

adult.mort

The value for 45q15 (if no value is supplied, defaults to NULL and uses just the single child mortality measure to classify)

Value

• train
  The output from the function MclustDA, the parameters and other summary information for the model best fitting each class according to BIC.
• out.dens
  The output from the function predict(MclustDA model name), A matrix in which the [i,j]th entry is the density for observation i in the model for class j.
• classification
  The classification for the test data

Details

The training models based on the known classification of the training data can be found in the workspace MLTobs in the "data" sub-directory of this package.

References

S. Clark and D. Sharrow (2011) "Contemporary Model Life Tables for Developed Countries: An Application of Model-Based Clustering." Paper presented at the annual meeting of the Population Association of America. http://paa2011.princeton.edu/abstracts/111808

C. Fraley and A. E. Raftery (2002) "Model-based clustering, discriminant analysis, and density estimation." Journal of the American Statistical Association 97:611-631.

C. Fraley, A. E. Raftery, T. B. Murphy, L. Scrucca (2012) "mclust Version 4 for R: Normal Mixture Modeling for Model-Based Clustering, Classification, and Density Estimation", Technical Report no. 597, Department of Statistics, University of Washington

Human Mortality Database University of California, Berkeley (USA), and Max Planck Institute for Demographic Research (Germany). Available at www.mortality.org or www.humanmortality.de (data downloaded February 2015).

See Also

MclustDA, Mclust, MLTobs

Examples

# classifying with 1m0
examp.out <- hmd.DA(x=.05, sex="male", child.mort=1)
examp.out$classification

# classifying with 5q0 and 45q15 
examp.out2 <- hmd.DA(x=.03, sex="male", child.mort=4, adult.mort=.18)
examp.out2$classification