2013-01-23 3 views
2

J'ai estimé quelques données avec 'densité' et défini le noyau comme Epanechnikov, y a-t-il une fonction inverse qui peut me récupérer une valeur estimée à sa valeur d'origine?Fonction inverse d'Epanechnikov (noyau)

set.seed(1) 
data<-rnorm(518,10,3) 
hist(data) 


data_E<-density(data,bw="sj",kernel="epanechnikov") 
hist(data_E$y) 

k<-quantile(data_E$y, 0.95) 
k 

Dans l'exemple ci-dessus exposé, j'ai une taille d'échantillon 518, de ce que j'ai estimé la densité par la méthode du noyau (Epanechnikov), après cela, je voudrais trouver la valeur d'un certain élément de cette "nouvelle donnée" en termes de "données originales".

+0

Vous devriez ajouter plus de détails, peut-être nous donner un petit exemple avec quelques valeurs. – Spacedman

Répondre

1

Non. Vous avez supprimé des informations dans le processus d'estimation de la densité. La fonction density ne stocke pas son entrée comme le fait, par exemple, ecdf.

Modifier. Je peux peut-être adoucir un peu cela. Vous pouvez utiliser l'information en densité pour reconstruire une estimation approximative des valeurs d'origine. Vous devez connaître la taille de l'échantillon, qui est enregistrée dans l'objet renvoyé, puis utiliser l'estimation pour générer de nouvelles données. Intégrez la densité, divisez la fonction de probabilité cumulative de manière égale le long de sa plage de support et vous pourriez construire une estimation quantile empirique.