MATLAB kommandoer til
eksempel C i avsnitt 11.2.1 i Rice
% F?rst leser vi inn dataene i tabellen p? side 429.
% (Kommandoen forutsetter at filen
'iron.txt' er
% lagret i en mappe der MATLAB kan finne den.)
iron=load('iron.txt')?
% Som i Rice vil vi n?ye
oss med ? se p? dataene for
% konsentrasjone 1.2 millimolar. Vi lagrer dataene for
% denne konsentrasjonen for 2-verdig jern i vektoren
X
% og dataene for 3-verdig jern i vektoren Y
X=iron(:,5);
Y=iron(:,2);
% Vi lager s? et boksplott og et normalfordelingsplott
% av dataene (jf. figurene 11.2 og
11.3 p? side 430 i Rice,
% men merk at normalfordelingsplottet i Rice og det MATLAB
% lager bytter om x- og y-aksene)
boxplot([X,Y])
probplot([X,Y])
% Endelig tester vi om det er forskjell
% p? metodene og beregner 95 % konfidensintervall
[h,p,ci,stats]=ttest2(X,Y)
% Studentene kan selv pr?ve tilsvarende kommandoer
% p? de log-transformerte
observasjonene
XX=log(X);
YY=log(Y);