On cluster analysis of complex and heterogeneous data

3rd SMTDA Conference Proceedings, 11-14 June 2014, Lisbon Portugal.

Some reflections on acquisition, processing and analysis of statistical data in forest soils

Beyond the classical statistical approaches (determination of basic statistics, regression analysis, ANOVA, etc.) a new set of applications of different statistical techniques has increasingly gained relevance in the analysis, processing and interpretation of data concerning the characteristics of forest soils. This is possible to be seen in some of the recent publications in the context of Multivariate Statistics. These new methods require additional care that is not always included or refered in some approaches. In the particular case of geostatistical data applications it is necessary, besides to geo-reference all the data acquisition, to collect the samples in regular grids and in sufficient quantity so that the variograms can reflect the spatial distribution of soil properties in a representative manner. In the case of the great majority of Multivariate Statistics techniques (Principal Component Analysis, Correspondence Analysis, Cluster Analysis, etc.) despite the fact they do not require in most cases the assumption of normal distribution, they however need a proper and rigorous strategy for its utilization. In this work, some reflections about these methodologies and, in particular, about the main constraints that often occur during the information collecting process and about the various linking possibilities of these different techniques will be presented. At the end, illustrations of some particular cases of the applications of these statistical methods will also be presented.

The development of the N1 and N2 components in auditory oddball paradigms: a systematic review with narrative analysis and suggested normative values

Auditory event-related potentials (AERPs) are widely used in diverse fields of today’s neuroscience, concerning auditory processing, speech perception, language acquisition, neurodevelopment, attention and cognition in normal aging, gender, developmental, neurologic and psychiatric disorders. However, its transposition to clinical practice has remained minimal. Mainly due to scarce literature on normative data across age, wide spectrumof results, variety of auditory stimuli used and to different neuropsychological meanings of AERPs components between authors. One of the most prominent AERP components studied in last decades was N1, which reflects auditory detection and discrimination. Subsequently, N2 indicates attention allocation and phonological analysis. The simultaneous analysis of N1 and N2 elicited by feasible novelty experimental paradigms, such as auditory oddball, seems an objective method to assess central auditory processing. The aim of this systematic review was to bring forward normative values for auditory oddball N1 and N2 components across age. EBSCO, PubMed, Web of Knowledge and Google Scholarwere systematically searched for studies that elicited N1 and/or N2 by auditory oddball paradigm. A total of 2,764 papers were initially identified in the database, of which 19 resulted from hand search and additional references, between 1988 and 2013, last 25 years. A final total of 68 studiesmet the eligibility criteria with a total of 2,406 participants from control groups for N1 (age range 6.6–85 years; mean 34.42) and 1,507 for N2 (age range 9–85 years; mean 36.13). Polynomial regression analysis revealed thatN1latency decreases with aging at Fz and Cz,N1 amplitude at Cz decreases from childhood to adolescence and stabilizes after 30–40 years and at Fz the decrement finishes by 60 years and highly increases after this age. Regarding N2, latency did not covary with age but amplitude showed a significant decrement for both Cz and Fz. Results suggested reliable normative values for Cz and Fz electrode locations; however, changes in brain development and components topography over age should be considered in clinical practice.

Da sobrevivência física à sobrevivência mental : aspectos do funcionamento psíquico em ex-combatentes da guerra colonial portuguesa

Tese de Doutoramento apresentada ao ISPA - Instituto Universitário

Morphometric analysis of Capitella capitata (Polychaeta, Capitellidae)

Accurate size measurements are fundamental in characterizing the population structure and secondary production of a species. The purpose of this study was to determine the best morphometric parameter to estimate the size of individuals of Capitella capitata (Fabricius, 1780). The morphometric analysis was applied to individuals collected in the intertidal zones of two beaches on the northern coast of the state of São Paulo, Brazil: São Francisco and Araçá. The following measurements were taken: the width and length (height) of the 4th, 5th and 7th setigers, and the length of the thoracic region (first nine setigers). The area and volume of these setigers were calculated and a linear regression analysis was applied to the data. The data were log-transformed to fit the allometric equation y = ax b into a straight line (log y = log a + b * log x). The measurements which best correlated with the thoracic length in individuals from both beaches were the length of setiger 5 (r² = 0.722; p<0.05 in São Francisco and r² = 0.795; p<0.05 in Araçá) and the area of setiger 7 (r² = 0.705; p<0.05 in São Francisco and r² = 0.634; p<0.05 in Araçá). According to these analyses, the length of setiger 5 and/or the area of setiger 7 are the best parameters to evaluate the growth of individuals of C. capitata.

