J'utilise la procédure NLMIXED tester un effet aléatoire arbre dans mon modèle le code est le suivant :
proc nlmixed data=GourHetr.gour qpoints=5 ;
parms beta0=7.2 beta1=-0.03 s2u=0.01;
eta=beta0 + beta1 * L + u*L;
expeta=exp(eta);
p=expeta/(1+expeta);
bounds s2u>0;
model dis ~ binary(p);
random u ~ normal(0,s2u) subject=arbre;
run;
Je pense que mon code est correct, la méthode converge mais SAS ne calcule pas l'erreur standard et la statistique de Student pour s2u, comment ça se fait ?
Merci