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. ^

Models of DNA sequence evolution

Models of DNA sequence evolution and methods for estimating evolutionary distances are needed for studying the rate and pattern of molecular evolution and for inferring the evolutionary relationships of organisms or genes. In this dissertation, several new models and methods are developed.^ The rate variation among nucleotide sites: To obtain unbiased estimates of evolutionary distances, the rate heterogeneity among nucleotide sites of a gene should be considered. Commonly, it is assumed that the substitution rate varies among sites according to a gamma distribution (gamma model) or, more generally, an invariant+gamma model which includes some invariable sites. A maximum likelihood (ML) approach was developed for estimating the shape parameter of the gamma distribution $(\alpha)$ and/or the proportion of invariable sites $(\theta).$ Computer simulation showed that (1) under the gamma model, $\alpha$ can be well estimated from 3 or 4 sequences if the sequence length is long; and (2) the distance estimate is unbiased and robust against violations of the assumptions of the invariant+gamma model.^ However, this ML method requires a huge amount of computational time and is useful only for less than 6 sequences. Therefore, I developed a fast method for estimating $\alpha,$ which is easy to implement and requires no knowledge of tree. A computer program was developed for estimating $\alpha$ and evolutionary distances, which can handle the number of sequences as large as 30.^ Evolutionary distances under the stationary, time-reversible (SR) model: The SR model is a general model of nucleotide substitution, which assumes (i) stationary nucleotide frequencies and (ii) time-reversibility. It can be extended to SRV model which allows rate variation among sites. I developed a method for estimating the distance under the SR or SRV model, as well as the variance-covariance matrix of distances. Computer simulation showed that the SR method is better than a simpler method when the sequence length $L>1,000$ bp and is robust against deviations from time-reversibility. As expected, when the rate varies among sites, the SRV method is much better than the SR method.^ The evolutionary distances under nonstationary nucleotide frequencies: The statistical properties of the paralinear and LogDet distances under nonstationary nucleotide frequencies were studied. First, I developed formulas for correcting the estimation biases of the paralinear and LogDet distances. The performances of these formulas and the formulas for sampling variances were examined by computer simulation. Second, I developed a method for estimating the variance-covariance matrix of the paralinear distance, so that statistical tests of phylogenies can be conducted when the nucleotide frequencies are nonstationary. Third, a new method for testing the molecular clock hypothesis was developed in the nonstationary case. ^

A comparison of Markov and generalized linear models of hospital mortality

This paper reports a comparison of three modeling strategies for the analysis of hospital mortality in a sample of general medicine inpatients in a Department of Veterans Affairs medical center. Logistic regression, a Markov chain model, and longitudinal logistic regression were evaluated on predictive performance as measured by the c-index and on accuracy of expected numbers of deaths compared to observed. The logistic regression used patient information collected at admission; the Markov model was comprised of two absorbing states for discharge and death and three transient states reflecting increasing severity of illness as measured by laboratory data collected during the hospital stay; longitudinal regression employed Generalized Estimating Equations (GEE) to model covariance structure for the repeated binary outcome. Results showed that the logistic regression predicted hospital mortality as well as the alternative methods but was limited in scope of application. The Markov chain provides insights into how day to day changes of illness severity lead to discharge or death. The longitudinal logistic regression showed that increasing illness trajectory is associated with hospital mortality. The conclusion is reached that for standard applications in modeling hospital mortality, logistic regression is adequate, but for new challenges facing health services research today, alternative methods are equally predictive, practical, and can provide new insights. ^

A nonparametric method of comparing the partial areas under two correlated receiver operating characteristic curves

