bar(Kvec)
hold on;
plot([0, 9], [mean(Kvec), mean(Kvec)]);
plot([0, 9], [mean(Kvec)-std(Kvec), mean(Kvec)-std(Kvec)], 'm--');
plot([0, 9], [mean(Kvec)-2*std(Kvec), mean(Kvec)-2*std(Kvec)], 'r:', 'LineWidth', 2);
plot([0, 9], [mean(Kvec)+std(Kvec), mean(Kvec)+std(Kvec)], 'm--');
plot([0, 9], [mean(Kvec)+2*std(Kvec), mean(Kvec)+2*std(Kvec)], 'r:', 'LineWidth', 2);
axis([0, 9, (min(Kvec)-5), (max(Kvec)+5)]);
xlabel('Index of data', 'FontSize', 18);
ylabel('Estimation of K', 'FontSize', 18);
legend('K', 'mean', 'mean-sd', 'mean-2sd');
orient landscape;
print('-dpng', 'Ksbarchart.png');
hold off;