Site and Age Condition the Growth Responses to Climate and Drought of Relict Pinus nigra subsp. salzmannii Populations in Southern Spain

To assess if tree age may modulate the main climatic drivers of radial growth, two relict Pinus nigra subsp. salzmannii populations (Maria, most xeric site; Magina, least xeric site) were sampled in southern Spain near the limits of the species range. Tree-ring width residual chronologies for two age groups (mature trees, age <= 100 years (minimum 40 years); old trees, age > 100 years) were built to evaluate their responses to climate by relating them to monthly precipitation and temperature and a drought index (DRI) using correlation and response functions. We found that drought is the main driver of growth of relict P. nigra populations, but differences between sites and age classes were also observed. First, growth in the most xeric site depends on the drought severity during the previous autumn and the spring of the year of tree-ring formation, whereas in the relatively more mesic site growth is mainly enhanced by warm and wet conditions in spring. Second, growth of mature trees responded more to drought severity than that of old trees. Our findings indicate that drought severity will mainly affect growth of relict P. nigra populations dominated by mature trees in xeric sites. This conclusion may also apply to similar mountain Mediterranean conifer relicts.

Impact of biatrial defragmentation in patients with paroxysmal atrial fibrillation: Results from a randomized prospective study

Background Single procedure success rates of pulmonary vein isolation (PVI) in patients with paroxysmal atrial fibrillation (PAF) are still unsatisfactory. In patients with persistent atrial fibrillation (AF), ablation of complex fractionated atrial electrograms (CFAEs) after PVI results in improved outcomes. Objective We aimed to investigate if PAF-patients with intraprocedurally sustained AF after PVI might benefit from additional CFAE ablation. Methods A total of 1134 consecutive patients underwent a first catheter ablation procedure of PAF between June 2008 and December 2012. In most patients, AF was either not inducible or terminated during PVI. In 68 patients (6%), AF sustained after successful PVI. These patients were randomized to either cardioversion (PVI-alone group; n = 33) or additional CFAE ablation (PVI+CFAE group; n = 35) and followed up every 1–3 months and serial Holter recordings were also obtained. The primary end point was the recurrence of AF/atrial tachycardia (AT) after a blanking period of 3 months. Results Procedure duration (127 ± 6 minutes vs 174 ± 10 minutes), radiofrequency application time (44 ± 3 minutes vs 74 ± 5 minutes), and fluoroscopy time (26 ± 2 minutes vs 41 ± 3 minutes) were longer in the PVI+CFAE group (all P < .001). In 30 of 35 patients (86%) in the PVI+CFAE group, ablation terminated AF. There was no significant group difference with respect to freedom from AF/AT (22 of 33 [67%] vs 22 of 35 [63%]; P = .66). Subsequently, 10 of 11 patients in the PVI-alone group (91%) and 11 of 13 patients in PVI+CFAE group (85%) underwent repeat ablation (P = 1.00). Overall, 29 of 33 [88%] vs 30 of 35 [86%] patients (P = 1.00) were free from AF/AT after 1.4 ± 0.1 vs 1.4 ± 0.2 (P = .87) procedures. Conclusion Patients with sustained AF after PVI in a PAF cohort are rare. Regarding AF/AT recurrence, these patients did not benefit from further CFAE ablation compared to PVI alone, but are exposed to longer procedure duration, fluoroscopy time, and radiofrequency application time.

Plotting regression coefficients and other estimates

Graphical display of regression results has become increasingly popular in presentations and in scientific literature because graphs are often much easier to read than tables. Such plots can be produced in Stata by the marginsplot command (see [R] marginsplot). However, while marginsplot is versatile and flexible, it has two major limitations: it can only process results left behind by margins (see [R] margins), and it can handle only one set of results at a time. In this article, I introduce a new command called coefplot that overcomes these limitations. It plots results from any estimation command and combines results from several models into one graph. The default behavior of coefplot is to plot markers for coefficients and horizontal spikes for confidence intervals. However, coefplot can also produce other types of graphs. I illustrate the capabilities of coefplot by using a series of examples.

The Blinder-Oaxaca decomposition for linear regression models

The counterfactual decomposition technique popularized by Blinder (1973, Journal of Human Resources, 436–455) and Oaxaca (1973, International Economic Review, 693–709) is widely used to study mean outcome differences between groups. For example, the technique is often used to analyze wage gaps by sex or race. This article summarizes the technique and addresses several complications, such as the identification of effects of categorical predictors in the detailed decomposition or the estimation of standard errors. A new command called oaxaca is introduced, and examples illustrating its usage are given.

