Monsieur Kernel n'existe pas

Aller en bas

Monsieur Kernel n'existe pas

Message  Fred le Mar 15 Nov - 12:30

C'est en cherchant à reproduire un graphique stata "Kernel Density" au format excel que je me suis aperçu avec stupeur que monsieur Kernel n'existe pas.
En fait, on dit "estimateur par noyau" en bon français...

Mais tout ça ne règle pas mon problème. Je n'ai aucune idée de la manière dont on calcule une densité de ce type sous excel. Quelqu'un s'est-il déjà risqué à ce genre d'aventures?
avatar
Fred
Admin

Messages : 370
Date d'inscription : 17/11/2009
Age : 34
Localisation : Doctorant

Voir le profil de l'utilisateur http://themanwhobloggedtheworld.blogspot.com/

Revenir en haut Aller en bas

Re: Monsieur Kernel n'existe pas

Message  flbresson le Mar 15 Nov - 13:23

En cherchant à reproduire une inversion de matrice avec un boulier chinois, je me suis aperçu avec stupeur....

... Je n'ai aucune idée de la manière de réaliser un pivot de Gauss avec un boulier chinois. Quelqu'un s'est-il déjà risqué à ce genre d'aventures?


tongue


Un estimation en kernel avec un tableur... tu t'ennuies actuellement ?
avatar
flbresson

Messages : 360
Date d'inscription : 19/11/2009
Age : 39
Localisation : Clermont-Ferrand

Voir le profil de l'utilisateur http://sites.google.com/site/flbresson2/

Revenir en haut Aller en bas

Re: Monsieur Kernel n'existe pas

Message  Fred le Mar 15 Nov - 13:54

Non, je ne m'ennuie pas vraiment mais tu sais, certaines organisations internationales ont des valeurs différentes de nous autre chercheurs... Rolling Eyes
avatar
Fred
Admin

Messages : 370
Date d'inscription : 17/11/2009
Age : 34
Localisation : Doctorant

Voir le profil de l'utilisateur http://themanwhobloggedtheworld.blogspot.com/

Revenir en haut Aller en bas

Re: Monsieur Kernel n'existe pas

Message  flbresson le Mar 15 Nov - 14:01

Désolé d'avoir chambré mais c'était trop tentant. Ne peux tu pas simplement récupérer une série de points à partir de ton estimation en kernel depuis Stata (ou autre) et tracer par la suite la courbe ou l'ajouter sur un histogramme ?
avatar
flbresson

Messages : 360
Date d'inscription : 19/11/2009
Age : 39
Localisation : Clermont-Ferrand

Voir le profil de l'utilisateur http://sites.google.com/site/flbresson2/

Revenir en haut Aller en bas

Re: Monsieur Kernel n'existe pas

Message  Fred le Mar 15 Nov - 14:08

Tu fais comment pour récupérer les points?
avatar
Fred
Admin

Messages : 370
Date d'inscription : 17/11/2009
Age : 34
Localisation : Doctorant

Voir le profil de l'utilisateur http://themanwhobloggedtheworld.blogspot.com/

Revenir en haut Aller en bas

Re: Monsieur Kernel n'existe pas

Message  flbresson le Mar 15 Nov - 14:19

Avec Stata, je ne sais pas. Sous R:

Code:

% générons aléatoirement une série X, par exemple 10000 valeurs d'une loi normale centrée réduite
X<-rnorm(10000)

%estimation en kernel (gaussien par défaut) maintenant et conservation des résultats dans un objet machin
machin <- density(X)

% Les coordonnées sont dans les items x et y de l'objet machin. Pour obtenir sous forme matricielle les points de la courbe:
cbind(machin$x,machin$y)

% Pour exporter au format csv (version française) ces coordonnées:
write.csv2(cbind(machin$x,machin$y), file="machin.csv")
avatar
flbresson

Messages : 360
Date d'inscription : 19/11/2009
Age : 39
Localisation : Clermont-Ferrand

Voir le profil de l'utilisateur http://sites.google.com/site/flbresson2/

Revenir en haut Aller en bas

Re: Monsieur Kernel n'existe pas

Message  Fred le Mar 15 Nov - 16:39

Sous R:


Spoiler:
Bon, en fait, c'est pas compliquer d'extraire les valeurs sous stata. Il suffit de rajouter une option gen(var_x var_y). Stata va stocker les points d'estimation (var_x) et les densités (var_y). Si on rajoute l'option at(var), les points d'estimations correspondent aux valeurs de la variable désirée.
avatar
Fred
Admin

Messages : 370
Date d'inscription : 17/11/2009
Age : 34
Localisation : Doctorant

Voir le profil de l'utilisateur http://themanwhobloggedtheworld.blogspot.com/

Revenir en haut Aller en bas

Re: Monsieur Kernel n'existe pas

Message  flbresson le Mar 15 Nov - 16:41

Je me disais bien aussi qu'un logiciel aussi cher que Stata ne pouvait pas ne proposer une telle option. Je pense que ce sera mieux ainsi qu'avec ton boulier... pardon ton tableur.
avatar
flbresson

Messages : 360
Date d'inscription : 19/11/2009
Age : 39
Localisation : Clermont-Ferrand

Voir le profil de l'utilisateur http://sites.google.com/site/flbresson2/

Revenir en haut Aller en bas

Re: Monsieur Kernel n'existe pas

Message  Fred le Mar 15 Nov - 19:28

Une journée pour faire des graphs que j'ai mis 5 minutes à faire sous stata. Chapeau excel.
avatar
Fred
Admin

Messages : 370
Date d'inscription : 17/11/2009
Age : 34
Localisation : Doctorant

Voir le profil de l'utilisateur http://themanwhobloggedtheworld.blogspot.com/

Revenir en haut Aller en bas

Re: Monsieur Kernel n'existe pas

Message  Contenu sponsorisé


Contenu sponsorisé


Revenir en haut Aller en bas

Revenir en haut

- Sujets similaires

 
Permission de ce forum:
Vous ne pouvez pas répondre aux sujets dans ce forum