Multivariate analysis applied to clinical analysis data

Dissertação para obtenção do Grau de Mestre em Engenharia e Gestão Industrial

An analysis of strategic paths through the decomposition of total factor productivity: the case of Spanish savings banks

The Spanish savings banks attracted quite a considerable amount of interest within the scientific arena, especially subsequent to the disappearance of the regulatory constraints during the second decade of the 1980s. Nonetheless, a lack of research identified with respect to mainstream paths given by strategic groups, and the analysis of the total factor productivity. Therefore, on the basis of the resource-based view of the firm and cluster analysis, we make use of changes in structure and performance ratios in order to identify the strategic groups extant in the sector. We attain a threeways division, which we link with different input-output specifications defining strategic paths. Consequently, on the basis of these three dissimilar approaches we compute and decompose a Hicks-Moorsteen total factor productivity index. Obtained results put forward an interesting interpretation under a multi-strategic approach, together with the setbacks of employing cluster analysis within a complex strategic environment. Moreover, we also propose an ex-post method of analysing the outcomes of the decomposed total factor productivity index that could be merged with non-traditional techniques of forming strategic groups, such as cognitive approaches.

Effect of school based physical activity programme (KISS) on fitness and adiposity in primary schoolchildren: cluster randomised controlled trial.

To assess the effectiveness of a school based physical activity programme during one school year on physical and psychological health in young schoolchildren. Cluster randomised controlled trial. 28 classes from 15 elementary schools in Switzerland randomly selected and assigned in a 4:3 ratio to an intervention (n=16) or control arm (n=12) after stratification for grade (first and fifth grade), from August 2005 to June 2006. 540 children, of whom 502 consented and presented at baseline. Children in the intervention arm (n=297) received a multi-component physical activity programme that included structuring the three existing physical education lessons each week and adding two additional lessons a week, daily short activity breaks, and physical activity homework. Children (n=205) and parents in the control group were not informed of an intervention group. For most outcome measures, the assessors were blinded. Primary outcome measures included body fat (sum of four skinfolds), aerobic fitness (shuttle run test), physical activity (accelerometry), and quality of life (questionnaires). Secondary outcome measures included body mass index and cardiovascular risk score (average z score of waist circumference, mean blood pressure, blood glucose, inverted high density lipoprotein cholesterol, and triglycerides). 498 children completed the baseline and follow-up assessments (mean age 6.9 (SD 0.3) years for first grade, 11.1 (0.5) years for fifth grade). After adjustment for grade, sex, baseline values, and clustering within classes, children in the intervention arm compared with controls showed more negative changes in the z score of the sum of four skinfolds (-0.12, 95 % confidence interval -0.21 to -0.03; P=0.009). Likewise, their z scores for aerobic fitness increased more favourably (0.17, 0.01 to 0.32; P=0.04), as did those for moderate-vigorous physical activity in school (1.19, 0.78 to 1.60; P<0.001), all day moderate-vigorous physical activity (0.44, 0.05 to 0.82; P=0.03), and total physical activity in school (0.92, 0.35 to 1.50; P=0.003). Z scores for overall daily physical activity (0.21, -0.21 to 0.63) and physical quality of life (0.42, -1.23 to 2.06) as well as psychological quality of life (0.59, -0.85 to 2.03) did not change significantly. A school based multi-component physical activity intervention including compulsory elements improved physical activity and fitness and reduced adiposity in children. Trial registration Current Controlled Trials ISRCTN15360785.

Differentiation and genetic analysis of Rhodnius prolixus and Rhodnius colombiensis by rDNA and RAPD amplification

Domiciliated Rhodnius prolixus and sylvatic R. colombiensis were analyzed in order to confirm their genetic divergence and verify the risk that the latter represents in the domiciliation process, and to provide tools for identifying the sources of possible reinfestation by triatomines in human dwellings allowing control programs to be undertaken. Comparison of random amplified polymorphic DNA amplification patterns and cluster analysis suggests reproductive discontinuity between the two species. The calculated statistical F value of 0.24 and effective migration rate of 0.6 individuals per generation are insufficient to maintain genetic homogeneity between them and confirm the absence of present genetic flow. R. colombiensis presents higher intrapopulation variability. Polymerase chain reaction of ribosomal DNA supports these findings. The low genetic flow between the two species implies that R. colombiensis do not represent an epidemiological risk for the domiciliary transmission of Trypanosoma cruzi in the Tolima Department. The lower variability of the domiciliated R. prolixus could result in greater susceptibility to the use of pesticides in control programs.

An R Library for Compositional Data Analysis in Archaeometry

