Monsieur Kernel n'existe pas
2 participants
Page 1 sur 1
Monsieur Kernel n'existe pas
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?
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?
Re: Monsieur Kernel n'existe pas
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?
Un estimation en kernel avec un tableur... tu t'ennuies actuellement ?
... 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?
Un estimation en kernel avec un tableur... tu t'ennuies actuellement ?
Re: Monsieur Kernel n'existe pas
Non, je ne m'ennuie pas vraiment mais tu sais, certaines organisations internationales ont des valeurs différentes de nous autre chercheurs...
Re: Monsieur Kernel n'existe pas
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 ?
Re: Monsieur Kernel n'existe pas
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")
Re: Monsieur Kernel n'existe pas
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.
Re: Monsieur Kernel n'existe pas
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.
Re: Monsieur Kernel n'existe pas
Une journée pour faire des graphs que j'ai mis 5 minutes à faire sous stata. Chapeau excel.
Page 1 sur 1
Permission de ce forum:
Vous ne pouvez pas répondre aux sujets dans ce forum
|
|