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