title: New silicon photonics integration platform enabled by novel micron-scale bends

abstract: Even though submicron silicon waveguides have been proposed for dense integration of photonic devices, to date the lightwave circuits on the market mainly rely on waveguides with micron-scale core dimensions. These larger waveguides feature easier fabrication, higher reliability and better interfacing to optical fibres. Single-mode operation with large core dimensions is obtained with low lateral refractive index contrast. Hence, the main limitation in increasing the level of integration and in reducing the cost of micron-scale waveguide circuits is their mm- to cm-scale minimum bending radius. Fortunately, single-mode rib waveguides with a micron-scale silicon core can be locally transformed into multi-mode strip waveguides that have very high lateral index contrast. Here we show how Euler spiral bends realized with these waveguides can have bending radii below 10 {\mu}m and losses below 0.02 dB/90{\deg} for the fundamental mode, paving way for a novel densely integrated platform based on micron-scale waveguides.


title: On the double-ridge effect at the LHC

abstract: We discuss a possible explanation of the double-ridge observation in pPb--collisions at the LHC emphasizing that this double structure in the two-particle correlation function can result from the rotation of the transient state of matter.


title: Truncated many-body dynamics of interacting bosons: A variational principle with error monitoring

abstract: We develop a method to describe the temporal evolution of an interacting system of bosons, for which the field operator expansion is truncated after a finite number $M$ of modes, in a rigorously controlled manner. Using McLachlan's principle of least error, we find a self-consistent set of equations for the many-body state. As a particular benefit, and in distinction to previously proposed approaches, the presently introduced method facilitates the dynamical increase of the number of orbitals during the temporal evolution, due to the fact that we can rigorously monitor the error made by increasing the truncation dimension $M$. The additional orbitals, determined by the condition of least error of the truncated evolution relative to the exact one, are obtained from an initial trial state by steepest $constrained$ descent.


title: A Visual Grammar Approach for TV Program Identification

abstract: Automatic identification of TV programs within TV streams is an important task for archive exploitation. This paper proposes a new spatial-temporal approach to identify programs in TV streams in two main steps: First, a reference catalogue for video grammars of visual jingles is constructed. We exploit visual grammars characterizing instances of the same program type in order to identify the various program types in the TV stream. The role of video grammar is to represent the visual invariants for each visual jingle using a set of descriptors appropriate for each TV program. Secondly, programs in TV streams are identified by examining the similarity of the video signal to the visual grammars in the catalogue. The main idea of identification process consists in comparing the visual similarity of the video signal signature in TV stream to the catalogue elements. After presenting the proposed approach, the paper overviews the encouraging experimental results on several streams extracted from different channels and composed of several programs.


title: Quantum computer with atomic logical qubits encoded on macroscopic three-level systems in common quantum electrodynamic cavity

abstract: We propose an effective realization of the universal set of elementary quantum gates in solid state quantum computer based on macroscopic (or mesoscopic) resonance systems - multi-atomic coherent ensembles, squids or quantum dots in quantum electrodynamic cavity. We exploit an encoding of logical qubits by the pairs of the macroscopic two- or three-level atoms that is working in a Hilbert subspace of all states inherent to these atomic systems. In this subspace, logical single qubit gates are realized by the controlled reversible transfer of single atomic excitation in the pair via the exchange of virtual photons and by the frequency shift of one of the atomic ensembles in a pair. In the case of two-level systems, the logical two - qubit gates are performed by the controlling of Lamb shift magnitude in one atomic ensemble, allowing/blocking the excitation transfer in a pair, respectively, that is controlled by the third atomic system of another pair. When using three-level systems, we describe the NOT-gate in the atomic pair controlled by the transfer of working atomic excitation to the additional third level caused by direct impact of the control pair excitation. Finally, we discuss advantages of the proposed physical system for accelerated computation of some useful quantum gates.


title: A Simple Formula for Scalar Curvature of Level Sets in Euclidean Spaces

abstract: A simple formula is derived for the Ricci scalar curvature of any smooth level set ${\psi(x_0,x_1,...,x_n)=C}$ embedded in the Euclidean space $ \mathbb R^{n+1}$, in terms of the gradient $ \nabla\psi$ and the Laplacian $ \Delta\psi$. Some applications are given to the geometry of low-dimensional $p$-harmonic functions and high-dimensional harmonic functions.


title: Low-Mass X-ray Binaries in Globular Clusters: Puzzles and Solutions

abstract: In dense stellar systems, dynamical interactions between objects lead to frequent formation of exotic stellar objects, unusual binaries, and systems of higher multiplicity. They are especially important for the formation of low mass X-ray binaries (LMXBs), which are not only formed 100 times more efficiently than in the field, but also have a puzzling dependence on metallicity. In this contribution we review how compact objects are formed and retained, the mechanisms of dynamical formation and the specifics of the evolution of mass-transferring binaries with neutron stars and black holes in globular clusters - those two kinds of compact objects have different favored paths to become luminous in X-ray. We describe how stellar evolution affects ostensibly purely dynamical formation, producing the observed metallicity dependence for LMXBs. We also discuss the next puzzle to be solved on our journey to understand the link between LMXBs and millisecond pulsars formation.


title: Highlights of Supersymmetric Hypercharge $\pm1$ Triplets

abstract: The discovery of a standard model (SM)-like Higgs boson with a relatively heavy mass $m_h$ and hints of di-photon excess has deep implication to supersymmetric standard models (SSMs). We consider the SSM extended with hypercharge $\pm1$ triplets, and investigate two scenarios of it: (A) Triplets significantly couple to the Higgs doublets, which can substantially raise $m_h$ and simultaneously enhance the Higgs to di-photon rate via light chargino loops; (B) Oppositely, these couplings are quite weak and thus $m_h$ can not be raised. But the doubly-charged Higgs bosons, owing to the gauge group structure, naturally interprets why there is an excess rather than a deficient of Higgs to di-photon rate. Additionally, the pseudo Dirac triplet fermion is an inelastic non-thermal dark matter candidate. Light doubly-charged particles, especially the doubly-charged Higgs boson around 100 GeV in scenario B, are predicted. We give a preliminary discussion on their search at the LHC.


title: On cycles and coverings associated to a knot

abstract: We consider the space of all representations of the commutator subgroup of a knot group into a finite abelian group {\Sigma}, together with a shift map {\sigma}_x. This is a finite dynamical system, introduced by D.Silver and S. Williams. We describe the lengths of its cycles in terms of the roots of the Alexander polynomial of the knot. This generalizes our previous result for {\Sigma}= Z/p, p is prime, and gives a complete classification of depth 2 solvable coverings of the knot complement.


title: |V_us| from hadronic \tau{} decays

abstract: We review the status and future prospects of the determination of the CKM matrix element |V_us| using inclusive strange hadronic \tau{} decay data. We also review the results for |V_us| extracted from experimental measurements of some exclusive \tau{} decay channels such as \tau-> K \nu{} and \tau-> \pi \nu.


title: Landau theory-based estimates for viscosity coefficients of uniaxial and biaxial nematic liquid crystals

abstract: Using Landau theory, it is shown that eight phenomenological parameters are needed to describe and distinguish the twelve viscosity coefficients of a biaxial nematic phase, or the five viscosity coefficients of a uniaxial nematic phase. The dependence of the coefficients on the macroscopic uniaxial and biaxial order parameters is established. Since these order parameters are determined by the anisotropies of the dielectric constant, we show that it should be possible to determine values for all eight of the phenomenological parameters of the theory from measurements of the temperature dependence of the five viscosities of a uniaxial phase.


title: Finch-Skea star in (2+1) dimensions

abstract: The Ba\~{n}ados, Teitelboim and Zanelli\cite{BTZ} solution corresponding to the exterior space-time of a black hole in $(2+1)$ dimensions has been found to be very useful to understand various aspects relating to the gravitational field of a black hole. We present here a class of interior solutions corresponding to the BTZ exterior by making use of a model presented by Finch and Skea\cite{Finch} which was earlier found to be relevant for the description realistic stars in $(3+1)$ dimensions. We show physical viability of the model in lower dimensions as well.


title: Hybrid f(R) theories, local constraints, and cosmic speedup

