24小时热门版块排行榜

返回列表

sulingerl

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 1458.3
散金: 63
红花: 1
帖子: 180
在线: 28.3小时
虫号: 2229424
注册: 2013-01-08
性别: MM
专业: 工业工程与管理

[求助] 求助SAS nlmixed程序问题

SAS程序如下(关于可靠性截尾数据的分析)：
data glass;
      input group temp volt time cen @@;
      datalines;
1 170 200 439  1
1 170 200 904  1
1 170 200 1092 1
1 170 200 1105 1
1 170 200 1105 0
1 170 200 1105 0
1 170 200 1105 0
1 170 200 1105 0
2 170 250 572  1
2 170 250 690  1
2 170 250 904  1
2 170 250 1090 1
2 170 250 1090 0
2 170 250 1090 0
2 170 250 1090 0
2 170 250 1090 0
3 170 300 315  1
3 170 300 315  1
3 170 300 439  1
3 170 300 628  1
3 170 300 628  0
3 170 300 628  0
3 170 300 628  0
3 170 300 628  0
4 170 350 258  1
4 170 350 258  1
4 170 350 347  1
4 170 350 588  1
4 170 350 588  0
4 170 350 588  0
4 170 350 588  0
4 170 350 588  0
5 180 200 959  1
5 180 200 1065 1
5 180 200 1065 1
5 180 200 1087 1
5 180 200 1087 0
5 180 200 1087 0
5 180 200 1087 0
5 180 200 1087 0
6 180 250 216  1
6 180 250 315  1
6 180 250 455  1
6 180 250 473  1
6 180 250 473  0
6 180 250 473  0
6 180 250 473  0
6 180 250 473  0
7 180 300 241  1
7 180 300 315  1
7 180 300 332  1
7 180 300 380  1
7 180 300 380  0
7 180 300 380  0
7 180 300 380  0
7 180 300 380  0
8 180 350 241  1
8 180 350 241  1
8 180 350 435  1
8 180 350 455  1
8 180 350 455  0
8 180 350 455  0
8 180 350 455  0
8 180 350 455  0
;

ods output parameterestimates=est;
proc nlmixed data=glass method=gauss noad noadscale QPOINTS=20 tech=quanew update=bfgs cov;
      parms beta=2.78 theta0=6.7 theta1=-0.44 logsigma=-3;
      z = beta*(log(time)-theta0-theta1*volt-u);
      f1=beta/time*exp(z-exp(z));
      F2=exp(-exp(z));
      ll=log(f1)*(cen=1)+log(F2)*(cen=0);
      sigma2=exp(logsigma)**2;
      model time ~ general(ll);
      random u ~ normal(0,sigma2) subject=group;
run;

可是运行的时候出不来结果，出现如下问题，请高手指教是哪儿出了问题？
NOTE: Execution error for observation 1.
WARNING: 未创建输出“parameterestimates”。请确保输出对象名称、标签或路径拼写正确。
      还需确保使用了相应的过程选项来生成请求的输出对象。如，确保没有使用 NOPRINT 选项。
NOTE: “PROCEDURE NLMIXED”所用时间（总处理时间）:
   实际时间       0.04 秒
   CPU 时间       0.03 秒

回复此楼