Discussion :

[AustinPower999]

Bonjour à tous,

Voilà, je souhaiterais tracer une courbe de Bézier (ou polynomiale) avec R.
Voici un lien sur ce que c'est pour les personnes comme moi qui n'en savait rien : http://fr.wikipedia.org/wiki/Courbe_de_Bézier

En gros, ce que je souhaite faire c'est, à partir de 3 points par exemple, tracer une courbe qui ne passe ni par tous les points, ni qui soit une droite, mais plutôt une courbe "moyenne".

Je n'ai pas réussi à trouver de solutions qui marchent sur internet...
Les density.plot ce n'est même pas la peine, j'ai aussi trouvé une fonction qui "lisse" une courbe, mais ça ne donne pas exactement ce que je souhaite...

Est-ce que quelqu'un a une idée ?


Merci :-)

[collif]

Bonjour,

j'ai aussi trouvé une fonction qui "lisse" une courbe

: Laquelle ?

As-tu essayé la fonction spline ?

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
 
# avant de l'utiliser va voir l'aide de la fonction
?spline
 
# ci-dessous l'exemple tiré de l'aide de la fonction
op <- par(mfrow = c(2,1), mgp = c(2,.8,0), mar = .1+c(3,3,3,1))
n <- 9
x <- 1:n
y <- rnorm(n)
plot(x, y, main = paste("spline[fun](.) through", n, "points"))
lines(spline(x, y))
lines(spline(x, y, n = 201), col = 2)
 
y <- (x-6)^2
plot(x, y, main = "spline(.) -- 3 methods")
lines(spline(x, y, n = 201), col = 2)
lines(spline(x, y, n = 201, method = "natural"), col = 3)
lines(spline(x, y, n = 201, method = "periodic"), col = 4)
legend(6,25, c("fmm","natural","periodic"), col=2:4, lty=1)

Est-ce que ça répond à ta question ?

[AustinPower999]

Merci pour ta réponse très rapide et claire !

A vrai dire, la fonction pour lisser n'ayant pas du tout correspondu à mes attentes j'ai fini par l'oublier, donc impossible de te répondre...

J'aurais une dernière question, le "n" dans : "lines(spline(x, y, n = 201), col = 2)" tu le choisis arbitrairement ou il existe un moyen pour le calculer en fonction de la courbe ? Comme on calculerait une fréquence d'échantillonnage par exemple...


En tout cas, merci beaucoup !

[collif]

Dans l'aide de la fonction :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
 
spline(x, y = NULL, n = 3*length(x), method = "fmm",
            xmin = min(x), xmax = max(x), xout, ties = mean)
 
...
 n: if ‘xout’ is left unspecified, interpolation takes place at
          ‘n’ equally spaced points spanning the interval [‘xmin’,
          ‘xmax’].
...

Le diaporama de la présentation d'un cours dispensé à Agrocampus Ouest par David Causeur (Rennes) est disponible sur Internet. Si les explications sont un peu hardues sans le speech, il a tout de même le mérite de dresser un éventail de techniques d'ajustement non paramètriques dont les méthodes de lissage font partie.

http://math.agrocampus-ouest.fr/infoglueDeliverLive/digitalAssets/31936_NonParam.pdf


Un deuxième lien dans le même objectif :
http://wikistat.fr/pdf/st-m-app-non-param.pdf

[AustinPower999]

Ok, merci !