A categorical data analysis approach to estimation of demographic characteristics for Texas county populations

The need for timely population data for health planning and Indicators of need has Increased the demand for population estimates. The data required to produce estimates is difficult to obtain and the process is time consuming. Estimation methods that require less effort and fewer data are needed. The structure preserving estimator (SPREE) is a promising technique not previously used to estimate county population characteristics. This study first uses traditional regression estimation techniques to produce estimates of county population totals. Then the structure preserving estimator, using the results produced in the first phase as constraints, is evaluated.^ Regression methods are among the most frequently used demographic methods for estimating populations. These methods use symptomatic indicators to predict population change. This research evaluates three regression methods to determine which will produce the best estimates based on the 1970 to 1980 indicators of population change. Strategies for stratifying data to improve the ability of the methods to predict change were tested. Difference-correlation using PMSA strata produced the equation which fit the data the best. Regression diagnostics were used to evaluate the residuals.^ The second phase of this study is to evaluate use of the structure preserving estimator in making estimates of population characteristics. The SPREE estimation approach uses existing data (the association structure) to establish the relationship between the variable of interest and the associated variable(s) at the county level. Marginals at the state level (the allocation structure) supply the current relationship between the variables. The full allocation structure model uses current estimates of county population totals to limit the magnitude of county estimates. The limited full allocation structure model has no constraints on county size. The 1970 county census age - gender population provides the association structure, the allocation structure is the 1980 state age - gender distribution.^ The full allocation model produces good estimates of the 1980 county age - gender populations. An unanticipated finding of this research is that the limited full allocation model produces estimates of county population totals that are superior to those produced by the regression methods. The full allocation model is used to produce estimates of 1986 county population characteristics. ^

NONLINEAR TIME SERIES/SYSTEM ANALYSIS (STATE-SPACE, DYNAMICS, INTERVENTION, RESILIENCE, KALMAN)

A discussion of nonlinear dynamics, demonstrated by the familiar automobile, is followed by the development of a systematic method of analysis of a possibly nonlinear time series using difference equations in the general state-space format. This format allows recursive state-dependent parameter estimation after each observation thereby revealing the dynamics inherent in the system in combination with random external perturbations.^ The one-step ahead prediction errors at each time period, transformed to have constant variance, and the estimated parametric sequences provide the information to (1) formally test whether time series observations y(,t) are some linear function of random errors (ELEM)(,s), for some t and s, or whether the series would more appropriately be described by a nonlinear model such as bilinear, exponential, threshold, etc., (2) formally test whether a statistically significant change has occurred in structure/level either historically or as it occurs, (3) forecast nonlinear system with a new and innovative (but very old numerical) technique utilizing rational functions to extrapolate individual parameters as smooth functions of time which are then combined to obtain the forecast of y and (4) suggest a measure of resilience, i.e. how much perturbation a structure/level can tolerate, whether internal or external to the system, and remain statistically unchanged. Although similar to one-step control, this provides a less rigid way to think about changes affecting social systems.^ Applications consisting of the analysis of some familiar and some simulated series demonstrate the procedure. Empirical results suggest that this state-space or modified augmented Kalman filter may provide interesting ways to identify particular kinds of nonlinearities as they occur in structural change via the state trajectory.^ A computational flow-chart detailing computations and software input and output is provided in the body of the text. IBM Advanced BASIC program listings to accomplish most of the analysis are provided in the appendix. ^

AN ESTIMATION OF IN UTERO DRUG EXPOSURE BY DETERMINANTS OF FREQUENCY AND DISTRIBUTION OF MATERNAL DRUG CONSUMPTION

