Mercurial > hg > ltpda
view m-toolbox/classes/+utils/@math/boxplot.m @ 44:409a22968d5e default
Add unit tests
author | Daniele Nicolodi <nicolodi@science.unitn.it> |
---|---|
date | Tue, 06 Dec 2011 18:42:11 +0100 |
parents | f0afece42f48 |
children |
line wrap: on
line source
% BOXPLOT draw box plot on data %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % % DESCRIPTION: Boxplot is a convenient way of graphically depicting groups % of numerical data. The bottom and top of the box are always the 25th and % 75th percentile (the lower and upper quartiles, respectively), and the % band near the middle of the box is always the 50th percentile (the % median). The ends of the whiskers are the percentiles corresponding to % the confidence level defined by the user. % % CALL: boxplot(d1,...,dn, ops) % % INPUT: - d1,...dn, are column vectors. Different data must be input % as different columns. % - ops, a struct containing input parameters % - PlotData. Decide to plot all data or only the data % outside the confidence levels. Default true % - ConfLevel. Confidence level for box whiskers. Default 95% % - FontSize. Font size for x and y labels. Default 22 % - LineWidth. Box line width. Default 2 % - BoxWidth. the width of the box. Default 1 % - Marker. Data marker type. Default '.' % - MarkerSize. Data marker size. Default 1 % - XTickLabel. A cell containing a cell aray of x tick % labels. Default {{''}} % % acknowledgements % A special thanks to Shane Lin who submitted the boxPlot function in % Matlab Central to which the present function is inspired % % % Luigi Ferraioli 13-02-2011 % % $Id: boxplot.m,v 1.4 2011/02/22 14:02:03 luigi Exp $ % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function boxplot(varargin) %%% check and set imput options % Default input struct defaultparams = struct(... 'PlotData',true,... 'ConfLevel',0.95,... 'FontSize',22,... 'LineWidth',2,... 'BoxWidth',1,... 'Marker','.',... 'MarkerSize',1,... 'XTickLabel',{{''}}); names = {'PlotData','ConfLevel','FontSize','LineWidth','BoxWidth',... 'Marker','MarkerSize','XTickLabel'}; N = numel(varargin); if isstruct(varargin{N}) ops = varargin{N}; Nloop = N-1; else ops = {}; Nloop = N; end % collecting input and default params if ~isempty(ops) for jj=1:length(names) if isfield(ops, names(jj)) defaultparams.(names{1,jj}) = ops.(names{1,jj}); end end end conf = defaultparams.ConfLevel; % confidence level for confidence bounds calculation if conf>1 conf = conf/100; end plotdata = defaultparams.PlotData; fontsize = defaultparams.FontSize; lwidth = defaultparams.LineWidth; bwidth = defaultparams.BoxWidth; mark = defaultparams.Marker; mksize = defaultparams.MarkerSize; xticklabel = defaultparams.XTickLabel; %%% draw the plot drawBox(varargin(1:Nloop),Nloop,conf,plotdata,lwidth,bwidth,mark,mksize) % set ylabel ylabel('Values','FontSize',fontsize) % set xlabel xlabel('Experiment','FontSize',fontsize) % set XTickLabel if ~isempty(xticklabel{1}) set(gca,'XTickLabel',xticklabel) end end function drawBox(data,Nloop,conf,plotdata,lwidth,bwidth,mark,mksize) % define box width unit = (1-1/(1+Nloop))/(1+9/(bwidth+3)); figure; hold on; v = zeros(5,1); for ii = 1:Nloop % sort data acsending sdata = sort(data{ii}); % get data median v(1) = median(sdata); % get cdf [F,x] = utils.math.ecdf(sdata); % set box lower and upper limit, corresponding to a 25% and 75% % probability v(2) = interp1(F,x,0.25); % 25% limit v(3) = interp1(F,x,0.75); % 75% limit % set wisker limits alpm = (1-conf)/2; v(4) = interp1(F,x,alpm); % lower limit v(5) = interp1(F,x,1-alpm); % upper limit % draw data if plotdata % plot all data stddata = std(sdata); % set x for data uunit = unit.*exp(-1.*((sdata-v(1)).^2)./(stddata).^2).*9./10; xdat = ii-uunit + 2*uunit.*rand(size(sdata)); plot(xdat,sdata,[mark 'b'], 'MarkerSize', mksize); else % plot only data outside confidence levels idx1 = find(sdata<v(4)); idx2 = find(sdata>v(5)); idx = [idx1;idx2]; xdat = ii.*ones(size(idx)); plot(xdat,sdata(idx),[mark 'b'], 'MarkerSize', mksize); end % draw the min line plot([ii-unit, ii+unit], [v(4), v(4)], 'k', 'LineWidth', lwidth); % draw the max line plot([ii-unit, ii+unit], [v(5), v(5)], 'k', 'LineWidth', lwidth); % draw middle line plot([ii-unit, ii+unit], [v(1), v(1)], 'r', 'LineWidth', lwidth); % draw vertical line plot([ii, ii], [v(4), v(2)], '--k', 'LineWidth', lwidth); plot([ii, ii], [v(3), v(5)], '--k', 'LineWidth', lwidth); % draw box plot([ii-unit, ii+unit, ii+unit, ii-unit, ii-unit], [v(3), v(3), v(2), v(2), v(3)], 'k', 'LineWidth', lwidth); end end