Ambiguous Encoding and Distorted Perception (Carlson and Kawasaki 2006)

 Download zip file 
Help downloading and running models
"... In the weakly electric fish Eigenmannia, P- and T-type primary afferent fibers are specialized for encoding the amplitude and phase, respectively, of electrosensory stimuli. We used a stimulus estimation technique to quantify the ability of P- and T-units to encode random modulations in amplitude and phase. As expected, P-units exhibited a clear preference for encoding amplitude modulations, whereas T-units exhibited a clear preference for encoding phase modulations. Surprisingly, both types of afferents also encoded their nonpreferred stimulus attribute when it was presented in isolation or when the preferred stimulus attribute was sufficiently weak. Because afferent activity can be affected by modulations in either amplitude or phase, it is not possible to unambiguously distinguish between these two stimulus attributes by observing the activity of a single afferent fiber. Simple model neurons with a preference for encoding either amplitude or phase also encoded their nonpreferred stimulus attribute when it was presented in isolation, suggesting that such ambiguity is unavoidable. ... " See paper for more and details.
1 . Carlson BA, Kawasaki M (2006) Ambiguous encoding of stimuli by primary sensory afferents causes a lack of independence in the perception of multiple stimulus attributes. J Neurosci 26:9173-83 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Neuron or other electrically excitable cell;
Brain Region(s)/Organism:
Cell Type(s): Electric fish P- and T-type primary afferent fibers;
Gap Junctions:
Simulation Environment: MATLAB;
Model Concept(s): Oscillations; Sensory coding;
%   genstim.m
%       Creates stimulus waveform for PrimAff_IntFire.m

global tim sig ram rpm

% Initialize time variables
s_rate = 1000/tstep;
dur = tstop/1000;

% Set the EOD frequency
EODfrq = 400;

% Prepare pre-stimulus to avoid edge effects
[pre_tim,pre_sig] = modsine(EODfrq,0,10,0,0,0,0,0,s_rate,0.1);

% Prepare stimulus
if stim=='No Modulation ',
    [tim,sig] = modsine(EODfrq,0,10,0,0,0,0,0,s_rate,dur);
elseif stim=='Sinusoidal AM ',
    [tim,sig] = modsine(EODfrq,0,10,freq,amdepth,0,0,0,s_rate,dur);
elseif stim=='Sinusoidal PM ',
    [tim,sig] = modsine(EODfrq,0,10,freq,0,pmdepth,0,0,s_rate,dur);
elseif stim=='+DF           ',
    [tim,sig] = modsine(EODfrq,0,10,freq,amdepth,pmdepth,0,90,s_rate,dur);
elseif stim=='-DF           ',
    [tim,sig] = modsine(EODfrq,0,10,freq,amdepth,pmdepth,0,-90,s_rate,dur);
elseif stim=='Random AM     ',
    [tim,sig,ram,rpm] = blim_whnoise(EODfrq,10,s_rate,dur,freq,amdepth,1,0);
elseif stim=='Random PM     ',
    [tim,sig,ram,rpm] = blim_whnoise(EODfrq,10,s_rate,dur,freq,pmdepth,2,0);
elseif stim=='Random AM + PM',
    [tim,sig,ram,rpm] = blim_whnoise(EODfrq,10,s_rate,dur,freq,[amdepth pmdepth],3,0);
tim = [1000*pre_tim'; 1000*(tim+0.1)'];
sig = [pre_sig'; sig'];

Loading data, please wait...