평균 및 표준 편차 만의 신뢰 구간

Question

데이터의 평균 및 표준 편차 (표준 편차)로부터 신뢰 구간을 계획하려고합니다. 나는 정규 분포를 플롯 경우

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 ]; 

가 그들의 cofidence의 inteval 내가 신뢰 구간에서 중복이 플롯을 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를 계산하고

말은 다른 사람의 도움이 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 인 간격이 정확합니다. this image (나는 이미지 인라인을 넣고 싶어하지만 충분한 담당자 부족) 생산

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 

:

는 MATLAB의 errorbar 기능 ( documentation)를 사용하여 다음을 시도, 오차 막대로이 값을 플롯합니다.

신뢰 구간이 중복되지 않는 것을 알 수 있습니다. 그러나, 우리는 n = 1, 즉 CI = mu +/- 1.96*sigma을 설정하여 밀도 플롯에서 결과를 얻을 수 있습니다. 두 번째 이미지를 게시 할 예정이지만 더 많은 담당자가 필요합니다!