Code R pour l'estimation de la densité d'une loi normal par des fonctions spine ( ca me rends fou!)

[invited5106e4f]

Bonjour ,
Je rencontre des difficultés monstres a coder et tracer l'estimation de la densité d'une loi normale(0,1) par des fonctions splines .
Tout ce que j'ai reussi a produire pour l'instant est ce bout de code suivant:

library(splines)
xn=rnorm(100)
xn = sort(xn)

b0=bs(xn)
str(b0)
bsp=data.frame(y=xn,x1=b0[,1],x2=b0[,2],x3=b0[,3])
reg=lm(y~x1+x2+x3,data=bsp)
summary(reg)
beta0=reg$coef[1]
beta1=reg$coef[2]
beta2=reg$coef[3]

f=function(xn){beta0*b0[,1]+beta1*b0[,2]+beta2*b0[,3]}

curve(f,0,100,col="red")
points(x,y,pch="+")
b0
g=function(x){-4.79*x+2.94*x^2+5.88*x^3}
curve(g,-20,20)


Ca ne fonctionne pas et je n'arrive toujours pas a savoir pourquoi
Aidez moi svp , je vous en serai tres reconnaissant
-----

[minushabens]

Ta fonction f a pour argument xn mais n'en fait pas usage, i.e. elle renvoie une valeur constante, indépendante de xn.

Par ailleurs je te déconseille d'utiliser "=" pour l'affectation d'une valeur à un objet.

[invited5106e4f]

Merci de m'avoir répondu, en fait tout fonctionne jusqu'à beta2 mais le problème c'est que je ne comprends pas comment passer de la base des b splines b0 a un graphe des points simulés selon la N(0,1) par lesquels sont supposés passé une spline ayant la forme d'une N(0,1)
En plus une fois la spline obtenue , on m'a dit qu'il fallait normaliser la fonction spline pour obtenir la densité.

Toute la partie apres la ligne beta2 est incorrecte.
Si tu as un code qui correspond au probleme je veux bien le tester, toutes les propositions sont les bienvenues , merci encore