Discussion :

[Eco_prince]

Bonjour,

J'aimerais faire un test anova sur mes 4 échantillons de tailles différentes (Je ne peux donc pas faire de matrice avec chaque variable dans une colonne, comme la prévoit la fonction arima1 de base..)
Le problème, c'est qu'ils ont des tailles assez conséquentes (15'900, 5300, 21'000 et 3800 observations).

Si j'ai bien compris la fonction anova1(X,group), je dois spécifier pour chaque observations dans quel groupe elle va... Et donc dans mon vecteur group, je devrais spécifier 45'000 fois si c'est le groupe 1, le 2,...

Cela risque de me prendre quelques années pour spécifier cela!

Je suis sur qu'il y a un moyen plus rationnel de faire, je ne dois pas être le premier avec ce type de problème!

Merci d'avance pour vos propositions/aides/...

[Jerome Briot]

Je ne connais pas grand chose à ces fonctions mais si je suis bien ce qui est dis dans la documentation, tu devrais pouvoir t'en sortir comme ceci :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
X1 = [...]; % Echantillon 1
X2 = [...]; % Echantillon 2
X3 = [...]; % Echantillon 3
X4 = [...]; % Echantillon 4
 
X = [X1(:);X2(:);X3(:);X4(:)];
 
group = [repmat(1,numel(X1(:)),1)
         repmat(2,numel(X2(:)),1)
         repmat(3,numel(X3(:)),1)
         repmat(4,numel(X4(:)),1)];
 
p = anova1(X,group);

C'est l'idée...

[Eco_prince]

Merci bien pour la réponse...

Je crois que c'est à peu près la même idée que ce que je suis entrain de faire:

Avec une boucle, je crée un vecteur de dimension n pour chaque variable:

for i=1:size(X_1)
X_1_group(i,1)={'1er groupe'};
end

après j'ai mis les 4 vecteurs ainsi trouvé ensemble dans le vecteur colonne group. Et mes X_1,..., X_4 dans un vecteur colonne stats.
Il faut que je controle encore mais il me semble que ça marche!

Merci quand même pour la réponse (peut-être plus pro que mes 4 boucles et mon bricolage!)