loadRepeatRecord_chunk('buster359'); loadRepeatRecord_chunk('buster364'); loadRepeatRecord_chunk('cinderella19307'); loadRepeatRecord_chunk('cinderella19309');
loadRepeatRecord_chunk('cinderella19307');
loadRepeatRecord_chunk('cinderella19309');
clc
SSprocess('cinderella19307');
edit loadRepeatRecord_chunk.m
loadRepeatRecord_chunk('cinderella19309');
SSprocess('cinderella19309');
load('C:\Users\USER\Documents\mingfaifong\buster364\buster364_FullRecNoWaveformNoClean.mat')
bins = 0:1:max(time);
asdr = hist(time,bins);
plot(bins,asdr);
load('C:\Users\USER\Documents\mingfaifong\buster359\buster359_FullRecNoWaveformNoClean.mat')
bins = 0:1:max(time);
asdr = hist(time,bins);
plot(bins,asdr);
buster359_time = time - 5.64e4;
buster359_channel = channel;
bins = 0:3600:max(time);
asdr = hist(time,bins);
plot(bins,asdr);
plot(bins,asdr,'.');
bins = min(time):3600:max(time);
asdr = hist(time,bins);
plot(bins,asdr,'.');
bins = min(buster359_time):3600:max(buster359_time);
asdr = hist(buster359_time,bins);
plot(bins,asdr,'.');
plot(bins,asdr,'.','MarkerSize',10);
plot(bins,asdr,'.','MarkerSize',15);
min(buster359_time)/3600
set(gca,[-16:2:
max(buster359_time)/3600
set(gca,'XTick',[-16:2:50]*3600,'XTickLabel',[-16:2:50])
xlim([-16:2:50]*3600)
xlim([-16 50]*3600)
grid on;
set(gca,'YTick',[0:20000:1800000]);
ylim([0:20000:1800000]);
ylim([0 1800000]);
ylim([0 180000]);
set(gca,'YTick',[0:20000:1800000]);
xlabel('Time (hours)'); ylabel('Firing Rate (Hz)');
plot(bins,asdr/3600,'.','MarkerSize',15);
xlim([-16 50]*3600)
set(gca,'XTick',[-16:2:50]*3600,'XTickLabel',[-16:2:50])
ylime([0 50]); set(gca,'YTick',[0:10:50]);
ylim([0 50]); set(gca,'YTick',[0:10:50]);
grid on;
ylim([0 50]); set(gca,'YTick',[0:5:50]);
xlabel('Time (hours)'); ylabel('Firing Rate (Hz)');
xlabel('Time (hours)','FontSize',18); ylabel('Firing Rate (Hz)','FontSize',18);
title('buster359 ASDR 1-hour bins','FontSize',18);
title('buster359 ASDR 1-hour bins','FontSize',14);
asdr'
-16:1:48
[ans' asdr']
x = -16:1:48
[x' asdr'/3600]
save buster359 asdr bins buster359*
load('C:\Users\USER\Documents\mingfaifong\buster364\buster364_FullRecNoWaveformNoClean.mat')
buster364_channel = channel;
plot(time,channel);
plot(time,channel,'.');
%-- 10/22/11 12:16 PM --%
cd buster359\
ls
cd ..
load buster359\
cd ..
cd mingfaifong\
min(buster359_time)
ISI_predrug = diff(buster359_time(buster359_time<=0))
ISI_predrug = diff(buster359_time(buster359_time<=0));
min(ISI_predrug)
max(ISI_predrug)
histc(ISI_predrug,100);
help histc
hist(ISI_predrug,100);
hist(ISI_predrug,10000);
min(ISI_predrug)
ISI_predrug = diff(sortrows(buster359_time(buster359_time<=0)));
hist(ISI_predrug,1000);
histc(ISI_predrug,0:0.01:1);
help hist
hist(ISI_predrug,0:0.01:1);
hist(ISI_predrug,0:0.001:1);
hist(ISI_predrug,0:0.0005:1);
hist(ISI_predrug,0:0.0005:max(ISI_predrug));
hist(ISI_predrug,0:0.0001:max(ISI_predrug));
hist(ISI_predrug,0:1/25000:max(ISI_predrug));
help histc
bins_ISI = [0:1/25000:max(time))];
bins_ISI = [0:1/25000:max(time)];
bins_ISI = [0:1/25000:max(buster359_time)];
bins_ISI = 0:1/25000:max(ISI_predrug);
bins_ISI = 0:1/25000:max(ISI_predrug)';
bins_ISI = [0:1/25000:max(ISI_predrug)]';
buster359_histISI = histc[buster359_time,bins_ISI);
buster359_histISI = histc(buster359_time,bins_ISI);
help loglog
semilog(bins_ISI,buster359_histISI);
semilogy(bins_ISI,buster359_histISI);
buster359_histISI = hist(buster359_time,bins_ISI);
%-- 10/22/11 12:53 PM --%
load buster359\
predrug.ISI = diff(sortrows(buster359_time(buster359_time<=0)));
bins_ISI = 0:1/25000:max(predrug.ISI)';
save buster359
help hist
help histc
[x b] = histc(predrug.ISI,bins_ISI);
help bar
bar(b,x);
help histc
x = histc(predrug.ISI,bins_ISI);
bar(bins_ISI,x);
get(gca)
set(gca,'YScale','log');
%-- 10/22/11  1:31 PM --%
load buster359\
x = histc(predrug.ISI,bins_ISI);
bar(bins_ISI,x);set(gca,'YScale','log');
drug.ISI = diff(sortrows(buster359_time(buster359_time>=0)));
save buster359; bins_ISI_drug = 0:1/25000:max(drug.ISI)';
y = histc(drug.ISI,bins_ISI_drug);
figure;
bar(bins_ISI_drug,y); set(gca,'YScale','log');
xlim([0 5]);
xlim([-0.1 5]);
xlabel('Inter-spike Interval (seconds)');
max(buster359_time)/3600
ylabel(['Frequency']);
xlabel('Inter-spike Interval (seconds)','FontSize',18);
ylabel(['Frequency'],'FontSize',18);
title(['buster359',10,'ISI Histogram during CNQX treatment (48.5 hours)'],'FontSize',18);
title(['buster359',10,'ISI Histogram during CNQX treatment (48.5 hours)'],'FontSize',14);
set(gca,'YScale','linear');
ylim([0 1]);
ylim([0 10]);
ylim([0 1000]);
save buster359;
bar(bins_ISI,x);
set(gca,'YScale','log');
xlim([-0.1 5]);
xlabel('Inter-spike Interval (seconds)','FontSize',18);
ylabel(['Frequency'],'FontSize',18);
title(['buster359',10,'ISI Histogram before CNQX treatment (x hours)'],'FontSize',14);
min(buster359_time)/3600
title(['buster359',10,'ISI Histogram before CNQX treatment (15.67 hours)'],'FontSize',14);
set(gca,'YScale','linear');
ylim([0 1000]);
median(drug.ISI);
median(drug.ISI)
help median
length(predrug.ISI)/10
predrug.ISI(length(predrug.ISI)/10)
predrug.ISI(round(length(predrug.ISI)/10))
predrug.ISI(round(length(predrug.ISI)/10*2))
predrug.ISI(round(length(predrug.ISI)/10*3))
predrug.ISI(round(length(predrug.ISI)/10*4))
predrug.ISI(round(length(predrug.ISI)/10*5))
predrug.ISI(round(length(predrug.ISI)/10*6))
predrug.ISI(round(length(predrug.ISI)/10*7))
predrug.ISI(round(length(predrug.ISI)/10*8))
predrug.ISI(round(length(predrug.ISI)/10*9))
median(predrug.ISI)
predrug.ISI(round(length(predrug.ISI)/2))
help median
size(predrug.ISI)
predrug.ISI = sortrows(predrug.ISI);
size(predrug.ISI)
median(predrug.ISI)
predrug.ISI(round(length(predrug.ISI)/2))
predrug.ISI(round(length(predrug.ISI)/10))
predrug.ISI(round(length(predrug.ISI)/4))
predrug.ISI(round(length(predrug.ISI)/4*3))
predrug.ISI(round(length(predrug.ISI)/10*9))
tmp = sortrows(drug.ISI);
tmp(round(length(tmp.ISI)/10))
tmp(round(length(tmp)/10))
tmp(round(length(tmp)/4))
tmp(round(length(tmp)/2))
median(tmp)
tmp(round(length(tmp)/4*3))
tmp(round(length(tmp)/10*9))
xlim[(-0.001 0.0046])
xlim([-0.001 0.0046])
xlim([-0.0001 0.0046])
xlim([-0.0001 0.138])
xlim([-0.001 0.138])
set(gca,'YScale','linear');
xlim([-0.0001 0.0046])
bar(bins_ISI_drug,y); set(gca,'YScale','log');
xlim([-0.001 0.138])
ylabel(['Frequency'],'FontSize',18);
xlabel('Inter-spike Interval (seconds)','FontSize',18);
title(['buster359',10,'ISI Histogram during CNQX treatment (48.5 hours)'],'FontSize',14);
close all;
clear all;
pack;
oct22
oct22_tmp
cd cinderella19307\
cd SSfiles\
load cinderella19307_15to20hr.SS
load('cinderella19307_10to10.5hr.SS','-mat')
SS
max(SS.unit)
min(SS.unit)
x = SS.channel(SS.unit==1);
mean(SS.channel(SS.unit==1))
max(SS.channel(SS.unit==1))
min(SS.channel(SS.unit==1))
cd ..
batch_plotAvgWaveform('cinderella363');
clear;
batch_plotAvgWaveform('cinderella363');
SS
dbcont
batch_plotAvgWaveform('cinderella363');
grid on
set(gca,'Color','r')
set(gca,'Color',[.1 .1 .1])
set(gca,'Color',[.8 .8 .8])
set(gca,'LineColor',[.8 .8 .8])
set(gca,'MarkerColor','r')
get(gca)
batch_plotAvgWaveform('cinderella363');