平均と標準偏差のみからの信頼区間を設定する

Question

データの平均と標準偏差（std）から信頼区間をプロットしようとしています。

MeanA=1.876; %mean of A 
STDA=0.018; % std of A 
MeanB=1.821; 
STDB=0.039; 

MeanC=1.735; 
STDC=0.023; 
MeanD=1.667; 
STDD=0.039; 
Y = [MeanA MeanB ; 
     MeanC MeanD ];  
errY= [STDA STDB; 
     STDC STDC ]; 

私は正規分布をプロットした場合、その後彼らの自信のINTEVAL私は重なりが信頼区間（CI）であったプロットをobtined

alpha = 0.05;   % significance level 


for tt=1:length(Y) 
figure 
mu = Y(tt,1);    % mean 
sigma = errY(tt,1); 
cutoff1n = norminv(alpha, mu, sigma); 
cutoff2n = norminv(1-alpha, mu, sigma); 
xn = [linspace(mu-4*sigma,cutoff1n), ... 
linspace(cutoff1n,cutoff2n), ... 
linspace(cutoff2n,mu+4*sigma)]; 
yn = normpdf(xn, mu, sigma); 
plot(xn,yn) 
mu = Y(tt,2);    % mean 
sigma = errY(tt,2); 
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); 
hold on, plot(x,y) 
plot(x,y,'r-',xn,yn,'g-' , 'LineWidth',3) 
set(gca,'Fontsize', 32) 
if tt==1 
hleg1=legend('A', 'B'); 
title('study1') 
elseif tt==2 
hleg1=legend('C', 'D'); 
title('Study2') 
end 
set(hleg1,'Location','NorthEastoutside') 
set(gca,'Fontsize',22) 
xlo = [x(x<=cutoff1) cutoff1]; 
ylo = [y(x<=cutoff1) 0]; 
patch(xlo, ylo, 'r') 
xhi = [cutoff2 x(x>=cutoff2)]; 
yhi = [0 y(x>=cutoff2)]; 
patch(xhi, yhi, 'r') 
xlon = [xn(xn<=cutoff1n) cutoff1n]; 
ylon = [yn(xn<=cutoff1n) 0]; 
    patch(xlon, ylon, 'g') 
    xhin = [cutoff2n xn(xn>=cutoff2n)]; 
    yhin = [0 yn(xn>=cutoff2n)]; 
    patch(xhin, yhin, 'g') 
end 

に重なっているようだ：ここで私が書いたコードの一部であります。

今私はの形でCIをプロットする必要があります。

誰かが周りのCIをプロットするために助けることができる、私は方法を以下にCIを計算しています

%%

を意味する：

SE2=errY/sqrt(10); 
    CI2n=Y-1.96*(SE2); 
    CI2p=Y+1.96*(SE2); 

場合はその正しい方法を教えてくださいとはい、どうすればそれらをプロットすることができます。おかげ

Answer 1

標準偏差sigmaが知られており、サンプルサイズがnが

CI = mu +/- 1.96*sigma/sqrt(n) 

mu = Yとそう

とsigma = errY、95％信頼のあなたの計算式で与えられている平均muの95％信頼区間サンプルサイズがn = 10の場合、間隔は正しいです。

n = 10; 
SE2=errY/sqrt(n); 
for i = 1:size(SE2,1) 
    figure 
    y = Y(i,:); 
    x = 1:length(y); 
    e = 1.96*SE2(i,:); 
    errorbar(1:length(y),y,e,e,'k.','markersize',18); 
    if i == 1 
     c = ['A';'B']; 
    elseif i == 2 
     c = ['C';'D']; 
    end 
    set(gca,'xtick',x,'xticklabels',c) 
end 

this imageを生成し（私は画像をインラインを入れたいのですが、十分な担当者が不足している）：エラーバー、MATLABのerrorbar機能（documentation）を使用して、以下を試して、これらの値をプロットし

。

信頼区間が重複して表示されないことがあります。しかし、n = 1、つまりCI = mu +/- 1.96*sigmaを設定することで、濃度プロットの結果を得ることができます。私は2番目の画像を投稿しますが、私はより多くの担当者が必要です！