%naive model
clc
clear all
%load('losangels.mat')
load('3y14-5y29.mat')
for j=1:7
for i=0:1
b(i+1)=a7510(288,42+7*i+j-1);
c(i+1)=a7510(1,43+7*i+j-1);
end
d=a7510(288,56+j-1);
pre(1,j)=(d/sum(b))*sum(c);
for j1=2:288
for i=0:1
b(i+1)=a7510(j1-1,43+7*i+j-1);
c(i+1)=a7510(j1,44+7*i+j-1);
end
d=a7510(j1-1,57+j-1);
pre(j1,j)=(d/sum(b))*sum(c);
end
end