Uncovering Dependence Clusters and Linchpin Functions

Dependence clusters are (maximal) collections of mutually dependent source code entities according to some dependence relation. Their presence in software complicates many maintenance activities including testing, refactoring, and feature extraction. Despite several studies finding them common in production code, their formation, identification, and overall structure are not well understood, partly because of challenges in approximating true dependences between program entities. Previous research has considered two approximate dependence relations: a fine-grained statement-level relation using control and data dependences from a program’s System Dependence Graph and a coarser relation based on function-level controlflow reachability. In principal, the first is more expensive and more precise than the second. Using a collection of twenty programs, we present an empirical investigation of the clusters identified by these two approaches. In support of the analysis, we consider hybrid cluster types that works at the coarser function-level but is based on the higher-precision statement-level dependences. The three types of clusters are compared based on their slice sets using two clustering metrics. We also perform extensive analysis of the programs to identify linchpin functions – functions primarily responsible for holding a cluster together. Results include evidence that the less expensive, coarser approaches can often be used as e�ective proxies for the more expensive, finer-grained approaches. Finally, the linchpin analysis shows that linchpin functions can be e�ectively and automatically identified.

Electrical characterization of metal-oxide-polymer devices for non-volatile memory applications

The objective of this thesis is to study the properties of resistive switching effect based on bistable resistive memory which is fabricated in the form of Al2O3/polymer diodes and to contribute to the elucidation of resistive switching mechanisms. Resistive memories were characterized using a variety of electrical techniques, including current-voltage measurements, small-signal impedance, and electrical noise based techniques. All the measurements were carried out over a large temperature range. Fast voltage ramps were used to elucidate the dynamic response of the memory to rapid varying electric fields. The temperature dependence of the current provided insight into the role of trapped charges in resistive switching. The analysis of fast current fluctuations using electric noise techniques contributed to the elucidation of the kinetics involved in filament formation/rupture, the filament size and correspondent current capabilities. The results reported in this thesis provide insight into a number of issues namely: (i) The fundamental limitations on the speed of operation of a bi-layer resistive memory are the time and voltage dependences of the switch-on mechanism. (ii) The results explain the wide spread in switching times reported in the literature and the apparently anomalous behaviour of the high conductance state namely the disappearance of the negative differential resistance region at high voltage scan rates which is commonly attributed to a “dead time” phenomenon which had remained elusive since it was first reported in the ‘60s. (iii) Assuming that the current is filamentary, Comsol simulations were performed and used to explain the observed dynamic properties of the current-voltage characteristics. Furthermore, the simulations suggest that filaments can interact with each other. (iv) The current-voltage characteristics have been studied as a function of temperature. The findings indicate that creation and annihilation of filaments is controlled by filling and neutralizing traps localized at the oxide/polymer interface. (v) Resistive switching was also studied in small-molecule OLEDs. It was shown that the degradation that leads to a loss of light output during operation is caused by the presence of a resistive switching layer. A diagnostic tool that predicts premature failure of OLEDs was devised and proposed. Resistive switching is a property of oxides. These layers can grow in a number of devices including, organic light emitting diodes (OLEDs), spin-valve transistors and photovoltaic devices fabricated in different types of material. Under strong electric fields the oxides can undergo dielectric breakdown and become resistive switching layers. Resistive switching strongly modifies the charge injection causing a number of deleterious effects and eventually device failure. In this respect the findings in this thesis are relevant to understand reliability issues in devices across a very broad field.

Synthesis, characterization, electrochemical behavior and in vitro protein tyrosine kinase inhibitory activity of the cymene-halogenobenzohydroxamato [Ru(eta(6)-cymene)(bha)Cl] complexes

The ruthenium(II)-cymene complexes [Ru(eta(6)-cymene)(bha)Cl] with substituted halogenobenzohydroxamato (bha) ligands (substituents = 4-F, 4-Cl, 4-Br, 2,4-F-2, 3,4-F-2, 2,5-F-2, 2,6-F-2) have been synthesized and characterized by elemental analysis, IR, H-1 NMR, C-13 NMR, cyclic voltammetry and controlled-potential electrolysis, and density functional theory (DFT) studies. The compositions of their frontier molecular orbitals (MOs) were established by DFT calculations, and the oxidation and reduction potentials are shown to follow the orders of the estimated vertical ionization potential and electron affinity, respectively. The electrochemical E-L Lever parameter is estimated for the first time for the various bha ligands, which can thus be ordered according to their electron-donor character. All complexes exhibit very strong protein tyrosine kinase (PTK) inhibitory activity, even much higher than that of genistein, the clinically used PTK inhibitory drug. The complex containing the 2,4-difluorobenzohydroxamato ligand is the most active one, and the dependences of the PTK activity of the complexes and of their redox potentials on the ring substituents are discussed. (C) 2012 Elsevier B.V. All rights reserved.

