1234567891011121314151617181920212223242526 |
- function [output] = median_ci(data,nBS,quantileBS)
- % MEDIAN_CI computes median and confidence interval of the median using bootstrap method
- % Syntax: [output] = median_ci(data,nBS,quantileBS)
- % Example:
- % data = [1 1 1 1 1; 2 2 2 2 2; 3 3 3 3 3];
- % nBS=1000;
- % quantileBS = 0.5;
- % [output] = median_ci(data,nBS,quantileBS)
- nsubj = size(data,1);
- nlevel = size(data,2);
- quantileBS=1-quantileBS;
- quantileBS = [quantileBS/2 1-quantileBS/2];
- output = NaN(3,nlevel);
- output(2,:) = median(data,1);
- for l=1:nlevel
- temp = data(:,l);
- medianBS = NaN(nBS,1);
- for b=1:nBS
- medianBS(b) = median( datasample(temp,length(temp)) );
- end
- output([1 3],l) = quantile(medianBS,quantileBS);
- end
|