abstract: We present an extension of general relativity in which an $f(R)$ term \`{a} la Palatini is added to the usual metric Einstein-Hilbert Lagrangian. Expressing the theory in a dynamically equivalent scalar-tensor form, we show that it can pass the Solar System observational tests even if the scalar field is very light or massless. Applications to cosmology and astrophysics, and some exact solutions are discussed.


title: Nanoscale Heat Transfer: from Computation to Experiment

abstract: Heat transfer can differ distinctly at the nanoscale from that at the macroscale. Recent advancement in computational and experimental techniques has enabled a large number of interesting observations and understanding of heat transfer processes at the nanoscale. In this review, we will first discuss recent advances in computational and experimental methods used in nanoscale thermal transport studies, followed by reviews of novel thermal transport phenomena at the nanoscale observed in both computational and experimental studies, and discussion on current understanding of these novel phenomena. Our perspectives on challenges and opportunities on computational and experimental methods are also presented.


title: Growth of Equally-Sized Insulin Crystals

abstract: Guidelines for growing insulin crystals of a uniform size are formulated and tested experimentally. A simple theoretical model based on the balance of matter predicts the time evolution of the crystal size and supersaturation. The time dependence of the size is checked experimentally. The experimental approach decouples crystal nucleation and growth processes according to the classical nucleation-growth-separation principle. Strict control over the nucleation process is exerted. Crystalline substance dispersity is predetermined during the nucleation stage of a batch crystallization process. To avert nutrition competition during the crystal growth stage, the number density of nucleated crystals is preset to be optimal.


title: On the identifiability of copulas in bivariate competing risks models

abstract: In competing risks models, the joint distribution of the event times is not identifiable even when the margins are fully known, which has been referred to as the "identifiability crisis in competing risks analysis" (Crowder, 1991). We model the dependence between the event times by an unknown copula and show that identification is actually possible within many frequently used families of copulas. The result is then extended to the case where one margin is unknown.


title: A real seminorm with square property is submultiplicative

abstract: A seminorm with square property on a real associative algebra is submultiplicative


title: Corrigendum: The Plebanski sectors of the EPRL vertex

abstract: We correct what amounts to a sign error in the proof of part (i.) of theorem 3 in Class.Quant.Grav.28 225003 (arXiv:1107.0709). The Plebanski sectors isolated by the linear simplicity constraints do not change --- they are still the three sectors (deg), (II+), and (II-). What changes is the characterization of the continuum Plebanski two-form corresponding to the first two terms in the asymptotics of the EPRL vertex amplitude for Regge-like boundary data. These two terms do not correspond to Plebanski sectors (II+) and (II-), but rather to the two possible signs of the product of the sign of the sector --- +1 for (II+) and -1 for (II-) --- and the sign of the orientation $\epsilon_{IJKL}B^{IJ} \wedge B^{KL}$ determined by $B^{IJ}$. This is consistent with what one would expect, as this is exactly the sign which classically relates the BF action, in Plebanski sectors (II+) and (II-), to the Einstein-Hilbert action, whose discretization is the Regge action appearing in the asymptotics.


title: Proceedings of Answer Set Programming and Other Computing Paradigms (ASPOCP 2012), 5th International Workshop, September 4, 2012, Budapest, Hungary

abstract: This volume contains the papers presented at the fifth workshop on Answer Set Programming and Other Computing Paradigms (ASPOCP 2012) held on September 4th, 2012 in Budapest, co-located with the 28th International Conference on Logic Programming (ICLP 2012). It thus continues a series of previous events co-located with ICLP, aiming at facilitating the discussion about crossing the boundaries of current ASP techniques in theory, solving, and applications, in combination with or inspired by other computing paradigms.


title: Searching for evidence of jet-cloud interaction in radio galaxies. First results for 3C 381

abstract: We present results of Gemini spectroscopy and Hubble Space Telescope imaging of the 3C~381 radio galaxy. Possible ionising mechanisms for the Extended Emission-Line Region were studied through state-of-the-art diagnostic analysis employing line-ratios. Photoionisation from the central engine as well as mixed-medium photoionisation models fail in reproducing both the strengths and the behaviour of the highest-excitation lines, such as [NeV]3424, HeII, and [OIII}]5007, which are measured at very large distances from the AGN. Shock-ionisation models provide a better fit to the observation. Expanding shocks with velocities higher than 500 km/s are capable of reaching the observed intensity ratios for lines with different ionisation states and excitation degrees. This model also provide a direct explanation of the mechanical energy input needed to explain the high-velocity line-splitting observed in the velocity field.


title: Elastic curves and surfaces under long-range forces: A geometric approach

abstract: Using classical differential geometry, the problem of elastic curves and surfaces in the presence of long-range interactions $\Phi$, is posed. Starting from a variational principle, the balance of elastic forces and the corresponding projections ${\bf n}_i\cdot \nabla\Phi$, are found. In the case of elastic surfaces, a force coupling the mean curvature with the external potential, $K\Phi$, appears; it is also present in the shape equation along the normal principal in the case of curves. The potential $\Phi$ contributes to the effective tension of curves and surfaces and also to the orbital torque. The confinement of a curve on a surface is also addressed, in such a case, the potential contributes to the normal force through the terms $-\kappa\Phi -{\bf n}\cdot \nabla\Phi$. In general, the equation of motion becomes integro-differential that must be numerically solved.


title: Estimation from Relative Measurements in Mobile Networks with Markovian Switching Topology: Clock Skew and Offset Estimation for Time Synchronization

abstract: We analyze a distributed algorithm for estimation of scalar parameters belonging to nodes in a mobile network from noisy relative measurements. The motivation comes from the problem of clock skew and offset estimation for the purpose of time synchronization. The time variation of the network was modeled as a Markov chain. The estimates are shown to be mean square convergent under fairly weak assumptions on the Markov chain, as long as the union of the graphs is connected. Expressions for the asymptotic mean and correlation are also provided. The Markovian switching topology model of mobile networks is justified for certain node mobility models through empirically estimated conditional entropy measures.


title: More Remote Globular Clusters in the Outer Halo of M31

abstract: We searched the Sloan Digital Sky Survey for outer halo globular clusters (GCs) around M31. Our search of non-stellar objects, within the limits 0.3 < (g-i)o < 1.5 and 14.0 < ro < 19.0 concentrated in some remote areas of the extended halo, to a maximum projected distance of 240 kpc, for a total of approximately 200 square degrees. Another roughly 50 square degrees, approximately 5-75 kpc from M31, were surveyed as test areas. In these areas, we identified 39 GCs and 2 GC candidates, 84% of the previously known GCs (93% of the "classical GCs" and 40% of the "halo extended clusters", on the cluster classification scheme of Huxor et al.). For the entire survey, we visually inspected 78,516 objects for morphological evidence of cluster status, and we identified 18 new clusters, and 75 candidate clusters. The new clusters include 15 classical globulars and three clusters of lower density. Six of the clusters reside in the remote areas of the outer halo, beyond projected distances of 100 kpc. Previously, only MGC1 was found beyond this limit at 117 kpc. The farthest cluster discovered in this survey lies at a projected radius of 158 kpc from M31, assuming that the M31 distance is 780 kpc.


title: Providing Probabilistic Guarantees on the Time of Information Spread in Opportunistic Networks

abstract: A variety of mathematical tools have been developed for predicting the spreading patterns in a number of varied environments including infectious diseases, computer viruses, and urgent messages broadcast to mobile agent (e.g., humans, vehicles, and mobile devices). These tools have mainly focused on estimating the average time for the spread to reach a fraction (e.g., $\alpha$) of the agents, i.e., the so-called average completion time $E(T_{\alpha})$. We claim that providing probabilistic guarantee on the time for the spread $T_{\alpha}$ rather than only its average gives a much better understanding of the spread, and hence could be used to design improved methods to prevent epidemics or devise accelerated methods for distributing data. To demonstrate the benefits, we introduce a new metric $G_{\alpha, \beta}$ that denotes the time required to guarantee $\alpha$ completion with probability $\beta$, and develop a new framework to characterize the distribution of $T_\alpha$ for various spread parameters such as number of seeds, level of contact rates, and heterogeneity in contact rates. We apply our technique to an experimental mobility trace of taxies in Shanghai and show that our framework enables us to allocate resources (i.e., to control spread parameters) for acceleration of spread in a far more efficient way than the state-of-the-art.


title: On determinants of integrable operators with shifts

abstract: Integrable integral operator can be studied by means of a matrix Riemann--Hilbert problem. However, in the case of so-called integrable operators with shifts, the associated Riemann--Hilbert problem becomes operator valued and this complicates strongly the analysis. In this note, we show how to circumvent, in a very simple way, the use of such a setting while still being able to characterize the large-$x$ asymptotic behavior of the determinant associated with the operator.


title: A graph discretization of the Laplace-Beltrami operator

abstract: We show that eigenvalues and eigenfunctions of the Laplace-Beltrami operator on a Riemannian manifold are approximated by eigenvalues and eigenvectors of a (suitably weighted) graph Laplace operator of a proximity graph on an epsilon-net.


title: Disruptions in the U.S. Airport Network

abstract: Our project analyzes the United States domestic airport network. We attempt to determine which airports are most vital in maintaining the underlying infrastructure for all domestic flights within the United States. To perform our analysis, we use data from the first quarter of 2010 and use several methods and algorithms that are frequently used in network science. Using these statistics, we identified the most important airports in the United States and investigate the role and significance that these airports play in maintaining the structure of the entire domestic airport network. Some of these airports include Denver International and Ted Stevens Anchorage International. We also identified any structural holes and suggested improvements that can be made to the network. Finally, through our analysis, we developed a disaster response algorithm that calculates flight path reroutes in emergency situations.


title: The University of Texas Millimeter Wave Observatory

abstract: This is an account of the Millimeter Wave Observatory, a 4.9 meter diameter antenna facility that pioneered continuum observations of planets and interstellar molecular spectroscopy from 1971 to 1988. The circumstances of its founding, development of its instrumentation, and major research contributions are discussed. The MWO role in training of personnel in this new field is illustrated by a listing of student and postdoctoral observers, with titles of PhD theses that included MWO data.


title: Motivic L-Function Identities from CFT and Arithmetic Mirror Symmetry

abstract: Exactly solvable mirror pairs of Calabi-Yau threefolds of hypersurface type exist in the class of Gepner models that include nondiagonal affine invariants. Motivated by the string modular interpretation established previously for models in this class it is natural to ask whether the arithmetic structure of mirror pairs of varieties reflects the fact that as conformal field theories they are isomorphic. Mirror symmetry in particular predicts that the L-functions of the \Omega-motives of such pairs are identical. In the present paper this prediction is confirmed by showing that the \Omega-motives of exactly solvable mirror pairs are isomorphic. This follows as a corollary from a more general result establishing an isomorphism between nondiagonally and diagonally induced motives in this class of varieties. The motivic approach formulated here circumvents the difficulty that no mirror construction of the Hasse-Weil zeta function is known.


title: Topological Signatures in the Electronic Structure of Graphene Spirals

abstract: Topology is familiar mostly from mathematics, but also natural sciences have found its concepts useful. Those concepts have been used to explain several natural phenomena in biology and physics, and they are particularly relevant for the electronic structure description of topological insulators and graphene systems. Here, we introduce topologically distinct graphene forms - graphene spirals - and employ density-functional theory to investigate their geometric and electronic properties. We found that the spiral topology gives rise to an intrinsic Rashba spin-orbit splitting. Through a Hamiltonian constrained by space curvature, graphene spirals have topologically protected states due to time-reversal symmetry. In addition, we argue that the synthesis of such graphene spirals is feasible and can be achieved through advanced bottom-up experimental routes that we indicate in this work.


title: A note on the "logarithmic-W_3" octuplet algebra and its Nichols algebra

abstract: We describe a Nichols-algebra-motivated construction of an octuplet chiral algebra that is a "W_3-counterpart" of the triplet algebra of (p,1) logarithmic models of two-dimensional conformal field theory.


title: Random polarization dynamics in a resonant optical medium

abstract: Random optical-pulse polarization switching along an active optical medium in the $\Lambda$-configuration with spatially disordered occupation numbers of its lower energy sub-level pair is described using the idealized integrable Maxwell-Bloch model. Analytical results describing the light polarization-switching statistics for the single self-induced transparency pulse are compared with statistics obtained from direct Monte-Carlo numerical simulations.


title: Normal modes of the d>1 fermion gas

abstract: In the present work we solve the many-fermion problem in dimensions greater than one by identifying its normal modes. We start by decomposing the Fermi sea into independent sectors labeled by a wavevector q. The product of exactly calculable eigenstates, one from each sector, is an exact eigenstate of the full many-body problem. We examine bosonization of fermions, exchange and pairing forces and quasiparticles, all from this point of view.

categories: cond-mat.str-el math-ph math.MP quant-ph


title: Off-Shell Scalar Supermultiplet in the Unfolded Dynamics Approach

abstract: We show how manifestly supersymmetric action for Wess-Zumino model can be constructed within the unfolded dynamics approach. The off-shell unfolded system for N = 1, D = 4 scalar supermultiplet is found. The action is presented in the form of integral of a closed 4-form over any (4, 0) surface in superspace as well as a superspace integral of an integral form or a chiral integral form. The proposed method is argued to provide a most general tool for the analysis of manifestly supersymmetric functionals.


title: Singular General Relativity

abstract: This work presents the foundations of Singular Semi-Riemannian Geometry and Singular General Relativity, based on the author's research. An extension of differential geometry and of Einstein's equation to singularities is reported. Singularities of the form studied here allow a smooth extension of the Einstein field equations, including matter. This applies to the Big-Bang singularity of the FLRW solution. It applies to stationary black holes, in appropriate coordinates (since the standard coordinates are singular at singularity, hiding the smoothness of the metric). In these coordinates, charged black holes have the electromagnetic potential regular everywhere. Implications on Penrose's Weyl curvature hypothesis are presented. In addition, these singularities exhibit a (geo)metric dimensional reduction, which might act as a regulator for the quantum fields, including for quantum gravity, in the UV regime. This opens the perspective of perturbative renormalizability of quantum gravity without modifying General Relativity.


title: A bright thermonuclear X-ray burst simultaneously observed with Chandra and RXTE

abstract: The prototypical accretion-powered millisecond pulsar SAX J1808.4-3658 was observed simultaneously with Chandra-LETGS and RXTE-PCA near the peak of a transient outburst in November 2011. A single thermonuclear (type-I) burst was detected, the brightest yet observed by Chandra from any source, and the second-brightest observed by RXTE. We found no evidence for discrete spectral features during the burst; absorption edges have been predicted to be present in such bursts, but may require a greater degree of photospheric expansion than the rather moderate expansion seen in this event (a factor of a few). These observations provide a unique data set to study an X-ray burst over a broad bandpass and at high spectral resolution (lambda/delta-lambda=200-400). We find a significant excess of photons at high and low energies compared to the standard black body spectrum. This excess is well described by a 20-fold increase of the persistent flux during the burst. We speculate that this results from burst photons being scattered in the accretion disk corona. These and other recent observations of X-ray bursts point out the need for detailed theoretical modeling of the radiative and hydrodynamical interaction between thermonuclear X-ray bursts and accretion disks.


title: Measurement of the ratio of differential cross sections {\sigma}(ppbar -> Z + b jet)/{\sigma}(ppbar -> Z + jet) in ppbar collisions at sqrt(s)=1.96 TeV

abstract: We measure the ratio of cross sections, {\sigma}(ppbar -> Z + b jet)/{\sigma}(ppbar -> Z + jet), for associated production of a Z boson with at least one jet. The ratio is also measured as a function of the jet transverse momentum, jet pseudorapidity, Z boson transverse momentum, and the azimuthal angle between the Z boson and the closest jet for events with at least one b jet. These measurements use data collected by the D0 experiment in Run II of Fermilab's Tevatron ppbar Collider at a center-of-mass energy of 1.96 TeV, and correspond to an integrated luminosity of 9.7 fb$^{-1}$. The results are compared to predictions from next-to-leading order calculations and various Monte Carlo event generators.


title: All-loop infrared-divergent behavior of most-subleading-color gauge-theory amplitudes

abstract: The infrared singularities of gravitational amplitudes are one-loop exact, in that higher-loop divergences are characterized by the exponential of the one-loop divergence. We show that the contributions to SU(N) gauge-theory amplitudes that are most-subleading in the 1/N expansion are also one-loop exact, provided that the dipole conjecture holds. Possible corrections to the dipole conjecture, beginning at three loops, could violate one-loop-exactness, though would still maintain the absence of collinear divergences. We also demonstrate a relation between L-loop four-point N=8 supergravity and most-subleading-color N=4 SYM amplitudes that holds for the two leading IR divergences, O(1/\epsilon^L) and O(1/\epsilon^{L-1}), but breaks down at O(1/\epsilon^{L-2}).


title: Logarithmic ^sl(2) CFT models from Nichols algebras. 1

abstract: We construct chiral algebras that centralize rank-two Nichols algebras with at least one fermionic generator. This gives "logarithmic" W-algebra extensions of a fractional-level ^sl(2) algebra. We discuss crucial aspects of the emerging general relation between Nichols algebras and logarithmic CFT models: (i) the extra input, beyond the Nichols algebra proper, needed to uniquely specify a conformal model; (ii) a relation between the CFT counterparts of Nichols algebras connected by Weyl groupoid maps; and (iii) the common double bosonization U(X) of such Nichols algebras. For an extended chiral algebra, candidates for its simple modules that are counterparts of the U(X) simple modules are proposed, as a first step toward a functorial relation between U(X) and W-algebra representation categories.


title: User Profile-Driven Data Warehouse Summary for Adaptive OLAP Queries

abstract: Data warehousing is an essential element of decision support systems. It aims at enabling the user knowledge to make better and faster daily business decisions. To improve this decision support system and to give more and more relevant information to the user, the need to integrate user's profiles into the data warehouse process becomes crucial. In this paper, we propose to exploit users' preferences as a basis for adapting OLAP (On-Line Analytical Processing) queries to the user. For this, we present a user profile-driven data warehouse approach that allows dening user's profile composed by his/her identifier and a set of his/her preferences. Our approach is based on a general data warehouse architecture and an adaptive OLAP analysis system. Our main idea consists in creating a data warehouse materialized view for each user with respect to his/her profile. This task is performed off-line when the user defines his/her profile for the first time. Then, when a user query is submitted to the data warehouse, the system deals with his/her data warehouse materialized view instead of the whole data warehouse. In other words, the data warehouse view summaries the data warehouse content for the user by taking into account his/her preferences. Moreover, we are implementing our data warehouse personalization approach under the SQL Server 2005 DBMS (DataBase Management System).


title: Wyner's Common Information: Generalizations and A New Lossy Source Coding Interpretation

abstract: Wyner's common information was originally defined for a pair of dependent discrete random variables. Its significance is largely reflected in, hence also confined to, several existing interpretations in various source coding problems. This paper attempts to both generalize its definition and to expand its practical significance by providing a new operational interpretation. The generalization is two-folded: the number of dependent variables can be arbitrary, so are the alphabet of those random variables. New properties are determined for the generalized Wyner's common information of N dependent variables. More importantly, a lossy source coding interpretation of Wyner's common information is developed using the Gray-Wyner network. In particular, it is established that the common information equals to the smallest common message rate when the total rate is arbitrarily close to the rate distortion function with joint decoding. A surprising observation is that such equality holds independent of the values of distortion constraints as long as the distortions are within some distortion region. Examples about the computation of common information are given, including that of a pair of dependent Gaussian random variables.


title: PM = EM: Partially Massless Duality Invariance

abstract: In d=4 de Sitter space, novel conformally invariant photon-like theories consistently couple to charged matter. We show that these higher spin, maximal depth, partially massless systems enjoy a Maxwellian, "electric-magnetic" duality.


title: Semilocal and Hybrid Meta-Generalized Gradient Approximations Based on the Understanding of the Kinetic-Energy-Density Dependence

abstract: We present a global hybrid meta-generalized gradient approximation (meta-GGA) with three empirical parameters, as well as its underlying semilocal meta-GGA and a meta-GGA with only one empirical parameter. All of them are based on the new meta-GGA resulting from the understanding of kinetic-energy-density dependence [J. Chem. Phys. 137, 051101 (2012)]. The obtained functionals show robust performances on the considered molecular systems for the properties of heats of formation, barrier heights, and noncovalent interactions. The pair-wise additive dispersion corrections to the functionals are also presented.


title: Real Spectral Triples over Noncommutative Bieberbach Manifolds

abstract: We classify and construct all real spectral triples over noncommutative Bieberbach manifolds, which are restrictions of irreducible real equivariant spectral triple over the noncommutative three-torus. We show that in the classical case the constructed geometries correspond exactly to spin structures over Bieberbach manifolds and the Dirac operators constructed for a flat metric.


title: Effect of Eu doping and partial oxygen isotope substitution on magnetic phase transitions in (Pr$_{1-y}$Eu$_y$)$_{0.7}$Ca$_{0.3}$CoO$_3$ cobaltites

abstract: We study experimentally and theoretically the effect of Eu doping and partial oxygen isotope substitution on the transport and magnetic characteristics and spin-state transitions in (Pr$_{1-y}$Eu$_y$)$_{0.7}$Ca$_{0.3}$CoO$_3$ cobaltites. The Eu doping level $y$ is chosen in the range of the phase diagram near the crossover between the ferromagnetic and spin-state transitions ($0.10 < y < 0.20$). We prepared a series of samples with different degrees of enrichment by the heavy oxygen isotope $^{18}$O, namely, containing 90%, 67%, 43%, 17%, and 0% of $^{18}$O. Based on the measurements of ac magnetic susceptibility $\chi(T)$ and electrical resistivity $\rho(T)$, we analyze the evolution of the sample properties with the change of Eu and $^{18}$O content. It is demonstrated that the effect of increasing $^{18}$O content on the system is similar to that of increasing the Eu content. The band structure calculations of the energy gap between $t_{2g}$ and $e_g$ bands including the renormalization of this gap due to the electron-phonon interaction reveal the physical mechanisms underlying such similarity.


title: Maximal gaps between prime k-tuples: a statistical approach

abstract: Combining the Hardy-Littlewood k-tuple conjecture with a heuristic application of extreme-value statistics, we propose a family of estimator formulas for predicting maximal gaps between prime k-tuples. Computations show that the estimator a(log(x/a)-b) satisfactorily predicts the maximal gaps below x, where a is the expected average gap between the same type of k-tuples, a=O(log^k x). Heuristics suggest that maximal gaps between prime k-tuples near x are approximately a*log(x/a), and thus have the order O(log^{k+1}x). The distribution of maximal gaps around the trend curve a*log(x/a) is close to the Gumbel distribution. We explore two implications of this model of gaps: record gaps between primes and Legendre-type conjectures for prime k-tuples.


title: Bernstein-Gelfand-Gelfand resolutions for basic classical Lie superalgebras

abstract: We study Kostant cohomology and Bernstein-Gelfand-Gelfand resolutions for finite dimensional representations of basic classical Lie superalgebras and reductive Lie superalgebras based on them. For each choice of parabolic subalgebra and irreducible representation of such a Lie superalgebra, there is a natural definition of the derivative and coderivative, which define the (co)homology groups. We prove that a necessary condition to have a resolution of an irreducible module in terms of Verma modules is complete reducibility of the cohomology groups. Essentially, if it exists, every such a resolution is then given by modules induced by these cohomology groups. We also prove that if these cohomology groups are completely reducible, a sufficient condition for the existence of such a resolution is that these groups are isomorphic to the kernel of the Kostant quabla operator, which is equivalent with disjointness of the derivative and coderivative. Then we use these results to derive very explicit criteria under which BGG resolutions exist, which are particularly useful for the superalgebras of type I. For the unitarisable representations of gl(m|n) and osp(2|2n) we derive conditions on the parabolic subalgebra under which the BGG resolutions exist. This extends the BGG resolutions for gl(m|n) previously obtained through superduality and leads to entirely new results for osp(2|2n). We also apply the obtained theory to construct specific examples of BGG resolutions for osp(m|2n).


title: CDF results on $b \to s \mu \mu$ decays

abstract: We summarize recent b -> s mu mu results from the CDF experiment based on the full 9.6 fb-1 dataset collected in p-pbar collisions at sqrt(s) = 1.96 TeV. The decays analyzed are; B+ -> K+ mu+ mu-, B0 -> K*0(892) mu+ mu-, B0 -> K0_s mu+ mu-, B+ -> K*+(892) mu+ mu-, B0_s -> phi mu+ mu-, and Lambda0_b -> Lambda mu+ mu-. The latter two decays were first observed by CDF in 2011. From an angular analysis of the B -> K* mu+ mu- decays we also present updated results on the transverse polarization and T-odd CP asymmetries reported earlier.


title: Quantum-Reduced Loop Gravity: Cosmology

abstract: We introduce a new framework for loop quantum gravity: mimicking the spinfoam quantization procedure we propose to study the symmetric sectors of the theory imposing the reduction weakly on the full kinematical Hilbert space of the canonical theory. As a first application of Quantum-Reduced Loop Gravity we study the inhomogeneous Bianchi I model. The emerging quantum cosmological model represents a simplified arena on which the complete canonical quantization program can be tested. The achievements of this analysis could elucidate the relationship between Loop Quantum Cosmology and the full theory.


title: Precise evaluation of thermal response functions by optimized density matrix renormalization group schemes

abstract: This paper provides a study and discussion of earlier as well as novel more efficient schemes for the precise evaluation of finite-temperature response functions of strongly correlated quantum systems in the framework of the time-dependent density matrix renormalization group (tDMRG). The computational costs and bond dimensions as functions of time and temperature are examined for the example of the spin-1/2 XXZ Heisenberg chain in the critical XY phase and the gapped N\'eel phase. The matrix product state purifications occurring in the algorithms are in one-to-one relation with corresponding matrix product operators. This notational simplification elucidates implications of quasi-locality on the computational costs. Based on the observation that there is considerable freedom in designing efficient tDMRG schemes for the calculation of dynamical correlators at finite temperatures, a new class of optimizable schemes, as recently suggested in arXiv:1212.3570, is explained and analyzed numerically. A specific novel near-optimal scheme that requires no additional optimization reaches maximum times that are typically increased by a factor of two, when compared against earlier approaches. These increased reachable times make many more physical applications accessible. For each of the described tDMRG schemes, one can devise a corresponding transfer matrix renormalization group (TMRG) variant.


title: Evolutionary dynamics of group interactions on structured populations: A review

abstract: Interactions among living organisms, from bacteria colonies to human societies, are inherently more complex than interactions among particles and nonliving matter. Group interactions are a particularly important and widespread class, representative of which is the public goods game. In addition, methods of statistical physics have proven valuable for studying pattern formation, equilibrium selection, and self-organisation in evolutionary games. Here we review recent advances in the study of evolutionary dynamics of group interactions on structured populations, including lattices, complex networks and coevolutionary models. We also compare these results with those obtained on well-mixed populations. The review particularly highlights that the study of the dynamics of group interactions, like several other important equilibrium and non-equilibrium dynamical processes in biological, economical and social sciences, benefits from the synergy between statistical physics, network science and evolutionary game theory.

categories: physics.soc-ph cond-mat.stat-mech cs.SI nlin.AO q-bio.PE


title: Spoke-like differential rotation in a convective dynamo with a coronal envelope

abstract: We report on the results of four convective dynamo simulations with an outer coronal layer. The magnetic field is self-consistently generated by the convective motions beneath the surface. Above the convection zone, we include a polytropic layer that extends to 1.6 solar radii. The temperature increases in this region to $\approx8$ times the value at the surface, corresponding to $\approx1.2$ times the value at the bottom of the spherical shell. We associate this region with the solar corona. We find solar-like differential rotation with radial contours of constant rotation rate, together with a near-surface shear layer. This non-cylindrical rotation profile is caused by a non-zero latitudinal entropy gradient that offsets the Taylor--Proudman balance through the baroclinic term. The meridional circulation is multi-cellular with a solar-like poleward flow near the surface at low latitudes. In most of the cases, the mean magnetic field is oscillatory with equatorward migration in two cases. In other cases, the equatorward migration is overlaid by stationary or even poleward migrating mean fields.


title: Superloop Space

abstract: In this paper will construct and analyse the superloop space formulation of a $\mathcal{N} =1$ supergauge theory in three dimensions. We will obtain expressions for the connection and the curvature in this superloop space in terms of ordinary supergauge fields. This curvature will vanish, unless there is a monopole in the spacetime. We will also construct a quantity which will give the monopole charge in this formalism. Finally, we will show how these results even hold for a deformed superspace.


title: Same-sign Tops: A Powerful Diagnostic Test for Models of New Physics

abstract: We study the connection between the same sign top (SST) and the top quark forward-backward asymmetry $A^t_{FB}$. We find that a large class of new physics models that have been proposed to account for the $A^t_{FB}$ lead to SST quark production rate much larger than the observed rate at the LHC and consequently are severely constrained or ruled out. Our model independent, general, operator analysis shows that none of the tree-level flavor-changing operators are able to explain $A^t_{FB}$ and simultaneously remain consistent with the same-sign top-quark production constraints from the LHC data.


title: Anisotropic Universe Models with Perfect Fluid and Scalar Field in f(R,T) Gravity

abstract: In this paper, we study the behavior of perfect fluid and massless scalar field for homogeneous and anisotropic Bianchi type I universe model in $f(R,T)$ gravity, where $R$ is the Ricci scalar and $T$ is the trace of the energy-momentum tensor. We assume the variation law of mean Hubble parameter to obtain exact solutions of the modified field equations. The physical and kinematical quantities are discussed for both models in future evolution of the universe. We check the validity of null energy condition and conclude that our perfect fluid solution can behave like phantom model. Finally, we find that perfect fluid solutions correspond to massless scalar field models.


title: A Factorized Variational Technique for Phase Unwrapping in Markov Random Fields

abstract: Some types of medical and topographic imaging device produce images in which the pixel values are "phase-wrapped", i.e. measured modulus a known scalar. Phase unwrapping can be viewed as the problem of inferring the number of shifts between each and every pair of neighboring pixels, subject to an a priori preference for smooth surfaces, and subject to a zero curl constraint, which requires that the shifts must sum to 0 around every loop. We formulate phase unwrapping as a mean field inference problem in a Markov network, where the prior favors the zero curl constraint. We compare our mean field technique with the least squares method on a synthetic 100x100 image, and give results on a 512x512 synthetic aperture radar image from Sandia National Laboratories.<Long Text>


title: Efficient Approximation for Triangulation of Minimum Treewidth

abstract: We present four novel approximation algorithms for finding triangulation of minimum treewidth. Two of the algorithms improve on the running times of algorithms by Robertson and Seymour, and Becker and Geiger that approximate the optimum by factors of 4 and 3 2/3, respectively. A third algorithm is faster than those but gives an approximation factor of 4 1/2. The last algorithm is yet faster, producing factor-O(lg/k) approximations in polynomial time. Finding triangulations of minimum treewidth for graphs is central to many problems in computer science. Real-world problems in artificial intelligence, VLSI design and databases are efficiently solvable if we have an efficient approximation algorithm for them. We report on experimental results confirming the effectiveness of our algorithms for large graphs associated with real-world problems.


title: Markov Chain Monte Carlo using Tree-Based Priors on Model Structure

abstract: We present a general framework for defining priors on model structure and sampling from the posterior using the Metropolis-Hastings algorithm. The key idea is that structure priors are defined via a probability tree and that the proposal mechanism for the Metropolis-Hastings algorithm operates by traversing this tree, thereby defining a cheaply computable acceptance probability. We have applied this approach to Bayesian net structure learning using a number of priors and tree traversal strategies. Our results show that these must be chosen appropriately for this approach to be successful.


title: Graphical readings of possibilistic logic bases

abstract: Possibility theory offers either a qualitive, or a numerical framework for representing uncertainty, in terms of dual measures of possibility and necessity. This leads to the existence of two kinds of possibilistic causal graphs where the conditioning is either based on the minimum, or the product operator. Benferhat et al. (1999) have investigated the connections between min-based graphs and possibilistic logic bases (made of classical formulas weighted in terms of certainty). This paper deals with a more difficult issue : the product-based graphical representations of possibilistic bases, which provides an easy structural reading of possibilistic bases. Moreover, this paper also provides another reading of possibilistic bases in terms of comparative preferences of the form "in the context p, q is preferred to not q". This enables us to explicit preferences underlying a set of goals with different levels of priority.


title: Pre-processing for Triangulation of Probabilistic Networks

abstract: The currently most efficient algorithm for inference with a probabilistic network builds upon a triangulation of a network's graph. In this paper, we show that pre-processing can help in finding good triangulations forprobabilistic networks, that is, triangulations with a minimal maximum clique size. We provide a set of rules for stepwise reducing a graph, without losing optimality. This reduction allows us to solve the triangulation problem on a smaller graph. From the smaller graph's triangulation, a triangulation of the original graph is obtained by reversing the reduction steps. Our experimental results show that the graphs of some well-known real-life probabilistic networks can be triangulated optimally just by preprocessing; for other networks, huge reductions in their graph's size are obtained.


title: A Calculus for Causal Relevance

abstract: This paper presents a sound and completecalculus for causal relevance, based onPearl's functional models semantics.The calculus consists of axioms and rulesof inference for reasoning about causalrelevance relationships.We extend the set of known axioms for causalrelevance with three new axioms, andintroduce two new rules of inference forreasoning about specific subclasses ofmodels.These subclasses give a more refinedcharacterization of causal models than the one given in Halpern's axiomatizationof counterfactual reasoning.Finally, we show how the calculus for causalrelevance can be used in the task ofidentifying causal structure from non-observational data.


title: Instrumentality Tests Revisited

abstract: An instrument is a random variable thatallows the identification of parameters inlinear models when the error terms arenot uncorrelated.It is a popular method used in economicsand the social sciences that reduces theproblem of identification to the problemof finding the appropriate instruments.Few years ago, Pearl introduced a necessarytest for instruments that allows the researcher to discard those candidatesthat fail the test.In this paper, we make a detailed study of Pearl's test and the general model forinstruments. The results of this studyinclude a novel interpretation of Pearl'stest, a general theory of instrumentaltests, and an affirmative answer to aprevious conjecture. We also presentnew instrumentality tests for the casesof discrete and continuous variables.


title: UCP-Networks: A Directed Graphical Representation of Conditional Utilities

abstract: We propose a new directed graphical representation of utility functions, called UCP-networks, that combines aspects of two existing graphical models: generalized additive models and CP-networks. The network decomposes a utility function into a number of additive factors, with the directionality of the arcs reflecting conditional dependence of preference statements - in the underlying (qualitative) preference ordering - under a {em ceteris paribus} (all else being equal) interpretation. This representation is arguably natural in many settings. Furthermore, the strong CP-semantics ensures that computation of optimization and dominance queries is very efficient. We also demonstrate the value of this representation in decision making. Finally, we describe an interactive elicitation procedure that takes advantage of the linear nature of the constraints on "`tradeoff weights" imposed by a UCP-network. This procedure allows the network to be refined until the regret of the decision with minimax regret (with respect to the incompletely specified utility function) falls below a specified threshold (e.g., the cost of further questioning.


