我底下的SAS程式有一个问题 用不同的方法try出问题点 但却解决不了 程式不长 请大家可以的话帮我看看 也请大家给我一些想法 谢谢 data inhaler; input time d x; datalines; 1 0 19 1 1 59 1 0 15 1 0 15 1 0 17 1 1 51 1 0 12 2 1 51 2 0 11 2 0 11 2 0 17 2 1 57 2 0 14 2 1 54 2 0 11 3 1 58 3 0 15 3 0 14 3 0 13 3 1 51 3 0 13 4 0 13 4 0 13 4 1 53 4 0 10 4 1 50 4 0 15 4 1 55 run; proc nlmixed data=inhaler corr ecorr; parms beta0=0 beta1=1 lo=0.1 sd=1; ff = ((beta1*x + beta0 - f*(lo)**(1/2)) / ((1-lo)**(1/2))); K=1-d; L=(-1)*probnorm(ff); /* L=1-probnorm(ff); */ N=L**K; model time ~ general(N); random f ~ normal(0,sd*sd) subject=time; estimate 'beta0' beta0; estimate 'beta1' beta1; estimate 'lo' lo; run; 如果用绿体字的L 我可以利用nlmixed的方法 估计出beta等参数 结果会显示出很多表格 得到我要的结果 但是因为我要跑得方程式中L 的地方要改成黄体字 在不太了解probnorm()的语法在SAS有什麽特性之下 如果 我在L这个式子(绿体字)里加入1 (变成注解的黄体字) 程式虽然没有bug但是却没有跑出结果（结果只跑开头就无疾而终） 学长跟我说跑到开头就停表示程式有问题 但我苦思不出原因... 请大家帮帮我这个SAS新手 感激感激 --

※ 发信站: 批踢踢实业坊(ptt.cc) ◆ From: 140.113.129.127 ※ 编辑: Eagler 来自: 140.113.129.127 (04/12 07:32) ※ 编辑: Eagler 来自: 140.113.129.127 (04/12 09:46)