Challenging recommended oral and intravenous voriconazole doses for improved efficacy and safety: population pharmacokinetics-based analysis of adult patients with invasive fungal infections.

BACKGROUND: Recommended oral voriconazole (VRC) doses are lower than intravenous doses. Because plasma concentrations impact efficacy and safety of therapy, optimizing individual drug exposure may improve these outcomes. METHODS: A population pharmacokinetic analysis (NONMEM) was performed on 505 plasma concentration measurements involving 55 patients with invasive mycoses who received recommended VRC doses. RESULTS: A 1-compartment model with first-order absorption and elimination best fitted the data. VRC clearance was 5.2 L/h, the volume of distribution was 92 L, the absorption rate constant was 1.1 hour(-1), and oral bioavailability was 0.63. Severe cholestasis decreased VRC elimination by 52%. A large interpatient variability was observed on clearance (coefficient of variation [CV], 40%) and bioavailability (CV 84%), and an interoccasion variability was observed on bioavailability (CV, 93%). Lack of response to therapy occurred in 12 of 55 patients (22%), and grade 3 neurotoxicity occurred in 5 of 55 patients (9%). A logistic multivariate regression analysis revealed an independent association between VRC trough concentrations and probability of response or neurotoxicity by identifying a therapeutic range of 1.5 mg/L (>85% probability of response) to 4.5 mg/L (<15% probability of neurotoxicity). Population-based simulations with the recommended 200 mg oral or 300 mg intravenous twice-daily regimens predicted probabilities of 49% and 87%, respectively, for achievement of 1.5 mg/L and of 8% and 37%, respectively, for achievement of 4.5 mg/L. With 300-400 mg twice-daily oral doses and 200-300 mg twice-daily intravenous doses, the predicted probabilities of achieving the lower target concentration were 68%-78% for the oral regimen and 70%-87% for the intravenous regimen, and the predicted probabilities of achieving the upper target concentration were 19%-29% for the oral regimen and 18%-37% for the intravenous regimen. CONCLUSIONS: Higher oral than intravenous VRC doses, followed by individualized adjustments based on measured plasma concentrations, improve achievement of the therapeutic target that maximizes the probability of therapeutic response and minimizes the probability of neurotoxicity. These findings challenge dose recommendations for VRC.

Okun’s law – A meta analysis

This paper seeks to identify whether there is a representative empirical Okun’s Law coefficient (OLC) and to measure its size. We carry out a meta regression analysis on a sample of 269 estimates of the OLC to uncover reasons for differences in empirical results and to estimate the ‘true’ OLC. On statistical (and other) grounds, we find it appropriate to investigate two separate subsamples, using respectively (some measure of) unemployment or output as dependent variable. Our results can be summarized as follows. First, there is evidence of type II publication bias in both sub-samples, but a type I bias is present only among the papers using some measure of unemployment as the dependent variable. Second, after correction for publication bias, authentic and statistically significant OLC effects are present in both sub-samples. Third, bias-corrected estimated true OLCs are significantly lower (in absolute value) with models using some measure of unemployment as the dependent variable. Using a bivariate MRA approach, the estimated true effects are -0.25 for the unemployment sub-sample and -0.61 for the output-sub sample; with a multivariate MRA methodology, the estimated true effects are -0.40 and -1.02 for the unemployment and the output-sub samples respectively.

The application of chemometrics on Infrared and Raman spectra as a tool for the forensic analysis of paints

The aim of this work is to evaluate the capabilities and limitations of chemometric methods and other mathematical treatments applied on spectroscopic data and more specifically on paint samples. The uniqueness of the spectroscopic data comes from the fact that they are multivariate - a few thousands variables - and highly correlated. Statistical methods are used to study and discriminate samples. A collection of 34 red paint samples was measured by Infrared and Raman spectroscopy. Data pretreatment and variable selection demonstrated that the use of Standard Normal Variate (SNV), together with removal of the noisy variables by a selection of the wavelengths from 650 to 1830 cm−1 and 2730-3600 cm−1, provided the optimal results for infrared analysis. Principal component analysis (PCA) and hierarchical clusters analysis (HCA) were then used as exploratory techniques to provide evidence of structure in the data, cluster, or detect outliers. With the FTIR spectra, the Principal Components (PCs) correspond to binder types and the presence/absence of calcium carbonate. 83% of the total variance is explained by the four first PCs. As for the Raman spectra, we observe six different clusters corresponding to the different pigment compositions when plotting the first two PCs, which account for 37% and 20% respectively of the total variance. In conclusion, the use of chemometrics for the forensic analysis of paints provides a valuable tool for objective decision-making, a reduction of the possible classification errors, and a better efficiency, having robust results with time saving data treatments.

