我底下的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)