2009-10-21 2 views

ответ

5

Я не уверен, что вы имели в виду доверительных интервалов графика, но это пример того, как построить двойку односторонний 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') 

plot

+0

Спасибо Amro ....... но мне было интересно, как я могу получить график, как показано на рисунке по следующей ссылке: http://psycnet.apa.org/journals/bne/115/4/ images/bne_115_4_764_fig5a.gif Я хотел иметь бары для вариаций в графике, т. е. если на графике есть большая вариация, полоса больше, а если изменение ниже, то график должен иметь меньшие бары. сделал u получить то, что я хочу сказать? с нетерпением жду .... – Haqan

+1

Я верю, что вы запрашиваете функцию 'errorbar': http://www.mathworks.com/access/helpdesk/help/techdoc/ref/errorbar.html – Amro

+0

Спасибо alooooot .......! я точно нашел этот сюжет .... – Haqan

3

См., Например, эти м-файлы на Matlab файлообменный:

+0

благодарит RCS !! Как я начинаю в Matlab, можете ли вы рассказать мне, как я могу запустить этот файл m? coz, я вижу ошибку? Должен ли я указывать некоторые входные параметры или файл, чтобы иметь сюжет? – Haqan

1

После прочтения многочисленных потоков, вот моя попытка.

enter image description here

% 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

Смежные вопросы