A non-parametric method was developed and tested to compare the partial areas under two correlated Receiver Operating Characteristic curves. Based on the theory of generalized U-statistics the mathematical formulas have been derived for computing ROC area, and the variance and covariance between the portions of two ROC curves. A practical SAS application also has been developed to facilitate the calculations. The accuracy of the non-parametric method was evaluated by comparing it to other methods. By applying our method to the data from a published ROC analysis of CT image, our results are very close to theirs. A hypothetical example was used to demonstrate the effects of two crossed ROC curves. The two ROC areas are the same. However each portion of the area between two ROC curves were found to be significantly different by the partial ROC curve analysis. For computation of ROC curves with large scales, such as a logistic regression model, we applied our method to the breast cancer study with Medicare claims data. It yielded the same ROC area computation as the SAS Logistic procedure. Our method also provides an alternative to the global summary of ROC area comparison by directly comparing the true-positive rates for two regression models and by determining the range of false-positive values where the models differ. ^

Interrelationships between the structural and functional aspects of social relationships, comorbidity and health related quality of life in breast cancer survivors

It is widely acknowledged in theoretical and empirical literature that social relationships, comprising of structural measures (social networks) and functional measures (perceived social support) have an undeniable effect on health outcomes. However, the actual mechanism of this effect has yet to be clearly understood or explicated. In addition, comorbidity is found to adversely affect social relationships and health related quality of life (a valued outcome measure in cancer patients and survivors). ^ This cross sectional study uses selected baseline data (N=3088) from the Women's Healthy Eating and Living (WHEL) study. Lisrel 8.72 was used for the latent variable structural equation modeling. Due to the ordinal nature of the data, Weighted Least Squares (WLS) method of estimation using Asymptotic Distribution Free covariance matrices was chosen for this analysis. The primary exogenous predictor variables are Social Networks and Comorbidity; Perceived Social Support is the endogenous predictor variable. Three dimensions of HRQoL, physical, mental and satisfaction with current quality of life were the outcome variables. ^ This study hypothesizes and tests the mechanism and pathways between comorbidity, social relationships and HRQoL using latent variable structural equation modeling. After testing the measurement models of social networks and perceived social support, a structural model hypothesizing associations between the latent exogenous and endogenous variables was tested. The results of the study after listwise deletion (N=2131) mostly confirmed the hypothesized relationships (TLI, CFI >0.95, RMSEA = 0.05, p=0.15). Comorbidity was adversely associated with all three HRQoL outcomes. Strong ties were negatively associated with perceived social support; social network had a strong positive association with perceived social support, which served as a mediator between social networks and HRQoL. Mental health quality of life was the most adversely affected by the predictor variables. ^ This study is a preliminary look at the integration of structural and functional measures of social relationships, comorbidity and three HRQoL indicators using LVSEM. Developing stronger social networks and forming supportive relationships is beneficial for health outcomes such as HRQoL of cancer survivors. Thus, the medical community treating cancer survivors as well as the survivor's social networks need to be informed and cognizant of these possible relationships. ^

Clinical factors associated with growth and mortality of pediatric patients on HAART at Mulago PIDC, Uganda, 2003--2006

The global social and economic burden of HIV/AIDS is great, with over forty million people reported to be living with HIV/AIDS at the end of 2005; two million of these are children from birth to 15 years of age. Antiretroviral therapy has been shown to improve growth and survival of HIV-infected individuals. The purpose of this study is to describe a cohort of HIV-infected pediatric patients and assess the association between clinical factors, with growth and mortality outcomes. ^ This was a historical cohort study. Medical records of infants and children receiving HIV care at Mulago Pediatric Infectious Disease Clinic (PIDC) in Uganda between July 2003 and March 2006 were analyzed. Height and weight measurements were age and sex standardized to Centers for Disease Control and prevention (CDC) 2000 reference. Descriptive and logistic regression analyses were performed to identify covariates associated with risk of stunting or being underweight, and mortality. Longitudinal regression analysis with a mixed model using autoregressive covariance structure was used to compare change in height and weight before and after initiation of highly active antiretroviral therapy (HAART). ^ The study population was comprised of 1059 patients 0-20 years of age, the majority of whom were aged thirteen years and below (74.6%). Mean height-for-age before initiation of HAART was in the 10th percentile, mean weight-for-age was in the 8th percentile, and the mean weight-for-height was in the 23rd percentile. Initiation of HAART resulted in improvement in both the mean standardized weight-for-age Z score and weight-for-age percentiles (p <0.001). Baseline age, and weight-for-age Z score were associated with stunting (p <0.001). A negative weight-for-age Z score was associated with stunting (OR 4.60, CI 3.04-5.49). Risk of death decreased from 84% in the >2-8 years age category to 21% in the >13 years age category respectively, compared to the 0-2 years of age (p <0.05). ^ This pediatric population gained weight significantly more rapidly than height after starting HAART. A low weight-for-age Z score was associated with poor survival in children. These findings suggest that age, weight, and height measurements be monitored closely at Mulago PIDC. ^

