directionselectivity4celldemo

  DIRECTIONSELECTIVITY4CELLDEMO - Demo of a 4 input, single layer direction-selective output cell

   [MODELEL, DI, R_UP, R_DOWN] = DIRECTIONSELECTIVITY4CELLDEMO

   Enter "TYPE DIRECTIONSELECTIVITY4CELLDEMO" to see the code.
   Input cells: 1-4
   Cells 1 and 2 are at position 1, cells 3 and 4 are at position 2. Cells 1 and 3 fire instantly,
   while cells 2 and 4 have a latency.

   One can modify the default parameters by passing name/value pairs to the function like this
   [MODELEL, DI, UP, DOWN] = DIRECTIONSELECTIVITY4CELLDEMO(PARAM1NAME, PARAM1VALUE, ...)

   Default parameters are as follows:
   ----------------------------------------------------------------------
   dt (1e-4)                      |  Step size of model
   latency (0.200)                |  Latency of between cell 1, 2 and 3, 4
   lag (0.200)                    |  Lag of stimulus arrival between cell 1, cell3
   Syn_Gmax_initial (1x4, 4e-9)   |  Initial synaptic weight (about 75% of threshold)
   plasticity_params ([])         |  Synapse plasticity params 
   plasticity_method ('')         |  Synapse plasticity method
   isi (0.75)                     |  Inter-stimulus-interval (seconds)
   plotit (1)                     |  0/1 should we plot the behavior in a new figure?
   simit (1)                      |  0/1 Actually do simulation (1) or just build model(0)?
   simup (1)                      |  0/1 Actually simulate upward direction (1)?
   simdown (1)                    |  0/1 Actually simulate downward direction (1)?
   intfireparams                  |  {'name1','value1'} parameter list to pass to intfireleakyel_init
   synapseparams                  |  {'name1','value1'} parameter list to pass to synapset1t2el_init


   Outputs:
   MODELEL is the model produced by the function.  
   DI is (R_UP-R_DOWN) / (R_UP+R_DOWN), note this is -1 to 1
   R_UP is response to upward direction
   R_DOWN is response to downward direction

   Examples:
      % a totally non-direction-selective model
      [mel,di,rup,rdown]=directionselectivity4celldemo;

      % a direction-selective model
      [mel,di,rup,rdown]=directionselectivity4celldemo('Syn_Gmax_initial',[0 4 4 0]*1e-9);

   See also: INTFIRELEAKY_SYNAPSEDEMO