% load('C:\Users\USER\Documents\mingfaifong\2012-10-28_jujuI5372_11-13DIV_MEA_CNQX\jujuI5372_11-13DIV.mat')
% load('C:\Users\USER\Documents\mingfaifong\2013-03-19_dojo22493_11-13DIV_MEA_CNQX\dojo22493_11-13DIV.mat')
% load('C:\Users\USER\Documents\mingfaifong\2013-01-24_usain22493_12-14DIV_MEA_CNQX\usain22493_12-14DIV.mat')
% load('C:\Users\USER\Documents\mingfaifong\2013-02-03_wingnut22468_10-12DIV_MEA_CNQX\wingnut22468_10-12DIV.mat')
% load('C:\Users\USER\Documents\mingfaifong\2013-02-08_xinhai22493_11-13DIV_MEA_CNQX\xinhai22493_11-13DIV.mat')
% load('C:\Users\USER\Documents\mingfaifong\2013-03-01_atlas22493_10-12DIV_MEA_CNQX\atlas22493_10-12DIV.mat')
load('C:\Users\USER\Documents\mingfaifong\2013-01-21_trouble22467_12-14DIV_MEA_CNQX\trouble22467_12-14DIV.mat')

excludechannels = [1 8 33 57 64];
clear dat FR


predrug3hr = length(time(ismember(channel,excludechannels)==0 ...
    & time>=-3*3600 & time<0))/(3*3600);

CNQXhr1 = length(time(ismember(channel,excludechannels)==0 ...
    & time>=(0*3600) & time<(1*3600)))/(1*3600);

CNQXhr24 = length(time(ismember(channel,excludechannels)==0 ...
    & time>=(23*3600) & time<(24*3600)))/(1*3600);

CNQXfull24 = length(time(ismember(channel,excludechannels)==0 ...
    & time>=(0*3600) & time<=(24*3600)))/(24*3600);

deltaFR = [CNQXhr1 CNQXhr24 CNQXfull24]/predrug3hr