STATISTICAL MECHANICS OF THE INVERSE ISING MODEL

Similar documents
An Ising model on 2-D image

Constructing Free Energy Approximations and GBP Algorithms

Advanced Machine Learning. An Ising model on 2-D image

Probabilistic Graphical Models

LOCAL FRACTIONAL LAPLACE SERIES EXPANSION METHOD FOR DIFFUSION EQUATION ARISING IN FRACTAL HEAT TRANSFER

Geometric Correction or Georeferencing

Partially Observable Systems. 1 Partially Observable Markov Decision Process (POMDP) Formalism

Lecture 36. Finite Element Methods

Rank One Update And the Google Matrix by Al Bernstein Signal Science, LLC

Fall 2012 Analysis of Experimental Measurements B. Eisenstein/rev. S. Errede. with respect to λ. 1. χ λ χ λ ( ) λ, and thus:

6 Roots of Equations: Open Methods

Lecture 3 Camera Models 2 & Camera Calibration. Professor Silvio Savarese Computational Vision and Geometry Lab

Least squares. Václav Hlaváč. Czech Technical University in Prague

Chapter Newton-Raphson Method of Solving a Nonlinear Equation

Why BP Works STAT 232B

Lecture 4: Piecewise Cubic Interpolation

Katholieke Universiteit Leuven Department of Computer Science

CISE 301: Numerical Methods Lecture 5, Topic 4 Least Squares, Curve Fitting

Solubilities and Thermodynamic Properties of SO 2 in Ionic

Chapter Newton-Raphson Method of Solving a Nonlinear Equation

CIS587 - Artificial Intelligence. Uncertainty CIS587 - AI. KB for medical diagnosis. Example.

6. Chemical Potential and the Grand Partition Function

Machine Learning. Support Vector Machines. Le Song. CSE6740/CS7641/ISYE6740, Fall Lecture 8, Sept. 13, 2012 Based on slides from Eric Xing, CMU

ICS 252 Introduction to Computer Design

DCDM BUSINESS SCHOOL NUMERICAL METHODS (COS 233-8) Solutions to Assignment 3. x f(x)

Definition of Tracking

Multiple view geometry

Chemistry 163B Absolute Entropies and Entropy of Mixing

Chapter 2 Introduction to Algebra. Dr. Chih-Peng Li ( 李 )

ITERATIVE METHODS FOR SOLVING SYSTEMS OF LINEAR ALGEBRAIC EQUATIONS

UNIVERSITY OF IOANNINA DEPARTMENT OF ECONOMICS. M.Sc. in Economics MICROECONOMIC THEORY I. Problem Set II

ORDINARY DIFFERENTIAL EQUATIONS

HAMILTON-JACOBI TREATMENT OF LAGRANGIAN WITH FERMIONIC AND SCALAR FIELD

CS Lecture 13. More Maximum Likelihood

Introduction to Numerical Integration Part II

Lecture 8: Camera Calibration

Multipoint Analysis for Sibling Pairs. Biostatistics 666 Lecture 18

CHI-SQUARE DIVERGENCE AND MINIMIZATION PROBLEM

Electrochemical Thermodynamics. Interfaces and Energy Conversion

ψ ij has the eigenvalue

Least Squares Fitting of Data

An Introduction to Support Vector Machines

3/6/00. Reading Assignments. Outline. Hidden Markov Models: Explanation and Model Learning

Our focus will be on linear systems. A system is linear if it obeys the principle of superposition and homogenity, i.e.

Chapter 5 Supplemental Text Material R S T. ij i j ij ijk

Decomposition of Boolean Function Sets for Boolean Neural Networks

Fitting a Polynomial to Heat Capacity as a Function of Temperature for Ag. Mathematical Background Document

Department of Mechanical Engineering, University of Bath. Mathematics ME Problem sheet 11 Least Squares Fitting of data

For now, let us focus on a specific model of neurons. These are simplified from reality but can achieve remarkable results.

1 Input-Output Mappings. 2 Hebbian Failure. 3 Delta Rule Success.

18.7 Artificial Neural Networks

ME 501A Seminar in Engineering Analysis Page 1

Dennis Bricker, 2001 Dept of Industrial Engineering The University of Iowa. MDP: Taxi page 1

90 S.S. Drgomr nd (t b)du(t) =u()(b ) u(t)dt: If we dd the bove two equltes, we get (.) u()(b ) u(t)dt = p(; t)du(t) where p(; t) := for ll ; t [; b]:

