
load('spe_d1r0c0sAlln');

x=mean(x2022n(:,1:40)); %x=mean(x1pd2022n);
y=mean(x2022nx(:,1:40));

mdlpred = fitlm(x,y); %returns a linear regression model of the responses y, fit to the data matrix X.
anova(mdlpred,'summary'); %test fit by anova
figure;plot(mdlpred);title('mdle2 pred')

[totalClark, percentClark] = clarke(x,y);
dif = x-y;
figure; plot(dif)

[N M] = size(x);
m=[1:N];
m=m';
figure; plot(m,x,'-ro',m,y,'-bo')