%maeusestatistik
%plottet Mittelwert und Standardfehler 
%in der ersten Abbildung fuer Messungen von Tag zu Tag
%und in der zweite Abbildung fuer Messungen von Tier zu Tier

%load pellets
load pellets2
pellets=pellets2;

%Groesse bestimmen:
[tiere,tage]=size(pellets);

%mittelung und Standardabweichung ueber die Tiere jeweils eines Tages:
m1=mean(pellets);
abw1=std(pellets);

%mittelung ueber die Tage jeweils eines Tieres:
m2=mean(pellets');
abw2=std(pellets');


%Standardfehler in beide Richtungen:
fehler_ueber_Tage=standardfehler(tiere,m1)
fehler_ueber_Tiere=standardfehler(tage,m2)


%Fuer gleiche Skalierung der Abbildungen Minima und Maxima bestimmen:
mini=min([m1 m2])-max([abw1 abw2]);
maxi=max([m1 m2])+max([abw1 abw2]);

%Abbildungen
figure(1)
errorbar(m1,abw1)
ylabel('Anzahl gefressener Pellets')
xlabel('Tag')
title('Mittelwert und Standardabweichung von Messungen einzelner Tage')
axis([1 tage mini maxi])

figure(2)
errorbar(m2,abw2)
ylabel('Anzahl gefressener Pellets')
xlabel('Maus')
title('Mittelwert und Standardabweichung von Messungen einzelner Tiere')
axis([1 tiere mini maxi])