Mixture of bivariate Poisson regression models with an application to insurance

In a recent paper Bermúdez [2009] used bivariate Poisson regression models for ratemaking in car insurance, and included zero-inflated models to account for the excess of zeros and the overdispersion in the data set. In the present paper, we revisit this model in order to consider alternatives. We propose a 2-finite mixture of bivariate Poisson regression models to demonstrate that the overdispersion in the data requires more structure if it is to be taken into account, and that a simple zero-inflated bivariate Poisson model does not suffice. At the same time, we show that a finite mixture of bivariate Poisson regression models embraces zero-inflated bivariate Poisson regression models as a special case. Additionally, we describe a model in which the mixing proportions are dependent on covariates when modelling the way in which each individual belongs to a separate cluster. Finally, an EM algorithm is provided in order to ensure the models’ ease-of-fit. These models are applied to the same automobile insurance claims data set as used in Bermúdez [2009] and it is shown that the modelling of the data set can be improved considerably.

A logistic regression approach to estimating customer profit loss due to lapses in insurance

This article focuses on business risk management in the insurance industry. A methodology for estimating the profit loss caused by each customer in the portfolio due to policy cancellation is proposed. Using data from a European insurance company, customer behaviour over time is analyzed in order to estimate the probability of policy cancelation and the resulting potential profit loss due to cancellation. Customers may have up to two different lines of business contracts: motor insurance and other diverse insurance (such as, home contents, life or accident insurance). Implications for understanding customer cancellation behaviour as the core of business risk management are outlined.

Die Säuglingssterblichkeit in der Schweiz: multivariate Betrachtung. [Perinatal mortality in Switzerland: a multivariate analysis]

A regression analysis using a linked file of all Swiss births und perinatal deaths 1979-1981 showed a significant relation between birthweight and canton. Sex of infant and multiplicity of birth were significant, too. For live births, marital and socio-economic status of mother and father relate to birthweight. Logistic regressions brought out relationships between the risk of stillbirth and occupation of father, nationality and marital status of mother, apart from birthweight. For live births, only sex and (weakly) marital status and rank of the child were influencial after correction for birthweight.

Robust regression in biological assay: application to the evaluation of alternative experimental techniques.

Robust Huber type regression and testing of linear hypotheses are adapted to statistical analysis of parallel line and slope ratio assays. They are applied in the evaluation of results of several experiments carried out in order to compare and validate alternatives to animal experimentation based on embryo and cell cultures. Computational procedures necessary for the application of robust methods of analysis used the conversational statistical package ROBSYS. Special commands for the analysis of parallel line and slope ratio assays have been added to ROBSYS.

Factors associated with use of community mental health services by schizophrenia patients using multilevel analysis

BACKGROUND Persons with schizophrenia and related disorders may be particularly sensitive to a number of determinants of service use, including those related with illness, socio-demographic characteristics and organizational factors. The objective of this study is to identify factors associated with outpatient contacts at community mental health services of patients with schizophrenia or related disorders. METHODS This cross-sectional study analyzed 1097 patients. The main outcome measure was the total number of outpatient consultations during one year. Independent variables were related to socio-demographic, clinical and use of service factors. Data were collected from clinical records. RESULTS The multilevel linear regression model explained 46.35% of the variance. Patients with significantly more contacts with ambulatory services were not working and were receiving welfare benefits (p = 0.02), had no formal education (p = 0.02), had a global level of severity of two or three (four being the most severe) (p < 0.001), with one or more inpatient admissions (p < 0.001), and in contact with both types of professional (nurses and psychiatrists) (p < 0.001). The patients with the fewest ambulatory contacts were those with diagnoses of persistent delusional disorders (p = 0.04) and those who were attended by four of the 13 psychiatrists (p < 0.001). CONCLUSIONS As expected, the variables that explained the use of community service could be viewed as proxies for severity of illness. The most surprising finding, however, was that a group of four psychiatrists was also independently associated with use of ambulatory services by patients with schizophrenia or related disorders. More research is needed to carefully examine how professional support networks interact to affect use of mental health.

