%Normal, vår 09 kvinner gj.snitt. Konfidens intervall

figure
AntallUtval=200;
niutval=36;
mu=169.4;
sigma=6.0;
subplot(1,3,2);
plot([mu,mu],[1,AntallUtval],'r-', 'linewidt',3);
hold
title('KI og sann \mu', 'FontSize',14);
alfa=0.05;
z_alfahalve=norminv(1-alfa/2, 0,1);
subplot(1,3,3);
plot([1,AntallUtval],[alfa,alfa],'m-', 'linewidt',3);
hold
title('Andel KI utanfor sann \mu', 'FontSize',14)
AntallUtanfor=0;
KI=zeros(AntallUtval,2);
AndelUtanfor=zeros(AntallUtval,1);
pauseint=0.05;
for(i=1:AntallUtval)	
	utval=normrnd(mu,sigma, niutval,1);
	gjsnitt(i)=mean(utval);
	subplot(1,3,1);	
	hist(utval);
	title('Eit tilfeldig utval', 'FontSize',14)	
	KI(i,1)=gjsnitt(i)-z_alfahalve*sigma/sqrt(niutval);
	KI(i,2)=gjsnitt(i)+z_alfahalve*sigma/sqrt(niutval);
	subplot(1,3,2);
	plot(KI(1:i,:)',1:1:i,'*')
	if(KI(i,1)>mu || KI(i,2)<mu)
		AntallUtanfor=AntallUtanfor+1;
	end
	subplot(1,3,3);
	AndelUtanfor(i)=AntallUtanfor/i;
	plot(AndelUtanfor(1:i));
	pause(pauseint)
	if(i==1)
		pause
	end
end


%Rask
AntallUtval=10000;
KI=zeros(AntallUtval,2);
AndelUtanfor=zeros(AntallUtval,1);
utvalM=normrnd(mu,sigma, niutval,AntallUtval);
gjsnitt=mean(utvalM);
KI(:,1)=gjsnitt-z_alfahalve*sigma/sqrt(niutval);
KI(:,2)=gjsnitt+z_alfahalve*sigma/sqrt(niutval);
Utanfor= KI(:,1)>mu | KI(:,2)<mu;
for(i=1:AntallUtval)
	AndelUtanfor(i)=sum(Utanfor(1:i))/i;
end
figure
plot([1,AntallUtval],[alfa,alfa],'m-', 'linewidt',3);
hold
plot(AndelUtanfor);
xlabel('Antall utval', 'FontSize',16)
ylabel('Andel KI utanfor sann \mu', 'FontSize',16)