SAS 9.4 I am getting a "Termination due to floating point exception" error with the following program. Any ideas on what I'm doing wrong? (note: could not cut and paste so had to type in program, hopefully no typos) blk = 4 blocks (random) trt = 5 treatments (coded as dummy: trt1 trt2 trt3 trt4) cmp = 20 sites, 1 per trt per block time = time since trt (1, 3, 5, 8, 12, 16) abun=count poisson distribution (and eventually ZIP) proc sort; by blk; proc nlmixed; parms intercept=0 b1=0 b2=0 b3=0 b4=0 b5=0 b6=0 b7=0 b8=0 b9=0 varblk=0 varcomp=0; array blkE {4}; array cmpE {20}; eta = intercept + varblk*blkE{blk} + varcmp*cmpE{cmp} + b1*time + (b2 + b6*time)*trt1 + (b3 + b7*time)*trt2 + (b4 + b8*time)*trt3 + (b5 + b9*time)*trt4; mu=exp(eta); model abun ~ poisson(mu); random blkE1 blkE2 blkE3 blkE4 cmpE1 cmpE2 cmpE3 cmpE4 cmpE5 cmpE6 cmpE7 cmpE8 cmpE9 cmpE10 cmpE11 cmpE12 cmpE13 cmpE14 cmpE15 cmpE16 cmpE17 cmpE18 cmpE19 cmpE20 ~ normal([0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0], [varblk, 0, varblk, 0, 0, varblk, 0, 0, 0, varblk, 0, 0, 0, 0, varcmp, 0, 0, 0, 0, 0, varcmp, 0, 0, 0, 0, 0, 0, varcmp, 0, 0, 0, 0, 0, 0, 0, varcmp, 0, 0, 0, 0, 0, 0, 0, 0, varcmp, 0, 0, 0, 0, 0, 0, 0, 0, 0, varcmp, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, varcmp, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, varcmp, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, varcmp, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, varcmp, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, varcmp, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, varcmp, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, varcmp, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, varcmp, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, varcmp, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, varcmp, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, varcmp, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, varcmp 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, varcmp, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, varcmp]) subject=blk; run; Thanks for the help!
... View more