Commit 279f2332 authored by Turnhout, M.C. van's avatar Turnhout, M.C. van
Browse files

ani kappa epsilon

parent ecfda613
This diff is collapsed.
clear all; close all;
% sp = load('paramspace/pspaceaniparam.txt');
sp = load('paramspace/pspaceaniparam.txt');
partex = {'\mu', '\sigma', 'P'};
sigma = linspace(0, 2.5e-2, 26);
meansp = zeros(3, 26);
medsp = meansp;
minsp = medsp;
maxsp = minsp;
for s = 1:100
sim = sp(s, 1);
l(s, 1) = 4;
results = load(sprintf('paramspace/pspaceani%i_results.txt', sim));
% projected length
l(s, 2:size(results, 1)+1) = results(:, 3);
end
for s = 1:26
meansp(:, s) = [mean(sp(l(s, :) > 0, 2));
mean(sp(l(s, :) > 0, 3));
mean(sp(l(s, :) > 0, 4))];
medsp(:, s) = [median(sp(l(s, :) > 0, 2));
median(sp(l(s, :) > 0, 3));
median(sp(l(s, :) > 0, 4))];
minsp(:, s) = [min(sp(l(s, :) > 0, 2));
min(sp(l(s, :) > 0, 3));
min(sp(l(s, :) > 0, 4))];
maxsp(:, s) = [max(sp(l(s, :) > 0, 2));
max(sp(l(s, :) > 0, 3));
max(sp(l(s, :) > 0, 4))];
end
figure
hold on
plot(sigma, meansp(1, :), 'r', sigma, meansp(2, :), 'y', sigma, meansp(3, :), 'g', 'linewidth', 2)
plot(sigma, medsp(1, :), 'r--', sigma, medsp(2, :), 'y--', sigma, medsp(3, :), 'g--', 'linewidth', 2)
h = legend('$\mean \mu$', '$\mean \sigma$', '$\mean P$', '$\median \mu$', '$\median \sigma$', '$\median P$');
set(h, 'box', 'off')
% patch([sigma, flipdim(sigma, 2)], [minsp(1, :) flipdim(maxsp(1, :), 2)], 'r', 'facealpha', 0.5)
% patch([sigma, flipdim(sigma, 2)], [minsp(2, :) flipdim(maxsp(2, :), 2)], 'y', 'facealpha', 0.5)
% patch([sigma, flipdim(sigma, 2)], [minsp(3, :) flipdim(maxsp(3, :), 2)], 'g', 'facealpha', 0.5)
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment