Statistical analysis of survey-based event history data with application to modeling of unemployment duration

Longitudinal surveys are increasingly used to collect event history data on person-specific processes such as transitions between labour market states. Surveybased event history data pose a number of challenges for statistical analysis. These challenges include survey errors due to sampling, non-response, attrition and measurement. This study deals with non-response, attrition and measurement errors in event history data and the bias caused by them in event history analysis. The study also discusses some choices faced by a researcher using longitudinal survey data for event history analysis and demonstrates their effects. These choices include, whether a design-based or a model-based approach is taken, which subset of data to use and, if a design-based approach is taken, which weights to use. The study takes advantage of the possibility to use combined longitudinal survey register data. The Finnish subset of European Community Household Panel (FI ECHP) survey for waves 1–5 were linked at person-level with longitudinal register data. Unemployment spells were used as study variables of interest. Lastly, a simulation study was conducted in order to assess the statistical properties of the Inverse Probability of Censoring Weighting (IPCW) method in a survey data context. The study shows how combined longitudinal survey register data can be used to analyse and compare the non-response and attrition processes, test the missingness mechanism type and estimate the size of bias due to non-response and attrition. In our empirical analysis, initial non-response turned out to be a more important source of bias than attrition. Reported unemployment spells were subject to seam effects, omissions, and, to a lesser extent, overreporting. The use of proxy interviews tended to cause spell omissions. An often-ignored phenomenon classification error in reported spell outcomes, was also found in the data. Neither the Missing At Random (MAR) assumption about non-response and attrition mechanisms, nor the classical assumptions about measurement errors, turned out to be valid. Both measurement errors in spell durations and spell outcomes were found to cause bias in estimates from event history models. Low measurement accuracy affected the estimates of baseline hazard most. The design-based estimates based on data from respondents to all waves of interest and weighted by the last wave weights displayed the largest bias. Using all the available data, including the spells by attriters until the time of attrition, helped to reduce attrition bias. Lastly, the simulation study showed that the IPCW correction to design weights reduces bias due to dependent censoring in design-based Kaplan-Meier and Cox proportional hazard model estimators. The study discusses implications of the results for survey organisations collecting event history data, researchers using surveys for event history analysis, and researchers who develop methods to correct for non-sampling biases in event history data.

Probabilistic position and orientation estimation of rigid objects using an RGB-D sensor

The recent emergence of low-cost RGB-D sensors has brought new opportunities for robotics by providing affordable devices that can provide synchronized images with both color and depth information. In this thesis, recent work on pose estimation utilizing RGBD sensors is reviewed. Also, a pose recognition system for rigid objects using RGB-D data is implemented. The implementation uses half-edge primitives extracted from the RGB-D images for pose estimation. The system is based on the probabilistic object representation framework by Detry et al., which utilizes Nonparametric Belief Propagation for pose inference. Experiments are performed on household objects to evaluate the performance and robustness of the system.

Quantitative features of Cedrela odorata L. wood (Meliaceae)

Cedrela odorata L. (Meliaceae) occurs in the Atlantic forest, Amazon rain forest, riparian forest of the country, and wetlands, which demand species adapted to their water conditions. Studies in ecological wood anatomy demonstrated that weather factors' variations have direct influence on the wood anatomical structure and that the fragmentation of the natural habitats is a direct cause of the edge effect which alters the abiotic aspects of the location, interfering consequently in its vegetation. A comparative analysis of 20 anatomical quantitative features of the wood structure was performed in populations of Cedrela odorata growing inside and on the edge of the swamp forest and granulometric analysis was made on the soil. The quantitative data were submitted to the Mann-Whitney's nonparametric test, presenting a statistically significant value decrease in the eleven wood features mean for the specimens growing in the edge of swamp forest.

Transforming growth factor beta activity in urine of patients with type 2 diabetes and diabetic nephropathy

