Templates for Convex Cone Problems with Applications to Sparse Signal Recovery

This paper develops a general framework for solving a variety of convex cone problems that frequently arise in signal processing, machine learning, statistics, and other fields. The approach works as follows: first, determine a conic formulation of the problem; second, determine its dual; third, apply smoothing; and fourth, solve using an optimal first-order method. A merit of this approach is its flexibility: for example, all compressed sensing problems can be solved via this approach. These include models with objective functionals such as the total-variation norm, ||Wx||_1 where W is arbitrary, or a combination thereof. In addition, the paper also introduces a number of technical contributions such as a novel continuation scheme, a novel approach for controlling the step size, and some new results showing that the smooth and unsmoothed problems are sometimes formally equivalent. Combined with our framework, these lead to novel, stable and computationally efficient algorithms. For instance, our general implementation is competitive with state-of-the-art methods for solving intensively studied problems such as the LASSO. Further, numerical experiments show that one can solve the Dantzig selector problem, for which no efficient large-scale solvers exist, in a few hundred iterations. Finally, the paper is accompanied with a software release. This software is not a single, monolithic solver; rather, it is a suite of programs and routines designed to serve as building blocks for constructing complete algorithms.Comment: The TFOCS software is available at http://tfocs.stanford.edu This version has updated reference

Turning to God in the Face of Ostracism: Effects of Social Exclusion on Religiousness

The present research proposes that individuals who are socially excluded can turn to religion to cope with the experience. Empirical studies conducted to test this hypothesis consistently found that socially excluded persons reported (a) significantly higher levels of religious affiliation (Studies 1, 2, and 4) and (b) stronger intentions to engage in religious behaviors (Study 2) than comparable, nonexcluded individuals. Direct support for the stress-buffering function of religiousness was also found, with a religious prime reducing the aggression-eliciting effects of consequent social rejection (Study 5). These effects were observed in both Christian and Muslim samples, revealing that turning to religion can be a powerful coping response when dealing with social rejection. Theoretical and practical implications of these findings are discussed

Implementation of an Optimal First-Order Method for Strongly Convex Total Variation Regularization

We present a practical implementation of an optimal first-order method, due to Nesterov, for large-scale total variation regularization in tomographic reconstruction, image deblurring, etc. The algorithm applies to $\mu$-strongly convex objective functions with $L$-Lipschitz continuous gradient. In the framework of Nesterov both $\mu$ and $L$ are assumed known -- an assumption that is seldom satisfied in practice. We propose to incorporate mechanisms to estimate locally sufficient $\mu$ and $L$ during the iterations. The mechanisms also allow for the application to non-strongly convex functions. We discuss the iteration complexity of several first-order methods, including the proposed algorithm, and we use a 3D tomography problem to compare the performance of these methods. The results show that for ill-conditioned problems solved to high accuracy, the proposed method significantly outperforms state-of-the-art first-order methods, as also suggested by theoretical results.Comment: 23 pages, 4 figure

Bayesian system identification of a nonlinear dynamical system using a novel variant of Simulated Annealing

This work details the Bayesian identification of a nonlinear dynamical system using a novel MCMC algorithm: 'Data Annealing'. Data Annealing is similar to Simulated Annealing in that it allows the Markov chain to easily clear 'local traps' in the target distribution. To achieve this, training data is fed into the likelihood such that its influence over the posterior is introduced gradually - this allows the annealing procedure to be conducted with reduced computational expense. Additionally, Data Annealing uses a proposal distribution which allows it to conduct a local search accompanied by occasional long jumps, reducing the chance that it will become stuck in local traps. Here it is used to identify an experimental nonlinear system. The resulting Markov chains are used to approximate the covariance matrices of the parameters in a set of competing models before the issue of model selection is tackled using the Deviance Information Criterion. © 2014

The N-methyl-d-aspartate receptor antagonist CPP alters synapse and spine structure and impairs long-term potentiation and long-term depression induced morphological plasticity in dentate gyrus of the awake rat

Long-term morphological synaptic changes associated with homosynaptic long-term potentiation (LTP) and heterosynaptic long-term depression (LTD) in vivo, in awake adult rats were analyzed using three-dimensional (3-D) reconstructions of electron microscope images of ultrathin serial sections from the molecular layer of the dentate gyrus. For the first time in morphological studies, the specificity of the effects of LTP and LTD on both spine and synapse ultrastructure was determined using an N-methyl-d-aspartate (NMDA) receptor antagonist CPP (3-[(R)-2-carboxypiperazin-4-yl]-propyl-1-phosphonic acid). There were no differences in synaptic density 24 h after LTP or LTD induction, and CPP alone had no effect on synaptic density. LTP increased significantly the proportion of mushroom spines, whereas LTD increased the proportion of thin spines, and both LTP and LTD decreased stubby spine number. Both LTP and LTD increased significantly spine head evaginations (spinules) into synaptic boutons and CPP blocked these changes. Synaptic boutons were smaller after LTD, indicating a pre-synaptic effect. Interestingly, CPP alone decreased bouton and mushroom spine volumes, as well as post-synaptic density (PSD) volume of mushroom spines.These data show similarities, but also some clear differences, between the effects of LTP and LTD on spine and synaptic morphology. Although CPP blocks both LTP and LTD, and impairs most morphological changes in spines and synapses, CPP alone was shown to exert effects on aspects of spine and synaptic structure

