12345678910111213141516171819202122232425262728293031323334353637 |
- function spm_barh(E,C,P)
- % density plotting function (c.f. bar - horizontal)
- % FORMAT spm_barh(E,C,[P])
- % E - (n x 1) expectation
- % C - (n x 1) variances
- % P - (n x 1) priors
- %___________________________________________________________________________
- % Copyright (C) 2008 Wellcome Trust Centre for Neuroimaging
- % Karl Friston
- % $Id: spm_barh.m 1143 2008-02-07 19:33:33Z spm $
- % densities
- %---------------------------------------------------------------------------
- n = length(E);
- H = zeros(2*n + 1,64);
- x = linspace(min(E - 4*sqrt(C)),max(E + 4*sqrt(C)),64);
- for i = 1:n
- H(2*i,:) = exp(-[x - E(i)].^2/(2*C(i)));
- end
- imagesc(x,[0:n]+ 0.5,1 - H)
- set(gca,'Ytick',[1:n])
- grid on
- % confidence intervals based on conditional variance
- %---------------------------------------------------------------------------
- for i = 1:n
- z = spm_invNcdf(0.05)*sqrt(C(i));
- line([-z z] + E(i),[i i],'LineWidth',4);
- if nargin == 3
- line([P(i) P(i)],[-.4 .4] + i,'LineWidth',2,'Color','r');
- end
- end
|