aa=load('..\Données\AMBR.txt');   % On charge les données
t=aa(:,1);                        % le temps est dans la 1e colonne
ns=aa(:,5);                       % la composante N-S dans la 5e
h=plot(t,ns,'k.','markersize',14);% On dessine les résultats 
xlabel('Temps (année)');
ylabel('Coordonnée N-S (par rapport à une origine arbitraire)')
A=[t'>1999+11/12 ; t'*0+1 ; t'];  % définition de la matrice normale
y=A'\ns                           % Résolution des moindres carrés
modele=A'*y;                      % Calcul du modele
plot(t,ns,'k.',t,modele,'r','linewidth',2); % et dessin
dfm=2;                            % Calcul de l'ANOVA
dfd=length(ns)-3
SSM=sum(modele.^2)
SSD=sum((ns-modele).^2)
MSM=SSM/dfm
MSD=SSD/dfd
F=MSM/MSD                         % et de la statistique de Fisher 
yf=fp(F,dfm,dfd)                  % et on compare à la table

