Как построить доверительные интервалы в MATLAB?
-
05-07-2019 - |
Вопрос
Я хочу построить графики доверительных интервалов в MATLAB, но понятия не имею, как это сделать.У меня есть данные в файле .xls.
Может ли кто-нибудь дать мне подсказку, или кто-нибудь знает команды для построения графиков CI?
Решение
Я не уверен, что вы подразумевали под графиком доверительных интервалов, но это пример того, как построить двусторонний 95% -й доверительный интервал нормального распределения:
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')
Другие советы
После прочтения многочисленных тем, вот моя попытка. Р>
% 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
В основном на основе этого вопроса: Печать с прозрачностью
См., например.эти m-файлы в Matlab File Exchange:
Не связан с StackOverflow