Compositional data naturally arises from the scientific analysis of the chemicalcomposition of archaeological material such as ceramic and glass artefacts. Data of thistype can be explored using a variety of techniques, from standard multivariate methodssuch as principal components analysis and cluster analysis, to methods based upon theuse of log-ratios. The general aim is to identify groups of chemically similar artefactsthat could potentially be used to answer questions of provenance.This paper will demonstrate work in progress on the development of a documentedlibrary of methods, implemented using the statistical package R, for the analysis ofcompositional data. R is an open source package that makes available very powerfulstatistical facilities at no cost. We aim to show how, with the aid of statistical softwaresuch as R, traditional exploratory multivariate analysis can easily be used alongside, orin combination with, specialist techniques of compositional data analysis.The library has been developed from a core of basic R functionality, together withpurpose-written routines arising from our own research (for example that reported atCoDaWork'03). In addition, we have included other appropriate publicly availabletechniques and libraries that have been implemented in R by other authors. Availablefunctions range from standard multivariate techniques through to various approaches tolog-ratio analysis and zero replacement. We also discuss and demonstrate a smallselection of relatively new techniques that have hitherto been little-used inarchaeometric applications involving compositional data. The application of the libraryto the analysis of data arising in archaeometry will be demonstrated; results fromdifferent analyses will be compared; and the utility of the various methods discussed

Cluster recognition in spatial-temporal sequences: The case of forest fires

Forest fire sequences can be modelled as a stochastic point process where events are characterized by their spatial locations and occurrence in time. Cluster analysis permits the detection of the space/time pattern distribution of forest fires. These analyses are useful to assist fire-managers in identifying risk areas, implementing preventive measures and conducting strategies for an efficient distribution of the firefighting resources. This paper aims to identify hot spots in forest fire sequences by means of the space-time scan statistics permutation model (STSSP) and a geographical information system (GIS) for data and results visualization. The scan statistical methodology uses a scanning window, which moves across space and time, detecting local excesses of events in specific areas over a certain period of time. Finally, the statistical significance of each cluster is evaluated through Monte Carlo hypothesis testing. The case study is the forest fires registered by the Forest Service in Canton Ticino (Switzerland) from 1969 to 2008. This dataset consists of geo-referenced single events including the location of the ignition points and additional information. The data were aggregated into three sub-periods (considering important preventive legal dispositions) and two main ignition-causes (lightning and anthropogenic causes). Results revealed that forest fire events in Ticino are mainly clustered in the southern region where most of the population is settled. Our analysis uncovered local hot spots arising from extemporaneous arson activities. Results regarding the naturally-caused fires (lightning fires) disclosed two clusters detected in the northern mountainous area.

Community cluster of meningococcal disease by Neisseria meningitidis serogroup C in Andalusia, Spain, March to May 2011.

Between March and May of 2011, a cluster of three fatal cases of meningococcal sepsis occurred in Andalusia, Spain, in a municipality with a population of around 20,000 inhabitants. The cases were in their mid-teens to early thirties and were notified to the epidemiological surveillance system of Andalusia (Sistema de Vigilancia Epidemiológica de Andalucía, SVEA) during a 68-day period from March through May 2011. All three were infected with the same strain of Neisseria meningitidis serogroup C genosubtype VR1:5-1;VR2:10-8. None of the cases had been previously vaccinated against N. meningitidis serogroup C. Antibiotic post-exposure chemoprophylaxis was administered to close contacts of every diagnosed case. Once the cluster was confirmed, the local population was informed through the media about the control measures taken by the health authorities. The vaccination history against N. meningitidis serogroup C of the population under 25 years-old in the municipality was checked. Vaccination was offered to unimmunised individuals younger than 25 years of age and an additional dose of vaccine was offered to those who had been vaccinated between 2000 and 2006 with a vaccination schedule of three doses before the first year of age. No further cases occurred since the beginning of these actions.

Polymorphisms in DNA-Repair Genes in a Cohort of Prostate Cancer Patients from Different Areas in Spain: Heterogeneity between Populations as a Confounding Factor in Association Studies.

BACKGROUND Differences in the distribution of genotypes between individuals of the same ethnicity are an important confounder factor commonly undervalued in typical association studies conducted in radiogenomics. OBJECTIVE To evaluate the genotypic distribution of SNPs in a wide set of Spanish prostate cancer patients for determine the homogeneity of the population and to disclose potential bias. DESIGN SETTING AND PARTICIPANTS A total of 601 prostate cancer patients from Andalusia, Basque Country, Canary and Catalonia were genotyped for 10 SNPs located in 6 different genes associated to DNA repair: XRCC1 (rs25487, rs25489, rs1799782), ERCC2 (rs13181), ERCC1 (rs11615), LIG4 (rs1805388, rs1805386), ATM (rs17503908, rs1800057) and P53 (rs1042522). The SNP genotyping was made in a Biotrove OpenArray® NT Cycler. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS Comparisons of genotypic and allelic frequencies among populations, as well as haplotype analyses were determined using the web-based environment SNPator. Principal component analysis was made using the SnpMatrix and XSnpMatrix classes and methods implemented as an R package. Non-supervised hierarchical cluster of SNP was made using MultiExperiment Viewer. RESULTS AND LIMITATIONS We observed that genotype distribution of 4 out 10 SNPs was statistically different among the studied populations, showing the greatest differences between Andalusia and Catalonia. These observations were confirmed in cluster analysis, principal component analysis and in the differential distribution of haplotypes among the populations. Because tumor characteristics have not been taken into account, it is possible that some polymorphisms may influence tumor characteristics in the same way that it may pose a risk factor for other disease characteristics. CONCLUSION Differences in distribution of genotypes within different populations of the same ethnicity could be an important confounding factor responsible for the lack of validation of SNPs associated with radiation-induced toxicity, especially when extensive meta-analysis with subjects from different countries are carried out.

Spatial analysis of the relationship between mortality from cardiovascular and cerebrovascular disease and drinking water hardness.

Previously published scientific papers have reported a negative correlation between drinking water hardness and cardiovascular mortality. Some ecologic and case-control studies suggest the protective effect of calcium and magnesium concentration in drinking water. In this article we present an analysis of this protective relationship in 538 municipalities of Comunidad Valenciana (Spain) from 1991-1998. We used the Spanish version of the Rapid Inquiry Facility (RIF) developed under the European Environment and Health Information System (EUROHEIS) research project. The strategy of analysis used in our study conforms to the exploratory nature of the RIF that is used as a tool to obtain quick and flexible insight into epidemiologic surveillance problems. This article describes the use of the RIF to explore possible associations between disease indicators and environmental factors. We used exposure analysis to assess the effect of both protective factors--calcium and magnesium--on mortality from cerebrovascular (ICD-9 430-438) and ischemic heart (ICD-9 410-414) diseases. This study provides statistical evidence of the relationship between mortality from cardiovascular diseases and hardness of drinking water. This relationship is stronger in cerebrovascular disease than in ischemic heart disease, is more pronounced for women than for men, and is more apparent with magnesium than with calcium concentration levels. Nevertheless, the protective nature of these two factors is not clearly established. Our results suggest the possibility of protectiveness but cannot be claimed as conclusive. The weak effects of these covariates make it difficult to separate them from the influence of socioeconomic and environmental factors. We have also performed disease mapping of standardized mortality ratios to detect clusters of municipalities with high risk. Further standardization by levels of calcium and magnesium in drinking water shows changes in the maps when we remove the effect of these covariates.

Diversity of the bacterial community in the surface soil of a pear orchard based on 16S rRNA gene analysis

A cultivation-independent approach based on polymerase chain reaction (PCR)-amplified partial small subunit rRNA genes was used to characterize bacterial populations in the surface soil of a commercial pear orchard consisting of different pear cultivars during two consecutive growing seasons. Pyrus communis L. cvs Blanquilla, Conference, and Williams are among the most widely cultivated cultivars in Europe and account for the majority of pear production in Northeastern Spain. To assess the heterogeneity of the community structure in response to environmental variables and tree phenology, bacterial populations were examined using PCR-denaturing gradient gel electrophoresis (DGGE) followed by cluster analysis of the 16S ribosomal DNA profiles by means of the unweighted pair group method with arithmetic means. Similarity analysis of the band patterns failed to identify characteristic fingerprints associated with the pear cultivars. Both environmentally and biologically based principal-component analyses showed that the microbial communities changed significantly throughout the year depending on temperature and, to a lesser extent, on tree phenology and rainfall. Prominent DGGE bands were excised and sequenced to gain insight into the identities of the predominant bacterial populations. Most DGGE band sequences were related to bacterial phyla, such as Bacteroidetes, Cyanobacteria, Acidobacteria, Proteobacteria, Nitrospirae, and Gemmatimonadetes, previously associated with typical agronomic crop environments

Proteomic analysis of peach fruit mesocarp softening and chilling injury using difference gel electrophoresis (DIGE)

Background: Peach fruit undergoes a rapid softening process that involves a number of metabolic changes. Storing fruit at low temperatures has been widely used to extend its postharvest life. However, this leads to undesired changes, such as mealiness and browning, which affect the quality of the fruit. In this study, a 2-D DIGE approach was designed to screen for differentially accumulated proteins in peach fruit during normal softening as well as under conditions that led to fruit chilling injury. Results:The analysis allowed us to identify 43 spots -representing about 18% of the total number analyzed- that show statistically significant changes. Thirty-nine of the proteins could be identified by mass spectrometry. Some of the proteins that changed during postharvest had been related to peach fruit ripening and cold stress in the past. However, we identified other proteins that had not been linked to these processes. A graphical display of the relationship between the differentially accumulated proteins was obtained using pairwise average-linkage cluster analysis and principal component analysis. Proteins such as endopolygalacturonase, catalase, NADP-dependent isocitrate dehydrogenase, pectin methylesterase and dehydrins were found to be very important for distinguishing between healthy and chill injured fruit. A categorization of the differentially accumulated proteins was performed using Gene Ontology annotation. The results showed that the 'response to stress', 'cellular homeostasis', 'metabolism of carbohydrates' and 'amino acid metabolism' biological processes were affected the most during the postharvest. Conclusions: Using a comparative proteomic approach with 2-D DIGE allowed us to identify proteins that showed stage-specific changes in their accumulation pattern. Several proteins that are related to response to stress, cellular homeostasis, cellular component organization and carbohydrate metabolism were detected as being differentially accumulated. Finally, a significant proportion of the proteins identified had not been associated with softening, cold storage or chilling injury-altered fruit before; thus, comparative proteomics has proven to be a valuable tool for understanding fruit softening and postharvest.

Designing phase III or IV trials for vaccines: choosing between individual or cluster randomised trial designs.

The choice of design between individual randomisation, cluster or pseudo-cluster randomisation is often made difficult. Clear methodological guidelines have been given for trials in general practice, but not for vaccine trials. This article proposes a decisional flow-chart to choose the most adapted design for evaluating the effectiveness of a vaccine in large-scale studies. Six criteria have been identified: importance of herd immunity or herd protection, ability to delimit epidemiological units, homogeneity of transmission probability across sub-populations, population's acceptability of randomisation, availability of logistical resources, and estimated sample size. This easy to use decisional method could help sponsors, trial steering committees and ethical committees adopt the most suitable design.

Analysis of the genetic diversity in Metopolophium dirhodum (Walker) (Hemiptera, Aphididae) by RAPD markers

The emergence of host-races within aphids may constitute an obstacle to pest management by means of plant resistance. There are examples of host-races within cereals aphids, but their occurrence in Rose Grain Aphid, Metopolophium dirhodum (Walker, 1849), has not been reported yet. In this work, RAPD markers were used to assess effects of the hosts and geographic distance on the genetic diversity of M. dirhodum lineages. Twenty-three clones were collected on oats and wheat in twelve localitites of southern Brazil. From twenty-seven primers tested, only four primers showed polymorphisms. Fourteen different genotypes were revealed by cluster analysis. Five genotypes were collected only on wheat; seven only on oats and two were collected in both hosts. Genetic and geographical distances among all clonal lineages were not correlated. Analysis of molecular variance showed that some molecular markers are not randomly distributed among clonal lineages collected on oats and on wheat. These results suggest the existence of host-races within M. dirhodum, which should be further investigated using a combination of ecological and genetic data.

Multivariate analysis and geostatistics of the fertility of a humic rhodic hapludox under coffee cultivation

The spatial variability of soil and plant properties exerts great influence on the yeld of agricultural crops. This study analyzed the spatial variability of the fertility of a Humic Rhodic Hapludox with Arabic coffee, using principal component analysis, cluster analysis and geostatistics in combination. The experiment was carried out in an area under Coffea arabica L., variety Catucai 20/15 - 479. The soil was sampled at a depth 0.20 m, at 50 points of a sampling grid. The following chemical properties were determined: P, K+, Ca2+, Mg2+, Na+, S, Al3+, pH, H + Al, SB, t, T, V, m, OM, Na saturation index (SSI), remaining phosphorus (P-rem), and micronutrients (Zn, Fe, Mn, Cu and B). The data were analyzed with descriptive statistics, followed by principal component and cluster analyses. Geostatistics were used to check and quantify the degree of spatial dependence of properties, represented by principal components. The principal component analysis allowed a dimensional reduction of the problem, providing interpretable components, with little information loss. Despite the characteristic information loss of principal component analysis, the combination of this technique with geostatistical analysis was efficient for the quantification and determination of the structure of spatial dependence of soil fertility. In general, the availability of soil mineral nutrients was low and the levels of acidity and exchangeable Al were high.

How Much Does It Cost? Optimization of Costs in Sequence Analysis of Social Science Data

One major methodological problem in analysis of sequence data is the determination of costs from which distances between sequences are derived. Although this problem is currently not optimally dealt with in the social sciences, it has some similarity with problems that have been solved in bioinformatics for three decades. In this article, the authors propose an optimization of substitution and deletion/insertion costs based on computational methods. The authors provide an empirical way of determining costs for cases, frequent in the social sciences, in which theory does not clearly promote one cost scheme over another. Using three distinct data sets, the authors tested the distances and cluster solutions produced by the new cost scheme in comparison with solutions based on cost schemes associated with other research strategies. The proposed method performs well compared with other cost-setting strategies, while it alleviates the justification problem of cost schemes.