The effectiveness of a promotora-led intervention for Mexican Americans with type 2 diabetes

Purpose. The purpose of this randomized control repeated measures trial was to determine the effectiveness of a self-management intervention led by community lay workers called promotoras on the health outcomes of Mexican Americans with type 2 diabetes living in a major city on the Texas - Mexico border. The specific aims of this study, in relation to the intervention group participants, were to: (1) decrease the glycosylated hemoglobin (A1c) blood levels at the six-month assessment, (2) increase diabetes knowledge at the three and six-month assessments, and (3) strengthen the participants' beliefs in their ability to manage diabetes at the three and six-month assessments.^ Methods. One hundred and fifty Mexican American participants were recruited at a Catholic faith-based clinic and randomized into an intervention group and a usual-care control group. Personal characteristics, acculturation and baseline A1c, diabetes knowledge and diabetes health beliefs were measured. The six-month, two-phase intervention was culturally specific and it was delivered entirely by promotoras. Phase One of the intervention consisted of sixteen hours of participative group education and bi-weekly telephone contact follow-up. Phase Two consisted of bi-weekly follow-up using inspirational faith-based health behavior change postcards. The A1c levels, diabetes knowledge and diabetes health beliefs were measured at baseline, and three and six months post-baseline. The mean changes between the groups were analyzed using analysis of covariance. ^ Results. The 80% female sample, with a mean age of 58 years, demonstrated very low: acculturation, income, education, health insurance coverage, and strong Catholicism. No significant changes were noted at the three-month assessment, but the mean change of the A1c levels (F (1, 148 = 10.28, p < .001) and the diabetes knowledge scores (F (1, 148 = 9.0, p < .002) of the intervention group improved significantly at six months, adjusting for health insurance coverage. The diabetes health belief scores decreased in both groups.^ Conclusions. This study demonstrated that an intervention led by promotoras could result in decreased A1c levels and increased diabetes knowledge in spite of the very low acculturation, educational level and insurance coverage of the intervention group participants. Clinical implications and recommendations for future research are suggested. ^

Bayesian generalized linear models for meta-analysis of diagnostic tests

With the recognition of the importance of evidence-based medicine, there is an emerging need for methods to systematically synthesize available data. Specifically, methods to provide accurate estimates of test characteristics for diagnostic tests are needed to help physicians make better clinical decisions. To provide more flexible approaches for meta-analysis of diagnostic tests, we developed three Bayesian generalized linear models. Two of these models, a bivariate normal and a binomial model, analyzed pairs of sensitivity and specificity values while incorporating the correlation between these two outcome variables. Noninformative independent uniform priors were used for the variance of sensitivity, specificity and correlation. We also applied an inverse Wishart prior to check the sensitivity of the results. The third model was a multinomial model where the test results were modeled as multinomial random variables. All three models can include specific imaging techniques as covariates in order to compare performance. Vague normal priors were assigned to the coefficients of the covariates. The computations were carried out using the 'Bayesian inference using Gibbs sampling' implementation of Markov chain Monte Carlo techniques. We investigated the properties of the three proposed models through extensive simulation studies. We also applied these models to a previously published meta-analysis dataset on cervical cancer as well as to an unpublished melanoma dataset. In general, our findings show that the point estimates of sensitivity and specificity were consistent among Bayesian and frequentist bivariate normal and binomial models. However, in the simulation studies, the estimates of the correlation coefficient from Bayesian bivariate models are not as good as those obtained from frequentist estimation regardless of which prior distribution was used for the covariance matrix. The Bayesian multinomial model consistently underestimated the sensitivity and specificity regardless of the sample size and correlation coefficient. In conclusion, the Bayesian bivariate binomial model provides the most flexible framework for future applications because of its following strengths: (1) it facilitates direct comparison between different tests; (2) it captures the variability in both sensitivity and specificity simultaneously as well as the intercorrelation between the two; and (3) it can be directly applied to sparse data without ad hoc correction. ^

