.. currentmodule:: brian

.. index::
   pair: example usage; raster_plot
   pair: example usage; NeuronGroup
   pair: example usage; run
   pair: example usage; SpikeMonitor
   pair: example usage; linked_var

.. _example-misc_reliability:

Example: reliability (misc)
===========================

Reliability of spike timing.
See e.g. Mainen & Sejnowski (1995) for experimental results in vitro.

R. Brette

::

    from brian import *
    
    # The common noisy input
    N = 25
    tau_input = 5 * ms
    input = NeuronGroup(1, model='dx/dt=-x/tau_input+(2./tau_input)**.5*xi:1')
    
    # The noisy neurons receiving the same input
    tau = 10 * ms
    sigma = .015
    eqs_neurons = '''
    dx/dt=(0.9+.5*I-x)/tau+sigma*(2./tau)**.5*xi:1
    I : 1
    '''
    neurons = NeuronGroup(N, model=eqs_neurons, threshold=1, reset=0, refractory=5 * ms)
    neurons.x = rand(N)
    neurons.I = linked_var(input,'x') # input.x is continuously fed into neurons.I
    spikes = SpikeMonitor(neurons)
    
    run(500 * ms)
    raster_plot(spikes)
    show()