Bayesian non‐response models for categorical data from small areas: an application to BMD and age

Abstract

We provide a Bayesian analysis of data categorized into two levels of age (younger than 50 years, at least 50 years) and three levels of bone mineral density (normal, osteopenia, osteoporosis) for white females at least 20 years old in the third National Health and Nutrition Examination Survey. For the sample, the age of each individual is known, but some individuals did not have their BMD measured. We use two types of models: In the ignorable non-response models the propensity to respond does not depend on BMD and age of an individual, while in the non-ignorable non-response models it does. These are the baseline models which are used to derive all models for testing. Our non-ignorable non-response models are 'close' to the ignorable non-response models, thereby reducing the effects of the assumptions about non-respondents that cannot be tested in non-response models. We have data from 35 counties, small areas, and therefore our models are hierarchical, a feature that allows a 'borrowing of strength' across the counties, and they provide a substantial reduction in variation. The non-ignorable non-response models are generalizations of the ignorable non-response models, and therefore, the non-ignorable non-response models allow broader inference. The joint posterior density of the parameters for each model is complex, and therefore, we fit each model using Markov chain Monte Carlo methods to obtain samples which are used to make inference about BMD and age. For each county we can estimate the proportion of individuals in each BMD and age cell of the categorical table, and we can assess the relation between BMD and age using the Bayes factor. A sensitivity analysis shows that there are differences (typically small) in inference that permits different levels of association between BMD and age. A simulation study shows that there is not much difference between the baseline ignorable and non-ignorable non-response models.