Decision boundary for discrete Bayesian network classifiers

Bayesian network classifiers are a powerful machine learning tool. In order to evaluate the expressive power of these models, we compute families of polynomials that sign-represent decision functions induced by Bayesian network classifiers. We prove that those families are linear combinations of products of Lagrange basis polynomials. In absence of V -structures in the predictor sub-graph, we are also able to prove that this family of polynomials does indeed characterize the specific classifier considered. We then use this representation to bound the number of decision functions representable by Bayesian network classifiers with a given structure.

Discrete linear stability and sensitivity analysis of fluid flows

Este trabajo presenta un método discreto para el cálculo de estabilidad hidrodinámica y análisis de sensibilidad a perturbaciones externas para ecuaciones diferenciales y en particular para las ecuaciones de Navier-Stokes compressible. Se utiliza una aproximación con variable compleja para obtener una precisión analítica en la evaluación de la matriz Jacobiana. Además, mapas de sensibilidad para la sensibilidad a las modificaciones del flujo de base y a una fuerza constante permiten identificar las regiones del campo fluido donde una modificacin (ej. fuerza puntual) tiene un efecto estabilizador del flujo. Se presentan cuatro casos de prueba: (1) un caso analítico para comprobar la derivación discreta, (2) una cavidad cerrada a bajo Reynolds para mostrar la mayor precisión en el cálculo de los valores propios con la aproximación de paso complejo, (3) flujo 2D en un cilindro circular para validar la metodología, y (4) flujo en un cavidad abierta, presentado para validar el método en casos de inestabilidades convectivamente inestables. Los tres últimos casos mencionados (2-4) se resolvieron con las ecuaciones de Navier-Stokes compresibles, utilizando un método Discontinuous Galerkin Spectral Element Method. Se obtuvo una buena concordancia para el caso de validación (3), cuando se comparó el nuevo método con resultados de la literatura. Además, este trabajo muestra que para el cálculo de los modos propios directos y adjuntos, así como para los mapas de sensibilidad, el uso de variables complejas es de suprema importancia para obtener una predicción precisa. El método descrito es aplicado al análisis para la estabilización de la estela generada por un disco actuador, que representa un modelo sencillo para hélices, rotores de helicópteros o turbinas eólicas. Se explora la primera bifurcación del flujo para un disco actuador, y se sugiere que está asociada a una inestabilidad de tipo Kelvin-Helmholtz, cuya estabilidad se controla con en el número de Reynolds y en la resistencia del disco actuador (o fuerza resistente). En primer lugar, se verifica que la disminución de la resistencia del disco tiene un efecto estabilizador parecido a una disminución del Reynolds. En segundo lugar, el análisis hidrodinmico discreto identifica dos regiones para la colocación de una fuerza puntual que controle las inestabilidades, una cerca del disco y otra en una zona aguas abajo. En tercer lugar, se muestra que la inclusión de un forzamiento localizado cerca del actuador produce una estabilización más eficiente que al forzar aguas abajo. El análisis de los campos de flujo controlados confirma que modificando el gradiente de velocidad cerca del actuador es más eficiente para estabilizar la estela. Estos resultados podrían proporcionar nuevas directrices para la estabilización de la estela de turbinas de viento o de marea cuando estén instaladas en un parque eólico y minimizar las interacciones no estacionarias entre turbinas. ABSTRACT A discrete framework for computing the global stability and sensitivity analysis to external perturbations for any set of partial differential equations is presented. In particular, a complex-step approximation is used to achieve near analytical accuracy for the evaluation of the Jacobian matrix. Sensitivity maps for the sensitivity to base flow modifications and to a steady force are computed to identify regions of the flow field where an input could have a stabilising effect. Four test cases are presented: (1) an analytical test case to prove the theory of the discrete framework, (2) a lid-driven cavity at low Reynolds case to show the improved accuracy in the calculation of the eigenvalues when using the complex-step approximation, (3) the 2D flow past a circular cylinder at just below the critical Reynolds number is used to validate the methodology, and finally, (4) the flow past an open cavity is presented to give an example of the discrete method applied to a convectively unstable case. The latter three (2–4) of the aforementioned cases were solved with the 2D compressible Navier–Stokes equations using a Discontinuous Galerkin Spectral Element Method. Good agreement was obtained for the validation test case, (3), with appropriate results in the literature. Furthermore, it is shown that for the calculation of the direct and adjoint eigenmodes and their sensitivity maps to external perturbations, the use of complex variables is paramount for obtaining an accurate prediction. An analysis for stabilising the wake past an actuator disc, which represents a simple model for propellers, helicopter rotors or wind turbines is also presented. We explore the first flow bifurcation for an actuator disc and it suggests that it is associated to a Kelvin- Helmholtz type instability whose stability relies on the Reynolds number and the flow resistance applied through the disc (or actuator forcing). First, we report that decreasing the disc resistance has a similar stabilising effect to an decrease in the Reynolds number. Second, a discrete sensitivity analysis identifies two regions for suitable placement of flow control forcing, one close to the disc and one far downstream where the instability originates. Third, we show that adding a localised forcing close to the actuator provides more stabilisation that forcing far downstream. The analysis of the controlled flow fields, confirms that modifying the velocity gradient close to the actuator is more efficient to stabilise the wake than controlling the sheared flow far downstream. An interesting application of these results is to provide guidelines for stabilising the wake of wind or tidal turbines when placed in an energy farm to minimise unsteady interactions.

