bland2007 AGREEMENT BETWEEN METHODS OF MEASUREMENT WITH MULTIPLE OBSERVATIONS PER INDIVIDUAL ---代码实现

% bland2007 AGREEMENT BETWEEN METHODS OF MEASUREMENT WITH MULTIPLE OBSERVATIONS PER INDIVIDUAL
clear, clc,close all;

load RV2IC.mat       % RV:反射性心室造影， IC:心阻抗法 ， 心室射血时间

%%%%%%%%%%%%section 3%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%ss = zeros(1,12);
N = 60;
k = 12;
m = zeros(1,12);
samples = nan(7,12);
for i = 1:N % 总样本数
    j = b(i,1);  % subject number
    m(j) = m(j) + 1;
    %ss(j) = ss(j) + (b(i,2) - b(i,3))^2;
    samples(m(j),j) = b(i,2)-b(i,3);
end
mus = mean(samples,'omitnan');     % 每个组各自均值    
s = samples - repmat(mus,7,1);     % 每个组减去各自的均值
SSE = sum(sum(s.^2,'omitnan'))   % 去掉各自均值后，平方，求和。 去掉了各个组的影响,组内离散度

gm = mean(samples(:),'omitnan')   % 总体均值

RSS = dot(m,(mus - gm).^2)       % 组间离散度
TSS = sum((samples(:) - gm).^2,'omitnan&