title: Confidence Inference in Bayesian Networks

abstract: We present two sampling algorithms for probabilistic confidence inference in Bayesian networks. These two algorithms (we call them AIS-BN-mu and AIS-BN-sigma algorithms) guarantee that estimates of posterior probabilities are with a given probability within a desired precision bound. Our algorithms are based on recent advances in sampling algorithms for (1) estimating the mean of bounded random variables and (2) adaptive importance sampling in Bayesian networks. In addition to a simple stopping rule for sampling that they provide, the AIS-BN-mu and AIS-BN-sigma algorithms are capable of guiding the learning process in the AIS-BN algorithm. An empirical evaluation of the proposed algorithms shows excellent performance, even for very unlikely evidence.


title: Semi-Instrumental Variables: A Test for Instrument Admissibility

abstract: In a causal graphical model, an instrument for a variable X and its effect Y is a random variable that is a cause of X and independent of all the causes of Y except X. (Pearl (1995), Spirtes et al (2000)). Instrumental variables can be used to estimate how the distribution of an effect will respond to a manipulation of its causes, even in the presence of unmeasured common causes (confounders). In typical instrumental variable estimation, instruments are chosen based on domain knowledge. There is currently no statistical test for validating a variable as an instrument. In this paper, we introduce the concept of semi-instrument, which generalizes the concept of instrument. We show that in the framework of additive models, under certain conditions, we can test whether a variable is semi-instrumental. Moreover, adding some distribution assumptions, we can test whether two semi-instruments are instrumental. We give algorithms to estimate the p-value that a random variable is semi-instrumental, and the p-value that two semi-instruments are both instrumental. These algorithms can be used to test the experts' choice of instruments, or to identify instruments automatically.


