NLMIXED 参数初始值怎么都设不好啊？
如题，修改了很多次参数初始值“   parms a=0 b=0.001; ”， 还是没有搞出来。

高手帮忙啊~

多谢！


###############################################

SAS程序如下：

data trees;
  input n_ingrow n_live n_dead year live_off yr_interval;
  cards;
1    600    116    1968    715    7
19    495    124    1979    600    11
21    480    36    1984    495    5
41    485    36    1989    480    5
24    469    40    1994    485    5
32    456    45    1999    469    5
;

proc nlmixed data=trees;
   parms a=0 b=0.001;                                                         
   eta = a + b*year;                              
   p_ann_mort = exp(eta)/(1+exp(eta));                        
   p_ann_surv = 1-p_ann_mort;                                                
   p_surv = p_ann_surv**yr_interval;                                 
   p_mort = 1-p_surv;
   mean_mort = live_off*p_mort;
   ll=lgamma(1/alpha + n_dead) - lgamma(1 + n_dead) - lgamma(1/alpha)   + n_dead*log(alpha*mean_mort) - (1/alpha + n_dead)*log(1 + alpha*mean_mort);        
         * log likelihood function for negative binomial regression ;
   model n_dead ~ general(ll);
   predict live_off*p_mort out=predobs1;
   ods output ParameterEstimates =PEs; *create an output table PEs you needed;
run;