Confidence-Interval

從匯總統計中繪製置信區間條

  • July 19, 2013

有點像箱線圖。我的意思不一定是標準的上置信區間、下置信區間、均值和顯示數據范圍的箱線圖,但我的意思是只有三個數據的箱線圖:95% 置信區間和均值

這是一篇期刊文章的截圖,它正是我想要的:在此處輸入圖像描述

我還想知道如何使用回答者提到的軟件來創建這樣的情節。

在 MATLAB 中,您可能想嘗試errorbar函數: http: //www.mathworks.de/de/help/matlab/ref/errorbar.html

或者,您可以以愚蠢和手動的方式進行操作。例如,給定一個數據點矩陣“a”,您可以使用函數 m = mean(a) 計算您的均值,計算您的 CI(取決於您需要的 CI),然後手動繪製結果。

演示您是否已經知道均值和 CI,假設 CI 在矩陣CI(第一列和第二列)中並且均值在矩陣a中:

plot(1:length(CI),a,'o','markersize', 10)           % plot the mean
hold on;
plot(1:length(CI),CI(1,:),'v','markersize', 6)              % plot lower CI boundary
hold on;
plot(1:length(CI),CI(2,:),'^','markersize', 6)              % plot upper CI boundary
hold on;

for I = 1:length(CI)                                        % connect upper and lower bound with a line
line([I I],[CI(1,I) CI(2,I)])
hold on;
end;

axis([0 length(CI)+1 min(CI(1,:))*0.75 max(CI(2,:))*1.25])  % scale axis

在您知道單個測量值的情況下進行演示,對於重複測量實驗,3+ 個條件,每列一個條件,矩陣 a 中每行一個對象,沒有缺失樣本,95% CI,如 MATLAB 的ttest() 所示

[H,P,CI] = ttest(a);                                        % calculate 95% CIs for every column in matrix a
                                                           % CIs are now in the matrix CI!

plot(1:length(CI),[mean(a)],'o','markersize', 10)           % plot the mean
hold on;
plot(1:length(CI),CI(1,:),'v','markersize', 6)              % plot lower CI boundary
hold on;
plot(1:length(CI),CI(2,:),'^','markersize', 6)              % plot upper CI boundary
hold on;

for I = 1:length(CI)                                        % connect upper and lower bound with a line
line([I I],[CI(1,I) CI(2,I)])
hold on;
end;

axis([0 length(CI)+1 min(CI(1,:))*0.75 max(CI(2,:))*1.25])  % scale axis

引用自:https://stats.stackexchange.com/questions/64880

comments powered by Disqus