title: Conditions Under Which Conditional Independence and Scoring Methods Lead to Identical Selection of Bayesian Network Models

abstract: It is often stated in papers tackling the task of inferring Bayesian network structures from data that there are these two distinct approaches: (i) Apply conditional independence tests when testing for the presence or otherwise of edges; (ii) Search the model space using a scoring metric. Here I argue that for complete data and a given node ordering this division is a myth, by showing that cross entropy methods for checking conditional independence are mathematically identical to methods based upon discriminating between models by their overall goodness-of-fit logarithmic scores.


title: Linearity Properties of Bayes Nets with Binary Variables

abstract: It is "well known" that in linear models: (1) testable constraints on the marginal distribution of observed variables distinguish certain cases in which an unobserved cause jointly influences several observed variables; (2) the technique of "instrumental variables" sometimes permits an estimation of the influence of one variable on another even when the association between the variables may be confounded by unobserved common causes; (3) the association (or conditional probability distribution of one variable given another) of two variables connected by a path or trek can be computed directly from the parameter values associated with each edge in the path or trek; (4) the association of two variables produced by multiple treks can be computed from the parameters associated with each trek; and (5) the independence of two variables conditional on a third implies the corresponding independence of the sums of the variables over all units conditional on the sums over all units of each of the original conditioning variables.These properties are exploited in search procedures. It is also known that properties (2)-(5) do not hold for all Bayes nets with binary variables. We show that (1) holds for all Bayes nets with binary variables and (5) holds for all singly trek-connected Bayes nets of that kind. We further show that all five properties hold for Bayes nets with any DAG and binary variables parameterized with noisy-or and noisy-and gates.


