zwnSAS121 发表于 2017-3-24 13:19
proc sql noprint;
create table table3 as
select
还有 我想问:
我按下面程序:
proc nlmixed data=A;
parms b0=0.3619 b1=9.2584 b2=0.1223 b3=-0.7684 b4=-0.8946 b5=-0.5084 b6=-0.1055 b7=0.5279 b8=0.5699 b9=0.9933 b10=1.1098 b11=0.211 b12=3.2812 b13=-0.3272 b14=0.466 b15=-1.0409 b16=-0.5441 b17=22.462 b18=-0.5337 b19=0.1345 b20=0.4873 b21=0.5315 b22=10.3754 b23=-1.1663 b24=1.3985 b25=-1.2495 b26=-0.2036 b27=-0.2045 b28=0.05201 b29=-0.03028 b30=-0.04162 b31=-0.02524 b32=0.5906 b33=-0.2547 b34=-0.2493 b35=-0.2544 sigma=1;
y=b0+u+b1*month+b2*motheragedegree+b3*fatheragedegree+b4*mothereducation+b5*fatherBMIdegree+b6*fathereducation+b7*income+b8*inhabitant+b9*miscarriage+b10*parity+b11*gravidity+b12*gestation+b13*gender+b14*gwgdegree+b15*delivery+b16*preparation+b17*drink+b17*exercise+b19*vomit+b20*vaginal+b21*smoke+b22*noise+b23*computer+b24*phone+b25*depressiondegree+b26*nightfeed+b27*feeding+b28*breastnumber+b29*sleep+b30*sleeptime+b31*activity+b32*cold+b33*diarrhea+b34*complementaryfeedingtime+b35*motherBMIdegree;
p=1/(1+exp(-y));
model lowweightWAZ~binary(p);
random u~normal(0,sigma) subject=people;
title 'only slope is random';
run;
结果出现ERROR: Quadrature accuracy of 0.000100 could not be achieved with 31 points. The achieved accuracy was 1.000000.
这一问题该如何解决?