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.

Snow shielding factors for cosmogenic nuclide dating inferred from long-term neutron detector monitoring

The depth-dependent attenuation of the secondary cosmic-ray particle flux due to snow cover and its effects on production rates of cosmogenic nuclides constitutes a potential source of uncertainty for studies conducted in regions characterized by frequent seasonal snow burial. Recent experimental and numerical modelling studies have yielded new constraints on the effect of hydrogen-rich media on the production rates of cosmogenic nuclides by low- and high-energy neutrons (<10(-3) MeV and >10(2) MeV, respectively). Here we present long-term neutron-detector monitoring data from a natural setting that we use to quantify the effect of snow cover on the attenuation of fast neutrons (0.1-10 MeV), which are responsible for the production of Ne-21 from Mg and Cl-36 from K. We use data measured between July 2001 and May 2008 at seven stations located throughout the Ecrins-Pelvoux massif (French Western Alps) and its surroundings, at elevations ranging from 200 to 2500 m a.s.l. From the cosmic-ray fluxes recorded during summer, when snow is absent, we infer an apparent attenuation length of 148 g cm(-2) in the atmosphere at a latitude of similar to 45 degrees N and for altitudes ranging from similar to 200 to 2500 m a.s.l. Using snow water-equivalent (SWE) values obtained through snow-coring campaigns that overlap in time the neutron monitoring for five stations, we show that fast neutrons are much more strongly attenuated in snow than predicted by a conventional mass-shielding formulation and the attenuation length estimated in the atmosphere. We suggest that such strong attenuation results from boundary effects at the atmosphere/snow interface induced by the high efficiency of water as a neutron moderator. Finally, we propose an empirical model that allows calculating snow-shielding correction factors as a function of SWE for studies using Ne-21 and Cl-36 analyses in Mg- and K-rich minerals, respectively. This empirical model is of interest for studies with a focus on cosmic-ray exposure dating, particularly if the target rocks are made up of mafic to ultramafic units where seasonal snow-cover is a common phenomenon.

A large animal model of spinal muscular atrophy and correction of phenotype.

OBJECTIVES Spinal muscular atrophy (SMA) is caused by reduced levels of survival motor neuron (SMN) protein, which results in motoneuron loss. Therapeutic strategies to increase SMN levels including drug compounds, antisense oligonucleotides, and scAAV9 gene therapy have proved effective in mice. We wished to determine whether reduction of SMN in postnatal motoneurons resulted in SMA in a large animal model, whether SMA could be corrected after development of muscle weakness, and the response of clinically relevant biomarkers. METHODS Using intrathecal delivery of scAAV9 expressing an shRNA targeting pig SMN1, SMN was knocked down in motoneurons postnatally to SMA levels. This resulted in an SMA phenotype representing the first large animal model of SMA. Restoration of SMN was performed at different time points with scAAV9 expressing human SMN (scAAV9-SMN), and electrophysiology measurements and pathology were performed. RESULTS Knockdown of SMN in postnatal motoneurons results in overt proximal weakness, fibrillations on electromyography indicating active denervation, and reduced compound muscle action potential (CMAP) and motor unit number estimation (MUNE), as in human SMA. Neuropathology showed loss of motoneurons and motor axons. Presymptomatic delivery of scAAV9-SMN prevented SMA symptoms, indicating that all changes are SMN dependent. Delivery of scAAV9-SMN after symptom onset had a marked impact on phenotype, electrophysiological measures, and pathology. INTERPRETATION High SMN levels are critical in postnatal motoneurons, and reduction of SMN results in an SMA phenotype that is SMN dependent. Importantly, clinically relevant biomarkers including CMAP and MUNE are responsive to SMN restoration, and abrogation of phenotype can be achieved even after symptom onset.

U7 snRNAs induce correction of mutated dystrophin pre-mRNA by exon skipping.

