Multivariate methods for correlated effects sizes

Current statistical methods for estimation of parametric effect sizes from a series of experiments are generally restricted to univariate comparisons of standardized mean differences between two treatments. Multivariate methods are presented for the case in which effect size is a vector of standardized multivariate mean differences and the number of treatment groups is two or more. The proposed methods employ a vector of independent sample means for each response variable that leads to a covariance structure which depends only on correlations among the $p$ responses on each subject. Using weighted least squares theory and the assumption that the observations are from normally distributed populations, multivariate hypotheses analogous to common hypotheses used for testing effect sizes were formulated and tested for treatment effects which are correlated through a common control group, through multiple response variables observed on each subject, or both conditions.^ The asymptotic multivariate distribution for correlated effect sizes is obtained by extending univariate methods for estimating effect sizes which are correlated through common control groups. The joint distribution of vectors of effect sizes (from $p$ responses on each subject) from one treatment and one control group and from several treatment groups sharing a common control group are derived. Methods are given for estimation of linear combinations of effect sizes when certain homogeneity conditions are met, and for estimation of vectors of effect sizes and confidence intervals from $p$ responses on each subject. Computational illustrations are provided using data from studies of effects of electric field exposure on small laboratory animals. ^

Airliner cabin air quality exposure assessment

The airliner cabin environment and its effects on occupant health have not been fully characterized. This dissertation is: (1) A review of airliner environmental control systems (ECSs) that modulate the ventilation, temperature, relative humidity (RH), and barometric pressure (PB) of the cabin environment---variables related to occupant comfort and health. (2) A review and assessment of the methods and findings of key cabin air quality (CAQ) investigations. Several significant deficiencies impede the drawing of inferences about CAQ, e.g., lack of detail about investigative methods, differences in methods between investigations, limited assessment of CAQ variables, small sample sizes, and technological deficiencies of data collection. (3) A comprehensive evaluation of the methods used in the subsequent NIOSH-FAA Airliner CAQ Exposure Assessment Feasibility Study (STUDY) in which this author participated. A number of problems were identified which limit the usefulness of the data. (4) An analysis of the reliable 10-flight STUDY data. Univariate and multivariate methods applied to CO2 (a surrogate for air contaminants), temperature, RH, and PB, in association with percent passenger load, ventilation system, flight duration, airliner body type, and measurement location within the cabin, revealed neither the measured values nor their variability exceeded established health-based exposure limits. Regression analyses suggest CO2, temperature, and RH were affected by percent passenger load. In-flight measurements of CO2 and RH were relatively independent of ventilation system type or flight duration. Cabin temperature was associated with percent passenger load, ventilation system type, and flight duration. (5) A synthesis of the implications of the airliner ECS and cabin O2 environment on occupant health. A model was developed to predict consequences of the airliner cabin pressure altitude 8,000 ft limit and resulting model-estimated PO2 on cardiopulmonary status. Based on the PB, altitude, and environmental data derived from the 10 STUDY flights, the predicted PaO2 of adults with COPD, or elderly adults with or without COPD, breathing ambient cabin air could be < 55 mm Hg (SaO2 < 88%). Reduction in cabin PB found in the STUDY flights could aggravate various medical conditions and require the use of in-flight supplemental O2. ^

Atlas of cancer mortality in the United Arab Emirates aided by application of geographic information system methods

The aim of this study was to determine cancer mortality rates for the United Arab Emirates (UAE) and to create an atlas of cancer mortality for the UAE. This atlas is the first of its kind in the Gulf country and the Middle East. Death certificates were reviewed for a period from January 1, 1990 to December 31, 1999 and cancer deaths were identified. Cancer mortality cases were verified by comparing with medical records. Age-adjusted cancer mortality rates were calculated by gender, emirate/medical district and nationality (UAE nationals and overall UAE population). Individual rates for each emirate were compared to the overall rate of the corresponding population for the same cancer site and gender. Age-adjusted rates were mapped using MapInfo software. High rates for liver, lung and stomach cancer were observed in Abu Dhabi, Dubai and the northern emirates, respectively. Rates for UAE nationals were greater compared to the overall UAE population. Several factors were suggested that may account for high rates of specific cancers observed in certain emirates. It is hoped that this atlas will provide leads that will guide further epidemiologic and public health activities aimed at preventing cancer. ^