On the Velocity and Acceleration Estimation from Discrete Time-Position Sensors

In this work is addressed the topic of estimation of velocity and acceleration from digital position data. It is presented a review of several classic methods and implemented with real position data from a low cost digital sensor of a hydraulic linear actuator. The results are analyzed and compared. It is shown that static methods have a limited bandwidth application, and that the performance of some methods may be enhanced by adapting its parameters according to the current state.

Estimation of Symmetric Channels for Discrete Cosine Transform Type-I Multicarrier Systems: A Compressed Sensing Approach

The problem of channel estimation for multicarrier communications is addressed. We focus on systems employing the Discrete Cosine Transform Type-I (DCT1) even at both the transmitter and the receiver, presenting an algorithm which achieves an accurate estimation of symmetric channel filters using only a small number of training symbols. The solution is obtained by using either matrix inversion or compressed sensing algorithms. We provide the theoretical results which guarantee the validity of the proposed technique for the DCT1. Numerical simulations illustrate the good behaviour of the proposed algorithm.

Economic considerations and patients’ preferences affect treatment selection for rheumatoid arthritis patients : A discrete choice experiment among European rheumatologists

Peer reviewed

Voltage dependence of the pattern and frequency of discrete Ca2+ release events after brief repriming in frog skeletal muscle

Applying a brief repolarizing pre-pulse to a depolarized frog skeletal muscle fiber restores a small fraction of the transverse tubule membrane voltage sensors from the inactivated state. During a subsequent depolarizing test pulse we detected brief, highly localized elevations of myoplasmic Ca2+ concentration (Ca2+ “sparks”) initiated by restored voltage sensors in individual triads at all test pulse voltages. The latency histogram of these events gives the gating pattern of the sarcoplasmic reticulum (SR) calcium release channels controlled by the restored voltage sensors. Both event frequency and clustering of events near the start of the test pulse increase with test pulse depolarization. The macroscopic SR calcium release waveform, obtained from the spark latency histogram and the estimated open time of the channel or channels underlying a spark, exhibits an early peak and rapid marked decline during large depolarizations. For smaller depolarizations, the release waveform exhibits a smaller peak and a slower decline. However, the mean use time and mean amplitude of the individual sparks are quite similar at all test depolarizations and at all times during a given depolarization, indicating that the channel open times and conductances underlying sparks are essentially independent of voltage. Thus, the voltage dependence of SR Ca2+ release is due to changes in the frequency and pattern of occurrence of individual, voltage-independent, discrete release events.

Cell cycle-dependent colocalization of BARD1 and BRCA1 proteins in discrete nuclear domains

Germ-line mutations of the BRCA1 gene predispose women to early-onset breast and ovarian cancer by compromising the gene’s presumptive function as a tumor suppressor. Although the biochemical properties of BRCA1 polypeptides are not understood, their expression pattern and subcellular localization suggest a role in cell-cycle regulation. When resting cells are induced to proliferate, the steady-state levels of BRCA1 increase in late G1 and reach a maximum during S phase. Moreover, in S phase cells, BRCA1 polypeptides are hyperphosphorylated and accumulate into discrete subnuclear foci termed “BRCA1 nuclear dots.” BRCA1 associates in vivo with a structurally related protein termed BARD1. Here we show that the steady-state levels of BARD1, unlike those of BRCA1, remain relatively constant during cell cycle progression. However, immunostaining revealed that BARD1 resides within BRCA1 nuclear dots during S phase of the cell cycle, but not during the G1 phase. Nevertheless, BARD1 polypeptides are found exclusively in the nuclear fractions of both G1- and S-phase cells. Therefore, progression to S phase is accompanied by the aggregation of nuclear BARD1 polypeptides into BRCA1 nuclear dots. This cell cycle-dependent colocalization of BARD1 and BRCA1 indicates a role for BARD1 in BRCA1-mediated tumor suppression.

A Discrete Stage of Baculovirus GP64-mediated Membrane Fusion

