added surface plotting code

Author: cliffckerr

code/csd.m

@@ -2,8 +2,10 @@
 % so label accordingly 
 
 % load the files 
-filename = 'C:/Users/Felicity/macaqueERP/data/epochedtemp.mat';
-errorfilename = 'C:/Users/Felicity/macaqueERP/data/errordata.mat';
+% filename = 'C:/Users/Felicity/macaqueERP/data/epochedtemp.mat';
+% errorfilename = 'C:/Users/Felicity/macaqueERP/data/errordata.mat';
+filename = '/u/cliffk/drive/usyd/macaqueERP/data/epochedtemp.mat';
+errorfilename = '/u/cliffk/drive/usyd/macaqueERP/data/errordata.mat';
 norm = load(filename);
 error = load(errorfilename);
 

code/plotcsd.m

@@ -0,0 +1,12 @@
+% errorfilename = '/u/cliffk/drive/usyd/macaqueERP/data/errordata.mat';
+errorfilename = 'C:/Users/Felicity/macaqueERP/data/errordata.mat';
+error = load(errorfilename);
+errordata = error.data{2,1}.noleverresponse(:, :, :);
+errorthingall = zeros(size(errordata,1)-2,size(errordata,2),size(errordata,3));
+for i=1:size(errordata,2)
+    tmp = squeeze(errordata(:, i, :));
+    ecsd = diff(diff(tmp));
+    errorthingall(:, i, :) = ecsd;
+end
+meanerror = squeeze(mean(errorthingall,2));
+pcolor(meanerror)