A multivariate frailty model for disease recurrences and survival

A multivariate frailty hazard model is developed for joint-modeling of three correlated time-to-event outcomes: (1) local recurrence, (2) distant recurrence, and (3) overall survival. The term frailty is introduced to model population heterogeneity. The dependence is modeled by conditioning on a shared frailty that is included in the three hazard functions. Independent variables can be included in the model as covariates. The Markov chain Monte Carlo methods are used to estimate the posterior distributions of model parameters. The algorithm used in present application is the hybrid Metropolis-Hastings algorithm, which simultaneously updates all parameters with evaluations of gradient of log posterior density. The performance of this approach is examined based on simulation studies using Exponential and Weibull distributions. We apply the proposed methods to a study of patients with soft tissue sarcoma, which motivated this research. Our results indicate that patients with chemotherapy had better overall survival with hazard ratio of 0.242 (95% CI: 0.094 - 0.564) and lower risk of distant recurrence with hazard ratio of 0.636 (95% CI: 0.487 - 0.860), but not significantly better in local recurrence with hazard ratio of 0.799 (95% CI: 0.575 - 1.054). The advantages and limitations of the proposed models, and future research directions are discussed. ^

Characterization of the inflammatory cells in ascending thoracic aortic aneurysms in patients with Marfan syndrome, familial thoracic aortic aneurysms, and sporadic aneurysms.

OBJECTIVE: This study sought to characterize the inflammatory infiltrate in ascending thoracic aortic aneurysm in patients with Marfan syndrome, familial thoracic aortic aneurysm, or nonfamilial thoracic aortic aneurysm. BACKGROUND: Thoracic aortic aneurysms are associated with a pathologic lesion termed "medial degeneration," which is described as a noninflammatory lesion. Thoracic aortic aneurysms are a complication of Marfan syndrome and can be inherited in an autosomal dominant manner of familial thoracic aortic aneurysm. METHODS: Full aortic segments were collected from patients undergoing elective repair with Marfan syndrome (n = 5), familial thoracic aortic aneurysm (n = 6), and thoracic aortic aneurysms (n = 9), along with control aortas (n = 5). Immunohistochemistry staining was performed using antibodies directed against markers of lymphocytes and macrophages. Real-time polymerase chain reaction analysis was performed to quantify the expression level of the T-cell receptor beta-chain variable region gene. RESULTS: Immunohistochemistry of thoracic aortic aneurysm aortas demonstrated that the media and adventitia from Marfan syndrome, familial thoracic aortic aneurysm, and sporadic cases had increased numbers of T lymphocytes and macrophages when compared with control aortas. The number of T cells and macrophages in the aortic media of the aneurysm correlated inversely with the patient's age at the time of prophylactic surgical repair of the aorta. T-cell receptor profiling indicated a similar clonal nature of the T cells in the aortic wall in a majority of aneurysms, whether the patient had Marfan syndrome, familial thoracic aortic aneurysm, or sporadic disease. CONCLUSION: These results indicate that the infiltration of inflammatory cells contributes to the pathogenesis of thoracic aortic aneurysms. Superantigen-driven stimulation of T lymphocytes in the aortic tissues of patients with thoracic aortic aneurysms may contribute to the initial immune response.

Evaluating statistical methods used to estimate the number of postsynaptic receptors.