The linear system. The problem: solve

LECTURE :FACTOR ANALYSIS

5.03, Inorganic Chemistry Prof. Daniel G. Nocera Lecture 2 May 11: Ligand Field Theory

THE ARIMOTO-BLAHUT ALGORITHM FOR COMPUTATION OF CHANNEL CAPACITY. William A. Pearlman. References: S. Arimoto - IEEE Trans. Inform. Thy., Jan.

PHYS 215C: Quantum Mechanics (Spring 2017) Problem Set 3 Solutions

Prof. Dr. I. Nasser Phys 630, T Aug-15 One_dimensional_Ising_Model

Principle Component Analysis

Mechanics Physics 151

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

FUNDAMENTALS ON ALGEBRA MATRICES AND DETERMINANTS

Exploiting Structure in Probability Distributions Irit Gat-Viks

The Number of Rows which Equal Certain Row

Mean Field / Variational Approximations

Introduction to Density Functional Theory. Jeremie Zaffran 2 nd year-msc. (Nanochemistry)

Jens Siebel (University of Applied Sciences Kaiserslautern) An Interactive Introduction to Complex Numbers

Probability-Theoretic Junction Trees

ESCI 342 Atmospheric Dynamics I Lesson 1 Vectors and Vector Calculus

1 Definition of Rademacher Complexity

Meshless Surfaces. presented by Niloy J. Mitra. An Nguyen

Beam based calibration for beam position monitor

Bayesian belief networks

Outline. Review Quadrilateral Equation. Review Linear φ i Quadrilateral. Review x and y Derivatives. Review φ Derivatives

Machine Learning Support Vector Machines SVM

Perceptual Organization (IV)

ÉCOLE POLYTECHNIQUE FÉDÉRALE DE LAUSANNE

Least Squares Fitting of Data

Solution of Tutorial 5 Drive dynamics & control

International Journal of Pure and Applied Sciences and Technology

Applied Mathematics Letters

Approximate Inference: Mean Field Methods

Logistic Regression Maximum Likelihood Estimation

BEAM BASED CALIBRATION FOR BEAM POSITION MONITORS

MACHINE APPLIED MACHINE LEARNING LEARNING. Gaussian Mixture Regression

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION

II The Z Transform. Topics to be covered. 1. Introduction. 2. The Z transform. 3. Z transforms of elementary functions

VECTORS VECTORS VECTORS VECTORS. 2. Vector Representation. 1. Definition. 3. Types of Vectors. 5. Vector Operations I. 4. Equal and Opposite Vectors

Quiz: Experimental Physics Lab-I

Resonant Slepton Production

INTERPOLATION(1) ELM1222 Numerical Analysis. ELM1222 Numerical Analysis Dr Muharrem Mercimek

Many-Body Calculations of the Isotope Shift

The Schur-Cohn Algorithm

Lecture 10: Euler s Equations for Multivariable

10) Activity analysis

Study of Trapezoidal Fuzzy Linear System of Equations S. M. Bargir 1, *, M. S. Bapat 2, J. D. Yadav 3 1

CHALMERS, GÖTEBORGS UNIVERSITET. SOLUTIONS to RE-EXAM for ARTIFICIAL NEURAL NETWORKS. COURSE CODES: FFR 135, FIM 720 GU, PhD

Transcription:

STATISTICAL MECHANICS OF THE INVESE ISING MODEL Muro Cro Supervsors: rof. Mchele Cselle rof. ccrdo Zecchn uly 2009

INTODUCTION SUMMAY OF THE ESENTATION Defnton of the drect nd nverse prole Approton ethods of the drect prole: vrtonl pproches Bethe Overvew of drect nd nverse lgorths Sultons GOALS OF THE THESIS Use of lgorths whch generlze Bethe pproton Gp n order to solve the nverse prole Upgrde of the Gp lgorth to get ore ccurte clculton of the correltons nd pplcton to the nverse prole