Efficiency Comparison of DFT/IDFT Algorithms by Evaluating Diverse Hardware Implementations

In this paper we investigate various algorithms for performing Fast Fourier Transformation (FFT)/Inverse Fast Fourier Transformation (IFFT), and proper techniques for maximizing the FFT/IFFT execution speed, such as pipelining or parallel processing, and use of memory structures with pre-computed values (look up tables -LUT) or other dedicated hardware components (usually multipliers). Furthermore, we discuss the optimal hardware architectures that best apply to various FFT/IFFT algorithms, along with their abilities to exploit parallel processing with minimal data dependences of the FFT/IFFT calculations. An interesting approach that is also considered in this paper is the application of the integrated processing-in-memory Intelligent RAM (IRAM) chip to high speed FFT/IFFT computing. The results of the assessment study emphasize that the execution speed of the FFT/IFFT algorithms is tightly connected to the capabilities of the FFT/IFFT hardware to support the provided parallelism of the given algorithm. Therefore, we suggest that the basic Discrete Fourier Transform (DFT)/Inverse Discrete Fourier Transform (IDFT) can also provide high performances, by utilizing a specialized FFT/IFFT hardware architecture that can exploit the provided parallelism of the DFT/IDF operations. The proposed improvements include simplified multiplications over symbols given in polar coordinate system, using sinе and cosine look up tables, and an approach for performing parallel addition of N input symbols.

Efficiency Comparison of DFT/IDFT Algorithms by Evaluating Diverse Hardware : Implementations,Parallelization Prospectsand Possible Improvements

In this paper we investigate various algorithms for performing Fast Fourier Transformation (FFT)/Inverse Fast Fourier Transformation (IFFT), and proper techniquesfor maximizing the FFT/IFFT execution speed, such as pipelining or parallel processing, and use of memory structures with pre-computed values (look up tables -LUT) or other dedicated hardware components (usually multipliers). Furthermore, we discuss the optimal hardware architectures that best apply to various FFT/IFFT algorithms, along with their abilities to exploit parallel processing with minimal data dependences of the FFT/IFFT calculations. An interesting approach that is also considered in this paper is the application of the integrated processing-in-memory Intelligent RAM (IRAM) chip to high speed FFT/IFFT computing. The results of the assessment study emphasize that the execution speed of the FFT/IFFT algorithms is tightly connected to the capabilities of the FFT/IFFT hardware to support the provided parallelism of the given algorithm. Therefore, we suggest that the basic Discrete Fourier Transform (DFT)/Inverse Discrete Fourier Transform (IDFT) can also provide high performances, by utilizing a specialized FFT/IFFT hardware architecture that can exploit the provided parallelism of the DFT/IDF operations. The proposed improvements include simplified multiplications over symbols given in polar coordinate system, using sinе and cosine look up tables,and an approach for performing parallel addition of N input symbols.

Shot noise anomalies in elastic nondegenerate diffusive conductors

We present a theoretical investigation of shot-noise properties in nondegenerate elastic diffusive conductors. Both Monte Carlo simulations and analytical approaches are used. Two interesting phenomena are found: (i) the display of enhanced shot noise for given energy dependences of the scattering time, and (ii) the recovery of full shot noise for asymptotic high applied bias. The first phenomenon is associated with the onset of negative differential conductivity in energy space that drives the system towards a dynamical electrical instability in excellent agreement with analytical predictions. The enhancement is found to be strongly amplified when the dimensionality in momentum space is lowered from three to two dimensions. The second phenomenon is due to the suppression of the effects of long-range Coulomb correlations that takes place when the transit time becomes the shortest time scale in the system, and is common to both elastic and inelastic nondegenerate diffusive conductors. These phenomena shed different light in the understanding of the anomalous behavior of shot noise in mesoscopic conductors, which is a signature of correlations among different current pulses.

Entropy of a correlated system of nucleons

