Homec4science
Diffusion myhug (master) Browse

myh/dataProcessing/bci/analyses/clinical_analyses/sb.m
ab8049f9d919master

Last Change
Raw File

sb.m
View Options

This document is not UTF8. It was detected as Shift JIS and converted to UTF8 for display.
%% Play around gramm (ggplot for matlab)
% https://fr.mathworks.com/products/matlab/plot-gallery.html
% https://blogs.mathworks.com/images/pick/Sean/maingramm/grammexamples.html
%
load example_data;
clear g
g(1,1)=gramm('x',cars.Origin_Region,'y',cars.Horsepower,'color',cars.Cylinders,'subset',cars.Cylinders~=3 & cars.Cylinders~=5);
g(1,2)=copy(g(1));
g(2,1)=copy(g(1));
g(2,2)=copy(g(1));
%Raw data as scatter plot
g(1,1).geom_point();
g(1,1).set_title('geom_point()');
%Jittered scatter plot
g(1,2).geom_jitter('width',0.4,'height',0);
g(1,2).set_title('geom_jitter()');
%Averages with confidence interval
g(2,1).stat_summary('geom',{'bar','black_errorbar'});
g(2,1).set_title('stat_summary()');
%Boxplots
g(2,2).stat_boxplot();
g(2,2).set_title('stat_boxplot()');
% g(2,2).axe_property('XLim',[-5 5]);
%These functions can be called on arrays of gramm objects
g.set_names('x','Origin','y','Horsepower','color','# Cyl');
g.set_title('Visualization of Y~X relationships with X as categorical variable');
figure('Position',[100 100 800 550]);
g.draw();
%% Trying to change the legend name
load example_data;
clear g
g=gramm('x',cars.Origin_Region,'y',cars.Horsepower,'color',cars.Cylinders,'subset',cars.Cylinders~=3 & cars.Cylinders~=5);
g.stat_boxplot();
g.set_title('stat_boxplot()');
g.legend_axe_handle.Children(1).String='C';
figure('Position',[100 100 800 550]);
g.draw();
%%
g.set_color_options('lightness_range',[40 80],'chroma_range',[80 40]);
g.set_names('column','','row','');
%g.axe_property('color',[0.9 0.9 0.9],'XGrid','on','YGrid','on','GridColor',[1 1 1],'GridAlpha',0.8,'TickLength',[0 0],'XColor',[0.3 0.3 0.3],'YColor',[0.3 0.3 0.3])
%% Play with stat_summary
clear g
g=gramm('x',cars.Origin_Region,'y',cars.Horsepower,'color',cars.Cylinders,'subset',cars.Cylinders~=3 & cars.Cylinders~=5);
%Jittered scatter plot
g.geom_jitter('width',0.4,'height',0);
%These functions can be called on arrays of gramm objects
g.set_names('x','Origin','y','Horsepower','color','# Cyl');
g.set_title('Visualization of Y~X relationships with X as categorical variable');
figure('Position',[100 100 800 550]);
g.draw();
clear g
g=gramm('x',cars.Origin_Region,'y',cars.Horsepower,'color',cars.Cylinders,'subset',cars.Cylinders~=3 & cars.Cylinders~=5);
g.stat_summary('geom',{'bar','black_errorbar'});
%These functions can be called on arrays of gramm objects
g.set_names('x','Origin','y','Horsepower','color','# Cyl');
g.set_title('Visualization of Y~X relationships with X as categorical variable');
figure('Position',[100 100 800 550]);
g.draw();
clear g
g=gramm('x',cars.Origin_Region,'y',cars.Horsepower,'color',cars.Cylinders,'subset',cars.Cylinders~=3 & cars.Cylinders~=5);
% g.stat_summary('type', 'std', 'geom', 'black_errorbar');
% g.stat_summary('geom', {'point', 'area'});
g.stat_summary('type', 'std');
%These functions can be called on arrays of gramm objects
g.set_names('x','Origin','y','Horsepower','color','# Cyl');
g.set_title('Visualization of Y~X relationships with X as categorical variable');
figure('Position',[100 100 800 550]);
g.draw();
%%Corner histogram
clear g
g=gramm('x',(cars.Horsepower-nanmean(cars.Horsepower))/nanstd(cars.Horsepower),'y',-(cars.Acceleration-nanmean(cars.Acceleration))/nanstd(cars.Acceleration),'color',cars.Cylinders,'subset',cars.Cylinders~=3 & cars.Cylinders~=5);
g.geom_point();
g.stat_cornerhist('edges',-4:0.2:4,'aspect',0.6);
g.geom_abline();
g.set_title('stat_cornerhist()');
g.set_names('x','z(Horsepower)','y','-z(Acceleration)');
figure('Position',[100 100 800 550]);
g.draw();
%% Dodge with boxplot
%Create data
clear
x=repmat(1:10,1,100);
catx=repmat({'A' 'B' 'C' 'F' 'E' 'D' 'G' 'H' 'I' 'J'},1,100);
y=randn(1,1000)*3;
c=repmat([1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 3 2 2 2 1 1 2 2 2 2 3 2 3 2 1 1 2 2 2 2 2 2 2],1,25);
y=2+y+x+c*0.5;
clear g
g(1,1)=gramm('x',catx,'y',y,'color',c);
g(2,1)=copy(g(1));
g(3,1)=copy(g(1));
g(4,1)=copy(g(1));
g(5,1)=copy(g(1));
g(1,1).stat_boxplot();
g(1,1).geom_vline('xintercept',0.5:1:10.5,'style','k-');
g(1,1).set_title('''width'',0.6,''dodge'',0.7 (Default)');
g(2,1).stat_boxplot('width',0.5,'dodge',0);
g(2,1).geom_vline('xintercept',0.5:1:10.5,'style','k-');
g(2,1).set_title('''width'',0.5,''dodge'',0');
g(3,1).stat_boxplot('width',1,'dodge',1);
g(3,1).geom_vline('xintercept',0.5:1:10.5,'style','k-');
g(3,1).set_title('''width'',1,''dodge'',1');
g(4,1).stat_boxplot('width',0.6,'dodge',0.4);
g(4,1).geom_vline('xintercept',0.5:1:10.5,'style','k-');
g(4,1).set_title('''width'',0.6,''dodge'',0.4');
g(5,1).facet_grid([],c);
g(5,1).stat_boxplot('width',0.5,'dodge',0);
g(5,1).set_title('''width'',0.5,''dodge'',0');
g.set_title('Dodge and spacing options for stat_boxplot()');
figure('Position',[100 100 800 1000]);
g.draw();
%% Facet_grid Generating fake data
clear
N=2000;
colval={'A' 'B' 'C'};
rowval={'I' 'II'};
cind=randi(3,N,1);
c=colval(cind);
rind=randi(2,N,1);
r=rowval(rind);
x=randn(N,1);
y=randn(N,1);
x(cind==1 & rind==1) = x(cind==1 & rind==1)*5;
x=x+cind*3;
y(cind==3 & rind==2)=y(cind==3 & rind==2)*3;
y=y-rind*4;
clear g
g(1,1)=gramm('x',x,'y',y,'color',c,'lightness',r);
g(1,2)=copy(g(1));
g(2,1)=copy(g(1));
g(2,2)=copy(g(1));
g(3,1)=copy(g(1));
g(3,2)=copy(g(1));
g(1,1).geom_point();
g(1,1).set_title('No facets');
g(1,2).facet_grid(r,c);
g(1,2).geom_point();
g(1,2).no_legend();
g(1,2).set_title('facet_grid()');
g(2,1).facet_grid(r,c,'scale','free');
g(2,1).geom_point();
g(2,1).no_legend();
g(2,1).set_title('facet_grid(''scale'',''free'')');
g(2,2).facet_grid(r,c,'scale','free','space','free');
g(2,2).geom_point();
g(2,2).no_legend();
g(2,2).set_title('facet_grid(''scale'',''free'',''space'',''free'')');
g(3,1).facet_grid(r,c,'scale','free_x');
g(3,1).geom_point();
g(3,1).no_legend();
g(3,1).set_title('facet_grid(''scale'',''free_x'')');
g(3,2).facet_grid(r,c,'scale','independent');
g(3,2).geom_point();
g(3,2).no_legend();
g(3,2).set_title('facet_grid(''scale'',''independent'')');
g.set_color_options('lightness_range',[40 80],'chroma_range',[80 40]);
g.set_names('column','','row','');
%g.axe_property('color',[0.9 0.9 0.9],'XGrid','on','YGrid','on','GridColor',[1 1 1],'GridAlpha',0.8,'TickLength',[0 0],'XColor',[0.3 0.3 0.3],'YColor',[0.3 0.3 0.3])
g.set_title('facet_grid() options');
figure('Position',[100 100 800 800]);
g.draw();
%% stat
help ranksum % jamar et les autres
p = ranksum(x,y)
% x : bci
% y : sham
% p : pvalue
% une valeur par patient. non param騁rique
[h,p] = ttest2(x,y) %FMA
c4science · Help