Data derived from 1,194 gravidas presenting at the observation unit of a city/county hospital between October 11, 1979 through December 7, 1979 were evaluated with respect to the proportion ingesting drugs during pregnancy. The mean age of the mother at the time of the interview was 22.0 years; 43.0 percent were Black; 34.0 percent Latin-American, 21.0 percent White and 2.0 percent other; mean gravida was 2.5 pregnancies; mean parity was 1.0; and mean number of previous abortions was 0.34. Completed interview data was available for 1,119 gravida, corresponding urinalyses for 997 subjects. Ninety and one-tenth percent (90.1 percent) of the subjects reported ingestion of one or more drug preparation(s) (prescription, OTC, or substances used for recreational purposes) during pregnancy with a range of 0 to 11 substances and a mean of 2.7. Dietary supplements (vitamins and minerals) were most frequently reported followed by non-narcotic analgesics. Seventy-six and one tenth percent (76.1 percent) of the population reported consumption of prescription medication, 42.5 percent reported consumption of over-the-counter medications, 45.7 percent reported consumption of a substance for recreational purposes and 4.3 percent reported illicit consumption of a substance. For selected substances, no measurable difference was found between obtaining the information from the interview method or from a urinalysis assay. ^

AN EMPIRICAL ESTIMATION OF MEDICAL AND NON-MEDICAL INPUTS IN THE PRODUCTION OF CARDIO-RESPIRATORY PHYSICAL HEALTH STATUS

The purpose of this study was to examine, in the context of an economic model of health production, the relationship between inputs (health influencing activities) and fitness.^ Primary data were collected from 204 employees of a large insurance company at the time of their enrollment in an industrially-based health promotion program. The inputs of production included medical care use, exercise, smoking, drinking, eating, coronary disease history, and obesity. The variables of age, gender and education known to affect the production process were also examined. Two estimates of fitness were used; self-report and a physiologic estimate based on exercise treadmill performance. Ordinary least squares and two-stage least squares regression analyses were used to estimate the fitness production functions.^ In the production of self-reported fitness status the coefficients for the exercise, smoking, eating, and drinking production inputs, and the control variable of gender were statistically significant and possessed theoretically correct signs. In the production of physiologic fitness exercise, smoking and gender were statistically significant. Exercise and gender were theoretically consistent while smoking was not. Results are compared with previous analyses of health production. ^

New methods for quantification and analysis of quantitative real-time polymerase chain reaction data

Quantitative real-time polymerase chain reaction (qPCR) is a sensitive gene quantitation method that has been widely used in the biological and biomedical fields. The currently used methods for PCR data analysis, including the threshold cycle (CT) method, linear and non-linear model fitting methods, all require subtracting background fluorescence. However, the removal of background fluorescence is usually inaccurate, and therefore can distort results. Here, we propose a new method, the taking-difference linear regression method, to overcome this limitation. Briefly, for each two consecutive PCR cycles, we subtracted the fluorescence in the former cycle from that in the later cycle, transforming the n cycle raw data into n-1 cycle data. Then linear regression was applied to the natural logarithm of the transformed data. Finally, amplification efficiencies and the initial DNA molecular numbers were calculated for each PCR run. To evaluate this new method, we compared it in terms of accuracy and precision with the original linear regression method with three background corrections, being the mean of cycles 1-3, the mean of cycles 3-7, and the minimum. Three criteria, including threshold identification, max R2, and max slope, were employed to search for target data points. Considering that PCR data are time series data, we also applied linear mixed models. Collectively, when the threshold identification criterion was applied and when the linear mixed model was adopted, the taking-difference linear regression method was superior as it gave an accurate estimation of initial DNA amount and a reasonable estimation of PCR amplification efficiencies. When the criteria of max R2 and max slope were used, the original linear regression method gave an accurate estimation of initial DNA amount. Overall, the taking-difference linear regression method avoids the error in subtracting an unknown background and thus it is theoretically more accurate and reliable. This method is easy to perform and the taking-difference strategy can be extended to all current methods for qPCR data analysis.^

Magnesium/Calcium ratios and sea surface temperature estimation for sediments of the Timor Sea

