Home > data-utility > granger_causality > gcause_libs > ksplotting.m

ksplotting

PURPOSE ^

SYNOPSIS ^

This is a script file.

DESCRIPTION ^

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 clear all;
0002 
0003 % Load  data
0004 load data_real_catM1.mat;
0005 load result_real_catM1.mat;
0006 
0007 % Dimension of X (# Channels x # Samples x # Trials)
0008 [CHN SMP TRL] = size(X);
0009 
0010 % Selected spiking history orders by AIC
0011 ht = 3*[1 11 1 5 18 17 1 1 15 17 1 11 1 9 5];
0012 % ht = [1 1 1 1 1 1 1 1 1 1 1 7 1 2 1];       % Non-movement
0013 
0014 for neuron = 1:CHN:
0015 
0016     NEURON = num2str(neuron);
0017     HHT = num2str(ht(neuron));
0018 
0019     [z_sorted,b,bp,bn] = ksplot_trial(X,bhat{ht(neuron),neuron},ht(neuron),neuron);
0020     % load(['KSplotting_neuron' NEURON '_ht' HHT]);
0021 
0022     figure(neuron);plot(b,z_sorted);ylim([0 1]);
0023     hold on
0024     plot(b,b,'r');
0025     plot(b,bp,':');
0026     plot(b,bn,':');
0027     xlabel('Qunatiles');
0028     ylabel('Cumulative Distribution Function');
0029     title(['Neuron-' NEURON ' order-' HHT]);
0030 
0031     % save(['KSplotting_neuron' NEURON '_ht' HHT],'z_sorted','b','bn','bp');
0032 
0033 end
0034

Generated on Wed 24-May-2017 00:00:56 by m2html © 2005