The impact of nutrition intervention on glycemic control in a pediatric type 1 diabetes mellitus population

Aim: The goal of this study was to evaluate the change in hemoglobin A1C and glycemic control after nutrition intervention among a population of type 1 diabetic pediatric patients. Methods: Data was collected from all type 1 diabetic patients who were scheduled for a consultation with the diabetes/endocrine RD from January 2006 through December 2006. Two groups were compared, those who kept their RD appointment and those who did not keep their appointment. The main outcome measure was HgbA1C. An independent samples t-test compared the two groups with respect to change in HbgA1C before and after the most recent scheduled appointment with the RD. Baseline characteristics were used as covariates and analyzed and controlled for using analysis of covariance (ANCOVA). Results: There was no difference in HgbA1c after either attending an RD appointment or not having attended an RD appointment. Those who arrived for and attended their RD appointment and those who did not arrive for and attend their RD appointment, had statistically different HgbA1C's before their scheduled appointment as well as after the RD appointment. However, the two groups were not equal at the beginning of the study period. Discussion: A study design with inclusion criteria of a specified range of HgbA1C values within which the study subjects needed to fall, would have potentially eliminated the difference between the two groups at the beginning of the study period. Conducting either another retrospective study that controlled for the initial HgbA1C value or conducting a prospective study that designated a range of HgbA1C values would be worth investigating to evaluate the impact of medical nutrition therapy intervention and the role of the RD in diabetes management. It is an interesting finding that there was a significant difference in the initial HgbA1c for those who came to the RD appointment compared to those who did not come. The fact that in this study those who did not arrive for their RD appointment had worse control of their diabetes suggests that this is a high-risk group. Targeting diabetes education toward this group of patients may prove to be beneficial. ^

"Teens talk healthy weight": A video technology project

Purpose. The purpose of this study was to investigate the impact of a motivational weight management DVD on knowledge of obesity related diseases, readiness, motivation, and self-efficacy to lose weight, connectedness to their care provider, and patients return to clinic. Design. A randomized control trial was conducted in which 40 overweight/obese adolescents and their parents/caregivers were randomly assigned to standard care alone or standard care plus DVD. Subjects completed a set of pre- and post-questionnaire measures. A group of 22 patients was also formed as a historical control group in order to account for the potential effect of extra attention given to subjects prospectively enrolled. Methods. The adolescents and their parent/caregiver were placed into a patient room. Consent was obtained and a set of written pre-questionnaires were given to both the parent and the adolescent. Standard care was provided to all patients by the Registered Dietitian and physician; the DVD was shown in addition to standard care among the intervention group. A set of post-questionnaires were given and compensation was provided. Analysis. Groups were compared to determine equivalence at baseline. Analysis of covariance was used to evaluate changes over time, while controlling for pre-test scores and race/ethnicity. Results. Parents who viewed the DVD experienced greater changes in correct knowledge as compared to parents who did not view the DVD. Conclusion. Our study found only one substantial benefit of the DVD beyond standard clinical practices. This is an important area for change as it increased awareness of obesity as a serious disease and has future clinical implications.^

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. ^

A GENETIC MARKER STUDY OF NON - INSULIN DEPENDENT DIABETICS IN THE MEXICAN-AMERICAN COMMUNITY FROM STARR COUNTY, TEXAS (PLASMA GLUCOSE, HEMOGLOBIN ANALYSIS)

