% New data
smallpath_rss = smallpath(:,3);
smallpath_accuracy = smallpath(:,4);

% Divide newly collected data into the correct squares
smallpath_squares = {...
  smallpath_rss(1:19) ...  
  smallpath_rss(20:29) ... 
  smallpath_rss(30:40) ... 
  smallpath_rss(41:51) ... 
  smallpath_rss(52:61) ... 
  smallpath_rss(62:71) ... 
  smallpath_rss(72:81) ... 
  smallpath_rss(82:92) ... 
  smallpath_rss(93:102) ... 
  smallpath_rss(103:115) ... 
  smallpath_rss(114:123) ... 
  smallpath_rss(124:134) ... 
  smallpath_rss(135:143) ... 
  smallpath_rss(144:154) ... 
  smallpath_rss(155:164) ... 
  smallpath_rss(155:174) ... 
  smallpath_rss(175:184) ... 
  smallpath_rss(185:194) ... 
  smallpath_rss(195:204) ... 
  smallpath_rss(205:214) ... 
  smallpath_rss(215:224) ... 
  smallpath_rss(225:234) ... 
  smallpath_rss(235:244) ... 
  smallpath_rss(245:254) ...
};

smallpath_squares_size = [...
  length(smallpath_rss(1:19)) ...  
  length(smallpath_rss(20:29)) ... 
  length(smallpath_rss(30:40)) ... 
  length(smallpath_rss(41:51)) ... 
  length(smallpath_rss(52:61)) ... 
  length(smallpath_rss(62:71)) ... 
  length(smallpath_rss(72:81)) ... 
  length(smallpath_rss(82:92)) ... 
  length(smallpath_rss(93:102)) ... 
  length(smallpath_rss(103:115)) ... 
  length(smallpath_rss(114:123)) ... 
  length(smallpath_rss(124:134)) ... 
  length(smallpath_rss(135:143)) ... 
  length(smallpath_rss(144:154)) ... 
  length(smallpath_rss(155:164)) ... 
  length(smallpath_rss(155:174)) ... 
  length(smallpath_rss(175:184)) ... 
  length(smallpath_rss(185:194)) ... 
  length(smallpath_rss(195:204)) ... 
  length(smallpath_rss(205:214)) ... 
  length(smallpath_rss(215:224)) ... 
  length(smallpath_rss(225:234)) ... 
  length(smallpath_rss(235:244)) ... 
  length(smallpath_rss(245:254)) ...
];

% Our collected data
collected_accuracy = all_in_squares_path_acc;
collected_rss = all_in_squares_path_rss;

high_square_rss = high_square_all(:,1);
% high_square_acc = high_square_all(:,2);
% 
max_square_rss = max_square_all(:,1);
% max_square_acc = max_square_all(:,2);


figure(1);
plot(collected_rss,collected_accuracy,'r*');
hold on;
plot(smallpath_rss,smallpath_accuracy,'bo');
hold off;

xlabel('RSS [dBm]');
ylabel('Precision [m]');
legend('Observations in signal map','New data');


max_square = weighted_cells{(cell_max_row-1)*174+cell_max_column}(:,3);
high_square = weighted_cells{(55-1)*174+143}(:,3);

%figure(2);
% counter = 0;
% for x=1:24
%     for i=1:length(smallpath_squares{x})
%         counter = counter + 1;
%         pos(i)= counter;
%     end          
% end
% 
% figure(3);
% nr_of_observations_in_squares = [];
% ppmcc_for_squares = [];
% rss_for_a_square = [];
% rss_in_path = [];
% for x=1:24
%    q = length(per_squares_path_rss{1,x}(:,1));
%    nr_of_observations_in_squares(x) = q;
%    rss_for_a_square = per_squares_path_rss{1,x};
%    rss_in_path = smallpath_squares{1,x};
%    t = length(rss_in_path);
%    % ppmcc for points with higher amount of points
%    if q > 2 && t > 2  
%        if nr_of_observations_in_squares(x) >= t
%             rss_for_a_square = rss_for_a_square(1:t);
%        else
%            rss_in_path = rss_in_path(1:q);
%        end
%        [r,p] = corr(rss_for_a_square,rss_in_path);
%        ppmcc_for_squares(x) = r;
%    else
%        ppmcc_for_squares(x) = 0;
%    end
% end
% bar(nr_of_observations_in_squares);
% xlabel('Square');
% ylabel('Observations');


figure(4);
rss_for_a_square = [];
rss_in_path = [];
hold on;
for x=1:24
   rss_for_a_square = per_squares_path_rss{1,x};
   rss_in_path = smallpath_squares{1,x};
   in_path_x = [];
   a_square_x = [];
   for y=1:length(rss_for_a_square)
       a_square_x(y) = x;
   end
   for i=1:length(rss_in_path)
       in_path_x(i) = x;
   end
   plot(a_square_x,rss_for_a_square,'bx');
   plot(in_path_x,rss_in_path,'rs');
end

hold off;
xlabel('Square');
ylabel('RSS [dBm]');
legend('Signal Map', 'New data');



square_labels = 1:1:24;
figure(5);
ax_square = subplot(1,1,1);
bar(ax_square,smallpath_squares_size,'r');
set(gca,'xticklabel',square_labels) 
xlabel('Square');
ylabel('Observations');


% 
% figure(2);
% plot(max_square_rss,max_square_acc,'yx');
% hold on;
% plot(high_square_rss,high_square_acc,'ps');
% hold off;
% 
% xlabel('RSS [dBm]');
% ylabel('Precision [m]');