Viral fusion protein trimers can play a critical role in limiting lipids in membrane fusion. Because the trimeric oligomer of many viral fusion proteins is often stabilized by hydrophobic 4-3 heptad repeats, higher-order oligomers might be stabilized by similar sequences. There is a hydrophobic 4-3 heptad repeat contiguous to a putative oligomerization domain of Autographa californica multicapsid nucleopolyhedrovirus envelope glycoprotein GP64. We performed mutagenesis and peptide inhibition studies to determine if this sequence might play a role in catalysis of membrane fusion. First, leucine-to-alanine mutants within and flanking the amino terminus of the hydrophobic 4-3 heptad repeat motif that oligomerize into trimers and traffic to insect Sf9 cell surfaces were identified. These mutants retained their wild-type conformation at neutral pH and changed conformation in acidic conditions, as judged by the reactivity of a conformationally sensitive mAb. These mutants, however, were defective for membrane fusion. Second, a peptide encoding the portion flanking the GP64 hydrophobic 4-3 heptad repeat was synthesized. Adding peptide led to inhibition of membrane fusion, which occurred only when the peptide was present during low pH application. The presence of peptide during low pH application did not prevent low pH–induced conformational changes, as determined by the loss of a conformationally sensitive epitope. In control experiments, a peptide of identical composition but different sequence, or a peptide encoding a portion of the Ebola GP heptad motif, had no effect on GP64-mediated fusion. Furthermore, when the hemagglutinin (X31 strain) fusion protein of influenza was functionally expressed in Sf9 cells, no effect on hemagglutinin-mediated fusion was observed, suggesting that the peptide does not exert nonspecific effects on other fusion proteins or cell membranes. Collectively, these studies suggest that the specific peptide sequences of GP64 that are adjacent to and include portions of the hydrophobic 4-3 heptad repeat play a dynamic role in membrane fusion at a stage that is downstream of the initiation of protein conformational changes but upstream of lipid mixing.

The importance of being discrete: Life always wins on the surface

Many systems in chemistry, biology, finance, and social sciences present emerging features that are not easy to guess from the elementary interactions of their microscopic individual components. In the past, the macroscopic behavior of such systems was modeled by assuming that the collective dynamics of microscopic components can be effectively described collectively by equations acting on spatially continuous density distributions. It turns out that, to the contrary, taking into account the actual individual/discrete character of the microscopic components of these systems is crucial for explaining their macroscopic behavior. In fact, we find that in conditions in which the continuum approach would predict the extinction of all of the population (respectively the vanishing of the invested capital or the concentration of a chemical substance, etc.), the microscopic granularity insures the emergence of macroscopic localized subpopulations with collective adaptive properties that allow their survival and development. In particular it is found that in two dimensions “life” (the localized proliferating phase) always prevails.

Identification of Discrete Classes of Endosome-derived Small Vesicles as a Major Cellular Pool for Recycling Membrane Proteins

Vesicles carrying recycling plasma membrane proteins from early endosomes have not yet been characterized. Using Chinese hamster ovary cells transfected with the facilitative glucose transporter, GLUT4, we identified two classes of discrete, yet similarly sized, small vesicles that are derived from early endosomes. We refer to these postendosomal vesicles as endocytic small vesicles or ESVs. One class of ESVs contains a sizable fraction of the pool of the transferrin receptor, and the other contains 40% of the total cellular pool of GLUT4 and is enriched in the insulin-responsive aminopeptidase (IRAP). The ESVs contain cellubrevin and Rab4 but are lacking other early endosomal markers, such as EEA1 or syntaxin13. The ATP-, temperature-, and cytosol-dependent formation of ESVs has been reconstituted in vitro from endosomal membranes. Guanosine 5′-[γ-thio]triphosphate and neomycin, but not brefeldin A, inhibit budding of the ESVs in vitro. A monoclonal antibody recognizing the GLUT4 cytoplasmic tail perturbs the in vitro targeting of GLUT4 to the ESVs without interfering with the incorporation of IRAP or TfR. We suggest that cytosolic proteins mediate the incorporation of recycling membrane proteins into discrete populations of ESVs that serve as carrier vesicles to store and then transport the cargo from early endosomes, either directly or indirectly, to the cell surface.

Specific Soybean Lipoxygenases Localize to Discrete Subcellular Compartments and Their mRNAs Are Differentially Regulated by Source-Sink Status1