Realistic nucleon-nucleon interactions induce correlations to the nuclear many-body system, which lead to a fragmentation of the single-particle strength over a wide range of energies and momenta. We address the question of how this fragmentation affects the thermodynamical properties of nuclear matter. In particular, we show that the entropy can be computed with the help of a spectral function, which can be evaluated in terms of the self-energy obtained in the self-consistent Green's function approach. Results for the density and temperature dependences of the entropy per particle for symmetric nuclear matter are presented and compared to the results of lowest order finite-temperature Brueckner-Hartree-Fock calculations. The effects of correlations on the calculated entropy are small, if the appropriate quasiparticle approximation is used. The results demonstrate the thermodynamical consistency of the self-consistent T-matrix approximation for the evaluation of the Green's functions.

Modelo integrador para las rutas de señalización diferencial del receptor ionotrópico de glutamato activado por el N-metil-D-aspartato

El receptor ionotrópico de glutamato activado por N-metil-D-aspartato (iGluR-NMDA) es un complejo macromolecular heteromultimérico constituido por entre 3 y 5 subunidades de tres diferentes tipos, a saber: NR1, NR2A-D y NR3A y B. Se ha demostrado su participación activa en prácticamente todos los procesos fisiológicos, patológicos e intermediarios de efectos farmacológicos que ocurren en las células de tejidos excitables, inclusive se ha reportado su presencia en otros tejidos no excitables. En el sistema nervioso central (SNC) participa en los procesos de aprendizaje, memoria, plasticidad, diferenciación, migración de la célula neural y apoptosis. Además, en los eventos de índole farmacológica se ha demostrado su intervención en excitotoxicidad, drogadicción y alcoholismo. Surge entonces la pregunta de cómo un mismo complejo macromolecular puede participar en tantos y tan diversos procesos. La revisión de literatura en la que se demuestra la interacción del iGluR-NMDA con proteínas de señalización, soporte, adaptadoras, moduladoras, de adhesión celular, de citoesqueleto y enzimas reporta un conjunto de más de 160 moléculas que participan en las cascadas que generan las señales a diferentes niveles de interacción y con diferentes sustratos. En este artículo se presenta un modelo predictivo estructural y funcional que permite distinguir, por lo menos, tres rutas diferenciadas de señalización.

Numerical model of solid phase transformations governed by nucleation and growth: microstructure development during isothermal crystallization

A simple numerical model which calculates the kinetics of crystallization involving randomly distributed nucleation and isotropic growth is presented. The model can be applied to different thermal histories and no restrictions are imposed on the time and the temperature dependences of the nucleation and growth rates. We also develop an algorithm which evaluates the corresponding emerging grain-size distribution. The algorithm is easy to implement and particularly flexible, making it possible to simulate several experimental conditions. Its simplicity and minimal computer requirements allow high accuracy for two- and three-dimensional growth simulations. The algorithm is applied to explore the grain morphology development during isothermal treatments for several nucleation regimes. In particular, thermal nucleation, preexisting nuclei, and the combination of both nucleation mechanisms are analyzed. For the first two cases, the universal grain-size distribution is obtained. The high accuracy of the model is stated from its comparison to analytical predictions. Finally, the validity of the Kolmogorov-Johnson-Mehl-Avrami model SSSR, is verified for all the cases studied

Estudio de la aplicación de sistemas basados en el conocimiento a la operación de una planta de tratamiento de residuos sólidos urbanos por valorización energética

