Bonjour à tous,
j'ai un souci d'utilisation de la fonction "lsq_splin" dans le cadre d'une régression des moindres carrés pondérés (weighted least squares cubic spline fitting), ma recherche n'a donné qu'une source anglophone, mais qui discute plutôt d'un souci de taille du vecteur "breakpoint" (http://comp.soft-sys.math.scilab.nar...plin-lsq-splin).
La fonction scilab est donc "lsq_splin". Dans l'exemple de l'aide, la pondération (vecteur de même dimension appelée wd) n'est pas utilisée, ce qui implique un poids égal pour toutes les données par défaut. Mais lorsque je veux rajouter un vecteur de pondération (qui est bien de la même taille que les abscisses et ordonnées expérimentales), j'ai toujours le message d'erreur:
malgré le fait qu'il ne manque pas de parenthèse (enfin, je crois) et que les vecteurs sont tous bien de mêmes dimensions:
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4 -->[y, d] = lsq_splin(xd, yd [, wd], x); !--error 3 En attente d'une parenthèse droite.
A titre de test, j'ai modifié l'exemple de l'aide Scilab pour les pondérations, j'ai utilisé les lignes suivantes, mais sans succès:
Ce qui signifierait que la régression serait de plus en plus sensible au bruit à xd croissant (mais ça, c'est toujours pour l'exemple test).
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2 wd = linspace(5,1000,m)'; [y, d] = lsq_splin(xd, yd [, wd], x);
Quelqu'un a déjà eu ce problème?
Merci de vos lectures et de vos réponses,
Jojo
Partager