下面是我写的,可是不知道为啥run了之后什么都不显示,就是个specification出来了。是什么问题? 多谢指教!
proc nlmixed data=mp.hsb2;
parms beta0 0 theta 0;
science=beta0*math+theta;
random theta~normal(0,sigma**2) subject=science;
prob=(1/sqrt(2*pi)*sigma)*exp(-(science-beta1*socst-theta)**2/2*(signma*2));
ll=log(prob);
model science~general(ll);
run;
[此贴子已经被作者于2008-6-16 9:42:43编辑过]