Discussion :

[jojux]

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 :

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.

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 :

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

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

[phryte]

Bonjour,

Code :

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

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

Il manque la virgule après yd.

[rouxph]

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

Code :

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

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

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

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

• appel avec le paramètre "wd"
  

  Code :

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

  [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] ..."

[jojux]

Merci à tous de vos réponses .