title: Using Bayesian Networks to Identify the Causal Effect of Speeding in Individual Vehicle/Pedestrian Collisions

abstract: On roads showing significant violations of posted speed limits, one measure of the safety effect of speeding is the difference between the road's actual accident count and the count that would have occurred if the posted speed limit had been strictly obeyed. An estimate of this accident reduction can be had by computing the probability that speeding was a necessary condition for each of set of accidents. This is an instance of assessing individual probabilities of causation, which is generally not possible absent prior knowledge of causal structure. For traffic accidents such prior knowledge is often available and this paper illustrates how, for a commonly occurring class of vehicle/pedestrian accidents, approaches to uncertainty and causal analyses appearing in the accident reconstruction literature can be unified using Bayesian networks. Measured skidmarks, pedestrian throw distances, and pedestrian injury severity are treated as evidence, and using the Gibbs Sampling routine BUGS, the posterior probability distribution over exogenous variables, such as the vehicle's initial speed, location, and driver reaction time, is computed. This posterior distribution is then used to compute the "probability of necessity" for speeding.


title: Hybrid Processing of Beliefs and Constraints

abstract: This paper explores algorithms for processing probabilistic and deterministic information when the former is represented as a belief network and the latter as a set of boolean clauses. The motivating tasks are 1. evaluating beliefs networks having a large number of deterministic relationships and2. evaluating probabilities of complex boolean querie over a belief network. We propose a parameterized family of variable elimination algorithms that exploit both types of information, and that allows varying levels of constraint propagation inferences. The complexity of the scheme is controlled by the induced-width of the graph {em augmented} by the dependencies introduced by the boolean constraints. Preliminary empirical evaluation demonstrate the effect of constraint propagation on probabilistic computation.


