Received 14 January 2016; accepted 16 February 2016; published 19 February 2016
The cardiovascular disease of heart disease is the number one cause of death in the World with its mortality rate representing one third of all global deaths. Its documented risk factors include smoking, physical activity, diet, obesity, cholesterol, and diabetes . As for heart attack mortality rates, research shows that despite the progress made in faster hospital care, its rates remain unchanged  .The third health outcome of interest is high blood pressure mortality rate which contributes to thirteen percent of all deaths in the US and is continually increasing  .The three health outcomes have been associated with air pollution as a consequence of oxidative stress leading to inflammation, which generates the physiological processes evolving as cardiovascular symptoms like narrowing of airways, shortness of breath, wheezing, cough, and the ability of particles to penetrate the lung wall accumulating in the pulmonary interstitium between the lung and the bloodstream  - . The outcomes have also been negatively associated with socioeconomic status (SES) through life style effects   .
The purpose of this study is to quantify the association between the three health outcomes and exposure to fine particulate matter (PM2.5) and the different aspects of SES in the US,assessing causal inferences between exposure to fine particulate matter and mortality rates of heart attacks, high blood pressure, and heart disease provides epidemiological evidence regarding the adverse health effects of air pollution. Studying this association at the county level can address the question of whether the current network of air quality monitoring stations adequately represent the populations and locations with highest risk for respiratory and cardiovascular thathas not been fully answered yet. An important underlying reason for this problem is the uniform approach that has been implemented to select monitoring sites locations, which ignores the dissimilarities among sites e.g., Houston, Texas versus Jackson, Wyoming   . In addition, the effect of local air pollution versus neighboring levels is not fully explored yet   , which points to the importance of studying community-related characteristics because of the previously-mentioned factors’ contribution to the causal inference between air pollution and its respiratory and cardiovascular health.
To achieve the study objective, data was collected from the Centers for Disease Control and Prevention (CDC)  and the US Environmental Protection Agency (EPA)  , for the year 2013. Table 1 gives the twenty four variables collected and analyzed.Statistical analyses were conducted using R version 3.1.1.
Table 1.Data collected from CDC and EPA.
Table 2 gives some of the descriptive statistics for the main variables. The high departures from skewness and kurtosis from 0 and 3, respectively, which are the characteristic values of normal distribution show non symmetric nature of heart attack mortality, high blood pressure mortality, heart disease mortality, educational attainment, poverty, unemployment rate, percentage of population aged 65 years and older, percentage of population within a half a mile of a park, annual average ambient PM2.5, and percentages of households living with severe housing problems.Figure 1 presents average mortality rates for heart attacks, heard disease and high blood pressure per type of county. Heart disease mortality rates are lowest in rural areas compared to urban and metro counties. High blood pressure mortality rates follow a somewhat opposite pattern whereas heart attack mortality rates are lowest in metro counties and highest in urban counties.As for pairwise associations (Table 3), the percentage of households living with severe housing problems was positively associated with poverty, percentage of food stamp/supplemental nutrition assistance program recipients, and percentages of Native Americans, Hispanics, and African Americans. It was negatively associated with percentage of population aged 65 years and older and represented about third of the population.Ease of access to nearby parks (Percentage of Population Living within Half a Mile of a Park) was negatively associated with both heart attack and heart disease mortalities. Heart attack mortality rate was positively associated with heart disease mortality rates, educational attainment, and poverty. High blood pressure mortality rate was positively associated with heart disease, poverty, educational attainment, and percentage of African Americans. Heard disease mortality rate was positively associated with educational attainment, poverty, and percentage of African Americans.Understandingly, percentage of population in poverty was positively associated with the percentage of food stamps recipients, median income and unemployment.
Multivariate analysis of covariance (MANOVA) tests with the three dependent variables heart attack, high blood pressure, and heart disease mortality rates using Pillai Trace, Wilk’s Lambda, and Hotelling Trace gave p-values of zeros, which reject the null hypothesis that the population has no effect on the three different health outcomes at the significance level of 0.05. Three univariate (one dependent variable) linear regressions were performed after applying the logarithmic transformation function on the independent variables of the most significant variables of Percentage without high school diploma, families with female head of household, percentage living in poverty, percentage of Asians, percentage of African Americans, percentage of Hispanics, percentage
Table 2.Some descriptive statistics.
Table 3.The most notable pairwise associations.
Figure 1.Rates of health outcomes in rural/urban counties.
of white population, percentage aged 65 years or older, percentage of population living within half a mile of a park, and annual average for PM2.5. The other variables were excluded because of their minimal association or because of their strong association with the other predictors. Table 4 presents the results of the three univariate multiple linear regressions for each of the health outcomes whereas Table 5 presents the results for the multivariate multiple linear regression. The three univariate models are statistically significant. Also, the significant predictors are the same between the three univariate and the multivariate multiple regression model. However, one of the advantages of conducting a multivariate multiple regression is the ability to test the coefficients across the different health outcome variables.In this case, ten different hypotheses were tested simultaneously: testing that the coefficient is zero for each of the ten predictors in all three equations.The residuals mean is near 0.0001 but its variance is as high as 769. Spatial distribution of the residuals is given in Figure 2. Positive residuals indicating over-prediction and negative residuals indicating under-prediction indicate the need to improve the
Table 4.Univariate multiple linear regression results.
Table 5.Multivariate multiple linear regression results.
Figure 2.Residuals of the multivariate multiple linear regression model.
model, which is supported by the low adjusted R square values.
Codispersion coefficients for the three possible combinations of dependent variables (heart attack mortality vs. high blood pressure mortality, heart attack mortality vs. heart disease mortality, and high blood pressure mortality vs. heart disease mortality) were not constant. Thus, the hypothesis of intrinsic correlation could not be accepted. This means that the correlation structure of the three dependent variables is influenced by the spatial scale, which calls for cokriging  .Figure 3 gives the direct variograms (diagonal) and cross variograms (off-diagonal) along with fitted linear model of coregionalization (straight line), with the distance maxed at 90 kilometers.The three mortality rates are predicted to keep rising with the persisting spatial concentrations (hotspots) shown in Figure 4.
Research has shown associations between air pollution and mortality of respiratory and cardiovascular systems - . Socioeconomic status (SES) has also been associated with these outcomes. Associations between heart attack, high blood pressure, and heart disease mortality rates and between them and SES have been docu-
Figure 3.Direct variograms (diagonal) and cross variograms (off-diagonal) along with fitted linear model of coregionalization (straight line).
Figure 4.Predicted heart attack mortality rate.
mented as well  - . Nevertheless, such associations are seldom investigated in the same context of exposure to air pollution. Application of multivariate geostatistical techniques can assist in this regards because it has the advantage of making use of the spatial correlations between observations pertinent to communities to predict values at unsampled places. These spatial associations can exist even with lack of pair-wise correlations, known as “intrinsic correlation”. Multivariate intrinsic correlations allow simplification of data modeling, which in turn enables us to address previously unidentified or seldom studied associations  - . Consequently, new public health findings can be sought that are relevant to each community.
“Hotspots” are geographic areas that form at specific times and places where the mortality or morbidity rate is consistently rising. This kind of clustering reflects high occurrence and frequency of disease, to the point that rates become highly predictable. However, the definition of a hotspot can be ambiguous due to lack of quantitative measures that can be incorporated or agreed upon among other issues. Nevertheless, studying the three mortality rates at the county level shows consistent spatial distribution where certain areas persistently have higher rates than the rest. Paying more attention to such hotspots can lead research to potentially stronger causal inferences  - .
This study investigated the three mortality rates of heart attack, high blood pressure, and heart disease in the United States, and their association with fine particulate matter and socioeconomic status represented by educational attainment, poverty, household income, unemployment rate, demographics, access to a nearby park, and urban-rural classification, in 2013.The multivariate geostatistical model used latitude, longitude, educational attainment, poverty, percentage Asian, percentage African American, percentage white population, and annual PM2.5 concentration as the mean covariates.These variables were also significant in univariate and multivariate multiple linear regressions. This study had the limitations of not exploring temporal variations, inherent issues with the data, non-consideration of other air pollutants and demographic factors like gender and age. Nevertheless, it pointed to significant findings including the county level which showed consistent spatial distribution where certain areas persistently had higher rates than the rest. Paying more attention to such hotspots could lead research to potentially stronger causal inferences.
 Go, A.S., Mozaffarian, D., Roger, V.L., Benjamin, E.J., Berry, J.D., Borden, W.B., Bravata, D.M., Dai, S., Ford, E.S., Fox, C.S., Franco, S., Fullerton, H.J., Gillespie, C., Hailpern, S.M., Heit, J.A., Howard, V.J., Huffman, M.D., Kissela, B.M., Kittner, S.J., Lackland, D.T., Lichtman, J.H., Lisabeth, L.D., Magid, D., Marcus, G.M., Marelli, A., Matchar, D.B., McGuire, D.K., Mohler, E.R., Moy, C.S., Mussolino, M.E., Nichol, G., Paynter, N.P., Schreiner, P.J., Sorlie, P.D., Stein, J., Turan, T.N., Virani, S.S., Wong, N.D., Woo, D. and Turner, M.B. (2013) Heart Disease and Stroke Statistics-2013 Update: A Report from the American Heart Association. Circulation, 127, e6-e245. http://dx.doi.org/10.1161/CIR.0b013e31828124ad
 Arbor, A. (2015) Heart Attack Mortality Rates Unchanged in Spite of Faster Care at Hospitals|University of Michigan Health System.
 Donaldson, K., Mills, N., MacNee, W., Robinson, S. and Newby, D. (2005) Role of Inflammation in Cardiopulmonary Health Effects of PM. Toxicology and Applied Pharmacology, 207, 483-488.
 Mills, N.L., Miller, M.R., Lucking, A.J., Beveridge, J., Flint, L., Boere, A.J.F., Fokkens, P.H., Boon, N.A., Sandstrom, T., Blomberg, A., Duffin, R., Donaldson, K., Hadoke, P.W.F., Cassee, F.R. and Newby, D.E. (2011) Combustion-Derived Nanoparticulate Induces the Adverse Vascular Effects of Diesel Exhaust Inhalation. European Heart Journal, 32, 2660-2671. http://dx.doi.org/10.1093/eurheartj/ehr195
 Zigler, C.M. and Dominici, F. (2014) Point: Clarifying Policy Evidence with Potential-Outcomes Thinking—Beyond Exposure-Response Estimation in Air Pollution Epidemiology. American Journal of Epidemiology, 180, 1133-1140. http://dx.doi.org/10.1093/aje/kwu263
 Vanderweele, T.J. and Vansteelandt, S. (2009) Conceptual Issues Concerning Mediation, Interventions and Composition. Statistics and Its Interface, 2, 457-468. http://dx.doi.org/10.4310/SII.2009.v2.n4.a7
 Atkinson, R.W., Bremner, S.A., Anderson, H.R., Strachan, D.P., Bland, J.M. and de Leon, A.P. (1999) Short-Term Associations between Emergency Hospital Admissions for Respiratory and Cardiovascular Disease and Outdoor Air Pollution in London. Archives of Environmental Health, 54, 398-411.
 Wong, T.W., Tam, W.S., Yu, T.S. and Wong, A.H.S. (2002) Associations between Daily Mortalities from Respiratory and Cardiovascular Diseases and Air Pollution in Hong Kong, China. Occupational and Environmental Medicine, 59, 30-35. http://dx.doi.org/10.1136/oem.59.1.30
 Crapo, J.D., Broaddus, V.C., Brody, A.R., Malindzak, G., Samet, J. and Wright, J.R. (2003) Workshop on Lung Disease and the Environment: Where Do We Go from Here? American Journal of Respiratory and Critical Care Medicine, 168, 250-254. http://dx.doi.org/10.1164/rccm.168.2.250
 Jerrett, M., Gale, S. and Kontgis, C. (2010) Spatial Modeling in Environmental and Public Health Research. International Journal of Environmental Research and Public Health, 7, 1302-1329.
 Glass, T.A., Goodman, S.N., Hernán, M.A. and Samet, J.M. (2013) Causal Inference in Public Health. Annual Review of Public Health, 34, 61-75. http://dx.doi.org/10.1146/annurev-publhealth-031811-124606
 Noffsinger, M.A., Pfefferbaum, B., Pfefferbaum, R.L., Sherrib, K. and Norris, F.H. (2012) The Burden of Disaster: Part I. Challenges and Opportunities within a Child’s Social Ecology. International Journal of Emergency Mental Health, 14, 3-13.
 Anderson, F., Carson, A., Whitehead, L. and Burau, K. (2014) Spatiotemporal Analysis of the Effect of Ozone and Fine Particulate on CVD Emergency Room Visits in Harris County, Texas. Open Journal of Air Pollution, 3, 87-99. http://dx.doi.org/10.4236/ojap.2014.34009
 Goovaerts, P. (2006) Geostatistical Analysis of Disease Data: Accounting for Spatial Support and Population Density in the Isopleth Mapping of Cancer Mortality Risk Using Area-to-Point Poisson Kriging. International Journal of Health Geographics, 5, 52. http://dx.doi.org/10.1186/1476-072X-5-52
 Kulldorff, M., Tango, T. and Park, P.J. (2003) Power Comparisons for Disease Clustering Tests. Computational Statistics & Data Analysis, 42, 665-684.