Optimal Concentration of Information Content For Log-Concave Densities

An elementary proof is provided of sharp bounds for the varentropy of random vectors with log-concave densities, as well as for deviations of the information content from its mean. These bounds significantly improve on the bounds obtained by Bobkov and Madiman ({\it Ann. Probab.}, 39(4):1528--1543, 2011).Comment: 15 pages. Changes in v2: Remark 2.5 (due to C. Saroglou) added with more general sufficient conditions for equality in Theorem 2.3. Also some minor corrections and added reference

Inconsistency of the MLE for the joint distribution of interval censored survival times and continuous marks

This paper considers the nonparametric maximum likelihood estimator (MLE) for the joint distribution function of an interval censored survival time and a continuous mark variable. We provide a new explicit formula for the MLE in this problem. We use this formula and the mark specific cumulative hazard function of Huang and Louis (1998) to obtain the almost sure limit of the MLE. This result leads to necessary and sufficient conditions for consistency of the MLE which imply that the MLE is inconsistent in general. We show that the inconsistency can be repaired by discretizing the marks. Our theoretical results are supported by simulations.Comment: 27 pages, 4 figure

D-optimal designs via a cocktail algorithm

A fast new algorithm is proposed for numerical computation of (approximate) D-optimal designs. This "cocktail algorithm" extends the well-known vertex direction method (VDM; Fedorov 1972) and the multiplicative algorithm (Silvey, Titterington and Torsney, 1978), and shares their simplicity and monotonic convergence properties. Numerical examples show that the cocktail algorithm can lead to dramatically improved speed, sometimes by orders of magnitude, relative to either the multiplicative algorithm or the vertex exchange method (a variant of VDM). Key to the improved speed is a new nearest neighbor exchange strategy, which acts locally and complements the global effect of the multiplicative algorithm. Possible extensions to related problems such as nonparametric maximum likelihood estimation are mentioned.Comment: A number of changes after accounting for the referees' comments including new examples in Section 4 and more detailed explanations throughou

HL-1 cells express an inwardly rectifying K+ current activated via muscarinic receptors comparable to that in mouse atrial myocytes

An inwardly rectifying K^+ current is present in atrial cardiac myocytes that is activated by acetylcholine (I_{KACh}). Physiologically, activation of the current in the SA node is important in slowing the heart rate with increased parasympathetic tone. It is a paradigm for the direct regulation of signaling effectors by the Gβγ G-protein subunit. Many questions have been addressed in heterologous expression systems with less focus on the behaviour in native myocytes partly because of the technical difficulties in undertaking comparable studies in native cells. In this study, we characterise a potassium current in the atrial-derived cell line HL-1. Using an electrophysiological approach, we compare the characteristics of the potassium current with those in native atrial cells and in a HEK cell line expressing the cloned Kir3.1/3.4 channel. The potassium current recorded in HL-1 is inwardly rectifying and activated by the muscarinic agonist carbachol. Carbachol-activated currents were inhibited by pertussis toxin and tertiapin-Q. The basal current was time-dependently increased when GTP was substituted in the patch-clamp pipette by the non-hydrolysable analogue GTPγS. We compared the kinetics of current modulation in HL-1 with those of freshly isolated atrial mouse cardiomyocytes. The current activation and deactivation kinetics in HL-1 cells are comparable to those measured in atrial cardiomyocytes. Using immunofluorescence, we found GIRK4 at the membrane in HL-1 cells. Real-time RT-PCR confirms the presence of mRNA for the main G-protein subunits, as well as for M2 muscarinic and A1 adenosine receptors. The data suggest HL-1 cells are a good model to study IKAch

Flow and retreat of the Late Quaternary Pine Island-Thwaites palaeo-ice stream, West Antarctica

Multibeam swath bathymetry and sub-bottom profiler data are used to establish constraints on the flow and retreat history of a major palaeo-ice stream that carried the combined discharge from the parts of the West Antarctic Ice Sheet now occupied by the Pine Island and Thwaites glacier basins. Sets of highly elongated bedforms show that, at the last glacial maximum, the route of the Pine Island-Thwaites palaeo-ice stream arced north-northeast following a prominent cross-shelf trough. In this area, the grounding line advanced to within similar to 68 km of, and probably reached, the shelf edge. Minimum ice thickness is estimated at 715 m on the outer shelf, and we estimate a minimum ice discharge of similar to 108 km(3) yr(-1) assuming velocities similar to today's Pine Island glacier (similar to 2.5 km yr(-1)). Additional bed forms observed in a trough northwest of Pine Island Bay likely formed via diachronous ice flows across the outer shelf and demonstrate switching ice stream behavior. The "style" of ice retreat is also evident in five grounding zone wedges, which suggest episodic deglaciation characterized by halts in grounding line migration up-trough. Stillstands occurred in association with changes in ice bed gradient, and phases of inferred rapid retreat correlate to higher bed slopes, supporting theoretical studies that show bed geometry as a control on ice margin recession. However, estimates that individual wedges could have formed within several centuries still imply a relatively rapid overall retreat. Our findings show that the ice stream channeled a substantial fraction of West Antarctica's discharge in the past, just as the Pine Island and Thwaites glaciers do today

