全部版块 我的主页
论坛 站务区 十一区 新手入门区
1730 1
2014-01-19
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 秒



二维码

扫码加我 拉你入群

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

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

全部回复
2014-1-19 17:40:16
求助啊,求助啊!脑袋都大了
二维码

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群