Faking like a woman? Towards an interpretative theorization of sexual pleasure.

This article explores the possibility of developing a feminist approach to gendered and sexual embodiment which is rooted in the pragmatist/interactionist tradition derived from G.H. Mead, but which in turn develops this perspective by inflecting it through more recent feminist thinking. In so doing we seek to rebalance some of the rather abstract work on gender and embodiment by focusing on an instance of 'heterosexual' everyday/night life - the production of the female orgasm. Through engaging with feminist and interactionist work, we develop an approach to embodied sexual pleasure that emphasizes the sociality of sexual practices and of reflexive sexual selves. We argue that sexual practices and experiences must be understood in social context, taking account of the situatedness of sex as well as wider socio-cultural processes the production of sexual desire and sexual pleasure (or their non-production) always entails interpretive, interactional processes

On the selection of AGN neutrino source candidates for a source stacking analysis with neutrino telescopes

The sensitivity of a search for sources of TeV neutrinos can be improved by grouping potential sources together into generic classes in a procedure that is known as source stacking. In this paper, we define catalogs of Active Galactic Nuclei (AGN) and use them to perform a source stacking analysis. The grouping of AGN into classes is done in two steps: first, AGN classes are defined, then, sources to be stacked are selected assuming that a potential neutrino flux is linearly correlated with the photon luminosity in a certain energy band (radio, IR, optical, keV, GeV, TeV). Lacking any secure detailed knowledge on neutrino production in AGN, this correlation is motivated by hadronic AGN models, as briefly reviewed in this paper. The source stacking search for neutrinos from generic AGN classes is illustrated using the data collected by the AMANDA-II high energy neutrino detector during the year 2000. No significant excess for any of the suggested groups was found.Comment: 43 pages, 12 figures, accepted by Astroparticle Physic

Measurement of the Ratio Gamma(KL -> pi+ pi-)/Gamma(KL -> pi e nu) and Extraction of the CP Violation Parameter |eta+-|

We present a measurement of the ratio of the decay rates Gamma(KL -> pi+ pi-)/Gamma(KL -> pi e nu), denoted as Gamma(K2pi)/Gamma(Ke3). The analysis is based on data taken during a dedicated run in 1999 by the NA48 experiment at the CERN SPS. Using a sample of 47000 K2pi and five million Ke3 decays, we find Gamma(K2pi)/Gamma(Ke3) = (4.835 +- 0.022(stat) +- 0.016(syst)) x 10^-3. From this we derive the branching ratio of the CP violating decay KL -> pi+ pi- and the CP violation parameter |eta+-|. Excluding the CP conserving direct photon emission component KL -> pi+ pi- gamma, we obtain the results BR(KL -> pi+ pi-) = (1.941 +- 0.019) x 10^-3 and |eta+-| = (2.223 +- 0.012) x 10^-3.Comment: 20 pages, 7 figures, accepted by Phys. Lett.

Measurement of K^0_e3 form factors

The semileptonic decay of the neutral K meson, KL -> pi e nu (Ke3), was used to study the strangeness-changing weak interaction of hadrons. A sample of 5.6 million reconstructed events recorded by the NA48 experiment was used to measure the Dalitz plot density. Admitting all possible Lorentz-covariant couplings, the form factors for vector (f_+(q^2)), scalar (f_S) and tensor (f_T) interactions were measured. The linear slope of the vector form factor lambda_+ = 0.0284+-0.0007+-0.0013 and values for the ratios |f_S/f_+(0)| = 0.015^{+0.007}_{-0.010}+-0.012 and |f_T/f_+(0)| = 0.05^{+0.03}_{-0.04}+-0.03 were obtained. The values for f_S and f_T are consistent with zero. Assuming only Vector-Axial vector couplings, lambda_+ = 0.0288+-0.0004+-0.0011 and a good fit consistent with pure V-A couplings were obtained. Alternatively, a fit to a dipole form factor yields a pole mass of M = 859+-18 MeV, consistent with the K^*(892) mass.Comment: 16 pages, 7 figures. submitted to Phys. Lett.

A new measurement of direct CP violation in two pion decays of the neutral kaon

The NA48 experiment at CERN has performed a new measurement of direct CP violation, based on data taken in 1997 by simultaneously collecting K_L and K_S decays into pi0pi0 and pi+pi-. The result for the CP violating parameter Re(epsilon'/epsilon) is (18.5 +/- 4.5(stat)} +/- 5.8 (syst))x10^{-4}.Comment: 18 pages, 6 figure