The Standard Biplot

Biplots are graphical displays of data matrices based on the decomposition of a matrix as the product of two matrices. Elements of these two matrices are used as coordinates for the rows and columns of the data matrix, with an interpretation of the joint presentation that relies on the properties of the scalar product. Because the decomposition is not unique, there are several alternative ways to scale the row and column points of the biplot, which can cause confusion amongst users, especially when software packages are not united in their approach to this issue. We propose a new scaling of the solution, called the standard biplot, which applies equally well 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. The standard biplot also handles data matrices with widely different levels of inherent variance. Two concepts taken from correspondence analysis are important to this idea: the weighting of row and column points, and the contributions made by the points to the solution. In the standard biplot one set of points, usually the rows of the data matrix, optimally represent the positions of the cases or sample units, which are weighted and usually standardized in some way unless the matrix contains values that are comparable in their raw form. The other set of points, usually the columns, is represented in accordance with their contributions to the low-dimensional solution. As for any biplot, the projections of the row points onto vectors defined by the column points approximate the centred and (optionally) standardized data. The method is illustrated with several examples to demonstrate how the standard biplot copes in different situations to give a joint map which needs only one common scale on the principal axes, thus avoiding the problem of enlarging or contracting the scale of one set of points to make the biplot readable. The proposal 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.

Meta-analysis of gene-level associations for rare variants based on single-variant statistics.

Meta-analysis of genome-wide association studies (GWASs) has led to the discoveries of many common variants associated with complex human diseases. There is a growing recognition that identifying "causal" rare variants also requires large-scale meta-analysis. The fact that association tests with rare variants are performed at the gene level rather than at the variant level poses unprecedented challenges in the meta-analysis. First, different studies may adopt different gene-level tests, so the results are not compatible. Second, gene-level tests require multivariate statistics (i.e., components of the test statistic and their covariance matrix), which are difficult to obtain. To overcome these challenges, we propose to perform gene-level tests for rare variants by combining the results of single-variant analysis (i.e., p values of association tests and effect estimates) from participating studies. This simple strategy is possible because of an insight that multivariate statistics can be recovered from single-variant statistics, together with the correlation matrix of the single-variant test statistics, which can be estimated from one of the participating studies or from a publicly available database. We show both theoretically and numerically that the proposed meta-analysis approach provides accurate control of the type I error and is as powerful as joint analysis of individual participant data. This approach accommodates any disease phenotype and any study design and produces all commonly used gene-level tests. An application to the GWAS summary results of the Genetic Investigation of ANthropometric Traits (GIANT) consortium reveals rare and low-frequency variants associated with human height. The relevant software is freely available.

The representations of adolescents about gynecological consultation

Objective: To analyze the social representation of adolescents about gynecological consultation and the influence of those in searching for consultations. Method: Qualitative descriptive study based on the Social Representations Theory, conducted with 50 adolescents in their last year of middle school. The data was collected between April and May of 2010 by Evocations and a Focal Group. The software EVOC and contextual analysis were used in the data treatment. Results: The elements fear and constraint, constant in the central nucleus, can justify the low frequency of adolescents in consultations. The term embarrassment in the peripheral system reinforce current sociocultural norms, while prevention, associated with learning about sex and clarifying doubts, allows to envision an educative function. Obtained testimonies in the focal groups exemplify and reinforce those findings. Conclusion: For an effective health education, professionals, including nurses, need to clarify the youth individually and collectively about their rights to privacy, secrecy, in addition to focus the gynecological consultation as a promotion measure to sexual and reproductive health.

The contributions of rare objects in correspondence analysis

