%diameter ved bryst høyde (i fot) og alder (antall ˚ar) av 26 kastanjetrær
fileID = fopen('DBHTREES.DAT','r');
formatSpec = '%d %f';
d = fscanf(fileID,formatSpec,[2 26]);
fclose(fileID);
[p n] = size(d);
corrmat = corr(d')
rhohat = corrmat(1,2)
B = 10000;
rhosim = zeros(1,B);
B = 10000
for b = 1:B
    id = randsample(1:n,n,true);
    dstar = d(1:2,id);
    corrmat = corr(dstar');
    rhosim(b) = corrmat(1,2);
end
hist(rhosim)
%Forventningsskjevhet
biasrho = mean(rhosim) - rhohat
%Varians
Vrho = var(rhosim)
%Standardfeil
srho = sqrt(Vrho)
%95% konfidensinvervall
KIrho = [rhohat-1.96*srho rhohat+1.96*srho]


%phi=log(1+rho)-log(1-rho)
phihat = log(1+rhohat)-log(1-rhohat)
phisim = log(1+rhosim)-log(1-rhosim)
hist(phisim)
biasphi = mean(phisim)-phihat
Vphi = var(phisim)
sphi = sqrt(Vphi)
KIphi = [phihat-1.96*sphi phihat+1.96*sphi]


KIrho2 = (exp(KIphi)-1)./(exp(KIphi)+1)












%Regresjonslinje
plot(d(1,1:26),d(2,1:26),'*')
%Tilpasning "best" mulig linje
[r,bhat,ahat] = regression(d(1,1:26),d(2,1:26));
refline(bhat,ahat);
B = 10000
bsim = zeros(1,B);
for b = 1:B
    id = randsample(1:n,n,true);
    dstar = d(1:2,id);
    [r,bstar,astar] = regression(dstar(1,1:26),dstar(2,1:26));
    bsim(b) = bstar;
end
hist(bsim)
%Forventningsskjevhet
biasb = mean(bsim) - bhat
%Varians
Vb = var(bsim)
%Standardfeil
sb = sqrt(Vb)
%95% konfidensinvervall
KIb = [bhat-1.96*sb bhat+1.96*sb]
