Discussion :

[opium]

Bonjour,

Ceci est un appel d'aide, afin de programmer une fonction mathématique.

J'ai besoin d'implémenté une fonction qui retourne l'estimation de la fonction de densité d'un échantillon (X1....Xn) de valeurs numérique.
En effet l'estimation de cette densité peut se faire par l'estimation par noyau "kernel density estimation"
En resumé j'ai besoin de coder ceci:
http://support.sas.com/onlinedoc/913...ist_sect24.htm
en C++ ou C...
Savez vous si il existe des codes déja implémenté qui réaliserai ceci, j'ai juste besoin de l'implémentation en univarié (un seul échantillon de valeur X1....Xn)

Cordialement.

[gbdivers]

Bonjour et bienvenue sur le forum

As tu fais une petite recherche google ? Il y a plusieurs liens et tu auras déjà un début de réponse

Merci

[opium]

Ui j'ai déjà effectué une recherche, mais en vain...

J'ai dejà l'implémentation R et Matlab de cette fonction mais j'en ai besoin en C ou C++.
Apres les recherches google je trouve des choses vraiment compliquées souvent du code qui fait "KDE" pour des problèmes multivariés, souvent des programmes assez lourd, qui intègre un tas de fichier et je me perds...

Je n'est pas vraiment d'expérience en C et C++...

[gbdivers]

Google me donne wikipédia en premier, avec 2 libs en référence
il y a également un exemple d'implémentation pour Boost.random
Regarde ces liens et leurs implémentations et si tu as un problème, on pourra t'aider

[opium]

parmi ces 2 librairies en référence une seule semble vraiment faire ce qui m'interesse:
http://www.ics.uci.edu/~ihler/code/kde.html
mais c'est pour du multivarié et quand je telecharge le fichier il ya un tas de routine .m et .cpp mélangés et je ne comprend clairement pas....
ce que je demande me semble beacoup plus simple....

[oodini]

C'est dans C++11 (bibliothèque random) ou Boost.

[opium]

perdu dans mes recherches je ne sais plus ou donner de la tête...
need help, qq'un qui me guide et qui me pose les bases...

[Invité]

slt,

voilà un exemple si tu souhaites réinventer la roue...ou éviter de te taper du boost.

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
 
#include <iostream>
#include <vector>
typedef std::vector<float> Data;
float Ko(const float &t){
  //implement here Ko(t)
  return t;
}
 
float estimator(const Data& v, const float &y, const float &lambda, float(*ko)(const float &t)){
  float sum=0;
  for(size_t i=0; i<v.size(); ++i){
    sum+=ko((y-v[i])/lambda);
  }
  return sum;
}
int main(int argc, char* argv[]){
  Data v;
  v.push_back(1.0);
  v.push_back(2.0);
  v.push_back(3.0);
  std::cout<<estimator(v, 5.0, 2.0, &Ko)<<std::endl;
}

je te laisse vérifier que pour Ko(t)=t, tu as bien le résultat escompté et corriger le code le cas échéant.

[opium]

Merci bien je vais essayer