Calcium levels in spines play a significant role in determining the sign and magnitude of synaptic plasticity. The magnitude of calcium influx into spines is highly dependent on influx through N-methyl D-aspartate (NMDA) receptors, and therefore depends on the number of postsynaptic NMDA receptors in each spine. We have calculated previously how the number of postsynaptic NMDA receptors determines the mean and variance of calcium transients in the postsynaptic density, and how this alters the shape of plasticity curves. However, the number of postsynaptic NMDA receptors in the postsynaptic density is not well known. Anatomical methods for estimating the number of NMDA receptors produce estimates that are very different than those produced by physiological techniques. The physiological techniques are based on the statistics of synaptic transmission and it is difficult to experimentally estimate their precision. In this paper we use stochastic simulations in order to test the validity of a physiological estimation technique based on failure analysis. We find that the method is likely to underestimate the number of postsynaptic NMDA receptors, explain the source of the error, and re-derive a more precise estimation technique. We also show that the original failure analysis as well as our improved formulas are not robust to small estimation errors in key parameters.

Association of part-time employment and sexual risk-taking behaviors among middle school students

Background. Previous research shows inconsistent results as to the association between part-time employment and sexual behavior among younger teens. Studies of older teens cannot be generalized to younger teens because of the wide differences in types of work performed, nature of work environments, and work intensity. Objective. Examine the relationship between part-time employment and sexual behavior in a cross-sectional sample of public middle school students in Houston, Texas. Methods . The study presents a secondary analysis of data from the It’s Your Game…Keep it Real baseline data collection (11/2004–1/2005). It’s Your Game… is an intervention program for middle school students designed to prevent Sexually Transmitted Infections. Statistical analysis. Univariate and multivariate logistic regression analyses were conducted to examine the association between part-time employment and vaginal intercourse: (a) ever had sex; and (b) current sexual activity. Results. Overall, 13.2% of students worked for pay; male students were 1.5 times as likely as females to be working. Of all the students, 11.0% had had sexual intercourse; students who worked were 3 times more likely to be sexually experienced than those who did not. Among students who were sexually experienced, 67.0% were currently sexually active. After adjusting for the other covariates, Hispanic students were almost 3.6 times more likely to report current sexual activity compared to students in other racial/ethnic groups. In univariate analysis, students who worked 1-5 hrs/week were more likely to be sexually experienced than those not currently employed, and the likelihood increased with number of hours worked. There is a similar pattern in the multivariate model, but the odds ratios are too close for the evidence to be more than suggestive. Of sexually experienced students, students working 1-5 hrs/week were 2.7 times more likely to report current sexual intercourse than those not working; those working >5 hrs/week were 4.7 times more likely. The multivariate model showed a similar increase in likelihood, and adjustment for covariates increased these associations: students who worked 1-5 hrs/week were 3.6 times more likely to report current sexual intercourse, and students who worked >5 hrs/week were 4.5 times more likely, than students not currently employed.^

Statistical and methodological challenges for disaster preparedness and medical needs assessment in Rio Grande Valley of Texas