Diabetes mellitus occurs in two forms, insulin-dependent (IDDM, formerly called juvenile type) and non-insulin dependent (NIDDM, formerly called adult type). Prevalence figures from around the world for NIDDM, show that all societies and all races are affected; although uncommon in some populations (.4%), it is common (10%) or very common (40%) in others (Tables 1 and 2).^ In Mexican-Americans in particular, the prevalence rates (7-10%) are intermediate to those in Caucasians (1-2%) and Amerindians (35%). Information about the distribution of the disease and identification of high risk groups for developing glucose intolerance or its vascular manifestations by the study of genetic markers will help to clarify and solve some of the problems from the public health and the genetic point of view.^ This research was designed to examine two general areas in relation to NIDDM. The first aims to determine the prevalence of polymorphic genetic markers in two groups distinguished by the presence or absence of diabetes and to observe if there are any genetic marker-disease association (univariate analysis using two by two tables and logistic regression to study the individual and joint effects of the different variables). The second deals with the effect of genetic differences on the variation in fasting plasma glucose and percent glycosylated hemoglobin (HbAl) (analysis of Covariance for each marker, using age and sex as covariates).^ The results from the first analysis were not statistically significant at the corrected p value of 0.003 given the number of tests that were performed. From the analysis of covariance of all the markers studied, only Duffy and Phosphoglucomutase were statistically significant but poor predictors, given that the amount they explain in terms of variation in glycosylated hemoglobin is very small.^ Trying to determine the polygenic component of chronic disease is not an easy task. This study confirms the fact that a larger and random or representative sample is needed to be able to detect differences in the prevalence of a marker for association studies and in the genetic contribution to the variation in glucose and glycosylated hemoglobin. The importance that ethnic homogeneity in the groups studied and standardization in the methodology will have on the results has been stressed. ^

THE EFFECT OF PERSON-ENVIRONMENT INTERACTION ON ANXIETY, DEPRESSION, AND SUBJECTIVE DISTRESS OF OPTOMETRY STUDENTS

The purpose of this study was to investigate whether an incongruence between personality characteristics of individuals and concomitant charcteristics of health professional training environments on salient dimensions contributes to aspects of mental health. The dimensions examined were practical-theoretical orientation and the degree of structure-unstructure. They were selected for study as they are particularly important attributes of students and of learning environments. It was proposed that when the demand of the environment is disparate from the proclivities of the individual, strain arises. This strain was hypothesized to contribute to anxiety, depression, and subjective distress.^ Select subscales on the Omnibus Personality Inventory (OPI) were the operationalized measures for the personality component of the dimensions studied. An environmental index was developed to assess students' perceptions of the learning environment on these same dimensions. The Beck Depression Inventory, State-Trait Anxiety Inventory and General Well-Being schedule measured the outcome variables.^ A congruence model was employed to determine person-environment (P-E) interaction. Scores on the scales of the OPI and the environmental index were divided into high, medium, and low based on the range of scores. Congruence was defined as a match between the level of personality need and the complementary level of the perception of the environment. Alternatively, incongruence was defined as a mismatch between the person and the environment. The consistent category was compared to the inconsistent categories by an analysis of variance procedure. Furthermore, analyses of covariance were conducted with perceived supportiveness of the learning environment and life events external to the learning environment as the covariates. These factors were considered critical influences affecting the outcome measures.^ One hundred and eighty-five students (49% of the population) at the College of Optometry at the University of Houston participated in the study. Students in all four years of the program were equally represented in the study. However, the sample differed from the total population on representation by sex, marital status, and undergraduate major.^ The results of the study did not support the hypotheses. Further, after having adjusted for perceived supportiveness and life events external to the learning environment, there were no statistically significant differences between the congruent category and incongruent categories. Means indicated than the study sample experienced significantly lower depression and subjective distress than the normative samples.^ Results are interpreted in light of their utility for future study design in the investigation of the effects of P-E interaction. Emphasized is the question of the feasibility of testing a P-E interaction model with extant groups. Recommendations for subsequent research are proposed in light of the exploratory nature of the methodology. ^

