xls = xlsread('data.xls'); time = 1:399; data = xls(time,1:22); oil = data(:,1); gold = data(:,2); iron = data(:,3); logs = data(:,4); % all_goods = [oil gold iron logs maize beef chicken gas liquid_gas tea tobacco wheat sugar soy silver rice platinum cotton copper coffee coal aluminum]; % , : ids = {'oil','gold','iron','logs','maize','beef','chicken','liquid_gas','gas','tea','tobacco','wheat','sugar','soy','silver','rice','platinum','cotton','copper','coffee','coal','aluminum'}; goods_count = size(all_goods, 2); geom_average = ones(size(time))'; %' for i = 1:goods_count geom_average = geom_average .* all_goods(:,i); end geom_average = geom_average .^ (1/goods_count); all_goods_rel = zeros(size(all_goods)); for i = 1:goods_count all_goods_rel(:,i) = all_goods(:,i) ./ geom_average; end
R = corrcoef(all_goods_rel);
% : threshold = 0.25; % 0.33 0.4 0.45 0.55 0.6 0.65 0.7 % : links = R>threshold; % , : bg = biograph(links, ids); view(bg);
Source: https://habr.com/ru/post/240323/