Fuzzy coding in constrained ordinations

Canonical correspondence analysis and redundancy analysis are two methods of constrained ordination regularly used in the analysis of ecological data when several response variables (for example, species abundances) are related linearly to several explanatory variables (for example, environmental variables, spatial positions of samples). In this report I demonstrate the advantages of the fuzzy coding of explanatory variables: first, nonlinear relationships can be diagnosed; second, more variance in the responses can be explained; and third, in the presence of categorical explanatory variables (for example, years, regions) the interpretation of the resulting triplot ordination is unified because all explanatory variables are measured at a categorical level.

Quetiapine affects neuropeptide Y and corticotropin-releasing hormone in cerebrospinal fluid from schizophrenia patients: relationship to depression and anxiety symptoms and to treatment response.

Cumulative evidence indicates that neuropeptides play a role in the pathophysiology of schizophrenia. Early data showed increased neuropeptide Y (NPY) in cerebrospinal fluid (CSF) from schizophrenia patients and data from rodents show that antipsychotic drugs modulate NPY levels in and release from selected rat brain regions. In view of these findings we investigated whether the atypical antipsychotic quetiapine, originally used as an antipsychotic but subsequently shown to be efficient also in major depressive disorder and in both poles of bipolar disorder, would affect NPY-like immunoreactivity (-LI), and corticotropin-releasing hormone (CRH)-LI levels in CSF of schizophrenia patients. NPY-LI and CRH-LI in CSF were determined in 22 patients with schizophrenia. Lumbar puncture was performed at baseline and again after 4 wk of quetiapine treatment (600 mg/d). Patients were assessed with the Positive and Negative Syndrome Scale (PANSS) at baseline and at weekly intervals. Quetiapine treatment was associated with a significant increase in NPY-LI (p<0.001) and decrease in CRH-LI (p<0.01). Stepwise multiple regression analysis revealed that ΔNPY-LI and ΔCRH-LI levels predicted 63% (p<0.001) of the variability of the ΔPANSS total score, ΔNPY-LI 42% (p<0.05) of the ΔPANSS anxiety items (G2) and ΔCRH-LI 40% (p=0.05) of the ΔPANSS depression items (G6). These results suggest that while quetiapine's effects on monoamines are probably related to its antipsychotic properties, the modulation of NPY and CRH accounts for its antidepressant and anxiolytic effects and can be markers of response.

Contribution biplots

In order to interpret the biplot it is necessary to know which points usually variables are the ones that are important contributors to the solution, and this information is available separately as part of the biplot s numerical results. We propose a new scaling of the display, called the contribution biplot, which incorporates this diagnostic directly into the graphical display, showing visually the important contributors and thus facilitating the biplot interpretation and often simplifying the graphical representation considerably. The contribution biplot can be applied to a wide variety of analyses such as correspondence analysis, principal component analysis, log-ratio analysis and the graphical results of a discriminant analysis/MANOVA, in fact to any method based on the singular-value decomposition. In the contribution biplot one set of points, usually the rows of the data matrix, optimally represent the spatial positions of the cases or sample units, according to some distance measure that usually incorporates some form of standardization unless all data are comparable in scale. The other set of points, usually the columns, is represented by vectors that are related to their contributions to the low-dimensional solution. A fringe benefit is that usually only one common scale for row and column points is needed on the principal axes, thus avoiding the problem of enlarging or contracting the scale of one set of points to make the biplot legible. Furthermore, this version of the biplot also solves the problem in correspondence analysis of low-frequency categories that are located on the periphery of the map, giving the false impression that they are important, when they are in fact contributing minimally to the solution.

Dynamic perceptual mapping