title: Variational MCMC

abstract: We propose a new class of learning algorithms that combines variational approximation and Markov chain Monte Carlo (MCMC) simulation. Naive algorithms that use the variational approximation as proposal distribution can perform poorly because this approximation tends to underestimate the true variance and other features of the data. We solve this problem by introducing more sophisticated MCMC algorithms. One of these algorithms is a mixture of two MCMC kernels: a random walk Metropolis kernel and a blockMetropolis-Hastings (MH) kernel with a variational approximation as proposaldistribution. The MH kernel allows one to locate regions of high probability efficiently. The Metropolis kernel allows us to explore the vicinity of these regions. This algorithm outperforms variationalapproximations because it yields slightly better estimates of the mean and considerably better estimates of higher moments, such as covariances. It also outperforms standard MCMC algorithms because it locates theregions of high probability quickly, thus speeding up convergence. We demonstrate this algorithm on the problem of Bayesian parameter estimation for logistic (sigmoid) belief networks.


title: Efficient Stepwise Selection in Decomposable Models

abstract: In this paper, we present an efficient way of performing stepwise selection in the class of decomposable models. The main contribution of the paper is a simple characterization of the edges that canbe added to a decomposable model while keeping the resulting model decomposable and an efficient algorithm for enumerating all such edges for a given model in essentially O(1) time per edge. We also discuss how backward selection can be performed efficiently using our data structures.We also analyze the complexity of the complete stepwise selection procedure, including the complexity of choosing which of the eligible dges to add to (or delete from) the current model, with the aim ofminimizing the Kullback-Leibler distance of the resulting model from the saturated model for the data.


title: Incorporating Expressive Graphical Models in Variational Approximations: Chain-Graphs and Hidden Variables

abstract: Global variational approximation methods in graphical models allow efficient approximate inference of complex posterior distributions by using a simpler model. The choice of the approximating model determines a tradeoff between the complexity of the approximation procedure and the quality of the approximation. In this paper, we consider variational approximations based on two classes of models that are richer than standard Bayesian networks, Markov networks or mixture models. As such, these classes allow to find better tradeoffs in the spectrum of approximations. The first class of models are chain graphs, which capture distributions that are partially directed. The second class of models are directed graphs (Bayesian networks) with additional latent variables. Both classes allow representation of multi-variable dependencies that cannot be easily represented within a Bayesian network.


title: Learning the Dimensionality of Hidden Variables

abstract: A serious problem in learning probabilistic models is the presence of hidden variables. These variables are not observed, yet interact with several of the observed variables. Detecting hidden variables poses two problems: determining the relations to other variables in the model and determining the number of states of the hidden variable. In this paper, we address the latter problem in the context of Bayesian networks. We describe an approach that utilizes a score-based agglomerative state-clustering. As we show, this approach allows us to efficiently evaluate models with a range of cardinalities for the hidden variable. We show how to extend this procedure to deal with multiple interacting hidden variables. We demonstrate the effectiveness of this approach by evaluating it on synthetic and real-life data. We show that our approach learns models with hidden variables that generalize better and have better structure than previous approaches.


title: Multivariate Information Bottleneck

abstract: The Information bottleneck method is an unsupervised non-parametric data organization technique. Given a joint distribution P(A,B), this method constructs a new variable T that extracts partitions, or clusters, over the values of A that are informative about B. The information bottleneck has already been applied to document classification, gene expression, neural code, and spectral analysis. In this paper, we introduce a general principled framework for multivariate extensions of the information bottleneck method. This allows us to consider multiple systems of data partitions that are inter-related. Our approach utilizes Bayesian networks for specifying the systems of clusters and what information each captures. We show that this construction provides insight about bottleneck variations and enables us to characterize solutions of these variations. We also present a general framework for iterative algorithms for constructing solutions, and apply it to several examples.


title: A Comparison of Axiomatic Approaches to Qualitative Decision Making Using Possibility Theory

abstract: In this paper we analyze two recent axiomatic approaches proposed by Dubois et al and by Giang and Shenoy to qualitative decision making where uncertainty is described by possibility theory. Both axiomtizations are inspired by von Neumann and Morgenstern's system of axioms for the case of probability theory. We show that our approach naturally unifies two axiomatic systems that correspond respectively to pessimistic and optimistic decision criteria proposed by Dubois et al. The simplifying unification is achieved by (i) replacing axioms that are supposed to reflect two informational attitudes (uncertainty aversion and uncertainty attraction) by an axiom that imposes order on set of standard lotteries and (ii) using a binary utility scale in which each utility level is represented by a pair of numbers.


title: Enumerating Markov Equivalence Classes of Acyclic Digraph Models

abstract: Graphical Markov models determined by acyclic digraphs (ADGs), also called directed acyclic graphs (DAGs), are widely studied in statistics, computer science (as Bayesian networks), operations research (as influence diagrams), and many related fields. Because different ADGs may determine the same Markov equivalence class, it long has been of interest to determine the efficiency gained in model specification and search by working directly with Markov equivalence classes of ADGs rather than with ADGs themselves. A computer program was written to enumerate the equivalence classes of ADG models as specified by Pearl & Verma's equivalence criterion. The program counted equivalence classes for models up to and including 10 vertices. The ratio of number of classes to ADGs appears to approach an asymptote of about 0.267. Classes were analyzed according to number of edges and class size. By edges, the distribution of number of classes approaches a Gaussian shape. By class size, classes of size 1 are most common, with the proportions for larger sizes initially decreasing but then following a more irregular pattern. The maximum number of classes generated by any undirected graph was found to increase approximately factorially. The program also includes a new variation of orderly algorithm for generating undirected graphs.


title: Robust Combination of Local Controllers

abstract: Planning problems are hard, motion planning, for example, isPSPACE-hard. Such problems are even more difficult in the presence of uncertainty. Although, Markov Decision Processes (MDPs) provide a formal framework for such problems, finding solutions to high dimensional continuous MDPs is usually difficult, especially when the actions and time measurements are continuous. Fortunately, problem-specific knowledge allows us to design controllers that are good locally, though having no global guarantees. We propose a method of nonparametrically combining local controllers to obtain globally good solutions. We apply this formulation to two types of problems : motion planning (stochastic shortest path) and discounted MDPs. For motion planning, we argue that usual MDP optimality criterion (expected cost) may not be practically relevant. Wepropose an alternative: finding the minimum cost path,subject to the constraint that the robot must reach the goal withhigh probability. For this problem, we prove that a polynomial number of samples is sufficient to obtain a high probability path. For discounted MDPs, we propose a formulation that explicitly deals with model uncertainty, i.e., the problem introduced when transition probabilities are not known exactly. We formulate the problem as a robust linear program which directly incorporates this type of uncertainty.


title: Similarity Measures on Preference Structures, Part II: Utility Functions

abstract: In previous work cite{Ha98:Towards} we presented a case-based approach to eliciting and reasoning with preferences. A key issue in this approach is the definition of similarity between user preferences. We introduced the probabilistic distance as a measure of similarity on user preferences, and provided an algorithm to compute the distance between two partially specified {em value} functions. This is for the case of decision making under {em certainty}. In this paper we address the more challenging issue of computing the probabilistic distance in the case of decision making under{em uncertainty}. We provide an algorithm to compute the probabilistic distance between two partially specified {em utility} functions. We demonstrate the use of this algorithm with a medical data set of partially specified patient preferences,where none of the other existing distancemeasures appear definable. Using this data set, we also demonstrate that the case-based approach to preference elicitation isapplicable in domains with uncertainty. Finally, we provide a comprehensive analytical comparison of the probabilistic distance with some existing distance measures on preferences.


title: Causes and Explanations: A Structural-Model Approach --- Part 1: Causes

abstract: We propose a new definition of actual causes, using structural equations to model counterfactuals.We show that the definitions yield a plausible and elegant account ofcausation that handles well examples which have caused problems forother definitions and resolves major difficulties in the traditionalaccount. In a companion paper, we show how the definition of causality can beused to give an elegant definition of (causal) explanation.


title: A Dynamic Programming Model for Determining Bidding Strategies in Sequential Auctions: Quasi-linear Utility and Budget Constraints

abstract: In this paper, we develop a new method for finding an optimal biddingstrategy in sequential auctions, using a dynamic programming technique. Theexisting method assumes that the utility of a user is represented in anadditive form. Thus, the remaining endowment of money must be explicitlyrepresented in each state, and the calculation of the optimal biddingstrategy becomes time-consuming when the initial endowment of money mbecomes large.In this paper, we develop a new problem formalization that avoids explicitlyrepresenting the remaining endowment, by assuming the utility of a user canbe represented in a quasi-linear form, and representing the payment as astate-transition cost. Experimental evaluations show that we can obtainmore than an m-fold speed-up in the computation time. Furthermore, we havedeveloped a method for obtaining a semi-optimal bidding strategy underbudget constraints, and have experimentally confirmed the efficacy of thismethod.


title: A Clustering Approach to Solving Large Stochastic Matching Problems

abstract: In this work we focus on efficient heuristics for solving a class of stochastic planning problems that arise in a variety of business, investment, and industrial applications. The problem is best described in terms of future buy and sell contracts. By buying less reliable, but less expensive, buy (supply) contracts, a company or a trader can cover a position of more reliable and more expensive sell contracts. The goal is to maximize the expected net gain (profit) by constructing a dose to optimum portfolio out of the available buy and sell contracts. This stochastic planning problem can be formulated as a two-stage stochastic linear programming problem with recourse. However, this formalization leads to solutions that are exponential in the number of possible failure combinations. Thus, this approach is not feasible for large scale problems. In this work we investigate heuristic approximation techniques alleviating the efficiency problem. We primarily focus on the clustering approach and devise heuristics for finding clusterings leading to good approximations. We illustrate the quality and feasibility of the approach through experimental data.