Correspondence analysis, when used to visualize relationships in a table of counts(for example, abundance data in ecology), has been frequently criticized as being too sensitiveto objects (for example, species) that occur with very low frequency or in very few samples. Inthis statistical report we show that this criticism is generally unfounded. We demonstrate this inseveral data sets by calculating the actual contributions of rare objects to the results ofcorrespondence analysis and canonical correspondence analysis, both to the determination ofthe principal axes and to the chi-square distance. It is a fact that rare objects are oftenpositioned as outliers in correspondence analysis maps, which gives the impression that theyare highly influential, but their low weight offsets their distant positions and reduces their effecton the results. An alternative scaling of the correspondence analysis solution, the contributionbiplot, is proposed as a way of mapping the results in order to avoid the problem of outlying andlow contributing rare objects.

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.

On the robust effects of technology shocks on hours worked and output

We analyze the effects of neutral and investment-specific technology shockson hours and output. Long cycles in hours are captured in a variety of ways.Hours robustly fall in response to neutral shocks and robustly increase inresponse to investment specific shocks. The percentage of the variance ofhours (output) explained by neutral shocks is small (large); the opposite istrue for investment specific shocks. News shocks are uncorrelated with theestimated technology shocks.

Relación entre anchoveta y ambiente a diferentes escalas temporales

Relación entre anchoveta y ambiente a diferentes escalas temporales. Bol Inst Mar Perú 25(1-2):13-21.- Se estudió la relación del ambiente y la anchoveta peruana (Engraulis ringens) en el Ecosistema de Afloramiento Peruano (EAP). Se hipotetiza que en la zona de afloramiento comprendida entre 5°S - 13°S se presenta una relación negativa entre la temperatura, concentración de oxígeno y desembarques de anchoveta en escalas de tiempo decadales. En series de tiempo mensuales observadas entre 1950 – 2008, se analizó la temperatura superficial, subsuperficial, oxígeno, clorofila-a y desembarques de anchoveta y sardina, aplicando métodos estadísticos y espectrales para obtener modos temporales decadales y caracterizar la modulación a baja frecuencia de sus ciclos estacionales. También se usaron datos de reanálisis para caracterizar los cambios estacionales a partir de los forzamientos (remoto y local) ambientales del EAP. Se encontró una fuerte asociación a escala decadal entre las series ambientales costeras y los desembarques de anchoveta y se sugiere que la covariabilidad ambiental en escalas temporales interanuales, intraestacionales, decadales, seculares frente al EAP se incrementará durante la próxima década.

Synchronous and baroceptor-sensitive oscillations in skin microcirculation: evidence for central autonomic control.

To determine whether skin blood flow is local or takes part in general regulatory mechanisms, we recorded laser-Doppler flowmetry (LDF; left and right index fingers), blood pressure, muscle sympathetic nerve activity (MSNA), R-R interval, and respiration in 10 healthy volunteers and 3 subjects after sympathectomy. We evaluated 1) the synchronism of LDF fluctuations in two index fingers, 2) the relationship with autonomically mediated fluctuations in other signals, and 3) the LDF ability to respond to arterial baroreflex stimulation (by neck suction at frequencies from 0.02 to 0.20 Hz), using spectral analysis (autoregressive uni- and bivariate, time-variant algorithms). Synchronous LDF fluctuations were observed in the index fingers of healthy subjects but not in sympathectomized patients. LDF fluctuations were coherent with those obtained for blood pressure, MSNA, and R-R interval. LDF fluctuations were leading blood pressure in the low-frequency (LF; 0.1 Hz) band and lagging in the respiratory, high-frequency (HF; approximately 0.25 Hz) band, suggesting passive "downstream" transmission only for HF and "upstream" transmission for LF from the microvessels. LDF fluctuations were responsive to sinusoidal neck suction up to 0.1 Hz, indicating response to sympathetic modulation. Skin blood flow thus reflects modifications determined by autonomic activity, detectable by frequency analysis of spontaneous fluctuations.

Numerical simulation of ambient seismic wavefield modification caused by pore-fluid effects in an oil reservoir

