Re: MATLAB code for Poisson spike generator (Erik Larsen )


Subject: Re: MATLAB code for Poisson spike generator
From:    Erik Larsen  <elarsen@xxxxxxxx>
Date:    Thu, 26 Oct 2006 10:26:31 -0400

------=_Part_21517_2338923.1161872791799 Content-Type: text/plain; charset=ISO-8859-1; format=flowed Content-Transfer-Encoding: quoted-printable Content-Disposition: inline Dear Ana, a Poisson spike generator with a fixed mean rate for a given interval is very easy - you just need to generate random spike times according to a uniform probability density distribution during your time interval (number of spikes =3D mean rate * interval duration). This of course doesn't take refractoriness into account. Perhaps a simple way to incorporate that is to look at each successive spike time that you generate and check whether it occurs within the refractory period of another spike, and if so, reject tha= t new spike and generate a new one. Erik -- Erik Larsen Speech and Hearing Bioscience and Technology Harvard-MIT Division of Health Sciences and Technology On 10/26/06, Ana Alves Pinto <aalvespinto@xxxxxxxx> wrote: > > Dear List, > > I'm looking for MATLAB code for a Poisson spike generator, given a mean > rate, a certain duration and a refractory period. Could anyone provide th= e > code to do this? > Thank you in advance for any help. > Best regards, > > Ana > > Ana Alves Pinto > > Instituto de Neurociencias de Castilla y Le=F3n > Universidad de Salamanca > Avda. Alfonso X "El Sabio" s/n Tel. +34-923 294500 ext. 1957 > 37007 Salamanca Fax. +34-923 294 750 > E-mail: aalvespinto@xxxxxxxx > > ------=_Part_21517_2338923.1161872791799 Content-Type: text/html; charset=ISO-8859-1 Content-Transfer-Encoding: quoted-printable Content-Disposition: inline Dear Ana,<br><br>a Poisson spike generator with a fixed mean rate for a giv= en interval is very easy - you just need to generate random spike times acc= ording to a uniform probability density distribution during your time inter= val (number of spikes =3D mean rate * interval duration). This of course do= esn't take refractoriness into account. Perhaps a simple way to incorporate= that is to look at each successive spike time that you generate and check = whether it occurs within the refractory period of another spike, and if so,= reject that new spike and generate a new one. <br><br>Erik<br><br>--<br>Erik Larsen<br>Speech and Hearing Bioscience and = Technology<br>Harvard-MIT Division of Health Sciences and Technology<br><br= ><div><span class=3D"gmail_quote">On 10/26/06, <b class=3D"gmail_sendername= "> Ana Alves Pinto</b> &lt;<a href=3D"mailto:aalvespinto@xxxxxxxx">aalvespinto@xxxxxxxx= usal.es</a>&gt; wrote:</span><blockquote class=3D"gmail_quote" style=3D"bor= der-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-= left: 1ex;"> Dear List,<br><br>I'm looking for MATLAB code for a Poisson spike generator= , given a mean<br>rate, a certain duration and a refractory period. Could a= nyone provide the<br>code to do this?<br>Thank you in advance for any help. <br>Best regards,<br><br>Ana<br><br>Ana Alves Pinto<br><br>Instituto de Neu= rociencias de Castilla y Le=F3n<br>Universidad de Salamanca<br>Avda. Alfons= o X &quot;El Sabio&quot; s/n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Tel. +34-92= 3 294500 ext. 1957<br>37007 Salamanca&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n= bsp;&nbsp; Fax. +34-923 294 750 <br>E-mail: <a href=3D"mailto:aalvespinto@xxxxxxxx">aalvespinto@xxxxxxxx</a><= br><br></blockquote></div><br> ------=_Part_21517_2338923.1161872791799--


This message came from the mail archive
http://www.auditory.org/postings/2006/
maintained by:
DAn Ellis <dpwe@ee.columbia.edu>
Electrical Engineering Dept., Columbia University