如何在MATLAB中绘制置信区间?

o7jaxewo  于 2023-10-23  发布在  Matlab
关注(0)|答案(3)|浏览(269)

我想在MATLAB中绘制一些置信区间图,但我根本不知道如何做到这一点。我有一个.xls文件中的数据。
谁能给予点提示,或者谁知道绘制CI的命令?

vlf7wbxs

vlf7wbxs1#

在阅读了大量的帖子后,这是我的尝试。

% Get some random data
x       = linspace(0.3, pi-0.3, 10);
Data    = sin(x) + randn(1, 10)/10;
Data_sd = 0.1+randn(1,10)/30;

% prepare it for the fill function
x_ax    = 1:10;
X_plot  = [x_ax, fliplr(x_ax)];
Y_plot  = [Data-1.96.*Data_sd, fliplr(Data+1.96.*Data_sd)];

% plot a line + confidence bands
hold on 
plot(x_ax, Data, 'blue', 'LineWidth', 1.2)
fill(X_plot, Y_plot , 1,....
        'facecolor','blue', ...
        'edgecolor','none', ...
        'facealpha', 0.3);
hold off

主要基于这个问题:Plotting with transparency

xxhby3vn

xxhby3vn2#

我不知道你说的置信区间图是什么意思,但这是一个如何绘制正态分布双侧95%CI的示例:

alpha = 0.05;          % significance level
mu = 10;               % mean
sigma = 2;             % std
cutoff1 = norminv(alpha, mu, sigma);
cutoff2 = norminv(1-alpha, mu, sigma);
x = [linspace(mu-4*sigma,cutoff1), ...
    linspace(cutoff1,cutoff2), ...
    linspace(cutoff2,mu+4*sigma)];
y = normpdf(x, mu, sigma);
plot(x,y)

xlo = [x(x<=cutoff1) cutoff1];
ylo = [y(x<=cutoff1) 0];
patch(xlo, ylo, 'b')

xhi = [cutoff2 x(x>=cutoff2)];
yhi = [0 y(x>=cutoff2)];
patch(xhi, yhi, 'b')

6kkfgxo0

6kkfgxo03#

例如,参见Matlab文件交换中的这些m文件:

相关问题