We have modeled numerically the seismic response of a poroelastic inclusion with properties applicable to an oil reservoir that interacts with an ambient wavefield. The model includes wave-induced fluid flow caused by pressure differences between mesoscopic-scale (i.e., in the order of centimeters to meters) heterogeneities. We used a viscoelastic approximation on the macroscopic scale to implement the attenuation and dispersion resulting from this mesoscopic-scale theory in numerical simulations of wave propagation on the kilometer scale. This upscaling method includes finite-element modeling of wave-induced fluid flow to determine effective seismic properties of the poroelastic media, such as attenuation of P- and S-waves. The fitted, equivalent, viscoelastic behavior is implemented in finite-difference wave propagation simulations. With this two-stage process, we model numerically the quasi-poroelastic wave-propagation on the kilometer scale and study the impact of fluid properties and fluid saturation on the modeled seismic amplitudes. In particular, we addressed the question of whether poroelastic effects within an oil reservoir may be a plausible explanation for low-frequency ambient wavefield modifications observed at oil fields in recent years. Our results indicate that ambient wavefield modification is expected to occur for oil reservoirs exhibiting high attenuation. Whether or not such modifications can be detected in surface recordings, however, will depend on acquisition design and noise mitigation processing as well as site-specific conditions, such as the geologic complexity of the subsurface, the nature of the ambient wavefield, and the amount of surface noise.

Assessment of changes in body water by bioimpedance in acutely ill surgical patients.

OBJECTIVE: To evaluate the relationship between changes in body bioelectrical impedance (BI) at 0.5, 50 and kHz and the changes in body weight, as an index of total body water changes, in acutely ill surgical patients during the rapid infusion of isotonic saline solution. DESIGN: Prospective clinical study. SETTING: Multidisciplinary surgical ICU in a university hospital. PATIENTS: Twelve male patients treated for acute surgical illness (multiple trauma n = 5, major surgery n = 7). Selection criteria: stable cardiovascular parameters, normal cardiac function, signs of hypovolemia (CVP < or = 5 mmHg, urine output < 1 ml/kg x h). INTERVENTIONS: After baseline measurements, a 60 min fluid challenge test was performed with normal saline solution, 0.25 ml/kg/min [corrected]. MEASUREMENTS AND RESULTS: Body weight (platform digital scale), total body impedance (four-surface electrode technique; measurements at 0.5, 50 and 100 kHz) and urine output. Fluid retention induced a progressive decrease in BI at 0.5, 50 and 100 kHz, but the changes were significant for BI 0.5 and BI 100 only, from 40 min after the beginning of the fluid therapy onwards. There was a significant negative correlation between changes in water retention and BI 0.5, with individual correlation coefficients ranging from -0.72 to 0.95 (p < 0.01-0.0001). The slopes of the regression lines indicated that for each kg of water change, there was a mean decrease in BI of 18 ohm, but a substantial inter-individual variability was noted. CONCLUSION: BI measured at low frequency can represent a valuable index of acute changes in body water in a group of surgical patients but not in a given individual.

VLA multifrequency observations of RSCVn binaries

We present multiepoch Very Large Array (VLA) observations at 1.4 GHz, 4.9 GHz, 8.5 GHz and 14.9 GHz for a sample of eight RS CVn binary systems. Circular polarization measurements of these systems are also reported. Most of the fluxes observed are consistent with incoherent emission from mildly relativistic electrons. Several systems show an increase of the degree of circular polarization with increasing frequency in the optically thin regime, in conflict with predictions by gyrosynchrotron models. We observed a reversal in the sense of circular polarization with increasing frequency in three non-eclipsing systems: EI Eri, DM Uma and HD 8358. We find clear evidence for coherent plasma emission at 1.4 GHz in the quiescent spectrum of HD 8358 during the helicity reversal. The degrees of polarization of the other two systems could also be accounted for by a coherent emission process. The observations of ER Vul revealed two U-shaped flux spectra at the highest frequencies. The U-shape of the spectra may be accounted for by an optically thin gyrosynchrotron source for the low frequency part whereas the high frequency part is dominated by a thermal emission component.

