mean_fitness_distribution.m

clear all

% plots fitness or N distribution for each time step

init_file=15000; %365 files = 365 days
num_files=22000; %79027; % fitness shifted1overd

for j=init_file:1:num_files % number of lines in the file "fitness"

filename=['line',int2str(j)]

a=load(filename);

%time(j)=a(1);

s=1;

clear nonzero_fitness

for i=3:1:length(a)
    if(a(i)~=0)    
        nonzero_fitness(s)=a(i);
        s=s+1;
    end
end

most_fit(j)=max( nonzero_fitness );
least_fit(j)=min( nonzero_fitness );

end

max_fitness=max(most_fit);
min_fitness=min(least_fit);

bins=min_fitness:(max_fitness-min_fitness)/100:max_fitness;
n=zeros(1,length(bins));

for j=init_file:1:num_files % number of lines in the file "fitness"

filename=['line',int2str(j)]

a=load(filename);

%time(j)=a(1);

s=1;

clear nonzero_fitness

for i=3:1:length(a)
    if(a(i)~=0)    
        nonzero_fitness(s)=a(i);
        s=s+1;
    end
end

n=n+hist(nonzero_fitness,bins);

%title(sprintf('R_0 = %g, 1/nu = %g days, i(0) = %g',R0,1/nu*365,i0))

%pause(0.25)

end

h=figure(1)
bar(bins,n);
title('fitness distribution averaged over 19 years');
xlabel('fitness values');
ylabel('counts');

%saveas(h,['mean_fitness_distribution','.eps']);
%saveas(h,['mean_fitness_distribution','.fig']);