Logiciel R

[invitec5af4239]

Bonjour

je travaille sous le logiciel R et je rencontre un probleme au moment de tracer une fonction, voici la totalité du code :

vecteur=rnorm(1000,0,1)

K<-function(x,h){
g=ifelse(x<(-(h/2)),0,ifelse(x>(h/2),0,(1/h)))
g
}

curve(K(x,6),(-5),5,1000,add=FALSE,type="l",c ol="red")

vecteur2=matrix(vecteur,ncol=1 )

F<-function(x,h) {
m=matrix(NA,1000,1)
n=matrix(x,1000,1)
for (i in 1:1000) {m[i]=K(((n[i]-vecteur2[i])/h),h)
}
p=(1/(1000*h))*sum(m)
p
}

Mes fonctions marchent et le graphe de K est le bon, mais impossible de dessiner la fonction F(x,6) (par exemple) avec la fonction "curve" (pourtant les calculs de F(3,6) ou F(2,6) ou autres donnent des résultats). En fait R me repond que x et y ont des tailles differentes et qu'il ne peut pas tracer...

Pouvez vous svp m'aider à resoudre le problème ?
Merci a tous
-----

Publicité

[kinette]

Bonjour,
J'ai pas trop le temps de regarder ton code en détail, mais je vois que tu as utilisé des matrices... vérifie donc comment tes résultats sont donnés... je pense que le problème vient du fait que tu dois demander le plot d'un truc qui n'a pas les mêmes dimensions (d'ailleurs que souhaites-tu représenter, je n'ai pas compris).
J'avoue éviter l'utilisation de matrices quand c'est possible (au profit d'arrays qui sont plus "souples").

K

[invitec5af4239]

Voila le problème :

Je génère un vecteur aléatoire qui suit une loi normale N(0,1), je l'appelle "vecteur"...
Ma fonction K est assez claire je pense, c'est une sorte d'indicatrice de sauf qu'elle vaut plutot que 1. Elle fonctionne et je peux la traçer.
Le but est de créer la fonction F telle que :
.


De la façon dont j'ai procédé je ne peux pas traçer F pour differentes valeurs de h, or c'est ce que je veux faire... Vois tu mieux comment m'en sortir ?

Merci d'avance