hmd.DA {LifeTables} | R Documentation |
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. |
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) |
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.
Value
train |
The output from the function |
out.dens |
The output from the function |
classification |
The classification for the test data |
Author(s)
David Sharrow
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
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