Discussion :

[rollo]

Bonjour à tous,

Dans un article qui traitait de signaux périodiques, une analyse spectrale donnait ceci :



On y voit les lignes de moyenne du spectre et de l'intervalle de confiance ("95% CI") permettant d'affirmer que c'était un pic "réel", et pas un artefact.
Question: comment calcule-t-on (sous Matlab) la valeur de cet intervalle de confiance ?

Merci de votre aide !

Eric

[dionisis]

Salut Eric,

Pour obtenir l intervalle de confiance il faut faire une hypothese sur la distribution des donnees observees. Ensuite tu estimes les parametres de la distribution a partir de ces observations.

En matlab, on pourrait avoir quelquechose qui ressemble a ca:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
% simulons quelques donnees suivant une loi normale N(mu, sigma^2)
n = 50;
sigma = 450; 
mu = 250;
x = mu + sigma*randn(n,1);
 
% estimation de la moyenne et de la volatilite
est_mu = mean(x);
est_sigma = std(x);
 
% intervalle de confiance
p = 0.95; % on cherche le quantile a 95%;
CI = norminv(p,est_mu,est_sigma);

J espere que ca peut aider.

Cordialement,
Denis

[rollo]

Bonjour Denis,
Je viens seulement de voir ta réponse +6 mois après que tu l'aies écrite...
Il n'est jamais trop tard pour dire : merci !

Eric

Salut Eric,

Pour obtenir l intervalle de confiance il faut faire une hypothese sur la distribution des donnees observees. Ensuite tu estimes les parametres de la distribution a partir de ces observations.

En matlab, on pourrait avoir quelquechose qui ressemble a ca:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
% simulons quelques donnees suivant une loi normale N(mu, sigma^2)
n = 50;
sigma = 450; 
mu = 250;
x = mu + sigma*randn(n,1);
 
% estimation de la moyenne et de la volatilite
est_mu = mean(x);
est_sigma = std(x);
 
% intervalle de confiance
p = 0.95; % on cherche le quantile a 95%;
CI = norminv(p,est_mu,est_sigma);

J espere que ca peut aider.

Cordialement,
Denis