TMA4240 Statistikk, høsten 2016

MTBYGG, MTING, MTMT, MTTEKGEO

Matlabscript: Estimerer forventingsverdi ved gjennomsnitt

%
% Simulerer trekning av k kuler fra urne med totalt n kuler, 
% nR er røde og nG er gule. Trekker kuler gjentatte ganger og 
% estimerer forventningsverdier til X (antall røde kuler trukket ut)
% og Y (antall gule kuler trukket ut) ved å regne ut gjennomsnitt.
% Plotter resultatet som funksjon av antall forsøk utført.
%
 
 
n = 12;
nR = 3;
nG = 4;
k = 3;
 
m = 2000;
 
meanx = [];
meany = [];
sumx = 0;
sumy = 0;
for i=1:m
  draw = randperm(n,k);
  xx = sum(draw <= nR);
  yy = sum(draw > nR & draw <= nR + nG);	
  sumx = sumx + xx;
  sumy = sumy + yy;
  meanx = [meanx sumx/i];
  meany = [meany sumy/i];
 
  figure(1);
  plot(meanx);
  figure(2);
  plot(meany);
 
  if (i == 1)
    pause
  end
 
  pause(0.05);
end
 
[meanx(length(meanx)) meany(length(meany))]
2016-09-13, Håkon Tjelmeland