Changement de régimes avec PROC MODEL

[Mark531]

Bonjour à tous,

J'aimerais implémenter un modèle multi-régimes avec fonction de transition logistique. Le modèle à 2 régimes converge bien :

Code :

1
2
3
4
5
6
7
proc model DATA=indic;
g = 1 / (1 + exp(-gamma * (ident - c)));
lngn1 = a1 + b11 * lncarn + b12 * lnpgpp;
lngn2 = a2 + b21 * lncarn + b22 * lnpgpp;
lngn = g * lngn1 + (1 - g) * lngn2;
fit lngn /method=marquardt;
run;

Pour passer en multi-régimes, je dois variabiliser le coefficient de transition entre tous les régimes. Mais pour une raison que j'ignore, la modification suivante provoque une erreur à l'exécution de PROC MODEL :

Code :

1
2
lngn = k1 * g * lngn1 + k2 * g * lngn2;
k1 + k2 = 1;

Pourtant, les coefficients k1 et k2 permettent bien de "mélanger" les deux modèles. Si vous avez une méthode pour implémenter plus de 2 régimes, je suis preneur.

Merci,
Mark