In this appendix, we evaluate the derivative of Eq. 9 in the main text, i.e., we need to calculate

Similar documents
MAGIC058 & MATH64062: Partial Differential Equations 1

arxiv: v1 [math.co] 5 Jun 2015

Statistical Physics. Solutions Sheet 5.

Chapter 6 Notes, Larson/Hostetler 3e

1. The vibrating string problem revisited.

Complete Description of the Thelen2003Muscle Model

( dg. ) 2 dt. + dt. dt j + dh. + dt. r(t) dt. Comparing this equation with the one listed above for the length of see that

Math Calculus with Analytic Geometry II

The Trapezoidal Rule

Math 124B January 24, 2012

99/105 Comparison of OrcaFlex with standard theoretical results

1.2. Linear Variable Coefficient Equations. y + b "! = a y + b " Remark: The case b = 0 and a non-constant can be solved with the same idea as above.

Polynomial Approximations for the Natural Logarithm and Arctangent Functions. Math 230

Energy Balance of Solar Collector

UNIT 1 FUNCTIONS AND THEIR INVERSES Lesson 1.4: Logarithmic Functions as Inverses Instruction

STRAND B: NUMBER THEORY

We partition C into n small arcs by forming a partition of [a, b] by picking s i as follows: a = s 0 < s 1 < < s n = b.

2 b. , a. area is S= 2π xds. Again, understand where these formulas came from (pages ).

6.2 CONCEPTS FOR ADVANCED MATHEMATICS, C2 (4752) AS

W. We shall do so one by one, starting with I 1, and we shall do it greedily, trying

Mutipy by r sin RT P to get sin R r r R + T sin (sin T )+ P P = (7) ffi So we hve P P ffi = m (8) choose m re so tht P is sinusoi. If we put this in b

Heat flux and total heat

Read section 3.3, 3.4 Announcements:

1 The Riemann Integral

APPROXIMATE INTEGRATION

Goals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite

4 7x =250; 5 3x =500; Read section 3.3, 3.4 Announcements: Bell Ringer: Use your calculator to solve

Fluid Flow through a Tube

Chapter 3 Single Random Variables and Probability Distributions (Part 2)

Review Exercises for Chapter 4

Scientific notation is a way of expressing really big numbers or really small numbers.

13: Diffusion in 2 Energy Groups

4. Calculus of Variations

approaches as n becomes larger and larger. Since e > 1, the graph of the natural exponential function is as below

Section 5.1 #7, 10, 16, 21, 25; Section 5.2 #8, 9, 15, 20, 27, 30; Section 5.3 #4, 6, 9, 13, 16, 28, 31; Section 5.4 #7, 18, 21, 23, 25, 29, 40

Partial Derivatives. Limits. For a single variable function f (x), the limit lim

Section 10.2 Angles and Triangles

Linear Systems COS 323

SUMMER KNOWHOW STUDY AND LEARNING CENTRE

Calculus - Activity 1 Rate of change of a function at a point.

Calculus 2: Integration. Differentiation. Integration

Conducting Ellipsoid and Circular Disk

ARITHMETIC OPERATIONS. The real numbers have the following properties: a b c ab ac

Fundamental Theorem of Calculus

Generalized Fano and non-fano networks

Designing Information Devices and Systems I Discussion 8B

2008 Mathematical Methods (CAS) GA 3: Examination 2

7.2 The Definite Integral

Lesson 1: Quadratic Equations

Coordinate Geometry. Coordinate Geometry. Curriculum Ready ACMNA: 178, 214, 294.

Math 116 Calculus II

Bypassing no-go theorems for consistent interactions in gauge theories

THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS.

1 Part II: Numerical Integration

Anonymous Math 361: Homework 5. x i = 1 (1 u i )

Euler-Maclaurin Summation Formula 1

1.1. Linear Constant Coefficient Equations. Remark: A differential equation is an equation

Review of Calculus, cont d

Section 14.3 Arc Length and Curvature

Lecture 3 Gaussian Probability Distribution

Math 115 ( ) Yum-Tong Siu 1. Lagrange Multipliers and Variational Problems with Constraints. F (x,y,y )dx

Chapter 3 Exponential and Logarithmic Functions Section 3.1

Math 131. Numerical Integration Larson Section 4.6

For convenience, we rewrite m2 s m2 = m m m ; where m is repeted m times. Since xyz = m m m nd jxyj»m, we hve tht the string y is substring of the fir

CAAM 453 NUMERICAL ANALYSIS I Examination There are four questions, plus a bonus. Do not look at them until you begin the exam.

5.1 How do we Measure Distance Traveled given Velocity? Student Notes

Multiple Integrals. Review of Single Integrals. Planar Area. Volume of Solid of Revolution