Diabetic nephropathy (DN) is characterized structurally by progressive mesangial deposition of extracellular matrix (ECM). Transforming growth factor-ß (TGF-ß) is considered to be one of the major cytokines involved in the regulation of ECM synthesis and degradation. Several studies suggest that an increase in urinary TGF-ß levels may reflect an enhanced production of this polypeptide by the kidney cells. We evaluated TGF-ß in occasional urine samples from 14 normal individuals and 23 patients with type 2 diabetes (13 with persistent proteinuria >500 mg/24 h, DN, 6 with microalbuminuria, DMMA, and 4 with normal urinary albumin excretion, DMN) by enzyme immunoassay. An increase in the rate of urinary TGF-ß excretion (pg/mg UCreat.) was observed in patients with DN (296.07 ± 330.77) (P<0.001) compared to normal individuals (17.04 ± 18.56) (Kruskal-Wallis nonparametric analysis of variance); however, this increase was not observed in patients with DMMA (25.13 ± 11.30) or in DMN (18.16 ± 11.82). There was a positive correlation between the rate of urinary TGF-ß excretion and proteinuria (r = 0.70, a = 0.05) (Pearson's analysis), one of the parameters of disease progression.

Air pollution and neonatal deaths in São Paulo, Brazil

Air pollution has been associated with health effects on different age groups. The present study was designed to assess the impact of daily changes in air pollutants (NO2, SO2, CO, O3, and particle matter (PM10)) on total number of daily neonatal deaths (those that occur between the first and the 28th days of life) in São Paulo, from January 1998 to December 2000, since adverse outcomes such as neonatal deaths associated with air pollution in Brazil have not been evaluated before. Generalized additive Poisson regression models were used and nonparametric smooth functions (loess) were adopted to control long-term trend, temperature, humidity, and short-term trends. A linear term was used for holidays. The association between air pollutants and neonatal deaths showed a short time lag. Interquartile range increases in PM10 (23.3 µg/m³) and SO2 (9.2 µg/m³) were associated with increases of 4% (95% CI, 2-6) and 6% (95% CI, 4-8), respectively. Instead of adopting a two-pollutant model we created an index to represent PM10 and SO2 effects. For an interquartile range increase in the index an increase of 6.3% (95% CI, 6.1-6.5) in neonatal deaths was observed. These results agree with previous studies performed by our group showing the deleterious effects of air pollutants during the perinatal period. The method reported here represents an alternative approach to analyze the relationship between highly correlated pollutants and public health problems, reinforcing the idea of the synergic effects of air pollutants in public health.

Temperature, air pollution, and mortality from myocardial infarction in São Paulo, Brazil

An increase in daily mortality from myocardial infarction has been observed in association with meteorological factors and air pollution in several cities in the world, mainly in the northern hemisphere. The objective of the present study was to analyze the independent effects of environmental variables on daily counts of death from myocardial infarction in a subtropical region in South America. We used the robust Poisson regression to investigate associations between weather (temperature, humidity and barometric pressure), air pollution (sulfur dioxide, carbon monoxide, and inhalable particulate), and the daily death counts attributed to myocardial infarction in the city of São Paulo in Brazil, where 12,007 fatal events were observed from 1996 to 1998. The model was adjusted in a linear fashion for relative humidity and day-of-week, while nonparametric smoothing factors were used for seasonal trend and temperature. We found a significant association of daily temperature with deaths due to myocardial infarction (P < 0.001), with the lowest mortality being observed at temperatures between 21.6 and 22.6ºC. Relative humidity appeared to exert a protective effect. Sulfur dioxide concentrations correlated linearly with myocardial infarction deaths, increasing the number of fatal events by 3.4% (relative risk of 1.03; 95% confidence interval = 1.02-1.05) for each 10 µg/m³ increase. In conclusion, this study provides evidence of important associations between daily temperature and air pollution and mortality from myocardial infarction in a subtropical region, even after a comprehensive control for confounding factors.

Comparison of anaerobic threshold determined by visual and mathematical methods in healthy women

Several methods are used to estimate anaerobic threshold (AT) during exercise. The aim of the present study was to compare AT obtained by a graphic visual method for the estimate of ventilatory and metabolic variables (gold standard), to a bi-segmental linear regression mathematical model of Hinkley's algorithm applied to heart rate (HR) and carbon dioxide output (VCO2) data. Thirteen young (24 ± 2.63 years old) and 16 postmenopausal (57 ± 4.79 years old) healthy and sedentary women were submitted to a continuous ergospirometric incremental test on an electromagnetic braking cycloergometer with 10 to 20 W/min increases until physical exhaustion. The ventilatory variables were recorded breath-to-breath and HR was obtained beat-to-beat over real time. Data were analyzed by the nonparametric Friedman test and Spearman correlation test with the level of significance set at 5%. Power output (W), HR (bpm), oxygen uptake (VO2; mL kg-1 min-1), VO2 (mL/min), VCO2 (mL/min), and minute ventilation (VE; L/min) data observed at the AT level were similar for both methods and groups studied (P > 0.05). The VO2 (mL kg-1 min-1) data showed significant correlation (P < 0.05) between the gold standard method and the mathematical model when applied to HR (r s = 0.75) and VCO2 (r s = 0.78) data for the subjects as a whole (N = 29). The proposed mathematical method for the detection of changes in response patterns of VCO2 and HR was adequate and promising for AT detection in young and middle-aged women, representing a semi-automatic, non-invasive and objective AT measurement.

Estimation procedures affect the center of pressure frequency analysis

Even though frequency analysis of body sway is widely applied in clinical studies, the lack of standardized procedures concerning power spectrum estimation may provide unreliable descriptors. Stabilometric tests were applied to 35 subjects (20-51 years, 54-95 kg, 1.6-1.9 m) and the power spectral density function was estimated for the anterior-posterior center of pressure time series. The median frequency was compared between power spectra estimated according to signal partitioning, sampling rate, test duration, and detrending methods. The median frequency reliability for different test durations was assessed using the intraclass correlation coefficient. When increasing number of segments, shortening test duration or applying linear detrending, the median frequency values increased significantly up to 137%. Even the shortest test duration provided reliable estimates as observed with the intraclass coefficient (0.74-0.89 confidence interval for a single 20-s test). Clinical assessment of balance may benefit from a standardized protocol for center of pressure spectral analysis that provides an adequate relationship between resolution and variance. An algorithm to estimate center of pressure power density spectrum is also proposed.

Prevalence of vascular-endothelial growth factor, matrix metalloproteinases and tissue inhibitors of metalloproteinases in primary breast cancer

Our objective was to determine the presence of vascular endothelial growth factor (VEGF), matrix metalloproteinase-2 (MMP-2) and MMP-9 and specific tissue inhibitors of matrix metalloproteinase (TIMP-1 and TIMP-2) in tumor samples obtained from patients with primary breast cancer. We attempted to correlate these findings with the status of the sentinel lymph node (SLN) and clinical-pathological characteristics such as age, tumor size, histological type, histological grade, and vascular invasion. Tumor samples from 88 patients with primary breast cancer were analyzed. The immunoreactivity of VEGF, MMP-2, MMP-9, TIMP-1, and TIMP-2 in tumors was correlated with clinical and pathological features, as well as SLN status. Nonparametric, Mann-Whittney, Kruskal-Wallis, and Spearmann tests were used. Categorical variables were analyzed by the Pearson test. No statistically significant correlation was found between the amount of VEGF, MMP-2, MMP-9, TIMP-1, and TIMP-2 and the presence of tumor cells in the SLN. However, larger tumor diameter (P < 0.01) and the presence of vascular invasion (P < 0.01) were correlated positively with a positive SLN. A significant correlation of higher VEGF levels (P = 0.04) and lower TIMP-1 levels (P = 0.04) with ductal histology was also observed. Furthermore, lower TIMP-2 levels showed a statistically significant correlation with younger age (<50 years) and larger tumor diameter (2.0-5.0 cm). A positive SLN correlated significantly with a larger tumor diameter and the presence of vascular invasion. Higher VEGF and lower TIMP-1 levels were observed in patients with ductal tumors, while higher TIMP-1 levels were observed in lobular tumors.

Serum metabolic profiling of human gastric cancer based on gas chromatography/mass spectrometry

Research on molecular mechanisms of carcinogenesis plays an important role in diagnosing and treating gastric cancer. Metabolic profiling may offer the opportunity to understand the molecular mechanism of carcinogenesis and help to non-invasively identify the potential biomarkers for the early diagnosis of human gastric cancer. The aims of this study were to explore the underlying metabolic mechanisms of gastric cancer and to identify biomarkers associated with morbidity. Gas chromatography/mass spectrometry (GC/MS) was used to analyze the serum metabolites of 30 Chinese gastric cancer patients and 30 healthy controls. Diagnostic models for gastric cancer were constructed using orthogonal partial least squares discriminant analysis (OPLS-DA). Acquired metabolomic data were analyzed by the nonparametric Wilcoxon test to find serum metabolic biomarkers for gastric cancer. The OPLS-DA model showed adequate discrimination between cancer and non-cancer cohorts while the model failed to discriminate different pathological stages (I-IV) of gastric cancer patients. A total of 44 endogenous metabolites such as amino acids, organic acids, carbohydrates, fatty acids, and steroids were detected, of which 18 differential metabolites were identified with significant differences. A total of 13 variables were obtained for their greatest contribution in the discriminating OPLS-DA model [variable importance in the projection (VIP) value >1.0], among which 11 metabolites were identified using both VIP values (VIP >1) and the Wilcoxon test. These metabolites potentially revealed perturbations of glycolysis and of amino acid, fatty acid, cholesterol, and nucleotide metabolism of gastric cancer patients. These results suggest that gastric cancer serum metabolic profiling has great potential in detecting this disease and helping to understand its metabolic mechanisms.

Use of the Wasserman equation in optimization of the duration of the power ramp in a cardiopulmonary exercise test: a study of Brazilian men

This study aimed to analyze the agreement between measurements of unloaded oxygen uptake and peak oxygen uptake based on equations proposed by Wasserman and on real measurements directly obtained with the ergospirometry system. We performed an incremental cardiopulmonary exercise test (CPET), which was applied to two groups of sedentary male subjects: one apparently healthy group (HG, n=12) and the other had stable coronary artery disease (n=16). The mean age in the HG was 47±4 years and that in the coronary artery disease group (CG) was 57±8 years. Both groups performed CPET on a cycle ergometer with a ramp-type protocol at an intensity that was calculated according to the Wasserman equation. In the HG, there was no significant difference between measurements predicted by the formula and real measurements obtained in CPET in the unloaded condition. However, at peak effort, a significant difference was observed between oxygen uptake (V˙O2)peak(predicted)and V˙O2peak(real)(nonparametric Wilcoxon test). In the CG, there was a significant difference of 116.26 mL/min between the predicted values by the formula and the real values obtained in the unloaded condition. A significant difference in peak effort was found, where V˙O2peak(real)was 40% lower than V˙O2peak(predicted)(nonparametric Wilcoxon test). There was no agreement between the real and predicted measurements as analyzed by Lin’s coefficient or the Bland and Altman model. The Wasserman formula does not appear to be appropriate for prediction of functional capacity of volunteers. Therefore, this formula cannot precisely predict the increase in power in incremental CPET on a cycle ergometer.

Effects of computer-mediated communication and face-to-face communication on the quantity and quality of discourse produced by English as a second language students /

The effects oftwo types of small-group communication, synchronous computer-mediated and face-to-face, on the quantity and quality of verbal output were con^ared. Quantity was deiSned as the number of turns taken per minute, the number of Analysis-of-Speech units (AS-units) produced per minute, and the number ofwords produced per minute. Quality was defined as the number of words produced per AS-unit. In addition, the interaction of gender and type of communication was explored for any differences that existed in the output produced. Questionnaires were also given to participants to determine attitudes toward computer-mediated and face-to-face communication. Thirty intermediate-level students fi-om the Intensive English Language Program (lELP) at Brock University participated in the study, including 15 females and 15 males. Nonparametric tests, including the Wilcoxon matched-pairs test, Mann-Whitney U test, and Friedman test were used to test for significance at the p < .05 level. No significant differences were found in the effects of computer-mediated and face-to-face communication on the output produced during follow-up speaking sessions. However, the quantity and quality of interaction was significantly higher during face-to-face sessions than computer-mediated sessions. No significant differences were found in the output produced by males and females in these 2 conditions. While participants felt that the use of computer-mediated communication may aid in the development of certain language skills, they generally preferred face-to-face communication. These results differed fi-om previous studies that found a greater quantity and quality of output in addition to a greater equality of interaction produced during computer-mediated sessions in comparison to face-to-face sessions (Kern, 1995; Warschauer, 1996).

Quantum Monte Carlo study of electrostatic polarizabilities of H and He atoms /

The infinitesimal differential quantum Monte Carlo (QMC) technique is used to estimate electrostatic polarizabilities of the H and He atoms up to the sixth order in the electric field perturbation. All 542 different QMC estimators of the nonzero atomic polarizabilities are derived and used in order to decrease the statistical error and to obtain the maximum efficiency of the simulations. We are confident that the estimates are "exact" (free of systematic error): the two atoms are nodeless systems, hence no fixed-node error is introduced. Furthermore, we develope and use techniques which eliminate systematic error inherent when extrapolating our results to zero time-step and large stack-size. The QMC results are consistent with published accurate values obtained using perturbation methods. The precision is found to be related to the number of perturbations, varying from 2 to 4 significant digits.

Enhancing personal effectiveness : the impacts of Costa Rica Outward Bound on students in an independent school in southwestern Ontario

This project examines students in a private school in southwestern Ontario on a 17 -day Costa Rica Outward Bound Rainforest multielement course. The study attempted to discover whether voluntary teenage participants could increase their self-perceptions of life effectiveness by participating in a 17-day expedition. A total of9 students participated in the study. The experimental design that was implemented was a mixed methods design. Participants filled in a Life Effectiveness Questionnaire (LEQ) at four predesignated times during the study. These time intervals occurred (a) before the trip commenced, (b) the first day of the trip, ( c) the last day of the trip, and (d) 1 month after the trip ended. Fieldnotes and recordings from informal group debriefing sessions were also used to gather information. Data collected in this study were analyzed in a variety of ways by the researcher. Analyses that were run on the data included the Friedman test for covariance, means, medians, and the Wilcoxon Pairs Test. The questionnaires were analyzed quantitatively, and the fieldnotes were analyzed qualitatively. Nonparametric statistical analysis was implemented as a result of the small group size of participants. Both sets of data were grouped and discussed according to similarities and differences. The data indicate that voluntary teenage participants experience significant changes over time in the areas of time management, social competency, emotional control, active initiative, and self-confidence. The types of outcomes from this study illustrate that Outward Bound-type opportunities should be offered to teenagers in Ontario schools as a means to bring about self-development.

Using a Bayesian model for bankruptcy prediction : a comparative approach

The purpose of this study is to examine the impact of the choice of cut-off points, sampling procedures, and the business cycle on the accuracy of bankruptcy prediction models. Misclassification can result in erroneous predictions leading to prohibitive costs to firms, investors and the economy. To test the impact of the choice of cut-off points and sampling procedures, three bankruptcy prediction models are assessed- Bayesian, Hazard and Mixed Logit. A salient feature of the study is that the analysis includes both parametric and nonparametric bankruptcy prediction models. A sample of firms from Lynn M. LoPucki Bankruptcy Research Database in the U. S. was used to evaluate the relative performance of the three models. The choice of a cut-off point and sampling procedures were found to affect the rankings of the various models. In general, the results indicate that the empirical cut-off point estimated from the training sample resulted in the lowest misclassification costs for all three models. Although the Hazard and Mixed Logit models resulted in lower costs of misclassification in the randomly selected samples, the Mixed Logit model did not perform as well across varying business-cycles. In general, the Hazard model has the highest predictive power. However, the higher predictive power of the Bayesian model, when the ratio of the cost of Type I errors to the cost of Type II errors is high, is relatively consistent across all sampling methods. Such an advantage of the Bayesian model may make it more attractive in the current economic environment. This study extends recent research comparing the performance of bankruptcy prediction models by identifying under what conditions a model performs better. It also allays a range of user groups, including auditors, shareholders, employees, suppliers, rating agencies, and creditors' concerns with respect to assessing failure risk.