Most cases of Duchenne muscular dystrophy are caused by dystrophin gene mutations that disrupt the mRNA reading frame. Artificial exclusion (skipping) of a single exon would often restore the reading frame, giving rise to a shorter, but still functional dystrophin protein. Here, we analyzed the ability of antisense U7 small nuclear (sn)RNA derivatives to alter dystrophin pre-mRNA splicing. As a proof of principle, we first targeted the splice sites flanking exon 23 of dystrophin pre-mRNA in the wild-type muscle cell line C2C12 and showed precise exon 23 skipping. The same strategy was then successfully adapted to dystrophic immortalized mdx muscle cells where exon-23-skipped dystrophin mRNA rescued dystrophin protein synthesis. Moreover, we observed a stimulation of antisense U7 snRNA expression by the murine muscle creatine kinase enhancer. These results demonstrate that alteration of dystrophin pre-mRNA splicing could correct dystrophin gene mutations by expression of specific U7 snRNA constructs.

Benchmark problems for reactive transport modeling of the generation and attenuation of acid rock drainage

Acid rock drainage (ARD) is a problem of international relevance with substantial environmental and economic implications. Reactive transport modeling has proven a powerful tool for the process-based assessment of metal release and attenuation at ARD sites. Although a variety of models has been used to investigate ARD, a systematic model intercomparison has not been conducted to date. This contribution presents such a model intercomparison involving three synthetic benchmark problems designed to evaluate model results for the most relevant processes at ARD sites. The first benchmark (ARD-B1) focuses on the oxidation of sulfide minerals in an unsaturated tailing impoundment, affected by the ingress of atmospheric oxygen. ARD-B2 extends the first problem to include pH buffering by primary mineral dissolution and secondary mineral precipitation. The third problem (ARD-B3) in addition considers the kinetic and pH-dependent dissolution of silicate minerals under low pH conditions. The set of benchmarks was solved by four reactive transport codes, namely CrunchFlow, Flotran, HP1, and MIN3P. The results comparison focused on spatial profiles of dissolved concentrations, pH and pE, pore gas composition, and mineral assemblages. In addition, results of transient profiles for selected elements and cumulative mass loadings were considered in the intercomparison. Despite substantial differences in model formulations, very good agreement was obtained between the various codes. Residual deviations between the results are analyzed and discussed in terms of their implications for capturing system evolution and long-term mass loading predictions.

Full correction for spatially distributed speed-of-sound in echo ultrasound based on measuring aberration delays via transmit beam steering

Aberrations of the acoustic wave front, caused by spatial variations of the speed-of-sound, are a main limiting factor to the diagnostic power of medical ultrasound imaging. If not accounted for, aberrations result in low resolution and increased side lobe level, over all reducing contrast in deep tissue imaging. Various techniques have been proposed for quantifying aberrations by analysing the arrival time of coherent echoes from so-called guide stars or beacons. In situations where a guide star is missing, aperture-based techniques may give ambiguous results. Moreover, they are conceptually focused on aberrators that can be approximated as a phase screen in front of the probe. We propose a novel technique, where the effect of aberration is detected in the reconstructed image as opposed to the aperture data. The varying local echo phase when changing the transmit beam steering angle directly reflects the varying arrival time of the transmit wave front. This allows sensing the angle-dependent aberration delay in a spatially resolved way, and thus aberration correction for a spatially distributed volume aberrator. In phantoms containing a cylindrical aberrator, we achieved location-independent diffraction-limited resolution as well as accurate display of echo location based on reconstructing the speed-of-sound spatially resolved. First successful volunteer results confirm the clinical potential of the proposed technique.

TIMS measurements of full range of natural Ca isotopes with internally consistent fractionation correction