Multinomial goodness-of-fit: large sample tests with survey design correction and exact tests for small samples

I introduce the new mgof command to compute distributional tests for discrete (categorical, multinomial) variables. The command supports largesample tests for complex survey designs and exact tests for small samples as well as classic large-sample x2-approximation tests based on Pearson’s X2, the likelihood ratio, or any other statistic from the power-divergence family (Cressie and Read, 1984, Journal of the Royal Statistical Society, Series B (Methodological) 46: 440–464). The complex survey correction is based on the approach by Rao and Scott (1981, Journal of the American Statistical Association 76: 221–230) and parallels the survey design correction used for independence tests in svy: tabulate. mgof computes the exact tests by using Monte Carlo methods or exhaustive enumeration. mgof also provides an exact one-sample Kolmogorov–Smirnov test for discrete data.

Making regression tables simplified

estout, introduced by Jann (Stata Journal 5: 288–308), is a useful tool for producing regression tables from stored estimates. However, its syntax is relatively complex and commands may turn out long even for simple tables. Furthermore, having to store the estimates beforehand can be cumbersome. To facilitate the production of regression tables, I therefore present here two new commands called eststo and esttab. eststo is a wrapper for offcial Stata’s estimates store and simplifies the storing of estimation results for tabulation. esttab, on the other hand, is a wrapper for estout and simplifies compiling nice-looking tables from the stored estimates without much typing. I also provide updates to estout and estadd.

Making regression tables from stored estimates

Organizing and archiving statistical results and processing a subset of those results for publication are important and often underestimated issues in conducting statistical analyses. Because automation of these tasks is often poor, processing results produced by statistical packages is quite laborious and vulnerable to error. I will therefore present a new package called estout that facilitates and automates some of these tasks. This new command can be used to produce regression tables for use with spreadsheets, LaTeX, HTML, or word processors. For example, the results for multiple models can be organized in spreadsheets and can thus be archived in an orderly manner. Alternatively, the results can be directly saved as a publication-ready table for inclusion in, for example, a LaTeX document. estout is implemented as a wrapper for estimates table but has many additional features, such as support for mfx. However, despite its flexibility, estout is—I believe—still very straightforward and easy to use. Furthermore, estout can be customized via so-called defaults files. A tool to make available supplementary statistics called estadd is also provided.

Tabulation of multiple responses

Although multiple-response questions are quite common in survey research, Stata’s official release does not provide much capability for an effective analysis of multiple-response variables. For example, in a study on drug addiction an interview question might be, “Which substances did you consume during the last four weeks?” The respondents just list all the drugs they took, if any; e.g., an answer could be “cannabis, cocaine, heroin” or “ecstasy, cannabis” or “none”, etc. Usually, the responses to such questions are stored as a set of variables and, therefore, cannot be easily tabulated. I will address this issue here and present a new module to compute one- and two-way tables of multiple responses. The module supports several types of data structure, provides significance tests, and offers various options to control the computation and display of the results. In addition, tools to create graphs of multiple-response distributions are presented.

Rhythmic synaptic activity induced by mechanical injury of rat CA3 hippocampal area.

Mechanical injury of the CNS frequently results from accidents but also occurs in the course of neurosurgical interventions. A great variety of anatomical and physiological changes have been described to evolve after a brain trauma yet only little is known about processes that occur during a trauma. In the present study, I obtained whole-cell patch clamp recordings from pyramidal cells in hippocampal slice cultures while mechanically lesioning the CA3 area. Electrophysiological analysis revealed that traumatic injury massively increased excitatory and inhibitory synaptic activity in the entire CA3 region. Cutting the CA3 region induced highly rhythmic excitatory postsynaptic currents (EPSCs) that reached frequencies of around 70 Hz. Blocking voltage-dependent sodium channels with tetrodotoxin prevented the increase in synaptic activity and injury-induced neurotransmitter release in CA3 remote from the lesion site. With fast synaptic transmission blocked only neurons in the immediate vicinity of a lesion depolarized and fired action potentials upon mechanical damage. I hence suggest that mechanical injury damages the membrane and induces action potential firing in only a small population of neurons. This activity is then propagated throughout the undamaged CA3 network inducing highly rhythmic discharges. Thus mechanical brain injury initiates immediate functional changes that exceed the lesion site.