SAS nlmixed 程序的问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

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 秒

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：nlmixed Mixed mix nlm proc nlmixed 程序

求助啊，求助啊！脑袋都大了