Members of the lipoxygenase multigene family, found widely in eukaryotes, have been proposed to function in nitrogen partitioning and storage in plants. Lipoxygenase gene responses to source-sink manipulations in mature soybean (Glycine max [L.] Merr.) leaves were examined using gene-specific riboprobes to the five vegetative lipoxygenases (vlxA–vlxE). Steady-state levels of all vlx mRNAs responded strongly to sink limitation, but specific transcripts exhibited differential patterns of response as well. During reproductive sink limitation, vlxA and vlxB messages accumulated to high levels, whereas vlxC and vlxD transcript levels were modest. Immunolocalization using peptide-specific antibodies demonstrated that under control conditions, VLXB was present in the cytosol of the paraveinal mesophyll and with pod removal accumulated additionally in the bundle-sheath and adjacent cells. With sink limitation VLXD accumulated to apparent high levels in the vacuoles of the same cells. Segregation of gene products at the cellular and subcellular levels may thus permit complex patterns of differential regulation within the same cell type. Specific lipoxygenase isoforms may have a role in short-term nitrogen storage (VLXC/D), whereas others may simultaneously function in assimilate partitioning as active enzymes (VLXA/B).

Huntingtin-associated protein (HAP1): discrete neuronal localizations in the brain resemble those of neuronal nitric oxide synthase.

Huntington disease stems from a mutation of the protein huntingtin and is characterized by selective loss of discrete neuronal populations in the brain. Despite a massive loss of neurons in the corpus striatum, NO-generating neurons are intact. We recently identified a brain-specific protein that associates with huntingtin and is designated huntingtin-associated protein (HAP1). We now describe selective neuronal localizations of HAP1. In situ hybridization studies reveal a resemblance of HAP1 and neuronal nitric oxide synthase (nNOS) mRNA localizations with dramatic enrichment of both in the pedunculopontine nuclei, the accessory olfactory bulb, and the supraoptic nucleus of the hypothalamus. Both nNOS and HAP1 are enriched in subcellular fractions containing synaptic vesicles. Immunocytochemical studies indicate colocalizations of HAP1 and nNOS in some neurons. The possible relationship of HAP1 and nNOS in the brain is reminiscent of the relationship of dystrophin and nNOS in skeletal muscle and suggests a role of NO in Huntington disease, analogous to its postulated role in Duchenne muscular dystrophy.

Rock friction and its implications for earthquake prediction examined via models of Parkfield earthquakes.

The friction of rocks in the laboratory is a function of time, velocity of sliding, and displacement. Although the processes responsible for these dependencies are unknown, constitutive equations have been developed that do a reasonable job of describing the laboratory behavior. These constitutive laws have been used to create a model of earthquakes at Parkfield, CA, by using boundary conditions appropriate for the section of the fault that slips in magnitude 6 earthquakes every 20-30 years. The behavior of this model prior to the earthquakes is investigated to determine whether or not the model earthquakes could be predicted in the real world by using realistic instruments and instrument locations. Premonitory slip does occur in the model, but it is relatively restricted in time and space and detecting it from the surface may be difficult. The magnitude of the strain rate at the earth's surface due to this accelerating slip seems lower than the detectability limit of instruments in the presence of earth noise. Although not specifically modeled, microseismicity related to the accelerating creep and to creep events in the model should be detectable. In fact the logarithm of the moment rate on the hypocentral cell of the fault due to slip increases linearly with minus the logarithm of the time to the earthquake. This could conceivably be used to determine when the earthquake was going to occur. An unresolved question is whether this pattern of accelerating slip could be recognized from the microseismicity, given the discrete nature of seismic events. Nevertheless, the model results suggest that the most likely solution to earthquake prediction is to look for a pattern of acceleration in microseismicity and thereby identify the microearthquakes as foreshocks.

Logarithm of odds (lods) for linkage in complex inheritance.

Lod scores provide a method to unify linkage tests based on identity by descent and identity in marker state while permitting selection of the most informative individuals through their disease-related phenotypes and markers in relatives. After parametric lods are reviewed, a nonparametric approach that depends on a single logistic parameter beta is introduced. Lods for parents tested or unknown are derived, multiple pairwise mapping is presented, and power is shown to be good even for moderately small values of beta. Comparison of parametric and nonparametric approaches (yet to be made) will provide for polygenes the efficiency and reliability that lod scores gave to mapping of major loci 40 years ago.

SeismicWaveTool: Continuous and discrete wavelet analysis and filtering for multichannel seismic data

A MATLAB-based computer code has been developed for the simultaneous wavelet analysis and filtering of multichannel seismic data. The considered time–frequency transforms include the continuous wavelet transform, the discrete wavelet transform and the discrete wavelet packet transform. The developed approaches provide a fast and precise time–frequency examination of the seismograms at different frequency bands. Moreover, filtering methods for noise, transients or even baseline removal, are implemented. The primary motivation is to support seismologists with a user-friendly and fast program for the wavelet analysis, providing practical and understandable results.