summaryrefslogtreecommitdiff
path: root/291/CH10/EX10.4b/eg10_4b.sce
blob: 10372d0444f37cf134663fb53c3c27a174d11d5f (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
X=[75 73 60 70 86; 78 71 64 72 90; 80 69 62 70 85; 73 67 63 80 92 ];
Xidot = zeros(4,1);
for i=1:4
    for j=1:5
        Xidot(i)=Xidot(i) + X(i,j);
    end
end
Xidot = Xidot/5;
Xjdot = zeros(5,1);
for j=1:5
    for i=1:4
        Xjdot(j)=Xjdot(j) + X(i,j);
    end
end
Xjdot = Xjdot/4;
Xdotdot = sum(Xidot)/4;
//disp(Xdotdot)
meanhat = Xdotdot;
alphahat = Xidot - meanhat;
betahat = Xjdot - meanhat; 
disp(meanhat, "The estimator of the mean is");
disp("The alphas are-")
disp(alphahat)
disp("The betas are-")
disp(betahat)