24小时热门版块排行榜    

查看: 1206  |  回复: 0

sulingerl

金虫 (小有名气)

[求助] 求助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 秒

» 猜你喜欢

comeon
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子,希望有用哦~

科研从小木虫开始,人人为我,我为人人
相关版块跳转 我要订阅楼主 sulingerl 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] +3 q1zg6x4lf3 2026-06-12 3/150 2026-06-13 12:29 by 80j16ujcu8
[考研] O(∩_∩)O +3 q1zg6x4lf3 2026-06-12 3/150 2026-06-13 12:13 by 80j16ujcu8
[硕博家园] O(∩_∩)O +4 jk31rb6vyz 2026-06-12 5/250 2026-06-13 12:07 by 80j16ujcu8
[公派出国] π_π +3 m584r2z4j3 2026-06-11 3/150 2026-06-13 12:04 by 80j16ujcu8
[考博] T^T +3 m584r2z4j3 2026-06-11 3/150 2026-06-13 12:04 by 80j16ujcu8
[教师之家] 售SCI文章,我:8.O.5.51O.5.4,科目齐全,可十急 +3 uaz7enpl3v 2026-06-12 4/200 2026-06-13 09:59 by ZYNi吖
[公派出国] 售SCI文章,我:8.O.5.51O.5.4,科目齐全,可十急 +3 uaz7enpl3v 2026-06-12 3/150 2026-06-13 09:22 by t4hz8376wy
[找工作] 售SCI文章,我:8.O.5.51O.5.4,科目齐全,可十急 +4 q1zg6x4lf3 2026-06-12 4/200 2026-06-13 08:53 by t4hz8376wy
[找工作] *囍 +3 m584r2z4j3 2026-06-12 3/150 2026-06-13 08:29 by t4hz8376wy
[论文投稿] →_→ +5 m584r2z4j3 2026-06-11 5/250 2026-06-13 08:13 by t4hz8376wy
[考博] 售SCI文章,我:8.O.5.51O.5.4,科目齐全,可十急 +3 uaz7enpl3v 2026-06-12 3/150 2026-06-13 05:24 by fjz4noejhx
[考研] ╭(╯ε╰)╮ +3 m584r2z4j3 2026-06-11 3/150 2026-06-13 04:07 by 17a9ac4dm7
[考研] o>_<o +4 m584r2z4j3 2026-06-11 4/200 2026-06-13 03:49 by 17a9ac4dm7
[硕博家园] T^T +4 m584r2z4j3 2026-06-11 5/250 2026-06-13 03:47 by 17a9ac4dm7
[基金申请] 面上项目没有好文章就没希望了吗? +19 1234567wang 2026-06-07 21/1050 2026-06-12 21:24 by littlezl
[考博] [急招]吉林大学招收26年9月入学考核制博士生1名,截止今晚24点前 +3 hehejin 2026-06-10 4/200 2026-06-12 18:24 by 冻雨天
[基金申请] 前几天时间戳更新了 +8 vito刘 2026-06-11 13/650 2026-06-12 14:17 by BlakeReary
[教师之家] AASSDDFF +3 m584r2z4j3 2026-06-11 3/150 2026-06-12 11:28 by zs48m4lz8g
[考博] 还有课题组有博士名额吗 +5 冻雨天 2026-06-07 6/300 2026-06-08 20:59 by RCS_tyro
[考博] 博士申请 +3 张城斐 2026-06-07 3/150 2026-06-08 09:15 by A大朋同学
信息提示
请填处理意见