Perceptual maps have been used for decades by market researchers to illuminatethem about the similarity between brands in terms of a set of attributes, to position consumersrelative to brands in terms of their preferences, or to study how demographic and psychometricvariables relate to consumer choice. Invariably these maps are two-dimensional and static. Aswe enter the era of electronic publishing, the possibilities for dynamic graphics are opening up.We demonstrate the usefulness of introducing motion into perceptual maps through fourexamples. The first example shows how a perceptual map can be viewed in three dimensions,and the second one moves between two analyses of the data that were collected according todifferent protocols. In a third example we move from the best view of the data at the individuallevel to one which focuses on between-group differences in aggregated data. A final exampleconsiders the case when several demographic variables or market segments are available foreach respondent, showing an animation with increasingly detailed demographic comparisons.These examples of dynamic maps use several data sets from marketing and social scienceresearch.

Are Americans' musical preferences more omnivores today?

Although we found a general trend favouring the omnivorousness thesis, as soon as we adjusted it to a set of structural factors and consumers tastes it was clear that this was caused by elitist inclusive omnivores who had increased the scope of their tastes. In general, younger cohorts were becoming less omnivorous, nevertheless, they were also becoming more educated and had greater to higher levels of inc ome, making the youth moreomnivorous. As expected, upscale consumers set limits on their popular taste: musicalgenres, whose audiences had educational levels below the mean profile were less preferredby upscale respondents. In spite of this, as time passed, some popular brows gained socialstatus.

Nestmate recognition and the genetic relatedness of nests in the ant Formica pratensis

