Fonction lsq_splin et pondération (weighted)

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:
Code:

1 2 3 4 -->[y, d] = lsq_splin(xd, yd [, wd], x); !--error 3 En attente d'une parenthèse droite.
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:

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:
Code:

1 2 wd = linspace(5,1000,m)'; [y, d] = lsq_splin(xd, yd [, wd], x);
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).

Quelqu'un a déjà eu ce problème?

Merci de vos lectures et de vos réponses,

Jojo

Bonjour,

Code:

[y, d] = lsq_splin(xd, yd, [, wd], x);

Il manque la virgule après yd.

quand dans l'aide en ligne la structure d'appel de la fonction est décrit sous la forme :

Code:

[y, d] = lsq_splin(xd, yd [, wd], x);

ça signifie que le paramètre "wd" est optionnel, il y a donc 2 manières d'écrire l'appel de la fonction :

appel sans le paramètre "wd"

Code:

[y, d] = lsq_splin(xd, yd, x);

appel avec le paramètre "wd"

Code:

[y, d] = lsq_splin(xd, yd , wd, x);

c'est ce qui explique que scilab considère qu'il manque une virgule entre "yd" et l'ouverture d'une matrice (par "[") quand il lit "...yd [, wd] ..."

Merci à tous de vos réponses ;).