This study presents static measurements of the Ca isotopic composition of standard reference materials SRM 915 a/b on a Triton Plus™ thermal ionization mass spectrometer with a specially developed Faraday cup array allowing simultaneous measurement of 40Ca and 48Ca. The total amount of Ca in all analyses was kept < 1 µg. With this setup the measurement uncertainties were 0.06 ‰ for 40Ca/44Ca and 0.12 ‰ for 48Ca/40Ca. Measuring all isotopes simultaneously better allows to test the internal consistency of different Ca isotope abundances reported in the literature. The exponential law was observed to correct incompletely instrumental mass fractionation. An improved fractionation correction based on the exponential law is proposed. It changes the 40Ca/44Ca ratio of SRM 915a (corrected relative to 42Ca/44Ca = 0.31221; 48Ca/44Ca = 0.08871) from 47.1635 ± 0.0028 to 47.1649 ± 0.0047. The measurements of SRM 915b were performed with different analytical conditions (runs were prolonged till complete filament load depletion). Even if the 40Ca/44Ca ratio of SRM 915b, when corrected with the simple exponential law, appears different (47.1532 ± 0.0038) from that of SRM 915a, it becomes coincident (47.1613 ± 0.0028) when corrected with a second-order refinement. This supports the use of the improved exponential law to obtain internally consistent Ca isotope ratio for natural samples.

Electrocardiographic characteristics of patients with funnel chest before and after surgical correction using pectus bar: A new association with precordial J wave pattern.

OBJECTIVE Abnormal ECG findings suggestive of cardiac disease are frequent in patients with funnel chest, although structural heart disease is rare. Electrocardiographic characteristics and changes following new surgical treatments in young adults are not described so far. The aim of the study was to analyze electrocardiographic characteristics of patients with funnel chest before and after minimally invasive funnel chest correction by the Nuss procedure. METHODS Twenty-six patients with surgical correction of funnel chest using pectus bar were included. Twelve-lead ECGs before and later than one year after surgery were analyzed. RESULTS In postoperative ECGs, amplitude of P wave in lead II and negative terminal amplitude of P wave in lead V1 decreased from 0.13 to 0.10mV (p=0.03), and from 0.10 to 0.04mV (p<0.001), respectively. Mean QRS duration decreased from 108ms to 98ms (p=0.003) after correction. A pathological left and right Sokolow-Lyon index was observed in 35% and 23% of patients before, versus 8% (p=0.04) and 0% (p=0.01) after correction, respectively. In contrast, the rate of patients with J wave pattern in precordial leads V4-V6 increased from 8% before to 42% after surgery (p=0.004). CONCLUSIONS ECG abnormalities in patients with funnel chest are frequent, and can normalize after surgical correction by the Nuss procedure. De novo J wave pattern in precordial leads V4-V6 is a frequent finding after surgical funnel chest correction using pectus bar.

Correction: Usage of Complementary Medicine in Switzerland: Results of the Swiss Health Survey 2012 and Development Since 2007.

[This corrects the article DOI: 10.1371/journal.pone.0141985.]

Observations of the Sensitivity of Beam Attenuation to Particle Size in a Coastal Bottom Boundary Layer

The goal of this study was to test the hypothesis that the aggregated state of natural marine particles constrains the sensitivity of optical beam attenuation to particle size. An instrumented bottom tripod was deployed at the 12-m node of the Martha's Vineyard Coastal Observatory to monitor particle size distributions, particle size-versus-settling-velocity relationships, and the beam attenuation coefficient (c(p)) in the bottom boundary layer in September 2007. An automated in situ filtration system on the tripod collected 24 direct estimates of suspended particulate mass (SPM) during each of five deployments. On a sampling interval of 5 min, data from a Sequoia Scientific LISST 100x Type B were merged with data from a digital floc camera to generate suspended particle volume size distributions spanning diameters from approximately 2 mu m to 4 cm. Diameter-dependent densities were calculated from size-versus-settling-velocity data, allowing conversion of the volume size distributions to mass distributions, which were used to estimate SPM every 5 min. Estimated SPM and measured c(p) from the LISST 100x were linearly correlated throughout the experiment, despite wide variations in particle size. The slope of the line, which is the ratio of c(p) to SPM, was 0.22 g m(-2). Individual estimates of c(p):SPM were between 0.2 and 0.4 g m(-2) for volumetric median particle diameters ranging from 10 to 150 mu m. The wide range of values in c(p):SPM in the literature likely results from three factors capable of producing factor-of-two variability in the ratio: particle size, particle composition, and the finite acceptance angle of commercial beam-transmissometers.