title: Discovering Multiple Constraints that are Frequently Approximately Satisfied

abstract: Some high-dimensional data.sets can be modelled by assuming that there are many different linear constraints, each of which is Frequently Approximately Satisfied (FAS) by the data. The probability of a data vector under the model is then proportional to the product of the probabilities of its constraint violations. We describe three methods of learning products of constraints using a heavy-tailed probability distribution for the violations.


title: A Bayesian Approach to Tackling Hard Computational Problems

abstract: We are developing a general framework for using learned Bayesian models for decision-theoretic control of search and reasoningalgorithms. We illustrate the approach on the specific task of controlling both general and domain-specific solvers on a hard class of structured constraint satisfaction problems. A successful strategyfor reducing the high (and even infinite) variance in running time typically exhibited by backtracking search algorithms is to cut off and restart the search if a solution is not found within a certainamount of time. Previous work on restart strategies have employed fixed cut off values. We show how to create a dynamic cut off strategy by learning a Bayesian model that predicts the ultimate length of a trial based on observing the early behavior of the search algorithm. Furthermore, we describe the general conditions under which a dynamic restart strategy can outperform the theoretically optimal fixed strategy.


title: Estimating Well-Performing Bayesian Networks using Bernoulli Mixtures

abstract: A novel method for estimating Bayesian network (BN) parameters from data is presented which provides improved performance on test data. Previous research has shown the value of representing conditional probability distributions (CPDs) via neural networks(Neal 1992), noisy-OR gates (Neal 1992, Diez 1993)and decision trees (Friedman and Goldszmidt 1996).The Bernoulli mixture network (BMN) explicitly represents the CPDs of discrete BN nodes as mixtures of local distributions,each having a different set of parents.This increases the space of possible structures which can be considered,enabling the CPDs to have finer-grained dependencies.The resulting estimation procedure induces a modelthat is better able to emulate the underlying interactions occurring in the data than conventional conditional Bernoulli network models.The results for artificially generated data indicate that overfitting is best reduced by restricting the complexity of candidate mixture substructures local to each node. Furthermore, mixtures of very simple substructures can perform almost as well as more complex ones.The BMN is also applied to data collected from an online adventure game with an application to keyhole plan recognition. The results show that the BMN-based model brings a dramatic improvement in performance over a conventional BN model.


title: Graphical Models for Game Theory

abstract: In this work, we introduce graphical modelsfor multi-player game theory, and give powerful algorithms for computing their Nash equilibria in certain cases. An n-player game is given by an undirected graph on n nodes and a set of n local matrices. The interpretation is that the payoff to player i is determined entirely by the actions of player i and his neighbors in the graph, and thus the payoff matrix to player i is indexed only by these players. We thus view the global n-player game as being composed of interacting local games, each involving many fewer players. Each player's action may have global impact, but it occurs through the propagation of local influences.Our main technical result is an efficient algorithm for computing Nash equilibria when the underlying graph is a tree (or can be turned into a tree with few node mergings). The algorithm runs in time polynomial in the size of the representation (the graph and theassociated local game matrices), and comes in two related but distinct flavors. The first version involves an approximation step, and computes a representation of all approximate Nash equilibria (of which there may be an exponential number in general). The second version allows the exact computation of Nash equilibria at the expense of weakened complexity bounds. The algorithm requires only local message-passing between nodes (and thus can be implemented by the players themselves in a distributed manner). Despite an analogy to inference in Bayes nets that we develop, the analysis of our algorithm is more involved than that for the polytree algorithm in, owing partially to the fact that we must either compute, or select from, an exponential number of potential solutions. We discuss a number of extensions, such as the computation of equilibria with desirable global properties (e.g. maximizing global return), and directions for further research.


title: On characterizing Inclusion of Bayesian Networks

abstract: Every directed acyclic graph (DAG) over a finite non-empty set of variables (= nodes) N induces an independence model over N, which is a list of conditional independence statements over N.The inclusion problem is how to characterize (in graphical terms) whether all independence statements in the model induced by a DAG K are in the model induced by a second DAG L. Meek (1997) conjectured that this inclusion holds iff there exists a sequence of DAGs from L to K such that only certain 'legal' arrow reversal and 'legal' arrow adding operations are performed to get the next DAG in the sequence.In this paper we give several characterizations of inclusion of DAG models and verify Meek's conjecture in the case that the DAGs K and L differ in at most one adjacency. As a warming up a rigorous proof of well-known graphical characterizations of equivalence of DAGs, which is a highly related problem, is given.


title: Improved learning of Bayesian networks

abstract: The search space of Bayesian Network structures is usually defined as Acyclic Directed Graphs (DAGs) and the search is done by local transformations of DAGs. But the space of Bayesian Networks is ordered by DAG Markov model inclusion and it is natural to consider that a good search policy should take this into account. First attempt to do this (Chickering 1996) was using equivalence classes of DAGs instead of DAGs itself. This approach produces better results but it is significantly slower. We present a compromise between these two approaches. It uses DAGs to search the space in such a way that the ordering by inclusion is taken into account. This is achieved by repetitive usage of local moves within the equivalence class of DAGs. We show that this new approach produces better results than the original DAGs approach without substantial change in time complexity. We present empirical results, within the framework of heuristic search and Markov Chain Monte Carlo, provided through the Alarm dataset.


title: Classifier Learning with Supervised Marginal Likelihood

abstract: It has been argued that in supervised classification tasks, in practice it may be more sensible to perform model selection with respect to some more focused model selection score, like the supervised (conditional) marginal likelihood, than with respect to the standard marginal likelihood criterion. However, for most Bayesian network models, computing the supervised marginal likelihood score takes exponential time with respect to the amount of observed data. In this paper, we consider diagnostic Bayesian network classifiers where the significant model parameters represent conditional distributions for the class variable, given the values of the predictor variables, in which case the supervised marginal likelihood can be computed in linear time with respect to the data. As the number of model parameters grows in this case exponentially with respect to the number of predictors, we focus on simple diagnostic models where the number of relevant predictors is small, and suggest two approaches for applying this type of models in classification. The first approach is based on mixtures of simple diagnostic models, while in the second approach we apply the small predictor sets of the simple diagnostic models for augmenting the Naive Bayes classifier.


title: Plausible reasoning from spatial observations

abstract: This article deals with plausible reasoning from incomplete knowledge about large-scale spatial properties. The availableinformation, consisting of a set of pointwise observations,is extrapolated to neighbour points. We make use of belief functions to represent the influence of the knowledge at a given point to another point; the quantitative strength of this influence decreases when the distance between both points increases. These influences arethen aggregated using a variant of Dempster's rule of combination which takes into account the relative dependence between observations.


title: Iterative Markov Chain Monte Carlo Computation of Reference Priors and Minimax Risk

abstract: We present an iterative Markov chainMonte Carlo algorithm for computingreference priors and minimax risk forgeneral parametric families. Ourapproach uses MCMC techniques based onthe Blahut-Arimoto algorithm forcomputing channel capacity ininformation theory. We give astatistical analysis of the algorithm,bounding the number of samples requiredfor the stochastic algorithm to closelyapproximate the deterministic algorithmin each iteration. Simulations arepresented for several examples fromexponential families. Although we focuson applications to reference priors andminimax risk, the methods and analysiswe develop are applicable to a muchbroader class of optimization problemsand iterative algorithms.


title: Hypothesis Management in Situation-Specific Network Construction

abstract: This paper considers the problem of knowledge-based model construction in the presence of uncertainty about the association of domain entities to random variables. Multi-entity Bayesian networks (MEBNs) are defined as a representation for knowledge in domains characterized by uncertainty in the number of relevant entities, their interrelationships, and their association with observables. An MEBN implicitly specifies a probability distribution in terms of a hierarchically structured collection of Bayesian network fragments that together encode a joint probability distribution over arbitrarily many interrelated hypotheses. Although a finite query-complete model can always be constructed, association uncertainty typically makes exact model construction and evaluation intractable. The objective of hypothesis management is to balance tractability against accuracy. We describe an application to the problem of using intelligence reports to infer the organization and activities of groups of military vehicles. Our approach is compared to related work in the tracking and fusion literature.


title: Inference in Hybrid Networks: Theoretical Limits and Practical Algorithms

abstract: An important subclass of hybrid Bayesian networks are those that represent Conditional Linear Gaussian (CLG) distributions --- a distribution with a multivariate Gaussian component for each instantiation of the discrete variables. In this paper we explore the problem of inference in CLGs. We show that inference in CLGs can be significantly harder than inference in Bayes Nets. In particular, we prove that even if the CLG is restricted to an extremely simple structure of a polytree in which every continuous node has at most one discrete ancestor, the inference task is NP-hard.To deal with the often prohibitive computational cost of the exact inference algorithm for CLGs, we explore several approximate inference algorithms. These algorithms try to find a small subset of Gaussians which are a good approximation to the full mixture distribution. We consider two Monte Carlo approaches and a novel approach that enumerates mixture components in order of prior probability. We compare these methods on a variety of problems and show that our novel algorithm is very promising for large, hybrid diagnosis problems.


title: Exact Inference in Networks with Discrete Children of Continuous Parents