Una de las actuaciones posibles para la gestión de los residuos sólidos urbanos es la valorización energética, es decir la incineración con recuperación de energía. Sin embargo es muy importante controlar adecuadamente el proceso de incineración para evitar en lo posible la liberación de sustancias contaminantes a la atmósfera que puedan ocasionar problemas de contaminación industrial.Conseguir que tanto el proceso de incineración como el tratamiento de los gases se realice en condiciones óptimas presupone tener un buen conocimiento de las dependencias entre las variables de proceso. Se precisan métodos adecuados de medida de las variables más importantes y tratar los valores medidos con modelos adecuados para transformarlos en magnitudes de mando. Un modelo clásico para el control parece poco prometedor en este caso debido a la complejidad de los procesos, la falta de descripción cuantitativa y la necesidad de hacer los cálculos en tiempo real. Esto sólo se puede conseguir con la ayuda de las modernas técnicas de proceso de datos y métodos informáticos, tales como el empleo de técnicas de simulación, modelos matemáticos, sistemas basados en el conocimiento e interfases inteligentes. En [Ono, 1989] se describe un sistema de control basado en la lógica difusa aplicado al campo de la incineración de residuos urbanos. En el centro de investigación FZK de Karslruhe se están desarrollando aplicaciones que combinan la lógica difusa con las redes neuronales [Jaeschke, Keller, 1994] para el control de la planta piloto de incineración de residuos TAMARA. En esta tesis se plantea la aplicación de un método de adquisición de conocimiento para el control de sistemas complejos inspirado en el comportamiento humano. Cuando nos encontramos ante una situación desconocida al principio no sabemos como actuar, salvo por la extrapolación de experiencias anteriores que puedan ser útiles. Aplicando procedimientos de prueba y error, refuerzo de hipótesis, etc., vamos adquiriendo y refinando el conocimiento, y elaborando un modelo mental. Podemos diseñar un método análogo, que pueda ser implementado en un sistema informático, mediante el empleo de técnicas de Inteligencia Artificial.Así, en un proceso complejo muchas veces disponemos de un conjunto de datos del proceso que a priori no nos dan información suficientemente estructurada para que nos sea útil. Para la adquisición de conocimiento pasamos por una serie de etapas: - Hacemos una primera selección de cuales son las variables que nos interesa conocer. - Estado del sistema. En primer lugar podemos empezar por aplicar técnicas de clasificación (aprendizaje no supervisado) para agrupar los datos y obtener una representación del estado de la planta. Es posible establecer una clasificación, pero normalmente casi todos los datos están en una sola clase, que corresponde a la operación normal. Hecho esto y para refinar el conocimiento utilizamos métodos estadísticos clásicos para buscar correlaciones entre variables (análisis de componentes principales) y así poder simplificar y reducir la lista de variables. - Análisis de las señales. Para analizar y clasificar las señales (por ejemplo la temperatura del horno) es posible utilizar métodos capaces de describir mejor el comportamiento no lineal del sistema, como las redes neuronales. Otro paso más consiste en establecer relaciones causales entre las variables. Para ello nos sirven de ayuda los modelos analíticos - Como resultado final del proceso se pasa al diseño del sistema basado en el conocimiento. El objetivo principal es aplicar el método al caso concreto del control de una planta de tratamiento de residuos sólidos urbanos por valorización energética. En primer lugar, en el capítulo 2 Los residuos sólidos urbanos, se trata el problema global de la gestión de los residuos, dando una visión general de las diferentes alternativas existentes, y de la situación nacional e internacional en la actualidad. Se analiza con mayor detalle la problemática de la incineración de los residuos, poniendo especial interés en aquellas características de los residuos que tienen mayor importancia de cara al proceso de combustión.En el capítulo 3, Descripción del proceso, se hace una descripción general del proceso de incineración y de los distintos elementos de una planta incineradora: desde la recepción y almacenamiento de los residuos, pasando por los distintos tipos de hornos y las exigencias de los códigos de buena práctica de combustión, el sistema de aire de combustión y el sistema de humos. Se presentan también los distintos sistemas de depuración de los gases de combustión, y finalmente el sistema de evacuación de cenizas y escorias.El capítulo 4, La planta de tratamiento de residuos sólidos urbanos de Girona, describe los principales sistemas de la planta incineradora de Girona: la alimentación de residuos, el tipo de horno, el sistema de recuperación de energía, y el sistema de depuración de los gases de combustión Se describe también el sistema de control, la operación, los datos de funcionamiento de la planta, la instrumentación y las variables que son de interés para el control del proceso de combustión.En el capítulo 5, Técnicas utilizadas, se proporciona una visión global de los sistemas basados en el conocimiento y de los sistemas expertos. Se explican las diferentes técnicas utilizadas: redes neuronales, sistemas de clasificación, modelos cualitativos, y sistemas expertos, ilustradas con algunos ejemplos de aplicación.Con respecto a los sistemas basados en el conocimiento se analizan en primer lugar las condiciones para su aplicabilidad, y las formas de representación del conocimiento. A continuación se describen las distintas formas de razonamiento: redes neuronales, sistemas expertos y lógica difusa, y se realiza una comparación entre ellas. Se presenta una aplicación de las redes neuronales al análisis de series temporales de temperatura.Se trata también la problemática del análisis de los datos de operación mediante técnicas estadísticas y el empleo de técnicas de clasificación. Otro apartado está dedicado a los distintos tipos de modelos, incluyendo una discusión de los modelos cualitativos.Se describe el sistema de diseño asistido por ordenador para el diseño de sistemas de supervisión CASSD que se utiliza en esta tesis, y las herramientas de análisis para obtener información cualitativa del comportamiento del proceso: Abstractores y ALCMEN. Se incluye un ejemplo de aplicación de estas técnicas para hallar las relaciones entre la temperatura y las acciones del operador. Finalmente se analizan las principales características de los sistemas expertos en general, y del sistema experto CEES 2.0 que también forma parte del sistema CASSD que se ha utilizado.El capítulo 6, Resultados, muestra los resultados obtenidos mediante la aplicación de las diferentes técnicas, redes neuronales, clasificación, el desarrollo de la modelización del proceso de combustión, y la generación de reglas. Dentro del apartado de análisis de datos se emplea una red neuronal para la clasificación de una señal de temperatura. También se describe la utilización del método LINNEO+ para la clasificación de los estados de operación de la planta.En el apartado dedicado a la modelización se desarrolla un modelo de combustión que sirve de base para analizar el comportamiento del horno en régimen estacionario y dinámico. Se define un parámetro, la superficie de llama, relacionado con la extensión del fuego en la parrilla. Mediante un modelo linealizado se analiza la respuesta dinámica del proceso de incineración. Luego se pasa a la definición de relaciones cualitativas entre las variables que se utilizan en la elaboración de un modelo cualitativo. A continuación se desarrolla un nuevo modelo cualitativo, tomando como base el modelo dinámico analítico.Finalmente se aborda el desarrollo de la base de conocimiento del sistema experto, mediante la generación de reglas En el capítulo 7, Sistema de control de una planta incineradora, se analizan los objetivos de un sistema de control de una planta incineradora, su diseño e implementación. Se describen los objetivos básicos del sistema de control de la combustión, su configuración y la implementación en Matlab/Simulink utilizando las distintas herramientas que se han desarrollado en el capítulo anterior.Por último para mostrar como pueden aplicarse los distintos métodos desarrollados en esta tesis se construye un sistema experto para mantener constante la temperatura del horno actuando sobre la alimentación de residuos.Finalmente en el capítulo Conclusiones, se presentan las conclusiones y resultados de esta tesis.