In recent years, disaster preparedness through assessment of medical and special needs persons (MSNP) has taken a center place in public eye in effect of frequent natural disasters such as hurricanes, storm surge or tsunami due to climate change and increased human activity on our planet. Statistical methods complex survey design and analysis have equally gained significance as a consequence. However, there exist many challenges still, to infer such assessments over the target population for policy level advocacy and implementation. ^ Objective. This study discusses the use of some of the statistical methods for disaster preparedness and medical needs assessment to facilitate local and state governments for its policy level decision making and logistic support to avoid any loss of life and property in future calamities. ^ Methods. In order to obtain precise and unbiased estimates for Medical Special Needs Persons (MSNP) and disaster preparedness for evacuation in Rio Grande Valley (RGV) of Texas, a stratified and cluster-randomized multi-stage sampling design was implemented. US School of Public Health, Brownsville surveyed 3088 households in three counties namely Cameron, Hidalgo, and Willacy. Multiple statistical methods were implemented and estimates were obtained taking into count probability of selection and clustering effects. Statistical methods for data analysis discussed were Multivariate Linear Regression (MLR), Survey Linear Regression (Svy-Reg), Generalized Estimation Equation (GEE) and Multilevel Mixed Models (MLM) all with and without sampling weights. ^ Results. Estimated population for RGV was 1,146,796. There were 51.5% female, 90% Hispanic, 73% married, 56% unemployed and 37% with their personal transport. 40% people attained education up to elementary school, another 42% reaching high school and only 18% went to college. Median household income is less than $15,000/year. MSNP estimated to be 44,196 (3.98%) [95% CI: 39,029; 51,123]. All statistical models are in concordance with MSNP estimates ranging from 44,000 to 48,000. MSNP estimates for statistical methods are: MLR (47,707; 95% CI: 42,462; 52,999), MLR with weights (45,882; 95% CI: 39,792; 51,972), Bootstrap Regression (47,730; 95% CI: 41,629; 53,785), GEE (47,649; 95% CI: 41,629; 53,670), GEE with weights (45,076; 95% CI: 39,029; 51,123), Svy-Reg (44,196; 95% CI: 40,004; 48,390) and MLM (46,513; 95% CI: 39,869; 53,157). ^ Conclusion. RGV is a flood zone, most susceptible to hurricanes and other natural disasters. People in the region are mostly Hispanic, under-educated with least income levels in the U.S. In case of any disaster people in large are incapacitated with only 37% have their personal transport to take care of MSNP. Local and state government’s intervention in terms of planning, preparation and support for evacuation is necessary in any such disaster to avoid loss of precious human life. ^ Key words: Complex Surveys, statistical methods, multilevel models, cluster randomized, sampling weights, raking, survey regression, generalized estimation equations (GEE), random effects, Intracluster correlation coefficient (ICC).^

Outcomes associated with chemotherapy in elderly patients with early stage operable breast cancer

Background. Various clinical trials have proved the efficacy of adjuvant chemotherapy in women with breast cancer. Chemotherapy efficacy and guidelines for its use differ by stage of tumor and age of the patient with no clear recommendations for patients aged 70 and above. Objective. To examine the clinical and economic outcomes associated with chemotherapy use in and to examine the disparities in treatment and survival in elderly patients with early stage operable breast cancer by age and axillary node status. Methods. We studied a cohort of 23,110 node positive and 31,572 node negative women aged 65 and over diagnosed with incident American Joint Committee on Cancer (AJCC) stage I, II or IIIa breast cancer between January 1, 1991 and December 31, 2002 using SEER-Medicare data. Total patient costs were estimated using the phase of care approach and adjusted cost estimates were obtained from regression analysis using a 3% discount rate. Cox proportional hazard ratio of mortality was used to determine the effectiveness of chemotherapy. Propensity score approach was also used to minimize the bias associated with receipt of chemotherapy. To assess disparity in treatment, multivariate logistic regression analyses were performed to assess the relative odds of receiving surgery, chemotherapy and radiation after BCS for African Americans compared to Whites. Results. Regression adjusted cost estimates for all node positive patients receiving chemotherapy was approximately $2,300 and was significantly higher (p<0.05) than for patients not receiving chemotherapy. Mortality was significantly lower in node positive and node negative women aged 65-74 years receiving chemotherapy. There was a significant difference between African American and White women in receiving BCS and radiation after BCS; however this difference was explained by patient demographics, tumor characteristics and socioeconomic status (SES). African American node positive women were 21% less likely to receive chemotherapy than White women (OR, 0.79; CI, 0.68-0.92) in multivariate analysis. Conclusion. Chemotherapy is associated with increased survival in patients aged 65-74 and total costs attributable to chemotherapy differ by phase and age of the patient. Underutilization of systemic adjuvant chemotherapy in African American women requires attention and may serve as potential areas for appropriate intervention.^

Assessment of the effect on statistical power of regression model misspecification by using techniques of mathematical statistics and simulation study

