我有一组随时间变化的曲线,它们存储在MATLAB矩阵中.矩阵的每一行都是这些曲线之一,随着时间的推移展开.那些是随机实验的重复.
我需要绘制这些曲线随时间的平均值,以及95%置信区间.
我对统计数据的理解相当差,但我建议使用MATLAB的bootci函数来使用bootstrap置信区间.
我在MATLAB中实现了一个最小的例子,但我有些疑惑.我希望你能帮助我更好地掌握这个并避免愚蠢的错误.
这是一个例子:
NVARIABLES = 200;
NOBSERVATIONS = 1000;
RESAMPLING = 10000;
DATA = rand(NOBSERVATIONS, NVARIABLES);
[CI, STAT] = bootci(RESAMPLING, @mean, DATA);
MEAN = mean(DATA); %
x = 1:NVARIABLES;
figure;
hold on;
plot(x, MEAN, 'LineWidth', 2);
plot(x, CI(1,:), '--', 'LineWidth', 2); % [2]
plot(x, CI(2,:), '--', 'LineWidth', 2);
% plot(x, MEAN-CI(1,:)); % ?
% plot(x, MEAN+CI(2,:)); % ?
hold off;
这是我的问题:
>我是否正确使用该功能?
>在报告/绘制均值时,绘制均值(DATA)是否正确(参见第1行)或者我应该绘制自举程序导出的均值?我看到STAT包含每个bootstrap示例的平均值,但我不知道是否应该使用此信息,以及如何使用
>以我的方式绘制置信区间是否正确(见第[2]行),或者我应该绘制MEAN-CI(1,:)和MEAN CI(2,:)?
请查找附带代码生成的图表.