Acceptance Sampling by Attributes

Chapter 6 Continuous Random Variables and Distributions

Support vector machines for regression

Combined RKDG and LDG Method for the Simulation of the Bipolar Charge Transport in Solid Dielectrics

Figure XX.1.1 Plane truss structure

University of Washington Department of Chemistry Chemistry 453 Winter Quarter 2010 Homework Assignment 4; Due at 5p.m. on 2/01/10

5.7 Improper Integrals

Energy Bands Energy Bands and Band Gap. Phys463.nb Phenomenon

INTRODUCTION TO INTEGRATION

USA Mathematical Talent Search Round 1 Solutions Year 21 Academic Year

Lecture 6 Notes, Electromagnetic Theory I Dr. Christopher S. Baird University of Massachusetts Lowell

New data structures to reduce data size and search time

Math 8 Winter 2015 Applications of Integration

Population Dynamics Definition Model A model is defined as a physical representation of any natural phenomena Example: 1. A miniature building model.

Development of the Sinc Method for Nonlinear Integro-Differential Eequations

Generalized Uncertainty Principle, Extra-dimensions and Holography

The area under the graph of f and above the x-axis between a and b is denoted by. f(x) dx. π O

Suggested Solution to Assignment 5

Indefinite Integral. Chapter Integration - reverse of differentiation

Solution for Assignment 1 : Intro to Probability and Statistics, PAC learning

1 Probability Density Functions

Definite integral. Mathematics FRDIS MENDELU

GEOMETRIC PROBABILITY MODELS TO ANALYZE STRATEGIES FOR FINDING A BURIED CABLE

4.6 Numerical Integration

AQA Further Pure 2. Hyperbolic Functions. Section 2: The inverse hyperbolic functions

Non-Linear & Logistic Regression

Definite integral. Mathematics FRDIS MENDELU. Simona Fišnarová (Mendel University) Definite integral MENDELU 1 / 30

MATH SS124 Sec 39 Concepts summary with examples

BEAM DIAGRAMS AND FORMULAS. Nomenclature

Topics Covered AP Calculus AB

Practice Problem Set 3

Sections 5.2: The Definite Integral

Math& 152 Section Integration by Parts

Transcription:

Supporting Tet Evoution of the Averge Synptic Updte Rue In this ppendi e evute the derivtive of Eq. 9 in the min tet i.e. e need to ccute Py ( ) Py ( Y ) og γ og. [] P( y Y ) P% ( y Y ) Before e strt et us rec some nottion. The verge of n rbitrry function f ith rguments nd y is by definition f ( y ) = p( y ) f( y ) [2] y y here p( y ) denotes the oint probbiity of the pir ( y) to occur nd the sum runs over configurtions of nd y. The subscript indictes tht both the probbiity distribution p nd the function f my depend on prmeter. By definition e hve p( y ) = p( y p ) ( ) here p( ) is given input distribution nd p( y ) the (prmeter-dependent) condition probbiity of generting n output y given. Hence Eq. 2 cn be trnsformed into f ( y ) = p ( ) p( y ) f( y ) = p( y ) f( y ) y y y If e no te the derivtive ith respect to the prmeter the product rue yieds to terms f( y ) = p( y ) f( y) [4] y y + p( y ) og p( y ) f( y) y The first term contins the derivtive of the function f heres the second term contins the derivtive of the condition probbiity p. We note tht Eq. 4 cn so be ritten in the form f( y ) = f( y ) + og p( y ) f( y ) [5] y y y i.e. s n verge over the oint distribution of nd y. This formution i be usefu for the probem t hnd. The grdient in Eq. contins sever terms nd for the moment e pic ony one of these. The others i then be treted nogousy. Let us focus on the term og Py ( Y X ) nd ppy steps competey nogous to those [3]

eding from Eqs. 2-5. og ( Py ) Y X = og Py ( Y X ) + og PY ( X ) og P( y Y X ) [6] We no evute the verges using the identity = y of Eq. 6 vnishes since og Py ( Y X ). We find tht the first term on the right-hnd side y = og Py ( ) Py ( ) y {0 } = Py ( ) = 0 [7] y {0 } becuse of the normiztion of probbiities. The sme rgument cn be repeted to sho tht 0= og Py ( Y ). The reference y distribution Py % ( Y ) is by definition independent of. Hence the ony term tht gives nontrivi contribution on the right-hnd side of Eq. 6 is the second term. With n nogous rgument for the other fctors in Eq. e hve Py ( ) Py ( Y ) og γ og P( y Y ) P% ( y Y ) og P( Y X ) P( y ) P( y Y ) = og γ og [8] P( y Y ) P% ( y Y ) Y X An identifiction of the fctors F nd G in the min tet is strightforrd. From Eq. 4 in the min tet e hve og ( Py ) = y og( ρ ) + ( y )og( ρ ) [9] Hence e cn evute the fctors P y Y X ρ F = og = y og + ( y ) og ( ) ρ P( y Y ) ρ ρ P( y Y ) ρ ρ = og og ( ) og P% = ( y Y ) % ρ + % ρ G y y Furthermore e cn ccute the derivtive needed in Eq. 8 using the chin rue from Eq. 6 of the min tet i.e.