A STUDY ON THE EPIDEMIOLOGY OF POSITIVE PHYSICAL HEALTH (NHANES I)

HANES 1 detailed sample data were used to operationalize a definition of health in the absence of disease and to describe and compare the characteristics of the normal (healthy) group versus an abnormal (unhealthy) group.^ Parallel screening gave a 3.8 percent prevalence proportion of physical health, with a female:male ratio of 2:1 and younger ages in the healthy group. Statistically significant Mantel-Haenszel gender-age-adjusted odds ratios (MHOR) were estimated among abnormal non-migrants (1.53), skilled workers/unemployed (1.76), annual family incomes of less than $10,000 (1.56), having ever smoked (1.58), and started smoking before 18 years of age (1.58). Significant MHOR were also found for abnormals for health promoting measures: non-iodized salt use (1.94), needed dental care (1.91); and for fair to poor perceived health (4.28), perceiving health problems (2.52), and low energy level (1.68). Significant protective effects for much to moderate recreational exercise (MHOR 0.42) and very active to moderate non-recreational activity (MHOR 0.49) were also obtained. Covariance analysis additive models detected statistically significant higher mean values for abnormals than normals for serum magnesium, hemoglobin, hematocrit, urinary creatinine, and systolic and diastolic blood pressures, and lower values for abnormals than normals for serum iron. No difference was detected for serum cholesterol. Significant non-additive joint effects were found for body mass index.^ The results suggest positive physical health can be measured with cross-sectional survey data. Gender differentials, and associations between ecologic, socioeconomic, hazardous risk factors, health promoting activities and physical health are in general agreement with published findings on studies of morbidity. Longitudinal prospective studies are suggested to establish the direction of the associations and to enhance present knowledge of health and its promoting factors. ^

DIABETIC RETINOPATHY IN NON-INSULIN-DEPENDENT DIABETES MELLITUS IN BLACK WOMEN: A DESCRIPTIVE STUDY OF RISK FACTORS

This cross-sectional study examines the prevalence of selected potential risk factors by stage of diabetic retinopathy (DR) among Black American women with non-insulin-dependent diabetes mellitus (NIDDM) followed at a university diabetes clinic. DR was assessed by ophthalmoscopy and five-field retinography, and graded on counts of microaneurysms, hemorrhages and/or exudates, and presence of proliferative DR. Prevalence of other vascular diseases was assessed from medical records. Potential risk factors included age, known duration of diabetes, type of hypoglycemic treatment, concentrations of random capillary blood glucose, glycosylated hemoglobin, urine protein and fibrinogen, body mass index, and blood pressure. Prevalence of these risk factors is reported for three categories: No DR, mild background DR, severe background or proliferative DR (including surgically treated DR). Duration, age at diagnosis and treatment of diabetes, concentration of urine protein and average blood glucose, hypertension and cardiovascular disease were significantly associated with DR in univariate analysis. The covariance analysis employed stratification on duration, age at diagnosis and therapy of diabetes. The highest DR scores were calculated for those diagnosed before age 45, regardless of duration, therapy, or average blood glucose. Only individuals diagnosed before age 45 had high blood glucose concentrations in all categories of duration. These findings suggest that in this clinic population of Black women, those diagnosed with NIDDm before age 45 who eventually required insulin treatment were at the greatest risk of developing DR and that longterm poor glucose control is a contributing factor. These results suggest that greater emphasis be placed on this subgroup in allocating the limited resources available to improve the quality of glucose regulation, particularly through measures affecting compliance behavior.^ Findings concerning the association of DR with concentration of blood glucose and urine protein, blood pressure/hypertension and weight were compared with those reported from American Indian and Mexican American populations of the Southwestern United States where prevalence of NIDDM, hypertension and obesity is also high. Additional comparative analyses are outlined to substantiate the preliminary finding that there are systematic differences between these ethnic populations. ^