Proc NLMIXED - erreur paramètre

[alejandro27]

Bonjour,

J'ai un petit problème avec ma proc NLMIXED :

Code :

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
 
**** Macro pour LOG-NORMAL;
 	%Macro M_nlm_lognormal(dsin, var1);
		proc nlmixed DATA=&dsin. qmax=780;  
		parms  intercept=0  p_period=0 p_totcharge_6b=1 sd=1; %*Initiation des parametres;
 
%* Force les ecart-type a etre positif;
			linp=intercept + p_totcharge_6b*log_totcharge_6b;* + eta;
			mu = exp(linp); %* Introduction du lien log;
		model
			&var1. ~ normal(mu,sd*sd); 
		run; 
		quit;
	%mend;
	%M_nlm_lognormal(dsin=temp1, var1=totcharge);

et il me renvoit

Code :

1
2
3
4
 
WARNING: The final Hessian matrix IS NOT positive definite, AND therefore the estimated
         covariance matrix IS NOT full rank AND may be unreliable.  The variance of some
         parameter estimates IS zero OR some parameters are linearly related TO other parameters.

ou

Code :

1
2
 
ERROR : No valid parameter point were found

est-ce un problème d'initialisation ou de définition des paramètres ??

Merci à vous !

[Manoutz]

Bonjour,

Si tu es sur de l'écriture de ton modèle cela ne vient probablement pas de la. En tout cas, le modèle proposé ne converge pas et donc il faut opérer qq changements. La gestion de la convergence est beaucoup plus ardue avec les modèles non linéaires en comparaison aux modèles linéaires.

Tu peux par exemple:

• Changer tes parmètres initiaux, voire mettre en place une grille de recherche
• Essayer différentes méthodes d'approximation (si applicable)
• Essayer différentes techniques d'optimisation
• Si vraiment tu n'y arrives pas, tu peux esssayer avec la proc NLIN (si ma compréhension est bonne tu n'as pas d'effets aléatoires)

Manoutz