We use a multiproxy approach to monitor changes in the vertical profile of the Indonesian Throughflow as well as monsoonal wind and precipitation patterns in the Timor Sea on glacial-interglacial, precessional, and suborbital timescales. We focus on an interval of extreme climate change and sea level variation: marine isotope (MIS) 6 to MIS 5e. Paleoproductivity fluctuations in the Timor Sea follow a precessional beat related to the intensity of the Australian (NW) monsoon. Paired Mg/Ca and d18O measurements of surface- and thermocline-dwelling planktonic foraminifers (G. ruber and P. obliquiloculata) indicate an increase of >4°C in both surface and thermocline water temperatures during Termination II. Tropical sea surface temperature changed synchronously with ice volume (benthic d18O) during deglaciation, implying a direct coupling of high- and low-latitude climate via atmospheric and/or upper ocean circulation. Substantial cooling and freshening of thermocline waters occurred toward the end of Termination II and during MIS 5e, indicating a change in the vertical profile of the Indonesian Throughflow from surface- to thermocline-dominated flow.

Cretaceous stable isotopic record and sea surface temperture estimation for ODP Site 207-1259

The Turonian (93.5 to 89.3 million years ago) was one of the warmest periods of the Phanerozoic eon, with tropical sea surface temperatures over 35°C. High-amplitude sea-level changes and positive d18O excursions in marine limestones suggest that glaciation events may have punctuated this episode of extreme warmth. New d18O data from the tropical Atlantic show synchronous shifts ~91.2 million years ago for both the surface and deep ocean that are consistent with an approximately 200,000-year period of glaciation, with ice sheets of about half the size of the modern Antarctic ice cap. Even the prevailing supergreenhouse climate was not a barrier to the formation of large ice sheets, calling into question the common assumption that the poles were always ice-free during past periods of intense global warming.

Stable isotope ratios of foraminifera and sea surface temperature estimation of sediment core W2804A-14

Carbon isotopically based estimates of CO2 levels have been generated from a record of the photosynthetic fractionation of 13C (epsilon p) in a central equatorial Pacific sediment core that spans the last ~255 ka. Contents of 13C in phytoplanktonic biomass were determined by analysis of C37 alkadienones. These compounds are exclusive products of Prymnesiophyte algae which at present grow most abundantly at depths of 70-90 m in the central equatorial Pacific. A record of the isotopic compostion of dissolved CO2 was constructed from isotopic analyses of the planktonic foraminifera Neogloboquadrina dutertrei, which calcifies at 70-90 m in the same region. Values of epsilon p, derived by comparison of the organic and inorganic delta values, were transformed to yield concentrations of dissolved CO2 (c e) based on a new, site-specific calibration of the relationship between epsilon p and c e. The calibration was based on reassessment of existing epsilon p versus c e data, which support a physiologically based model in which epsilon p is inversely related to c e. Values of PCO2, the partial pressure of CO2 that would be in equilibrium with the estimated concentrations of dissolved CO2, were calculated using Henry's law and the temperature determined from the alkenone-unsaturation index UK 37. Uncertainties in these values arise mainly from uncertainties about the appropriateness (particularly over time) of the site-specific relationship between epsilon p and 1/c e. These are discussed in detail and it is concluded that the observed record of epsilon p most probably reflects significant variations in Delta pCO2, the ocean-atmosphere disequilibrium, which appears to have ranged from ~110 µatm during glacial intervals (ocean > atmosphere) to ~60 µatm during interglacials. Fluxes of CO2 to the atmosphere would thus have been significantly larger during glacial intervals. If this were characteristic of large areas of the equatorial Pacific, then greater glacial sinks for the equatorially evaded CO2 must have existed elsewhere. Statistical analysis of air-sea pCO2 differences and other parameters revealed significant (p < 0.01) inverse correlations of Delta pCO2 with sea surface temperature and with the mass accumulation rate of opal. The former suggests response to the strength of upwelling, the latter may indicate either drawdown of CO2 by siliceous phytoplankton or variation of [CO2]/[Si(OH)4] ratios in upwelling waters.