Metodi di ricampionamento: jacknife, bootstrap e test di ipotesi

t=[94 197 16 38 99 141
23];
c=[52 104 146 10 50 …
31 40 27 46];
[t c]
dati=[t c];
Delta=mean(t)-mean(c);
Codice
% ricampionamento
for i=1:100,
dati_r=dati(randperm(length(dati(1,:))));
delta(i)=mean(dati_r(1:length(t)))-mean(dati_r(length(t)+1:end));
end
hist(delta)
perc=length(find(delta>Delta))/100
for i=1:100,
bootstrap=dati(randsample(length(dati(1,:)),length(dati(1,:)),true))
dati_t=bootstrap(1:length(t))
dati_c=bootstrap(length(t)+1:end)
deltab(i)=mean(dati_t)-mean(dati_c);
end
figure
hist(deltab)
percb=length(find(deltab>Delta))/100