z = [length(pre.time>=-3*3600 & pre.time<0),...
    length(drug.time(drug.time>=drugstarttime & drug.time<drugstarttime+3600)),...
    length(drug.time(drug.time>=drugstarttime+23*3600 & drug.time<drugstarttime+24*3600)),...
    length(drug.time(drug.time>=drugstarttime & drug.time<drugstarttime+24*3600)),...
    length(pre.burst.time),...
    length(drug.burst.time(drug.burst.time>=drugstarttime & drug.burst.time<drugstarttime+3600)),...
    length(drug.burst.time(drug.burst.time>=drugstarttime+23*3600 & drug.burst.time<drugstarttime+24*3600)),...
    length(drug.burst.time(drug.burst.time>=drugstarttime & drug.burst.time<drugstarttime+24*3600)),...
    mean(pre.IBI),...
    std(pre.IBI),...
    median(pre.IBI),...
    quantile(pre.IBI,.25),...
    quantile(pre.IBI,.75),...
    mean(drug.IBI),...
    std(drug.IBI),...
    median(drug.IBI),...
    quantile(drug.IBI,.25),...
    quantile(drug.IBI,.75),...
    length(pre.time(~ismember(pre.time,pre.burst.snippet.spiketimes))),...
    length(drug.time(~ismember(drug.time,drug.burst.snippet.spiketimes))),...
    (3*3600-pre.burst.snippet.sniptime),...
    (24*3600-drug.burst.snippet.sniptime)];