abstract: Many real life domains contain a mixture of discrete and continuous variables and can be modeled as hybrid Bayesian Networks. Animportant subclass of hybrid BNs are conditional linear Gaussian (CLG) networks, where the conditional distribution of the continuous variables given an assignment to the discrete variables is a multivariate Gaussian. Lauritzen's extension to the clique tree algorithm can be used for exact inference in CLG networks. However, many domains also include discrete variables that depend on continuous ones, and CLG networks do not allow such dependencies to berepresented. No exact inference algorithm has been proposed for these enhanced CLG networks. In this paper, we generalize Lauritzen's algorithm, providing the first "exact" inference algorithm for augmented CLG networks - networks where continuous nodes are conditional linear Gaussians but that also allow discrete children ofcontinuous parents. Our algorithm is exact in the sense that it computes the exact distributions over the discrete nodes, and the exact first and second moments of the continuous ones, up to the accuracy obtained by numerical integration used within thealgorithm. When the discrete children are modeled with softmax CPDs (as is the case in many real world domains) the approximation of the continuous distributions using the first two moments is particularly accurate. Our algorithm is simple to implement and often comparable in its complexity to Lauritzen's algorithm. We show empirically that it achieves substantially higher accuracy than previous approximate algorithms.


title: Probabilistic Logic Programming under Inheritance with Overriding

abstract: We present probabilistic logic programming under inheritance with overriding. This approach is based on new notions of entailment for reasoning with conditional constraints, which are obtained from the classical notion of logical entailment by adding the principle of inheritance with overriding. This is done by using recent approaches to probabilistic default reasoning with conditional constraints. We analyze the semantic properties of the new entailment relations. We also present algorithms for probabilistic logic programming under inheritance with overriding, and program transformations for an increased efficiency.


title: Solving Influence Diagrams using HUGIN, Shafer-Shenoy and Lazy Propagation

abstract: In this paper we compare three different architectures for the evaluation of influence diagrams: HUGIN, Shafer-Shenoy, and Lazy Evaluation architecture. The computational complexity of the architectures are compared on the LImited Memory Influence Diagram (LIMID): a diagram where only the requiste information for the computation of the optimal policies are depicted. Because the requsite information is explicitly represented in the LIMID the evaluation can take advantage of it, and significant savings in computational can be obtained. In this paper we show how the obtained savings is considerably increased when the computations performed on the LIMID is according to the Lazy Evaluation scheme.


title: A Bayesian Multiresolution Independence Test for Continuous Variables

abstract: In this paper we present a method ofcomputing the posterior probability ofconditional independence of two or morecontinuous variables from data,examined at several resolutions. Ourapproach is motivated by theobservation that the appearance ofcontinuous data varies widely atvarious resolutions, producing verydifferent independence estimatesbetween the variablesinvolved. Therefore, it is difficultto ascertain independence withoutexamining data at several carefullyselected resolutions. In our paper, weaccomplish this using the exactcomputation of the posteriorprobability of independence, calculatedanalytically given a resolution. Ateach examined resolution, we assume amultinomial distribution with Dirichletpriors for the discretized tableparameters, and compute the posteriorusing Bayesian integration. Acrossresolutions, we use a search procedureto approximate the Bayesian integral ofprobability over an exponential numberof possible histograms. Our methodgeneralizes to an arbitrary numbervariables in a straightforward manner.The test is suitable for Bayesiannetwork learning algorithms that useindependence tests to infer the networkstructure, in domains that contain anymix of continuous, ordinal andcategorical variables.


title: Aggregating Learned Probabilistic Beliefs

abstract: We consider the task of aggregating beliefs of severalexperts. We assume that these beliefs are represented as probabilitydistributions. We argue that the evaluation of any aggregationtechnique depends on the semantic context of this task. We propose aframework, in which we assume that nature generates samples from a`true' distribution and different experts form their beliefs based onthe subsets of the data they have a chance to observe. Naturally, theideal aggregate distribution would be the one learned from thecombined sample sets. Such a formulation leads to a natural way tomeasure the accuracy of the aggregation mechanism.We show that the well-known aggregation operator LinOP is ideallysuited for that task. We propose a LinOP-based learning algorithm,inspired by the techniques developed for Bayesian learning, whichaggregates the experts' distributions represented as Bayesiannetworks. Our preliminary experiments show that this algorithmperforms well in practice.


title: Expectation Propagation for approximate Bayesian inference

abstract: This paper presents a new deterministic approximation technique in Bayesian networks. This method, "Expectation Propagation", unifies two previous techniques: assumed-density filtering, an extension of the Kalman filter, and loopy belief propagation, an extension of belief propagation in Bayesian networks. All three algorithms try to recover an approximate distribution which is close in KL divergence to the true distribution. Loopy belief propagation, because it propagates exact belief states, is useful for a limited class of belief networks, such as those which are purely discrete. Expectation Propagation approximates the belief states by only retaining certain expectations, such as mean and variance, and iterates until these expectations are consistent throughout the network. This makes it applicable to hybrid networks with discrete and continuous nodes. Expectation Propagation also extends belief propagation in the opposite direction - it can propagate richer belief states that incorporate correlations between nodes. Experiments with Gaussian mixture models show Expectation Propagation to be convincingly better than methods with similar computational cost: Laplace's method, variational Bayes, and Monte Carlo. Expectation Propagation also provides an efficient algorithm for training Bayes point machine classifiers.


title: Recognition Networks for Approximate Inference in BN20 Networks

abstract: We propose using recognition networks for approximate inference inBayesian networks (BNs). A recognition network is a multilayerperception (MLP) trained to predict posterior marginals given observedevidence in a particular BN. The input to the MLP is a vector of thestates of the evidential nodes. The activity of an output unit isinterpreted as a prediction of the posterior marginal of thecorresponding variable. The MLP is trained using samples generated fromthe corresponding BN.We evaluate a recognition network that was trained to do inference ina large Bayesian network, similar in structure and complexity to theQuick Medical Reference, Decision Theoretic (QMR-DT). Our networkis a binary, two-layer, noisy-OR network containing over 4000 potentially observable nodes and over 600 unobservable, hidden nodes. Inreal medical diagnosis, most observables are unavailable, and there isa complex and unknown bias that selects which ones are provided. Weincorporate a very basic type of selection bias in our network: a knownpreference that available observables are positive rather than negative.Even this simple bias has a significant effect on the posterior. We compare the performance of our recognition network tostate-of-the-art approximate inference algorithms on a large set oftest cases. In order to evaluate the effect of our simplistic modelof the selection bias, we evaluate algorithms using a variety ofincorrectly modeled observation biases. Recognition networks performwell using both correct and incorrect observation biases.


title: The Factored Frontier Algorithm for Approximate Inference in DBNs

abstract: The Factored Frontier (FF) algorithm is a simple approximate inferencealgorithm for Dynamic Bayesian Networks (DBNs). It is very similar tothe fully factorized version of the Boyen-Koller (BK) algorithm, butinstead of doing an exact update at every step followed bymarginalisation (projection), it always works with factoreddistributions. Hence it can be applied to models for which the exactupdate step is intractable. We show that FF is equivalent to (oneiteration of) loopy belief propagation (LBP) on the original DBN, andthat BK is equivalent (to one iteration of) LBP on a DBN where wecluster some of the nodes. We then show empirically that byiterating, LBP can improve on the accuracy of both FF and BK. Wecompare these algorithms on two real-world DBNs: the first is a modelof a water treatment plant, and the second is a coupled HMM, used tomodel freeway traffic.


title: A Case Study in Knowledge Discovery and Elicitation in an Intelligent Tutoring Application

abstract: Most successful Bayesian network (BN) applications to datehave been built through knowledge elicitation from experts.This is difficult and time consuming, which has lead to recentinterest in automated methods for learning BNs from data. We present a case study in the construction of a BN in anintelligent tutoring application, specifically decimal misconceptions. Wedescribe the BN construction using expert elicitation and then investigate how certainexisting automated knowledge discovery methods might support the BN knowledge engineering process.


title: Lattice Particle Filters

abstract: A standard approach to approximate inference in state-space models isto apply a particle filter, e.g., the Condensation Algorithm.However, the performance of particle filters often varies significantlydue to their stochastic nature.We present a class of algorithms, called lattice particle filters, thatcircumvent this difficulty by placing the particles deterministicallyaccording to a Quasi-Monte Carlo integration rule.We describe a practical realization of this idea, discuss itstheoretical properties, and its efficiency.Experimental results with a synthetic 2D tracking problem show that thelattice particle filter is equivalent to a conventional particle filterthat has between 10 and 60% more particles, depending ontheir "sparsity" in the state-space.We also present results on inferring 3D human motion frommoving light displays.


title: Approximating MAP using Local Search

abstract: MAP is the problem of finding a most probable instantiation of a set of variables in a Bayesian network, given evidence. Unlike computing marginals, posteriors, and MPE (a special case of MAP), the time and space complexity of MAP is not only exponential in the network treewidth, but also in a larger parameter known as the "constrained" treewidth. In practice, this means that computing MAP can be orders of magnitude more expensive than computingposteriors or MPE. Thus, practitioners generally avoid MAP computations, resorting instead to approximating them by the most likely value for each MAP variableseparately, or by MPE.We present a method for approximating MAP using local search. This method has space complexity which is exponential onlyin the treewidth, as is the complexity of each search step. We investigate the effectiveness of different local searchmethods and several initialization strategies and compare them to otherapproximation schemes.Experimental results show that local search provides a much more accurate approximation of MAP, while requiring few search steps.Practically, this means that the complexity of local search is often exponential only in treewidth as opposed to the constrained treewidth, making approximating MAP as efficient as other computations.


title: Direct and Indirect Effects

abstract: The direct effect of one eventon another can be defined and measured byholding constant all intermediate variables between the two.Indirect effects present conceptual andpractical difficulties (in nonlinear models), because they cannot be isolated by holding certain variablesconstant. This paper shows a way of defining any path-specific effectthat does not invoke blocking the remainingpaths.This permits the assessment of a more naturaltype of direct and indirect effects, one thatis applicable in both linear and nonlinear models. The paper establishesconditions under which such assessments can be estimated consistentlyfrom experimental and nonexperimental data,and thus extends path-analytic techniques tononlinear and nonparametric models.