| 查看: 1232 | 回复: 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 秒 |
回复此楼» 猜你喜欢
各位道友，我要去昆明玩几天，回来见。
已经有8人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有7人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有8人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有22人回复
-
是这周出结果还是下周出结果？
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有49人回复
找到一些相关的精华帖子，希望有用哦~
SAS编程的系统聚类分析
已经有4人回复
- 一些常用SAS的程序【转载】
已经有30人回复
- 如何用SAS 进行回归系数的显著性检验。。。。。
已经有3人回复
- SAS性别和体重血糖中的回归分析
已经有4人回复
- 求JMP 7.0 (SAS Institute, 2007)程序按照序列号
已经有3人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
