立即打开
nlmixed 过程出现错误,达人看看怎么解决?急求!
模型背景:重复测量计数数据。用泊松过程来模拟。数据集介绍:excel文件导入;有8个个体,但observation一共有432个,每个个体重复测量多次。
代码如下:
proc nlmixed data=Lfm.tfy2 qtol=0.2 fd=central;
parms beta1=0.7 beta2=-3 ,epthta1=0.07,epthta2=-0.1,epthta12=1.4;
x=log(t);
gamma=beta1+b1;
thta=exp(-(beta2+b2)/gamma);
eta=x*(beta1+b1)+beta2+b2;
mean=exp(eta);
model y ~ poisson(mean);
random b1 b2 ~ normal([0,0],[epthta1,epthta12,epthta2]) subject=group;
predict b1 out=p;
predict b2 out=d;
predict gamma out=f;
predict thta out=g;
run;
错误: Quadrature accuracy of 0.200000 could not be achieved with 31 points. The achieved accuracy was 1.000000.
模型并不复杂,数据集也不复杂。这个错误实在搞不定。哪位高手对NLMIXED过程中使用的参数估计方法原理熟的,讲讲原因及如何处理?
扫码加我 拉你入群
请注明:姓名-公司-职位
以便审核进群资格,未注明则拒绝
全部回复
扫码加我 拉你入群
请注明:姓名-公司-职位
以便审核进群资格,未注明则拒绝
扫码加我 拉你入群
请注明:姓名-公司-职位
以便审核进群资格,未注明则拒绝
相关推荐
栏目导航
热门文章
推荐文章
扫码加好友,拉您进群