Worksheet size: 40000000 cells MTB > Set c1 DATA> 1( 1 : 25 / .5 )1 DATA> End. MTB > let c2 = exp(1 + .1*c1) MTB > describe c1-c2 Descriptive Statistics Variable N Mean Median TrMean StDev SE Mean C1 49 13.00 13.00 13.00 7.14 1.02 C2 49 12.66 9.97 12.21 8.64 1.23 Variable Minimum Maximum Q1 Q3 C1 1.00 25.00 6.75 19.25 C2 3.00 33.12 5.34 18.64 ====================== GMACRO poissonregr DO K1 = 1 : 50 let k2 = c2(k1) random 1 c10; poisson k2. let c11(k1) = c10 erase c10 erase k2 ENDDO ENDMACRO ======================== MTB > %F:\drr99\ed257\web\exs\poissonregr.mac Executing from file: F:\drr99\ed257\web\exs\poissonregr.mac * ERROR * One or more arguments too large. * Macro exiting... MTB > print c11 Data Display C11 1 3 4 4 2 8 4 7 5 8 5 5 10 4 10 10 7 8 12 13 9 9 9 6 13 13 7 12 10 16 18 18 12 12 18 11 21 22 15 23 19 16 31 27 23 30 22 28 36 MTB > describe c11 Descriptive Statistics Variable N Mean Median TrMean StDev SE Mean C11 49 12.98 11.00 12.58 8.33 1.19 Variable Minimum Maximum Q1 Q3 C11 1.00 36.00 7.00 18.00 MTB > print c1 c2 c11 Data Display Row C1 C2 C11 1 1.0 3.0042 1 2 1.5 3.1582 3 3 2.0 3.3201 4 4 2.5 3.4903 4 5 3.0 3.6693 2 6 3.5 3.8574 8 7 4.0 4.0552 4 8 4.5 4.2631 7 9 5.0 4.4817 5 10 5.5 4.7115 8 11 6.0 4.9530 5 12 6.5 5.2070 5 13 7.0 5.4739 10 14 7.5 5.7546 4 15 8.0 6.0496 10 16 8.5 6.3598 10 17 9.0 6.6859 7 18 9.5 7.0287 8 19 10.0 7.3891 12 20 10.5 7.7679 13 21 11.0 8.1662 9 22 11.5 8.5849 9 23 12.0 9.0250 9 24 12.5 9.4877 6 25 13.0 9.9742 13 26 13.5 10.4856 13 27 14.0 11.0232 7 28 14.5 11.5883 12 29 15.0 12.1825 10 30 15.5 12.8071 16 31 16.0 13.4637 18 32 16.5 14.1540 18 33 17.0 14.8797 12 34 17.5 15.6426 12 35 18.0 16.4446 18 36 18.5 17.2878 11 37 19.0 18.1741 21 38 19.5 19.1060 22 39 20.0 20.0855 15 40 20.5 21.1153 23 41 21.0 22.1980 19 42 21.5 23.3361 16 43 22.0 24.5325 31 44 22.5 25.7903 27 45 23.0 27.1126 23 46 23.5 28.5027 30 47 24.0 29.9641 22 48 24.5 31.5004 28 49 25.0 33.1155 36 options linesize = 80 pagesize=80; data poisdat; infile 'F:\drr99\ed257\web\exs\poissonregr.dat'; input pred mean outcom; run; proc genmod data=poisdat; model outcom = pred / dist = poi link = log; title 'SAS poisson regression artificial data'; run; SAS poisson regression artificial data 1 11:27 Monday, May 17, 1999 The GENMOD Procedure Model Information Description Value Data Set WORK.POISDAT Distribution POISSON Link Function LOG Dependent Variable OUTCOM Observations Used 49 Criteria For Assessing Goodness Of Fit Criterion DF Value Value/DF Deviance 47 40.2547 0.8565 Scaled Deviance 47 40.2547 0.8565 Pearson Chi-Square 47 38.9313 0.8283 Scaled Pearson X2 47 38.9313 0.8283 Log Likelihood . 1098.8393 . Analysis Of Parameter Estimates Parameter DF Estimate Std Err ChiSquare Pr>Chi INTERCEPT 1 1.2775 0.1128 128.3121 0.0001 PRED 1 0.0854 0.0062 188.6911 0.0001 SCALE 0 1.0000 0.0000 . . NOTE: The scale parameter was held fixed.