An improved trap design for decoupling multinuclear RF coils.

PURPOSE: Multinuclear magnetic resonance spectroscopy and imaging require a radiofrequency probe capable of transmitting and receiving at the proton and non-proton frequencies. To minimize coupling between probe elements tuned to different frequencies, LC (inductor-capacitor) traps blocking current at the (1) H frequency can be inserted in non-proton elements. This work compares LC traps with LCC traps, a modified design incorporating an additional capacitor, enabling control of the trap reactance at the low frequency while maintaining (1) H blocking. METHODS: Losses introduced by both types of trap were analysed using circuit models. Radiofrequency coils incorporating a series of LC and LCC traps were then built and evaluated at the bench. LCC trap performance was then confirmed using (1) H and (13) C measurements in a 7T human scanner. RESULTS: LC and LCC traps both effectively block interaction between non-proton and proton coils at the proton frequency. LCC traps were found to introduce a sensitivity reduction of 5±2%, which was less than half of that caused by LC traps. CONCLUSION: Sensitivity of non-proton coils is critical. The improved trap design, incorporating one extra capacitor, significantly reduces losses introduced by the trap in the non-proton coil. Magn Reson Med 72:584-590, 2014. © 2013 Wiley Periodicals, Inc.

Hexokinase 2, Glycogen Synthase and Phosphorylase Play a Key Role in Muscle Glycogen Supercompensation

Background: Glycogen-depleting exercise can lead to supercompensation of muscle glycogen stores, but the biochemical mechanisms of this phenomenon are still not completely understood. Methods: Using chronic low-frequency stimulation (CLFS) as an exercise model, the tibialis anterior muscle of rabbits was stimulated for either 1 or 24 hours, inducing a reduction in glycogen of 90% and 50% respectively. Glycogen recovery was subsequently monitored during 24 hours of rest. Results: In muscles stimulated for 1 hour, glycogen recovered basal levels during the rest period. However, in those stimulated for 24 hours, glycogen was supercompensated and its levels remained 50% higher than basal levels after 6 hours of rest, although the newly synthesized glycogen had fewer branches. This increase in glycogen correlated with an increase in hexokinase-2 expression and activity, a reduction in the glycogen phosphorylase activity ratio and an increase in the glycogen synthase activity ratio, due to dephosphorylation of site 3a, even in the presence of elevated glycogen stores. During supercompensation there was also an increase in 59-AMP-activated protein kinase phosphorylation, correlating with a stable reduction in ATP and total purine nucleotide levels. Conclusions: Glycogen supercompensation requires a coordinated chain of events at two levels in the context of decreased cell energy balance: First, an increase in the glucose phosphorylation capacity of the muscle and secondly, control of the enzymes directly involved in the synthesis and degradation of the glycogen molecule. However, supercompensated glycogen has fewer branches.

Anomalous crossover between thermal and shot noise in macroscopic diffusive conductors

We predict the existence of an anomalous crossover between thermal and shot noise in macroscopic diffusive conductors. We first show that, besides thermal noise, these systems may also exhibit shot noise due to fluctuations of the total number of carriers in the system. Then we show that at increasing currents the crossover between the two noise behaviors is anomalous, in the sense that the low-frequency current spectral density displays a region with a superlinear dependence on the current up to a cubic law. The anomaly is due to the nontrivial coupling in the presence of the long-range Coulomb interaction among the three time scales relevant to the phenomenon, namely, diffusion, transit, and dielectric relaxation time.

The complete Raman spectrum of nanometric SnO2 particles

The complete Raman spectrum of SnO2 nanoparticles in presented and analyzed. In addition to the "classical" modes observed in the rutile structure, two other regions shown Raman activity for nanoparticles. The Raman bands in the low-frequency region are attributed to acoustic modes associated with the vibration of the individual nanoparticle as a whole. The high-frequency region is activated by surface disorder. A detailed analysis of these regions and the changes in the normal modes of SnO2 are presented as a function nanoparticle size.