Discussion :

[Anna76]

Bonjour,

J'ai une nouvelle fois besoin d'aide...
Le titre n'est pas très explicite donc je vais détaillé un peu plus !

J'ai en fait un fichier de données que je charge dans matlab.
Ce que je souhaite faire c'est pour chaque ligne, calculer la moyenne et l'écart type de celle-ci (jusque là pas de souci!) puis recalculer ces 2 paramètres mais en enlevant 1 donnée à chaque fois :
c'est à dire si j'ai 1 2 3 4 5 en première ligne, je veux ensuite calculer la moyenne et l'écart type de 2 3 4 5 puis de 1 3 4 5 puis de 1 2 4 5 etc...!
J'espère que je suis assez claire ^^!
Je me demande donc s'il existe un moyen de calculer l'écart type (ou la variance) en "sautant" une colonne à chaque fois.
Pour la moyenne je me suis débrouiller en bidouillant : multiplier par n soustraire la valeur puis rediviser par (n-1) mais je ne peux pas faire cela pour l'écart type!

Je précise que le but est de faire un test statistique pour comparer l'ensemble de la ligne, à chaque échantillon!

Merci par avance de votre aide

[shaiHulud]

Bonjour,

Une fois que tu as une matrice X ou chaque ligne est une variable, moyenne et variance sont donnés par

Code :

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

mean(X'), var(X')

Si tu veux enlever la 1ere colonne, il suffit de

Code :

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

mean(X(:,2:end)'), var(X(:,2:end)')

[Anna76]

Merci de cette réponse on ne peut plus rapide !

Pour la première colonne je sais effectivement qu'on peut faire comme cela mais par exemple si on veut calculer la variance en excluant uniquement la 2e colonne ou l'une des suivantes je ne vois pas comment procéder car

Code :

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

var(X(:,3:end)')

ne calculera que la variance des dernières colonnes et donc exclura 2 valeurs au lieu d'une.

Je ne sais pas si je me fais bien comprendre...

Au passage, j'ai des 0 dans certaines cases que je ne souhaite pas prendre dans le calcul donc c'est vrai que pour le moment, je calculais "manuellement" les moyennes et variances car je ne voyais pas non plus comment calculer cela en excluant les 0.

[Gooby]

Dans ce cas tu peux faire un truc de ce genre:

Code :

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

var(cat(1,X(1:i),X(i+1:end))

Avec une boucle sur le i, qui correspondra à la case qui va sauter.

Et pour les 0, tu peux trouver les indices où ils apparaissent en effectuant un find, puis de la même façon que pour exclure chacune des valeurs, tu peux exclure les 0 également.

[shaiHulud]

le cat des valeurs sera un poil plus gourmand que de cat les indices

Code :

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

1
2
3
4
 
var(X(:,setdiff(1:size(X,2),2))');% enleve la colonne 2
% ou
var(X(:,[1:5, 7:end])') % enleve la colonne 6

en ce qui concerne les 0, tu peux les remplacer par des nan et utiliser nanmean et nanstd qui excluent ces valeurs nan du calcul

Code :

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

1
2
3
 
X(abs(X)<eps) = nan;
nanstd(X(:,[1:5, 7:end])')  .^2 % variances en enlevant colonne 6 ET en excluant les 0

[Anna76]

Merci à vous deux pour vos idées !
c'est exactement ce qu'il me faut!

Pour les 0, j'avais déjà pensé au NaN mais malheureusement je n'ai pas la toolbox pour utiliser ces fonctions .
Je vais donc essayé avec un find !

Encore merci pour vos réponses aussi rapides !