You can use for loops to do the sum or use the sum function:
clc; clear; close all force;
M = 1; % Test day
N = 14; % Number of runs (or repetitions)
T = 101; % Number of time points
Y = rand(M,N,T); % random data to illustrate
% using for loops
Yi_for = zeros(M,1);
for i=1:M
for j=1:N
for t=1:T
Yi_for(i) = Yi_for(i) + Y(i,j,t);
end
end
Yi_for(i) = Yi_for(i)/N/T;
end
% using sums
Yi_sum = zeros(M,1);
for i=1:M
Yi_sum(i) = sum(sum(Y,2),3)/N/T;
end
You can expand this to calculate the average at time point t on the ith test day:
% Yit should be a matrix
Yit = zeros(M,T);
for i=1:M
for t=1:T
for j=1:N
Yit(i,t) = Yit(i,t)+Y(i,j,t);
end
Yit(i,t) = Yit(i,t)/N;
end
end
With Yi and Yit previously calculated, you can define Ra (I splitted the calculations of the numerator and denominator of the fraction to make it clearer):
% Ra
Yi = Yi_sum;
Ra = 1;
% partial sums
num = 0; % numerator
for i=1:M
for j=1:N
for t=1:T
num = num + (Y(i,j,t)-Yit(i,t))^2/M/T/(N-1);
end
end
end
den = 0; % denominator
for i=1:M
for j=1:N
for t=1:T
den = den + (Y(i,j,t)-Yi(i))^2/M/(N*T-1);
end
end
end
Ra2 = 1-num/den
