[fh_country]=textread('fh_countries','%q');
load fh_timedata;
n=size(fh_timedata,2);

pr=fh_timedata(:,1:3:n);
cl=fh_timedata(:,2:3:n);
t=1972:2002;

clf

for i=1:size(pr,1)
  for j=1:size(pr,2)
    if (pr(i,j)==-1) pr(i,j)=NaN; end
    if (cl(i,j)==-1) cl(i,j)=NaN; end
  end
end

dp=diff(pr')';
dc=diff(cl')';
sum1=zeros(1,size(dp,2));
sum2=zeros(1,size(dp,2));
ngov=zeros(1,size(dp,2));
for i=1:size(dp,1)
  for j=1:size(dp,2)
    if (~isnan(dp(i,j))) 
      ngov(j)=ngov(j)+1;
      sum1(j)=sum1(j)-(dp(i,j)); 
      sum2(j)=sum2(j)-(dc(i,j)); 
    end
  end
end
clf
plot(1+(1972:2001),sum1./ngov,'r')
hold on
plot(1+(1972:2001),sum2./ngov,'g')
plot(1+(1972:2001),sum2*0,'b:')
legend('Political freedom','Civil liberties');
	xlabel('Year');
	ylabel('Sum change from previous year per country')