Ajouter des barres d'erreur a mon histogramme

**keke0** · 30/07/2013, 23h22

Bonsoir tout le monde,

Je souhaiterais ajouter des errorbar a mon bar chart. Voici la partie de mon code qui permet d'afficher mon resultat.

Code :

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

1
2
3
4
5
6
7
8
9
subplot(length(dog),1,d)
bar(mean_m')
xlabel('month')
ylabel('snr')
legend(speed)
titre=dog(d,:);
title(['Dog ',num2str(titre)])
hold on
errorbar(mean_m', error','xr')

Et vous trouverez en pièce jointe, ce que j'obtiens.
Comment puis-je modifier mon code pour avoir une barre d'erreur exactement au-dessus de chaque colonne de mon histogramme?
Merci pour votre aide

**keke0** · 31/07/2013, 17h11

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
subplot(length(dog),1,d)
bar(mean_m')
xlabel('month')
ylabel('snr')
legend(speed)
titre=dog(d,:);
title(['Dog ',num2str(titre)])
 
hold on
X(1,:)=0.745:11.745;
X(2,:)=0.92:11.92;
X(3,:)=01.1:12.1;
X(4,:)=01.255:12.255;
for i=1:length(speed)
<div style="margin-left:40px">mean_m(i,:);
error(i,:);
X(i,:);
errorbar(X(i,:),mean_m(i,:),error(i,:),'xr')</div>end

Avec ce code je me rapproche un peu plus de ce que je veux, mais ce n'est pas encore ça

**Jerome Briot** · 04/08/2013, 15h38

Essaie avec cette contribution : Bar Chart with Error Bars by Martina Callaghan

**keke0** · 05/08/2013, 19h01

Merci pour ton aide. Mais lorsque j'applique ce code, j'ai une erreur que je ne comprends pas:

Attempted to access errors(:,2); index out of bounds because size(errors)=[6,1,1].

Error in barwitherr (line 97)
    upperErrors = errors(:,2);

Error in average (line 70)
    barwitherr(error',mean_m')

**Jerome Briot** · 05/08/2013, 21h16

Pour commencer, renomme la variable error car c'est déjà le nom d'une fonction MATLAB.

Ensuit, il faut que les deux arguments de barwitherr aient les mêmes dimensions.

**keke0** · 05/08/2013, 21h57

J'ai bien renomme la variable error comme tu me l'as conseille par ecart.

Ensuit, je pensais justement que les deux arguments de barwitherr avaient les mêmes dimensions.

voici comment je declare mes arguments:

Code :

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

1
2
mean_m(v,m)=nanmean(nanmean(valeursok));
ecart(v,m)=nanstd(nanstd(valeursok));

et voici comment j'appelle la fonction barwitherr:

Code :

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

barwitherr(ecart',mean_m')

**Jerome Briot** · 05/08/2013, 22h00

Que renvoie ceci ?

Code :

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

whos mean_m ecart

**keke0** · 05/08/2013, 22h07

Name        Size                    Bytes         Class     Attributes

  ecart         1x6                48  double              
  mean_m     4x12              384  double

Effectivement ils n'ont pas les memes dimensions. Comment puis-je corriger cela?

**Jerome Briot** · 05/08/2013, 22h19

Difficile de t'aider sans voir plus de code et sans savoir ce que contient valeursok

**keke0** · 05/08/2013, 22h29

Voici mon code:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
clear all
close all
speed=['2';'3';'4';'5'];
dog=['3196';'3197';'3198'];
month=['01';'02';'03';'04';'05';'06';'07';'08';'09';'10';'11';'12'];
index=0;
dimdog=size(dog);
counter=0;
Cvide=zeros(1,16);
C(1,:)=Cvide;
 
for s_chosen=1:1:length(speed)
    for d_chosen=1:1:dimdog(1)
        for m_chosen=1:1:length(month)
            index=index+1;
            s_index=speed(s_chosen,:);
            d_index=dog(d_chosen,:);
            m_index=month(m_chosen,:);
 
            for num=0:100
                if num == 0
                    nom=strcat(dog(d_chosen,:),'_',speed(s_chosen,:),'_',month(m_chosen,:),'.xls');
                end
                if num~=0
                    nom=strcat(dog(d_chosen,:),'_',speed(s_chosen,:),'_',month(m_chosen,:),'_',num2str(num),'.xls');
                end
                fileopen=fopen(nom);
 
                if fileopen~=-1
                    counter=counter+1;
                    C(counter,:)=Cvide;
                    C2=xlsread(nom);
                    C(counter,:)=C2(1,:);
                    dogstudied(counter)=d_chosen;
                    speedstudied(counter)=s_chosen;
                    monthstudied(counter)=m_chosen;         
                end   
 
            end      
            resultat(index)=1;
 
        end
    end
end
 
for d=1:1:dimdog(1)
    for v=1:1:length(speed)
 
        for m=1:1:length(month)
            mean_m(v,m)=0;
            indexchien=find(dogstudied==d);
            indexvitesse=find(speedstudied==v);
            indexmois=find(monthstudied==m);
            indexchienvitesse=intersect(indexchien,indexvitesse);
            indexchienmois=intersect(indexchien,indexmois);
            indexcomplet=intersect(indexchienvitesse,indexchienmois);
            valeursok=C(indexcomplet,:);
 
%             for i=1:size(C(indexcomplet,:))
                mean_m(v,m)=nanmean(nanmean(valeursok));
                ecart(v,m)=nanstd(nanstd(valeursok));
%             end 
 
        end
    end
 
    figure (1)
    subplot(length(dog),1,d)
    barwitherr(ecart',mean_m')
    xlabel('month')
    ylabel('snr')
    l=legend(speed,'Location','SouthEastOutside');
    v = get(l,'title');
    set(v,'string','Speed (mph)');
    titre=dog(d,:);
    title(['Dog ',num2str(titre)])

Desolee, c'est un peu indigeste vu qu'il n'y a pas de commentaires.

**keke0** · 06/08/2013, 18h28

Je pense avoir trouve mon problème.
C'est au niveau de :

Code :

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

C2=xlsread(nom);

Mes fichiers Excel sont tous de la forme une colonne a 16 valeurs. Et lorsqu'une valeur ne me convenait pas je la remplaçait par un 0. Mais après j'avais un problème au niveau du calcul de la moyenne parce qu'il prenait en compte mais 0. Du coup, j'ai remplace mes 0 par des Nan. Et xlsread ne recopie pas tous les NaN.

Ex:

Code :

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

Ce dernier fichier test est compose de 14 NaN et 2 données. Mais xlsread me recopie que mais 2 données.
Du coup j'ai utilise,

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
 [n,txt,raw] = xlsread(nom);
C2=raw
C2 = 
 
    'NaN'
    'NaN'
    'NaN'
    'NaN'
    'NaN'
    [2]
    [2]
    'NaN'
    'NaN'
    'NaN'
    'NaN'
    'NaN'
    'NaN'
    'NaN'
    'NaN'
    'NaN'

Donc je récupère bien mais 16 données mais la j'ai un problème pour calculer ma moyenne vu qu C2 est de classe cell alors que C est de classe double. Et je n'ai pas trouve de fonction cell2double

Comment je peux me dépatouiller?

Invité · 06/08/2013, 18h43

Bonsoir,

Une solution incluant la fonction cell2mat que tu cherchais

Code :

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

1
2
3
isNum = cellfun(@isnumeric,C2);
nums = nan(size(C2));
nums(isNum) = cell2mat(C2(isNum));

Voir

**keke0** · 06/08/2013, 19h39

Merci Winjerome, c'est parfait !

Sinon pour revenir a la fonction barwitherr, j'ai une nouvelle erreur:

Undefined function or variable "lowerErrors".

Error in barwitherr (line 103)
if any(size(values) ~= size(lowerErrors))

Error in average (line 76)
    barwitherr(ecart',mean_m')

**Jerome Briot** · 06/08/2013, 20h57

Encore une fois :

Envoyé par Dut

Que renvoie ceci ?

Code :

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

whos mean_m ecart

**keke0** · 06/08/2013, 21h01

whos mean_m ecart
  Name        Size            Bytes  Class     Attributes

  ecart       4x12              384  double              
  mean_m      4x12              384  double

Invité · 08/08/2013, 20h38

Que représentent ces deux variables mean_m et ecart, et que souhaites-tu tracer précisément par rapport à elles ?

Ceci peut-être ?

Code :

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

1
2
3
4
mean_m = rand(4,12);
err = cat(3, ecart.', ecart.');
 
barwitherr(err, mean_m.')

**keke0** · 08/08/2013, 21h40

J'ai fait un petit schema pour expliquer ce que je souhaite tracer.

Donc ma variable mean_m correspond aux colonnes bleues, vertes et rouges, c'est-a-dire a la moyenne de mes SNR en fonction du mois et de la vitesse. Et ma variable ecart correspond au "segment" bleu clair qui me permet de voir si ma moyenne est aberante ou non.

Lorsque j'utilise la fonction barwitherr, j'obtiens ce message d'erreur:

Undefined function or variable "lowerErrors".

Error in barwitherr (line 103)
if any(size(values) ~= size(lowerErrors))

Error in average (line 76)
    barwitherr(ecart',mean_m')

**keke0** · 08/08/2013, 21h47

Je viens d'essayer ton code et ca marche. Mais j'aimerai bien comprendre...

Invité · 08/08/2013, 22h21

Pour l'explication :

Code :

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

err = cat(3, ecart_inferieur, ecart_superieur);

ecart_inferieur étant la distance avec la barre de dessous, et ecart_superieur la distance avec celle du dessus.

**keke0** · 09/08/2013, 00h08

Merci beaucoup!!!

Et as-tu un lien pour expliquer ce que signifie le (.') que tu as utiliser apres mes variables?

Ajouter des barres d'erreur a mon histogramme [Débutant]

MATLAB

Discussions similaires

Partager

Partager