ISING MODEL We cn represent the usul Boltznn prolty dstruton wth fctor grph: > < h H e e e } { } { β β β > < h H } {

ISING MODEL H { } < > h h > < > < Inverse role!

EXAMLES OF ALICATIONS Neuron networks reconstructon [E.Schnedn M..Berry.Segev W.Blek 2006] Genes networks reconstructon [A.Brusten A.gn M.Wegt.Zecchn 2008] roten networks recontructons [G.Tkck 2007] Why Isng odel? Mu entropy prncple

VAIATIONAL AOACHES Men feld pprotons G[ ] < E > T S G[ ] ~ F Gs Free Energy Boltznn Dstr. Heloltz Free Energy Men feld pproton Mnzton of G In sudon We chose for for nd we nze the functonl G

Locl consstency: BELIEFS: 1 1 VAIATIONAL AOACHES Bethe pproton [ ] > < q p 1 } { [ ] [ ] log 1 log E q E G > < Mnu equtons Constrnts equtons 1 Couplngs f 2 Couplngs f Messges equtons t the f pont

B Drect Isng B equtons t t \ 1 η ν... k k t k t \ \ ν ψ η..k..

B Drect Isng * η t t s s \ * \ * η η ψ

GB [.S.Yedd W.T. Freen Y.Wess 2001] Drect Isng Let s defne: c c U E S Log : vrles close to functon node c 1 c U U A Generlzton of regons {} : ng[{}] Bethe pproton New for of G GB equtons

GB Drect Isng \ \ D I I N I I F ψ

GB Drect Isng \ \ D I I N I I F ψ Ψ \ ' ' D D E D D D A

GB [.S.Yedd W.T. Freen Y.Wess 2001] Drect Isng Is G vld? depends on regons c Condton: c I c I 1 I 1 f We wnt to count every node only once

INVESE ISING Solvng the nverse prole through n tertve ethod Coplete Grph χ Itertve process Updte rules for h Grph wth: < σ > < σσ > χ

B/GB Inverse Isng Self-consstent equtons n the essges Self consstent equtons n the nputs χ f M h χ g M h Messges fed h f 1 M χ g M 1 χ

B/GB Inverse Isng χ Otnng eperentl vlues of Intlze coplete grph wth rndo h For t 1T Itertons Bp/Gp sve essges Functon nodes updte n the grph h f M 1 g M 1 χ χ End

B Inverse Isng EXAMLE ~ h h Tnh h ATnh h ~ ν ν ν ν ~ ~ ~ ~ h h h h h h h h e e e e

B Inverse Isng EXAMLE χ e ν ν 2 1 ν 2 1 ν 1 Tnh Tnh χ ATnh χ χ 1

Sus.rop. [M.Mezrd T.Mor 2008] Drect nd Inverse Isng B lttons χ only f < > Fluctuton-esponse theore < σ σ > < σ >< σ > h h h wth η We hve to know dervtves of essges! B equtons ν η B dervtves ν η ν η χ f h Invertng f New updte rule for

GS Drect Isng A: prove ccurcy n χ n the GB schee B Sus.rop. Flutt.- esp. T. GB GS Fnd equtons for dervtves of essges n GB Gsp equtons Do the dervtves of the 1-elefs equtons Correltons

GS Isng Inverso Etrctng couplngs fro GS wth rtrry regons? Contrnt: t lest ll the 1 nd 2 regons 1-regons see Bethe ppro. rtl reuse of Sus.rop. forls

Sll correltons epnson [.Monsson V.Sessk 2008] Lkelhood zton Let s suppose to hve n copes of the syste Let s ze the prolty to hve these copes gven the theory: 1 [ Log { σ} n h ] n n where : S Log Z h ns h c Let s solve the prole for c 0 c β c Anltcl epressons for e h

Sultons eples 2 Coplete 2 2 N N 1 < Grph N 20 Bp Gp Nve M.F. Sus.rop. Gsp S.C.E.

Sultons eples 2 2 2 N N 1 < Tree N 20 Bp Gp Nve M.F. Sus.rop. Gsp S.C.E.

Conclusons nd future B Sus.rop. GB GS OBLEMS Bg correltons Glssy phse Loopy B Sll Correlton Epnson Necessty to redefne the prole Iterton ethods n g correltons rege Cvty ethods for GB n 1SB [T.zzo et l. 2009] New pprotons for the free energy

Acronys & Notton B GB Sus.rop. GS S.C.E. Nve M.F. GB 3 GS 3 Belef ropgton Generlzed Belef ropgton Susceptlty ropgton Generlzed Susceptlty ropgton Sll Correlton Epnson Men Feld Theory Fluctuton esponse Theore GB wth ll regons wth 1 or 2 or 3 spns GS wth ll regons wth 1 or 2 or 3 spns \ Identty true ecept for norlzton fctor Set of nodes whch re connected to wth lnk Set wthout the node node representng vrle leled y k node representng functon leled y c