Statistikuppgift

result.csv

result.m

M=csvread('result.csv');
K=M(:,[4 6]);

n=8132;

stdXut=std(K(:,1))
stdYut=std(K(:,2))
meanXut=mean(K(:,1))
meanYut=mean(K(:,2))

X=log10(K(:,1));
Y=log10(K(:,2));

X(X == -Inf) = 0;
Y(Y == -Inf) = 0;

stdX=std(X)
stdY=std(Y)

a=sum(X);
b=sum(power(X, 2));
c=sum(Y);
d=sum(X.*Y);
e=sum(power(Y, 2));

b1=(n*d-a*c)/(n*b-power(a,2));
b0=mean(Y)-b1*mean(X);

Sxx=(n*b-power(a,2))/n;
Syy=(n*e-power(c,2))/n;
Sxy=(n*d-a*c)/n;

SSE=Syy-b1*Sxy;
s2=SSE/(n-2);
s=sqrt(s2);
t=b1/(s/sqrt(Sxx));

t005=2.326; % From table
t05=1.282;
confInt=[b1-t005*s/sqrt(Sxx) b1+t005*s/sqrt(Sxx)];
t005*s/sqrt(Sxx)

confInt2=[b0-t05*s*sqrt(b)/sqrt(n*Sxx) b0+t05*s*sqrt(b)/sqrt(n*Sxx)]