Objectives. This paper seeks to assess the effect on statistical power of regression model misspecification in a variety of situations. ^ Methods and results. The effect of misspecification in regression can be approximated by evaluating the correlation between the correct specification and the misspecification of the outcome variable (Harris 2010).In this paper, three misspecified models (linear, categorical and fractional polynomial) were considered. In the first section, the mathematical method of calculating the correlation between correct and misspecified models with simple mathematical forms was derived and demonstrated. In the second section, data from the National Health and Nutrition Examination Survey (NHANES 2007-2008) were used to examine such correlations. Our study shows that comparing to linear or categorical models, the fractional polynomial models, with the higher correlations, provided a better approximation of the true relationship, which was illustrated by LOESS regression. In the third section, we present the results of simulation studies that demonstrate overall misspecification in regression can produce marked decreases in power with small sample sizes. However, the categorical model had greatest power, ranging from 0.877 to 0.936 depending on sample size and outcome variable used. The power of fractional polynomial model was close to that of linear model, which ranged from 0.69 to 0.83, and appeared to be affected by the increased degrees of freedom of this model.^ Conclusion. Correlations between alternative model specifications can be used to provide a good approximation of the effect on statistical power of misspecification when the sample size is large. When model specifications have known simple mathematical forms, such correlations can be calculated mathematically. Actual public health data from NHANES 2007-2008 were used as examples to demonstrate the situations with unknown or complex correct model specification. Simulation of power for misspecified models confirmed the results based on correlation methods but also illustrated the effect of model degrees of freedom on power.^

Statistical Methods for Differential Expressions of Genes Detected in Multiple-Condition Experiment of Microarray

Most studies of differential gene-expressions have been conducted between two given conditions. The two-condition experimental (TCE) approach is simple in that all genes detected display a common differential expression pattern responsive to a common two-condition difference. Therefore, the genes that are differentially expressed under the other conditions other than the given two conditions are undetectable with the TCE approach. In order to address the problem, we propose a new approach called multiple-condition experiment (MCE) without replication and develop corresponding statistical methods including inference of pairs of conditions for genes, new t-statistics, and a generalized multiple-testing method for any multiple-testing procedure via a control parameter C. We applied these statistical methods to analyze our real MCE data from breast cancer cell lines and found that 85 percent of gene-expression variations were caused by genotypic effects and genotype-ANAX1 overexpression interactions, which agrees well with our expected results. We also applied our methods to the adenoma dataset of Notterman et al. and identified 93 differentially expressed genes that could not be found in TCE. The MCE approach is a conceptual breakthrough in many aspects: (a) many conditions of interests can be conducted simultaneously; (b) study of association between differential expressions of genes and conditions becomes easy; (c) it can provide more precise information for molecular classification and diagnosis of tumors; (d) it can save lot of experimental resources and time for investigators.^

Extension of k-ratio methods to rank-based analyses

An extension of k-ratio multiple comparison methods to rank-based analyses is described. The new method is analogous to the Duncan-Godbold approximate k-ratio procedure for unequal sample sizes or correlated means. The close parallel of the new methods to the Duncan-Godbold approach is shown by demonstrating that they are based upon different parameterizations as starting points.^ A semi-parametric basis for the new methods is shown by starting from the Cox proportional hazards model, using Wald statistics. From there the log-rank and Gehan-Breslow-Wilcoxon methods may be seen as score statistic based methods.^ Simulations and analysis of a published data set are used to show the performance of the new methods. ^

Selection of body fat distribution indices in adults and associations with plasma lipid metabolism variables