Genetic relatedness of the mound-building ant Formica pratensis was determined by means of microsatellite DNA polymorphism, and its impact on nestmate recognition was tested in a population in Southern Sweden (Oeland). Recognition between nests was measured by testing aggression levels between single pairs of workers. The genetic distances of nests (Nei's genetic distance) and the spatial distance of nests were correlated and both showed a strong relation to the aggression behavior. Multiple regression analysis revealed a stronger impact of genetic relatedness rather than spatial distances on aggression behavior. Neighbouring nests were more closely related than distant nests, which may reflect budding as a possible spreading mechanism. The genetic distance data showed that nestmate recognition was strongly genetically influenced in F. pratensis.

Weighted metric multidimensional scaling

This paper establishes a general framework for metric scaling of any distance measure between individuals based on a rectangular individuals-by-variables data matrix. The method allows visualization of both individuals and variables as well as preserving all the good properties of principal axis methods such as principal components and correspondence analysis, based on the singular-value decomposition, including the decomposition of variance into components along principal axes which provide the numerical diagnostics known as contributions. The idea is inspired from the chi-square distance in correspondence analysis which weights each coordinate by an amount calculated from the margins of the data table. In weighted metric multidimensional scaling (WMDS) we allow these weights to be unknown parameters which are estimated from the data to maximize the fit to the original distances. Once this extra weight-estimation step is accomplished, the procedure follows the classical path in decomposing a matrix and displaying its rows and columns in biplots.

Forest fires: Evaluation of government intervention measures

The purpose of this paper is to examine the relation between government measures, volunteer participation, climate variables and forest fires. A number of studies have related forest fires to causes of ignition, to fire history in one area, to the type of vegetation and weathercharacteristics or to community institutions, but there is little research on the relation between fire production and government prevention and extinction measures from a policy evaluation perspective.An observational approach is first applied to select forest fires in the north east of Spain. Taking a selection of fires with a certain size, a multiple regression analysis is conducted to find significant relations between policy instruments under the control of the government and the number of hectares burn in each case, controlling at the same time the effect of weather conditions and other context variables. The paper brings evidence on the effects of simultaneity and the relevance of recurring to army soldiers in specific days with extraordinary high simultaneity. The analysis also brings light on the effectiveness of twopreventive policies and of helicopters for extinction tasks.

Resting metabolic rate in six- to ten-year-old obese and nonobese children.

The resting metabolic rate (RMR) and body composition of 130 obese and nonobese prepubertal children, aged 6 to 10 years, were assessed by indirect calorimetry and skin-fold thickness, respectively. The mean (+/- SD) RMR was 4619 +/- 449 kJ.day-1 (164 +/- 31 kJ.kg body weight-1 x day-1) in the 62 boys and 4449 +/- 520 kJ.day-1 (147 +/- 32 kJ.kg body weight-1 x day-1) in the 68 girls. Fat-free mass was the best single predictor of RMR (R2 = 0.64; p < 0.001). Step-down multiple regression analysis, with independent variables such as age, gender, weight, and height, allowed several RMR predictive equations to be developed. An equation for boys is as follows: RMR (kJ.day-1) = 1287 + 28.6 x Weight(kg) + 23.6 x Height(cm) - 69.1 x Age(yr) (R2 = 0.58; p < 0.001). An equation for girls is as follows: RMR (kJ.day-1 = 1552 + 35.8 x Weight (kg) + 15.6 x Height (cm) - 36.3 x Age (yr) (R2 = 0.69; p < 0.001). Comparison between the measured RMR and that predicted by currently used formulas showed that most of these equations tended to overestimate the RMR of both genders, especially in overweight children.

Measuring subcompositional incoherence

Subcompositional coherence is a fundamental property of Aitchison s approach to compositional data analysis, and is the principal justification for using ratios of components. We maintain, however, that lack of subcompositional coherence, that is incoherence, can be measured in an attempt to evaluate whether any given technique is close enough, for all practical purposes, to being subcompositionally coherent. This opens up the field to alternative methods, which might be better suited to cope with problems such as data zeros and outliers, while being only slightly incoherent. The measure that we propose is based on the distance measure between components. We show that the two-part subcompositions, which appear to be the most sensitive to subcompositional incoherence, can be used to establish a distance matrix which can be directly compared with the pairwise distances in the full composition. The closeness of these two matrices can be quantified using a stress measure that is common in multidimensional scaling, providing a measure of subcompositional incoherence. The approach is illustrated using power-transformed correspondence analysis, which has already been shown to converge to log-ratio analysis as the power transform tends to zero.

Country effects in ISSP-1993 environmental data: Comparison of SEM approaches

Structural equation models (SEM) are commonly used to analyze the relationship between variables some of which may be latent, such as individual ``attitude'' to and ``behavior'' concerning specific issues. A number of difficulties arise when we want to compare a large number of groups, each with large sample size, and the manifest variables are distinctly non-normally distributed. Using an specific data set, we evaluate the appropriateness of the following alternative SEM approaches: multiple group versus MIMIC models, continuous versus ordinal variables estimation methods, and normal theory versus non-normal estimation methods. The approaches are applied to the ISSP-1993 Environmental data set, with the purpose of exploring variation in the mean level of variables of ``attitude'' to and ``behavior''concerning environmental issues and their mutual relationship across countries. Issues of both theoretical and practical relevance arise in the course of this application.

Weighted Euclidean biplots

We construct a weighted Euclidean distance that approximates any distance or dissimilarity measure between individuals that is based on a rectangular cases-by-variables data matrix. In contrast to regular multidimensional scaling methods for dissimilarity data, the method leads to biplots of individuals and variables while preserving all the good properties of dimension-reduction methods that are based on the singular-value decomposition. The main benefits are the decomposition of variance into components along principal axes, which provide the numerical diagnostics known as contributions, and the estimation of nonnegative weights for each variable. The idea is inspired by the distance functions used in correspondence analysis and in principal component analysis of standardized data, where the normalizations inherent in the distances can be considered as differential weighting of the variables. In weighted Euclidean biplots we allow these weights to be unknown parameters, which are estimated from the data to maximize the fit to the chosen distances or dissimilarities. These weights are estimated using a majorization algorithm. Once this extra weight-estimation step is accomplished, the procedure follows the classical path in decomposing the matrix and displaying its rows and columns in biplots.

Interaction between widening of diameter of abdominal aorta and cardiovascular risk factors and atherosclerosis burden.

The aim of this study was to investigate influence of traditional cardiovascular risk factors (CVRF) and subclinical atherosclerosis (ATS) burden on early stages of abdominal aortic diameter (AAD) widening among adults. 2,052 consecutive patients (P) (39 % women), mean age 52 ± 13 years, were prospectively screened for CVRF, ATS, and AAD. B-mode ultrasound was used to evaluate the largest AAD and to detect carotid and femoral atherosclerotic plaques. Mean AAD was 15.2 ± 2.8 mm. Atherosclerotic plaques were detected in 71 % of patients. Significant univariate correlation between AAD, traditional CVRF, and ABS was found. However, multiple regression analysis showed that only seven of them were significantly and weakly correlated with AAD (R² = 0.27, p < 0.001). On the other hand, a multivariate logistic analysis was used to evaluate CVRF impact on enlarged AAD ≥25 mm (EAAD) as compared to those with AAD <25 mm. These factors did not account for more than 30 % of interaction (R² = 0.30, p = 0.001). Furthermore, despite a large proportion of patients with high number of CVRF, and subclinical ATS, rate of patients with AAD ≥25 mm was low (1 %) and scattered regardless their CHD risk score or ATS burden. In conclusion, these results suggest that although some traditional CVRF and presence of ATS are associated with early stages of EAAD, other determinants still need to be identified for a better understanding of abdominal aortic aneurysm pathogenesis.

The prediction of speed and incline in outdoor running in humans using accelerometry.

PURPOSE: To explore whether triaxial accelerometric measurements can be utilized to accurately assess speed and incline of running in free-living conditions. METHODS: Body accelerations during running were recorded at the lower back and at the heel by a portable data logger in 20 human subjects, 10 men, and 10 women. After parameterizing body accelerations, two neural networks were designed to recognize each running pattern and calculate speed and incline. Each subject ran 18 times on outdoor roads at various speeds and inclines; 12 runs were used to calibrate the neural networks whereas the 6 other runs were used to validate the model. RESULTS: A small difference between the estimated and the actual values was observed: the square root of the mean square error (RMSE) was 0.12 m x s(-1) for speed and 0.014 radiant (rad) (or 1.4% in absolute value) for incline. Multiple regression analysis allowed accurate prediction of speed (RMSE = 0.14 m x s(-1)) but not of incline (RMSE = 0.026 rad or 2.6% slope). CONCLUSION: Triaxial accelerometric measurements allows an accurate estimation of speed of running and incline of terrain (the latter with more uncertainty). This will permit the validation of the energetic results generated on the treadmill as applied to more physiological unconstrained running conditions.

Remaining phosphorus estimated by pedotransfer function

Although the determination of remaining phosphorus (Prem) is simple, accurate values could also be estimated with a pedotransfer function (PTF) aiming at the additional use of soil analysis data and/or Prem replacement by an even simpler determination. The purpose of this paper was to develop a pedotransfer function to estimate Prem values of soils of the State of São Paulo based on properties with easier or routine laboratory determination. A pedotransfer function was developed by artificial neural networks (ANN) from a database of Prem values, pH values measured in 1 mol L-1 NaF solution (pH NaF) and soil chemical and physical properties of samples collected during soil classification activities carried out in the State of São Paulo by the Agronomic Institute of Campinas (IAC). Furthermore, a pedotransfer function was developed by regressing Prem values against the same predictor variables of the ANN-based PTF. Results showed that Prem values can be calculated more accurately with the ANN-based pedotransfer function with the input variables pH NaF values along with the sum of exchangeable bases (SB) and the exchangeable aluminum (Al3+) soil content. In addition, the accuracy of the Prem estimates by ANN-based PTF were more sensitive to increases in the experimental database size. Although the database used in this study was not comprehensive enough for the establishment of a definitive pedotrasnfer function for Prem estimation, results indicated the inclusion of Prem and pH NaF measurements among the soil testing evaluations as promising ind order to provide a greater database for the development of an ANN-based pedotransfer function for accurate Prem estimates from pH NaF, SB, and Al3+ values.