Effectiveness of two interventions based on improving patient-practitioner communication on diabetes self-management in patients with low educational level: study protocol of a clustered randomized trial in primary care.

BACKGROUND In the last decades the presence of social inequalities in diabetes care has been observed in multiple countries, including Spain. These inequalities have been at least partially attributed to differences in diabetes self-management behaviours. Communication problems during medical consultations occur more frequently to patients with a lower educational level. The purpose of this cluster randomized trial is to determine whether an intervention implemented in a General Surgery, based in improving patient-provider communication, results in a better diabetes self-management in patients with lower educational level. A secondary objective is to assess whether telephone reinforcement enhances the effect of such intervention. We report the design and implementation of this on-going study. METHODS/DESIGN The study is being conducted in a General Practice located in a deprived neighbourhood of Granada, Spain. Diabetic patients 18 years old or older with a low educational level and inadequate glycaemic control (HbA1c > 7%) were recruited. General Practitioners (GPs) were randomised to three groups: intervention A, intervention B and control group. GPs allocated to intervention groups A and B received training in communication skills and are providing graphic feedback about glycosylated haemoglobin levels. Patients whose GPs were allocated to group B are additionally receiving telephone reinforcement whereas patients from the control group are receiving usual care. The described interventions are being conducted during 7 consecutive medical visits which are scheduled every three months. The main outcome measure will be HbA1c; blood pressure, lipidemia, body mass index and waist circumference will be considered as secondary outcome measures. Statistical analysis to evaluate the effectiveness of the interventions will include multilevel regression analysis with three hierarchical levels: medical visit level, patient level and GP level. DISCUSSION The results of this study will provide new knowledge about possible strategies to promote a better diabetes self-management in a particularly vulnerable group. If effective, this low cost intervention will have the potential to be easily incorporated into routine clinical practice, contributing to decrease health inequalities in diabetic patients. TRIAL REGISTRATION Clinical Trials U.S. National Institutes of Health, NCT01849731.

Trabecular bone score (TBS) the new parameter of 2D texture analysis for the evaluation of 3D bone micro architecture status

X-ray is a technology that is used for numerous applications in the medical field. The process of X-ray projection gives a 2-dimension (2D) grey-level texture from a 3- dimension (3D) object. Until now no clear demonstration or correlation has positioned the 2D texture analysis as a valid indirect evaluation of the 3D microarchitecture. TBS is a new texture parameter based on the measure of the experimental variogram. TBS evaluates the variation between 2D image grey-levels. The aim of this study was to evaluate existing correlations between 3D bone microarchitecture parameters - evaluated from μCT reconstructions - and the TBS value, calculated on 2D projected images. 30 dried human cadaveric vertebrae were acquired on a micro-scanner (eXplorer Locus, GE) at isotropic resolution of 93 μm. 3D vertebral body models were used. The following 3D microarchitecture parameters were used: Bone volume fraction (BV/TV), Trabecular thickness (TbTh), trabecular space (TbSp), trabecular number (TbN) and connectivity density (ConnD). 3D/2D projections has been done by taking into account the Beer-Lambert Law at X-ray energy of 50, 100, 150 KeV. TBS was assessed on 2D projected images. Correlations between TBS and the 3D microarchitecture parameters were evaluated using a linear regression analysis. Paired T-test is used to assess the X-ray energy effects on TBS. Multiple linear regressions (backward) were used to evaluate relationships between TBS and 3D microarchitecture parameters using a bootstrap process. BV/TV of the sample ranged from 18.5 to 37.6% with an average value at 28.8%. Correlations' analysis showedthat TBSwere strongly correlatedwith ConnD(0.856≤r≤0.862; p<0.001),with TbN (0.805≤r≤0.810; p<0.001) and negatively with TbSp (−0.714≤r≤−0.726; p<0.001), regardless X-ray energy. Results show that lower TBS values are related to "degraded" microarchitecture, with low ConnD, low TbN and a high TbSp. The opposite is also true. X-ray energy has no effect onTBS neither on the correlations betweenTBS and the 3Dmicroarchitecture parameters. In this study, we demonstrated that TBS was significantly correlated with 3D microarchitecture parameters ConnD and TbN, and negatively with TbSp, no matter what X-ray energy has been used. This article is part of a Special Issue entitled ECTS 2011. Disclosure of interest: None declared.