( ) = ( ) PY X P y Y X [0] = hich yieds og PY ( X ) = og Py ( ) [] = y y = ρ ε ( t t ) = ρ ρ n n n [2] We note tht in Eq. 8 the fctor og PY ( X ) hs to be mutipied ith F or ith G before ting the verge. Mutipiction genertes terms of the form y y = y y Y X y y ith Y X X For ny given input X the utocorretion < of the postsynptic neuron i hve trivi vue y y = y y for > Y X Y X Y X [3] here t is the idth of the utocorretion. As consequence y y F γ G 0 for ρ ρ = > Y X [4] Hence for > e cn truncte the sum over in Eq. 2 i.e. = hich yieds ecty the coincidence mesure = introduced in the min tet; cf. Eq. in the min tet nd hich e repet here for convenience y y n n = ρ ε ( t t ) = ρ ρ [5] n From Averges to n Onine Rue The coincidence mesure counts coincidences in rectngur time indo. If e repce the rectngur time indo by n eponenti one ith time constnt nd go to continuous time the summtion in = t Eq. 5 turns into n integr dt ep[ ( t t )/ ] hich cn be trnsformed into differenti eqution d () t ( t δ ) ( f ) = + ε( tt ) () ( ˆ St δ tt δ) gut ( ()) Rt () ; dt [6] f cf. Eq. 5 in the min tet. Bsed on the considertions in the previous prgrph the time constnt shoud best be chosen in the rnge

t 0 t. Simiry the verge firing rte ρ () t = g() t R() t cn be estimted using running verge dg() t g = gt () + gut ( ()) [7] dt ith time constnt g. In Fig. 6 e compre the performnce of three different updte schemes in numeric simutions. In prticur e sho tht (i) the ect vue of the trunction of the sum in Eq. 5 is not reevnt s ong s t is rger thn the idth of the utocorretion; nd (ii) tht the onine rue is good pproimtion to the ect soution. To do so e te the scenrio from Fig. 3 of the min tet. For ech segment of s e simute 00 pirs of input nd output spie trins. We evute numericy Eq. 8 by verging over the 00 smpes. After ech segment of second (=000 time steps) e updte the eights using rue ithout trunction in the sum of Eq. 5. We c this the fu btch updte; compre Fig. 6 (Top). Second e use the definition of ith the truncted sum nd repet the bove steps; Fig. 6 (Midde). The trunction is set to t = 200 ms hich is e bove the epected idth of the utocorretion function of the postsynptic neuron. We c this the truncted btch rue. Third e use the onine rue discussed in the min body of the pper ith = s; Fig. 6 (Bottom). omprison of top nd center grphs of Fig. 6 shos tht there is no difference in the evoution of men synptic efficcies i.e. the trunction of the sum is oed s epected from the theoretic rguments. A further comprison ith Fig. 6 Bottom shos tht updtes bsed on the onine rue dd some fuctutions to the resuts but its trend cptures nicey the evoution of the btch rues. Suppement to the Pttern Detection Prdigm In Fig. 3 e presented pttern detection prdigm here ptterns defined by input rtes ere chosen rndomy nd ppied for one second. After erning the spie count over one second is sensitive to the inde of the pttern. Fig. 7A shos the histogrm of spie counts for ech pttern. Optim cssifiction is chieved by choosing for ech spie count the pttern hich is most iey. With this criterion 8 percent of the ptterns i be cssified correcty. The updte of synptic efficcies depends on the choice of the prmeter γ in the erning rue. According the the optimity criterion in Eq. 8 of the min tet high eve of γ impies strong homeosttic contro of the firing rte of the postsynptic neuron heres o eve of γ induces ony e homeosttic contro. In order to study the roe of γ e repeted the numeric eperiments for the bove pttern detection prdigm ith vue of γ = 00 insted of our stndrd vue of γ =. Fig. 7B shos tht the output firing rte is sti moduted by the pttern inde the modution t

γ = 00 is hoever eer thn tht t γ =. As resut pttern detection is ess reiby ith 45 percent correct cssifiction ony. We note tht this is sti significnty higher thn the chnce eve of 25 percent.