Information bounds for Gaussian copulas

Often of primary interest in the analysis of multivariate data are the copula parameters describing the dependence among the variables, rather than the univariate marginal distributions. Since the ranks of a multivariate dataset are invariant to changes in the univariate marginal distributions, rank-based estimators are natural candidates for semiparametric copula estimation. Asymptotic information bounds for such estimators can be obtained from an asymptotic analysis of the rank likelihood, that is, the probability of the multivariate ranks. In this article, we obtain limiting normal distributions of the rank likelihood for Gaussian copula models. Our results cover models with structured correlation matrices, such as exchangeable or circular correlation models, as well as unstructured correlation matrices. For all Gaussian copula models, the limiting distribution of the rank likelihood ratio is shown to be equal to that of a parametric likelihood ratio for an appropriately chosen multivariate normal model. This implies that the semiparametric information bounds for rank-based estimators are the same as the information bounds for estimators based on the full data, and that the multivariate normal distributions are least favorable.Comment: Published in at http://dx.doi.org/10.3150/12-BEJ499 the Bernoulli (http://isi.cbs.nl/bernoulli/) by the International Statistical Institute/Bernoulli Society (http://isi.cbs.nl/BS/bshome.htm

A specific case in the classification of woods by FTIR and chemometric: discrimination of Fagales from Malpighiales

Fourier transform infrared (FTIR) spectroscopic data was used to classify wood samples from nine species within the Fagales and Malpighiales using a range of multivariate statistical methods. Taxonomic classification of the family Fagaceae and Betulaceae from Angiosperm Phylogenetic System Classification (APG II System) was successfully performed using supervised pattern recognition techniques. A methodology for wood sample discrimination was developed using both sapwood and heartwood samples. Ten and eight biomarkers emerged from the dataset to discriminate order and family, respectively. In the species studied FTIR in combination with multivariate analysis highlighted significant chemical differences in hemicelluloses, cellulose and guaiacyl (lignin) and shows promise as a suitable approach for wood sample classification

Characteristics of the Private Nuisance Wildlife Control Industry in New York

The nuisance wildlife control industry is rapidly expanding in New York State. To gain additional insight about this industry and the number of animals handled, we reviewed the 1989-90 annual logs submitted by Nuisance Wildlife Control Operators (NWCOs) to the New York State Department of Environmental Conservatlon (DEC). The specific objectives of this study were to determine: (1) the number and species of different wildlife responsible for damage incidents, (2) the cause of damage complaints, (3) the disposition of animals handled, (4) the location of damage events (i.e., urban, suburban, rural), and (5) an estimate of the economic impact of the nuisance wildlife industry in Upstate New York. The Nuisance Wildlife Logs (NWLs) were examined for 7 urban and 7 rural counties (25.5% of Upstate counties), and these data were used to estimate total NWCO activity in DEC Regions 3 through 9 (excludes Long Island). Approximately 75% of NWCOs licensed by DEC were active during 1989-90, and nearly 2,800 complaints were handled in the 14 counties sampled. More than 90% of complaints came from urban counties, and we estimated that NWCOs responded to more than 11,000 calls in Upstate New York. At a conservative estimate of $35/call, revenue generated by this industry exceeded$385,000 annually . Six wildlife species accounted for 85% of the nuisance complaints in urban and rural counties. During 1986 to 1993, the number of NWCOs licensed by DEC nearly quadrupled, and there is no indication that this trend will change in the near future

TGF-ß induces miR-100 and miR-125b but blocks let-7a through LIN28B controlling PDAC progression.

Abstract TGF-ß/Activin induces epithelial-to-mesenchymal transition (EMT) and stemness in pancreatic ductal adenocarcinoma (PDAC). However, the microRNAs (miRNAs) regulated during this response have remained yet undetermined. Here, we show that TGF-ß transcriptionally induces MIR100HG lncRNA, containing miR-100, miR-125b and let-7a in its intron, via SMAD2/3. Interestingly, we find that although the pro-tumourigenic miR-100 and miR-125b accordingly increase, the amount of anti-tumourigenic let-7a is unchanged, as TGF-ß also induces LIN28B inhibiting its maturation. Notably, we demonstrate that inactivation of miR-125b or miR-100 affects the TGF-ß-mediated response indicating that these miRNAs are important TGF-ß effectors. We integrated AGO2-RIP-seq with RNA-seq to identify the global regulation exerted by these miRNAs in PDAC cells. Transcripts targeted by miR-125b and miR-100 significantly overlap and mainly inhibit p53 and cell-cell junctions’ pathways. Together, we uncover that TGF-ß induces an lncRNA, whose encoded miRNAs, miR-100, let-7a and miR-125b, play opposing roles in controlling PDAC tumourigenesis