Hello.(adsbygoogle = window.adsbygoogle || []).push({});

I am having some trouble trying to reconstruct my firing rates using a Volterra expansion.

Basically it is known that , if :

s(t) = the value of a given stimulus at time t

r(t) = the firing rate of a neuron at time t

then we assume that a possible estimate of the firing rate can be achieved by

rest(t) = r_0 + ∫dτ D(τ)s(t-τ) where D(τ) is the Wiener Kernel.

Then the condition for the best estimate rest(t) is achieved when

FD(ω)*FQss(ω) = FQrs(-ω) . Where FA denotes the fourier transform of a function A.

While I understand all this mathematically I can not apply it using Matlab.

I have s as a vector of positions of a given stimulus and rg as a vector of

rates of a neuron .

Then I make

QSS = xcorr(S,S);

QrS = xcorr(rg,S);

FQSS = fft(QSS);

FQrS = fft(QrS);

G = FQrS(length(FQrS):-1:1); % G(ω) = FQrS(-ω) because xcorr seems to put the t = 0

%at the middle of the vector

FD = G./FQSS'

D = ifft(FD);

When I make 'rest' proportional to 'conv(D, S)' I obtain not my firing rates but rather

a very similar version of my own stimulus S.

What am I doing wrong? Thanks!

Beserra

**Physics Forums | Science Articles, Homework Help, Discussion**

Join Physics Forums Today!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# Stimulus / rate reconstruction with Wiener Kernels.

Loading...

Similar Threads for Stimulus rate reconstruction |
---|

High dose gamma ray exposure vs Low dose gamma ray exposure |

Why is the replacement fertility rate 2.1 instead of 2.0? |

Medical Vitamin B megadoses, smokers, and lung cancer rates |

**Physics Forums | Science Articles, Homework Help, Discussion**