Revision 26:b03ef38fe497 multithreshold 1.46/testAN.m
| multithreshold 1.46/testAN.m | ||
|---|---|---|
| 108 | 108 |
numHSRfibers=numLSRfibers; |
| 109 | 109 |
|
| 110 | 110 |
LSRspikes=ANoutput(1:numLSRfibers,:); |
| 111 |
PSTH=UTIL_makePSTH(LSRspikes, ANdt, localPSTHbinwidth);
|
|
| 111 |
PSTH=UTIL_PSTHmaker(LSRspikes, ANdt, localPSTHbinwidth);
|
|
| 112 | 112 |
PSTHLSR=mean(PSTH,1)/localPSTHbinwidth; % across fibers rates |
| 113 | 113 |
PSTHtime=localPSTHbinwidth:localPSTHbinwidth:... |
| 114 | 114 |
localPSTHbinwidth*length(PSTH); |
| ... | ... | |
| 117 | 117 |
|
| 118 | 118 |
% HSR |
| 119 | 119 |
HSRspikes= ANoutput(end- numHSRfibers+1:end, :); |
| 120 |
PSTH=UTIL_makePSTH(HSRspikes, ANdt, localPSTHbinwidth);
|
|
| 120 |
PSTH=UTIL_PSTHmaker(HSRspikes, ANdt, localPSTHbinwidth);
|
|
| 121 | 121 |
PSTH=mean(PSTH,1)/localPSTHbinwidth; % sum across fibers (HSR only) |
| 122 | 122 |
AN_HSRonset(levelNo)= max(PSTH); |
| 123 | 123 |
AN_HSRsaturated(levelNo)= mean(PSTH(round(length(PSTH)/2): end)); |
| ... | ... | |
| 150 | 150 |
[nCNneurons c]=size(CNoutput); |
| 151 | 151 |
nLSRneurons=round(nCNneurons/nTaus); |
| 152 | 152 |
CNLSRspikes=CNoutput(1:nLSRneurons,:); |
| 153 |
PSTH=UTIL_makePSTH(CNLSRspikes, ANdt, localPSTHbinwidth);
|
|
| 153 |
PSTH=UTIL_PSTHmaker(CNLSRspikes, ANdt, localPSTHbinwidth);
|
|
| 154 | 154 |
PSTH=sum(PSTH)/nLSRneurons; |
| 155 | 155 |
CNLSRrate(levelNo)=mean(PSTH(round(length(PSTH)/2):end))/localPSTHbinwidth; |
| 156 | 156 |
|
| 157 | 157 |
%CN HSR |
| 158 | 158 |
MacGregorMultiHSRspikes=... |
| 159 | 159 |
CNoutput(end-nLSRneurons:end,:); |
| 160 |
PSTH=UTIL_makePSTH(MacGregorMultiHSRspikes, ANdt, localPSTHbinwidth);
|
|
| 160 |
PSTH=UTIL_PSTHmaker(MacGregorMultiHSRspikes, ANdt, localPSTHbinwidth);
|
|
| 161 | 161 |
PSTH=sum(PSTH)/nLSRneurons; |
| 162 | 162 |
PSTH=mean(PSTH,1)/localPSTHbinwidth; % sum across fibers (HSR only) |
| 163 | 163 |
|
| ... | ... | |
| 174 | 174 |
[nICneurons c]=size(ICoutput); |
| 175 | 175 |
nLSRneurons=round(nICneurons/nTaus); |
| 176 | 176 |
ICLSRspikes=ICoutput(1:nLSRneurons,:); |
| 177 |
PSTH=UTIL_makePSTH(ICLSRspikes, ANdt, localPSTHbinwidth);
|
|
| 177 |
PSTH=UTIL_PSTHmaker(ICLSRspikes, ANdt, localPSTHbinwidth);
|
|
| 178 | 178 |
ICLSRsaturated(levelNo)=mean(PSTH(round(length(PSTH)/2):end))/localPSTHbinwidth; |
| 179 | 179 |
|
| 180 | 180 |
%IC HSR |
| 181 | 181 |
MacGregorMultiHSRspikes=... |
| 182 | 182 |
ICoutput(end-nLSRneurons:end,:); |
| 183 |
PSTH=UTIL_makePSTH(MacGregorMultiHSRspikes, ANdt, localPSTHbinwidth);
|
|
| 183 |
PSTH=UTIL_PSTHmaker(MacGregorMultiHSRspikes, ANdt, localPSTHbinwidth);
|
|
| 184 | 184 |
PSTH=sum(PSTH)/nLSRneurons; |
| 185 | 185 |
PSTH=mean(PSTH,1)/localPSTHbinwidth; % sum across fibers (HSR only) |
| 186 | 186 |
|
Also available in: Unified diff