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