Latitude and longitude vertical disparities

The literature on vertical disparity is complicated by the fact that several different definitions of the term “vertical disparity” are in common use, often without a clear statement about which is intended or a widespread appreciation of the properties of the different definitions. Here, we examine two definitions of retinal vertical disparity: elevation-latitude and elevation-longitude disparities. Near the fixation point, these definitions become equivalent, but in general, they have quite different dependences on object distance and binocular eye posture, which have not previously been spelt out. We present analytical approximations for each type of vertical disparity, valid for more general conditions than previous derivations in the literature: we do not restrict ourselves to objects near the fixation point or near the plane of regard, and we allow for non-zero torsion, cyclovergence, and vertical misalignments of the eyes. We use these expressions to derive estimates of the latitude and longitude vertical disparities expected at each point in the visual field, averaged over all natural viewing. Finally, we present analytical expressions showing how binocular eye position—gaze direction, convergence, torsion, cyclovergence, and vertical misalignment—can be derived from the vertical disparity field and its derivatives at the fovea.

A comparison of the reactivity of germylene and dimethylgermylene with some methylgermanes. Direct kinetic and quantum chemical studies

Time-resolved studies of germylene, GeH2, and dimethygermylene, GeMe2, generated by the 193 nm laser flash photolysis of appropriate precursor molecules have been carried out to try to obtain rate coefficients for their bimolecular reactions with dimethylgermane, Me2GeH2, in the gas-phase. GeH2 + Me2GeH2 was studied over the pressure range 1-100 Torr with SF6 as bath gas and at five temperatures in the range 296-553 K. Only slight pressure dependences were found (at 386, 447 and 553 K). RRKM modelling was carried out to fit these pressure dependences. The high pressure rate coefficients gave the Arrhenius parameters: log(A/cm(3) molecule(-1)s(-1)) = -10.99 +/- 0.07 and E-a = -(7.35 +/- 0.48) kJ mol(-1). No reaction could be found between GeMe2 + Me2GeH2 at any temperature up to 549 K, and upper limits of ca. 10(-14) cm(3) molecule(-1)s(-1) were set for the rate coefficients. A rate coefficient of (1.33 +/- 0.04) x 10(-11)cm(3) molecule(-1)s(-1) was also obtained for GeH2 + MeGeH3 at 296 K. No reaction was found between GeMe2 and MeGeH3. Rate coefficient comparisons showed, inter alia, that in the substrate germane Me-for-H substitution increased the magnitudes of rate coefficients significantly, while in the germylene Me-for-H substitution decreased the magnitudes of rate coefficients by at least four orders of magnitude. Quantum chemical calculations (G2(MP2,SVP)// B3LYP level) supported these findings and showed that the lack of reactivity of GeMe2 is caused by a positive energy barrier for rearrangement of the initially formed complexes. Full details of the structures of intermediate complexes and the discussion of their stabilities are given in the paper.

The addition reaction between silylene and ethyne: further isotope studies, pressure dependence studies, and quantum chemical calculations

Time-resolved kinetic studies of the reaction of dideutero-silylene, SiD2, generated by laser flash photolysis of phenylsilane-d(3), have been carried out to obtain rate constants for its bimolecular reaction with C2H2. The reaction was studied in the gas phase over the pressure range 1-100 Torr in SF6 bath gas, at five temperatures in the range 297-600 K. The second-order rate constants obtained by extrapolation to the high-pressure limits at each temperature fitted the Arrhenius equation log(k(infinity)/cm(3) molecule(-1) s(-1)) = (-10.05 +/- 0.05) + (3.43 +/- 0.36 kJ mol(-1))/RT ln 10. The rate constants were used to obtain a comprehensive set of isotope effects by comparison with earlier obtained rate constants for the reactions of SiH2 with C2H2 and C2D2. Additionally, pressure-dependent rate constants for the reaction of SiH2 with C2H2 in the presence of He (1-100 Tort) were obtained at 300, 399, and 613 K. Quantum chemical (ab initio) calculations of the SiC2H4 reaction system at the G3 level support the initial formation of silirene, which rapidly isomerizes to ethynylsilane as the major pathway. Reversible formation of vinylsilylene is also an important process. The calculations also indicate the involvement of several other intermediates, not previously suggested in the mechanism. RRKM calculations are in semiquantitative agreement with the pressure dependences and isotope effects suggested by the ab initio calculations, but residual discrepancies suggest the possible involvement of the minor reaction channel, SiH2 + C2H2 - SWPO + C2H4. The results are compared and contrasted with previous studies of this reaction system.

Aqueous gels of mixtures of ionic surfactant SDS with pluronic copolymers P123 or F127

Gel diagrams based on tube inversion and oscillatory rheometry are reported for Pluronic copolymers F127 (E98P67E98) and P123 (E21P67E21) in mixtures with anionic surfactant sodium dodecyl sulfate (SDS). Total concentrations (e, SDS+copolymer) were as high as 50 wt% with mole ratios SDS/copolymer (mr) in the ranges 1-5 (F127) a lid 1-7 (PI 23). Temperatures were its high as 90 degrees C. Determination of the temperature dependences of the dynamic moduli served to confirm the gel boundaries from tube inversion and to reveal the high elastic moduli of the gels, e.g., compared at corn parable positions in the gel phase, a 50 wt% SDS/P123 wit h mr = 7 had G' three times that of a corresponding gel of P123 alone. Sin all-angle X-ray scattering (SAX S) was used to show that the structures of all the SDS/F127 gels were bee and that the structures of the SDS/P123 gels with mr = I were either fcc(c = 30 wt%) or hex (c = 40 wt%). Assignment of structures to SDS/P123 gels with values of mr in the range 3-7 was more difficult, as high-order scattering peaks Could be very weak, and at the higher values of c and mr, the SAXS peaks included multiple reflections.

Mining gene networks with application to GAW15 problem 1

The genetic analysis workshop 15 (GAW15) problem 1 contained baseline expression levels of 8793 genes in immortalised B cells from 194 individuals in 14 Centre d’Etude du Polymorphisme Humane (CEPH) Utah pedigrees. Previous analysis of the data showed linkage and association and evidence of substantial individual variations. In particular, correlation was examined on expression levels of 31 genes and 25 target genes corresponding to two master regulatory regions. In this analysis, we apply Bayesian network analysis to gain further insight into these findings. We identify strong dependences and therefore provide additional insight into the underlying relationships between the genes involved. More generally, the approach is expected to be applicable for integrated analysis of genes on biological pathways.