12345678910111213141516171819202122232425262728293031323334353637383940414243444546 |
- function pval = circ_medtest(alpha,md)
- %
- % [pval, z] = circ_medtest(alpha,w)
- % Tests for significance of the median.
- % H0: the population has median angle md
- % HA: the population has not median angle md
- %
- % Input:
- % alpha sample of angles in radians
- % md median to test for
- %
- % Output:
- % pval p-value
- %
- % PHB 3/19/2009
- %
- % References:
- % Biostatistical Analysis, J. H. Zar, 27.4
- %
- % Circular Statistics Toolbox for Matlab
- % By Philipp Berens, 2009
- % berens@tuebingen.mpg.de - www.kyb.mpg.de/~berens/circStat.html
- if size(alpha,2) > size(alpha,1)
- alpha = alpha';
- end
- if length(md)>1
- error('Median can only be a single value.')
- end
- n = length(alpha);
- % compute deviations from median
- d = circ_dist(alpha,md);
- n1 = sum(d<0);
- n2 = sum(d>0);
- % compute p-value with binomial test
- pval = sum(binopdf([0:min(n1,n2) max(n1,n2):n],n,0.5));
|