Body fat distribution is a cardiovascular health risk factor in adults. Body fat distribution can be measured through various methods including anthropometry. It is not clear which anthropometric index is suitable for epidemiologic studies of fat distribution and cardiovascular disease. The purpose of the present study was to select a measure of body fat distribution from among a series of indices (those traditionally used in the literature and others constructed from the analysis) that is most highly correlated with lipid-related variables and is independent of overall fatness. Subjects were Mexican-American men and women (N = 1004) from a study of gallbladder disease in Starr County, Texas. Multivariate associations were sought between lipid profile measures (lipids, lipoproteins, and apolipoproteins) and two sets of anthropometric variables (4 circumferences and 6 skinfolds). This was done to assess the association between lipid-related measures and the two sets of anthropometric variables and guide the construction of indices.^ Two indices emerged from the analysis that seemed to be highly correlated with lipid profile measures independent of obesity. These indices are: 2*arm circumference-thigh skinfold in pre- and post-menopausal women and arm/thigh circumference ratio in men. Next, using the sum of all skinfolds to represent obesity and the selected body fat distribution indices, the following hypotheses were tested: (1) state of obesity and centrally/upper distributed body fat are equally predictive of lipids, lipoproteins and apolipoproteins, and (2) the correlation among the lipid-related measures is not altered by obesity and body fat distribution.^ With respect to the first hypothesis, the present study found that most lipids, lipoproteins and apolipoproteins were significantly associated with both overall fatness and anatomical location of body fat in both sex and menopausal groups. However, within men and post-menopausal women, certain lipid profile measures (triglyceride and HDLT among post-menopausal women and apos C-II, CIII, and E among men) had substantially higher correlation with body fat distribution as compared with overall fatness.^ With respect to the second hypothesis, both obesity and body fat distribution were found to alter the association among plasma lipid variables in men and women. There was a suggestion from the data that the pattern of correlations among men and post-menopausal women are more comparable. Among men correlations involving apo A-I, HDLT, and HDL$\sb2$ seemed greatly influenced by obesity, and A-II by fat distribution; among post-menopausal women correlations involving apos A-I and A-II were highly affected by the location of body fat.^ Thus, these data point out that not only can obesity and fat distribution affect levels of single measures, they also can markedly influence the pattern of relationship among measures. The fact that such changes are seen for both obesity and fat distribution is significant, since the indices employed were chosen because they were independent of one another. ^

THE USE OF MULTIVARIATE ANALYSIS OF COVARIANCE IN EVALUATING THE COMPARABILITY OF LABORATORY DETERMINATIONS IN COOPERATIVE STUDIES

The role of clinical chemistry has traditionally been to evaluate acutely ill or hospitalized patients. Traditional statistical methods have serious drawbacks in that they use univariate techniques. To demonstrate alternative methodology, a multivariate analysis of covariance model was developed and applied to the data from the Cooperative Study of Sickle Cell Disease.^ The purpose of developing the model for the laboratory data from the CSSCD was to evaluate the comparability of the results from the different clinics. Several variables were incorporated into the model in order to control for possible differences among the clinics that might confound any real laboratory differences.^ Differences for LDH, alkaline phosphatase and SGOT were identified which will necessitate adjustments by clinic whenever these data are used. In addition, aberrant clinic values for LDH, creatinine and BUN were also identified.^ The use of any statistical technique including multivariate analysis without thoughtful consideration may lead to spurious conclusions that may not be corrected for some time, if ever. However, the advantages of multivariate analysis far outweigh its potential problems. If its use increases as it should, the applicability to the analysis of laboratory data in prospective patient monitoring, quality control programs, and interpretation of data from cooperative studies could well have a major impact on the health and well being of a large number of individuals. ^

FACTORS INFLUENCING CHOICE OF CONTRACEPTIVE METHOD AMONG MARRIED FECUND WOMEN WHO INTEND NO ADDITIONAL BIRTHS: HEALTH BELIEF MODEL AND ECONOMIC PERSPECTIVES

This research examined to what extent Health Belief Model (HBM) and socioeconomic variables were useful in explaining the choice whether or not more effective contraceptive methods were used among married fecund women intending no additional births. The source of the data was the 1976 National Survey of Family Growth conducted under the auspices of the National Center for Health Statistics. Using the HBM as a framework for multivariate analyses limited support was found (using available measures) that the HBM components of motivation and perceived efficacy influence the likelihood of more effective contraceptive method use. Support was also found that modifying variables suggested by the HBM can influence the effects of HBM components on the likelihood of more effective method use. Socioeconomic variables were found, using all cases and some subgroups, to have a significant additional influence on the likelihood of use of more effective methods. Limited support was found for the concept that the greater the opportunity costs of an unwanted birth the greater the likelihood of use of more effective contraceptive methods. This research supports the use of HBM and socioeconomic variables to explain the likelihood of a protective health behavior, use of more effective contraception if no additional births are intended.^