hustlfm
data hust;
167 pump=_N_;
168 set Lfm.Tfy;
169
170 proc nlmixed data=hust;
171 parms beta1=0.7 beta2=-3 ,epthta1=0.07,epthta2=-0.1,epthta12=1.4;
172 x=log(fail_time);
173 gamma=beta1+b1;
174 thta=exp(-(beta2+b2)/gamma);
175 eta=x*(beta1+b1)+beta2+b2;
176 mean=exp(eta);
177 model y ~ poisson(mean);
178 random b1 b2 ~ normal([0,0],[epthta1,epthta12,epthta2]) subject=pump;
179 estimate 'gamma-thta' gamma-thta;
180 run;
ERROR: ESTIMATE statement expressions are not allowed to depend on random effects.
NOTE: SAS 系统由于错误而停止了该步的处理。
cran及各位达人,看看我这个错误是什么原因？实在是不懂到底什么地方出错了？是我的数据集出错还是模型构建的问题？
chen2015
[未知用户]
你好 我也做了这个程序 但是怎么都不出结果 麻烦你帮我看一看 谢谢 代码附在下面
proc nlmixed data=new;
parms b0=0 b1=0 b2=0 b3=0 k=0.001;
eta=b0+b1*age+b2*model+b3*year;
bounds 0<=k;
mu=exp(eta);
loglike=nclaims*log(mu/(1+k*mu))-log(fact(nclaims))+(nclaims-1)*log(1+k*nclaims)-mu*(1+k*nclaims)/(1+k*mu);
prob=exp(loglike);
model nclaims~general(loglike);
run;
