Resources

Example 3 for PROC GAMPL

/****************************************************************/
/*          S A S   S A M P L E   L I B R A R Y                 */
/*                                                              */
/*    NAME: hpgamex3                                            */
/*   TITLE: Example 3 for PROC GAMPL                            */
/*    DESC: Mackerel egg survey data set                        */
/*     REF: Bowman and Azzalini (1997)                          */
/* PRODUCT: STAT                                                */
/*  SYSTEM: ALL                                                 */
/*    KEYS:                                                     */
/*   PROCS: GAMPL                                               */
/*                                                              */
/* SUPPORT: Weijie Cai                                          */
/****************************************************************/

title 'Mackerel Egg Density Study';
data Mackerel;
   input Egg_Count Longitude Latitude  Net_Area Depth Distance;
   datalines;
   0     -4.65     44.57     0.242    4342    0.8395141177
   0     -4.48     44.57     0.242    4334    0.8591926336
   0      -4.3     44.57     0.242    4286    0.8930152895
   1     -2.87     44.02     0.242    1438    0.3956408691
   4     -2.07     44.02     0.242     166    0.0400088237
   3     -2.13     44.02     0.242     460    0.0974234463
   0     -2.27     44.02     0.242     810    0.2362566569
   1     -2.35     44.02     0.242     887    0.3160524314
   0     -2.42     44.02     0.242     946    0.3859431039
   0     -2.48     44.02     0.242    1012    0.4254223762
   0     -2.48     44.12     0.242    1011    0.4612031784
   0     -2.33     44.12     0.242     969    0.3185483454
   0     -2.27     44.12     0.242     929    0.2622465816
   2      -2.2     44.12     0.242     847    0.1965611905
   0     -2.13     44.12     0.242     704    0.1308757994
   1     -2.07     44.12     0.242     266    0.0745740357
   3     -2.07     44.37     0.242     181    0.0737069183
   1     -2.13     44.37     0.242     467    0.1308511425
   0      -2.2     44.37     0.242     696    0.1979967228
   0     -2.35     44.37     0.242    1012    0.3454873183
   0     -2.35     44.62     0.242    1335    0.3423067093
   0     -2.27     44.62     0.242     922    0.2623067093
   0     -2.18     44.62     0.242     606    0.1723067093
   0     -2.13     44.62     0.242     539    0.1223067093
   2     -2.07     44.62     0.242     883    0.0623067093
   1     -2.07     44.87     0.242     135    0.0268741384
   0     -2.13     44.88     0.242     147    0.0395908028
   0      -2.2     44.87     0.242     268    0.0806910854
   1     -2.27     44.87     0.242     783    0.1304256852
   0     -2.35     44.87     0.242    1892    0.1937029682
   0     -2.42     44.87     0.242    2034    0.2299645979
   0     -2.48     44.87     0.242    2580    0.2492961559
   0     -2.55     44.87     0.242    2709     0.271376999
   0     -2.62     44.87     0.242    2681    0.2823344738
   0      -2.7     44.87     0.242    2879     0.311043345
   0     -2.77     44.87     0.242    2852    0.3199188326
   0     -2.83     44.87     0.242    3158    0.3239188078
   0      -2.9     44.87     0.242    2943    0.3376793441
   0     -2.97     44.87     0.242    3097    0.3526674424
   0     -3.12     44.98     0.242    2249    0.3063093929
   0     -3.27     45.12     0.242    1202    0.2996330778
   0      -3.2     45.12     0.242     964    0.2504312196
   0     -3.13     45.12     0.242     161    0.2088774779
   0     -3.05     45.12     0.242    1039    0.1534121987
   1     -2.98     45.12     0.242     242    0.1157074109
   2     -2.92     45.12     0.242    1159     0.098536924
   0     -2.85     45.12     0.242     155    0.0866991148
   4     -2.77     45.12     0.242     156    0.0704395178
   0      -2.7     45.12     0.242     306    0.0614082461
   1     -2.63     45.12     0.242     126    0.0568025763
   0     -2.55     45.12     0.242     127    0.0232507398
   0     -2.48     45.12     0.242     125    0.0146964787
   0     -2.42     45.12     0.242     126    0.0003367127
   0     -2.35     45.12     0.242     120    0.0294471811
   0     -2.27     45.12     0.242     120    0.0562913573
   1      -2.2     45.12     0.242     116    0.0907005137
   0     -2.13     45.12     0.242     111    0.1334280664
   0     -2.07     45.12     0.242     104    0.1807325361
   0     -1.98     45.12     0.242     100    0.2388023189
   0     -1.92     45.12     0.242      95     0.262181456
   0     -1.78     45.12     0.242      79    0.3414250436
   0      -1.7     45.12     0.242      71    0.3992517839
   2     -1.57     45.12     0.242      60    0.5062630171
   0     -1.48     45.12     0.242      52    0.5858216749
   0     -1.42     45.12     0.242      40    0.6403984929
   0     -1.35     45.12     0.242      31     0.705185997
   0     -1.28     45.12     0.242      26    0.7708854332
   0     -1.28     45.18     0.242      27    0.7927131334
   0     -1.28     45.25     0.242      27    0.8229897728
   0     -1.28      45.3     0.242      29    0.8474984563
   0     -1.28     45.37     0.242      29    0.8854217569
   0     -1.37     45.37     0.242      35    0.8130725031
   0     -1.42     45.37     0.242      41    0.7744788998
   0     -1.48     45.37     0.242      46    0.7299961447
   0     -1.57     45.37     0.242      51    0.6678394859
   0     -1.63     45.37     0.242      51    0.6301479065
   0     -1.72     45.37     0.242      56    0.5797984059
   0     -1.78     45.37     0.242      66    0.5486850531
   0     -1.85     45.37     0.242      74    0.5087696952
   0     -1.92     45.37     0.242      78    0.4586245543
   1     -1.98     45.37     0.242      82    0.4201861232
   0     -2.07     45.37     0.242      92    0.3721646287
   0     -2.13     45.37     0.242      97     0.342784142
   0      -2.2     45.37     0.242     103    0.3135403699
   0     -2.28     45.52     0.242     102    0.4233335102
   0     -2.28     45.62     0.242      99    0.5184856328
   0      -2.2     45.62     0.242      92    0.5450129543
   1     -2.12     45.62     0.242      83    0.5759822965
   4     -2.07     45.62     0.242      79    0.5937739873
   0     -1.35     45.43     0.242      31    0.8666743262
   4     -1.35      45.3     0.242      34     0.788210995
   0     -1.35     45.23     0.242      33     0.751903264
   0     -1.42     45.23     0.242      41    0.6915066578
   1     -1.48     45.23     0.242      50    0.6412942093
   0     -1.57     45.25     0.242      57    0.5825376018
   2     -1.63     45.23     0.242      63    0.5248335695
   3      -1.7     45.25     0.242      68     0.492391919
   3     -1.78     45.25     0.242      74    0.4464470693
   0     -1.85     45.25     0.242      77    0.4098291238
   0     -1.92     45.23     0.242      84    0.3628981859
   0     -1.98     45.25     0.242      90    0.3319200273
   1     -2.13     45.25     0.242     105    0.2382686636
   0     -2.28     45.37     0.242     108    0.2844304813
   0     -2.35     45.37     0.242     108    0.2590744601
   1     -2.48     45.37     0.242     112    0.2307815007
   0     -2.63     45.37     0.242     120    0.1892521557
   0      -2.7     45.37     0.242     124    0.1835552819
   0     -2.77     45.37     0.242     133    0.1740475352
   1     -2.85     45.37     0.242     130    0.1562860501
   1     -2.92     45.37     0.242     139    0.1168201072
   0     -2.98     45.37     0.242     136    0.0815351079
   1     -3.05     45.37     0.242     147    0.0288824077
   5     -3.13     45.37     0.242     158     0.016790237
   0      -3.2     45.37     0.242     178    0.0520579023
   0     -3.27     45.37     0.242     592    0.1161080566
   9     -3.35     45.37     0.242     941    0.1943987105
   2     -3.42     45.37     0.242    2570    0.2512595743
   0     -3.48     45.37     0.242    3200     0.291093129
   0     -3.55     45.37     0.242    3527    0.3417890342
   0     -3.62     45.37     0.242    3611    0.3886139984
   0      -3.7     45.37     0.242    3723    0.4367806449
   0     -3.83     45.37     0.242    3691    0.5179946065
   0     -3.98     45.37     0.242    3960    0.6130237854
   2     -4.12     45.52     0.242    3789    0.5976801872
   5     -4.27     45.62     0.242    3964    0.6042301961
   1      -4.2     45.62     0.242    4015    0.5667162652
   1     -4.13     45.62     0.242    3498    0.5295171301
   2     -4.05     45.62     0.242    2602    0.4804855719
   0     -3.98     45.62     0.242    2183    0.4257784423
   0     -3.92     45.62     0.242    1879    0.3805299042
   1     -3.85     45.62     0.242    2354    0.3336735257
   1     -3.78     45.62     0.242    1794    0.2951277304
   1      -3.7     45.62     0.242    2101    0.2532669769
   2     -3.63     45.62     0.242    1294    0.2003720056
   1     -3.55     45.62     0.242    1144    0.1468598219
   0     -3.48     45.62     0.242     419    0.1047140062
   0     -3.42     45.62     0.242     169    0.0685890213
   0     -3.35     45.62     0.242     148    0.0245552715
   1     -3.27     45.62     0.242     136     0.033385505
   1      -3.2     45.62     0.242     135    0.0716968847
   0     -3.13     45.62     0.242     139     0.102588002
   1     -3.07     45.62     0.242     133     0.147611337
   0     -2.98     45.62     0.242     133    0.2241105112
   1     -2.92     45.62     0.242     125    0.2762483085
   0     -2.77     45.62     0.242     121     0.403859887
   0      -2.7     45.62     0.242     121    0.4321757844
   0     -2.63     45.62     0.242     117    0.4391348184
   0     -2.57     45.62     0.242     114    0.4444086724
   0     -2.48     45.62     0.242     112    0.4648272796
   0     -2.42     45.62     0.242     106    0.4830717505
   0      -2.2     45.48     0.242      90    0.4181302423
   0     -2.07     45.52     0.242      82    0.5033368395
   1     -2.07     45.87     0.242      73    0.8261550822
   0     -2.13     45.87     0.242      77    0.8032848849
   5      -2.2     45.87     0.242      79    0.7809725901
   6     -2.28     45.87     0.242      84    0.7587625502
   0     -2.35     45.87     0.242      91    0.7419150329
   0     -2.42     45.87     0.242      98    0.7218872066
   0     -2.48     45.87     0.242      98    0.7071149818
   0     -2.57     45.87     0.242     101    0.6924753269
   0      -2.7     45.87     0.242     106    0.5941331548
   0     -2.85     45.87     0.242     114    0.4767936492
   1     -2.98     45.87     0.242     120    0.3934782059
   0     -3.13     45.87     0.242     123    0.3075686928
   0     -3.28     45.87     0.242     135    0.2153109709
   0     -3.35     45.87     0.242     138    0.1731651552
   0     -3.42     45.87     0.242     138     0.117052398
   0     -3.48     45.87     0.242     139    0.0776678955
   0     -3.57     45.87     0.242     137    0.0278001088
   1     -3.63     45.87     0.242     151    0.0035954563
   1      -3.7     45.87     0.242     180    0.0533279819
   0     -3.78     45.87     0.242     175    0.1171976208
   1     -3.85     45.87     0.242     934    0.1603654846
   1     -3.92     45.87     0.242     598     0.203268098
   1     -3.98     45.87     0.242    1003    0.2387012687
   1     -4.07     45.87     0.242    1972    0.2857047672
   0     -4.13     45.87     0.242    2443    0.3188313644
   0      -4.2     45.87     0.242    2323    0.3644993417
   3     -4.28     45.87     0.242    3024    0.4196773185
   0     -4.35     45.87     0.242    3463    0.4618545078
   0     -4.42     45.87     0.242    4097    0.5036019328
   0     -4.57     45.87     0.242    4142     0.579097532
   0     -4.63     45.87     0.242    3640    0.6019552966
   0      -4.8     46.02     0.242    3496    0.5597152007
   1     -4.95     46.12     0.242    3883    0.5075978088
   0     -4.88     46.12     0.242    2765    0.5060715332
   0     -4.73     46.12     0.242    2254    0.4377232045
   1     -4.58     46.12     0.242    2208    0.3536401867
   2     -4.52     46.12     0.242    3518    0.3292787693
   0     -4.45     46.12     0.242    1581    0.3017949176
  18     -4.37     46.12     0.242    1346    0.2754879266
   8      -4.3     46.12     0.242    1698    0.2441731581
   6     -4.23     46.12     0.242    1200    0.1903899909
   1     -4.15     46.12     0.242    1220    0.1425503662
   0     -4.08     46.12     0.242     189     0.100714452
   1     -4.02     46.12     0.242     154    0.0593309693
   0     -3.93     46.12     0.242     151    0.0002530401
   0     -3.87     46.12     0.242     148    0.0319942153
   7      -3.8     46.12     0.242     152    0.0684702687
   0     -3.72     46.12     0.242     147    0.1168519893
   0     -3.65     46.12     0.242     141    0.1597546027
   0     -3.58     46.12     0.242     138    0.2029447134
   0     -3.48     46.12     0.242     137    0.2691384918
   0     -3.43     46.12     0.242     133     0.308403784
   0     -3.37     46.12     0.242     133    0.3455591358
   0     -3.28     46.12     0.242     128    0.3973143043
   0     -3.15     46.12     0.242     121    0.4847632693
   0     -2.98     46.12     0.242     117    0.5955562612
   0     -2.92     46.12     0.242     111    0.6337804213
   0     -2.72     46.12     0.242     102    0.7491849225
   0     -2.57     46.12     0.242      94    0.8520211674
   0     -2.48     46.12     0.242      84    0.9199659139
   0     -2.42     46.08     0.242      71    0.9252948749
   0     -2.35     46.08     0.242      64     0.942760965
   2     -2.28     46.08     0.242      51    0.9632590156
   2     -2.22     46.12     0.242      56    1.0158612029
   0     -2.13     46.12     0.242      58    1.0403837629
   0     -2.07     46.12     0.242      54     1.058953803
   0     -2.07     46.37     0.242      36    1.2967971616
   0     -2.22     46.37     0.242      41    1.2575435305
   0     -2.35     46.37     0.242      54    1.1838425346
   0     -2.48     46.37     0.242      67    1.0949776534
   0     -2.57     46.37     0.242      73    1.0385465076
   0     -2.63     46.37     0.242      83    1.0036503573
   0     -2.72     46.37     0.242      83    0.9536377808
   0     -2.87     46.37     0.242     100    0.8601297509
   0     -2.98     46.37     0.242     108      0.78150727
   0     -3.15     46.37     0.242     113    0.6763897421
   0     -3.22     46.37     0.242     117    0.6296923646
   0     -3.28     46.37     0.242     119    0.5871112624
   0     -3.37     46.37     0.242     119    0.5292610219
   0     -3.43     46.37     0.242     124    0.4922171758
   1     -3.52     46.37     0.242     124    0.4369720766
   0     -3.58     46.37     0.242     127     0.400198408
   0     -3.65     46.37     0.242     133    0.3600092117
   0     -3.73     46.37     0.242     136    0.3183222936
   0      -3.8     46.37     0.242     138    0.2734263915
   0     -3.87     46.37     0.242     139    0.2251456617
   0     -3.95     46.37     0.242     139    0.1645040131
   0     -4.02     46.37     0.242     141    0.0999022176
   0     -4.08     46.37     0.242     145    0.0557048402
   0     -4.17     46.37     0.242     148     0.009411632
   0     -4.23     46.37     0.242     154    0.0024671392
   0      -4.3     46.37     0.242     151    0.0163257056
   0     -4.38     46.37     0.242     153    0.0468258082
   1     -4.45     46.37     0.242     166    0.0720718139
   1     -4.53     46.37     0.242     360    0.1178848997
   5      -4.6     46.37     0.242     690    0.1670144382
  12     -4.67     46.37     0.242     914    0.2013499642
  13     -4.73     46.37     0.242    1317    0.2339516795
   0     -4.82     46.37     0.242    1330     0.249429604
   0     -4.88     46.37     0.242    1699    0.2575978088
   0     -4.95     46.37     0.242    2198    0.2575978088
   0     -5.03     46.37     0.242    2402    0.2664339843
   0      -5.1     46.37     0.242    2580    0.2921852719
   0     -5.17     46.37     0.242    2866    0.3019951109
   0     -5.25     46.37     0.242    3291    0.3186299755
   0     -5.32     46.37     0.242    3370    0.3446143006
  16     -5.47     46.52     0.242    3198    0.3327686065
   1     -5.62     46.62     0.242    3793    0.4284737606
   1     -5.55     46.62     0.242    3225    0.3668810962
   8     -5.47     46.62     0.242    2776    0.2900961236
  14      -5.4     46.62     0.242    2523    0.2247566026
   2     -5.33     46.62     0.242    2867    0.1621795694
   0     -5.25     46.62     0.242    1721    0.0949067619
   6     -5.18     46.62     0.242    1985    0.0602167335
   4      -5.1     46.62     0.242    1841    0.0466900789
   2     -5.03     46.62     0.242    1152    0.0375872025
  11     -4.97     46.62     0.242     667    0.0104181925
   0     -4.67     46.62     0.242     158    0.0188133483
   1     -4.88     46.62     0.242     388    0.0075978088
  16     -4.82     46.62     0.242     271    0.0075978088
   0     -4.73     46.62     0.242     189    0.0157147106
   2      -4.6     46.62     0.242     154    0.0242108657
   0     -4.53     46.62     0.242     150    0.0526485479
   0     -4.45     46.62     0.242     145     0.114848174
   2     -4.38     46.62     0.242     145    0.1597245396
   0     -4.32     46.62     0.242     143    0.2061424875
   0     -4.23     46.62     0.242     144    0.2393055029
   0     -4.17     46.62     0.242     141    0.2544631712
   0      -4.1     46.62     0.242     138    0.2684337243
   0     -4.02     46.62     0.242     136    0.2969379548
   1     -3.95     46.62     0.242     135    0.3350550567
   0     -3.88     46.62     0.242     132    0.3754848747
   0      -3.8     46.62     0.242     127    0.4286142415
   0     -3.73     46.62     0.242     121    0.4812393222
   0     -3.67     46.62     0.242     125    0.5295585358
   0     -3.58     46.62     0.242     120     0.605965081
   0     -3.52     46.62     0.242     113    0.6411264252
   0     -3.45     46.62     0.242     117    0.6785672801
   0     -3.37     46.62     0.242     118    0.7264474401
   0      -3.3     46.62     0.242     115    0.7693500535
   0     -3.22     46.62     0.242     108    0.8185334843
   0     -3.15     46.62     0.242     108    0.8617513047
   0     -3.08     46.62     0.242     102    0.9066221601
   0     -2.98     46.62     0.242     100    0.9758680617
   0     -2.93     46.63     0.242      92    1.0157715232
   0     -2.87     46.62     0.242      86    1.0434127718
   0     -2.78     46.62     0.242      81    1.1010506764
   0     -2.65     46.62     0.242      68    1.1914064642
   0     -2.57     46.62     0.242      61     1.241979781
   0     -2.48     46.62     0.242      58    1.2948416274
   0     -2.42     46.62     0.242      51    1.3288640445
   0     -2.48     46.75     0.242      39    1.3997992726
   0     -2.48     46.87     0.242      27    1.4883573323
   0     -2.65     46.87     0.242      44    1.3759856056
   0     -2.72     46.87     0.242      51     1.333267116
   0      -2.8     46.87     0.242      53    1.2801998414
   0     -2.87     46.87     0.242      57     1.233030138
   0     -2.93     46.87     0.242      50    1.1943828874
   0     -2.97     44.37     0.242    1420    0.7456408691
   0      -2.9     44.37     0.242    1427    0.7456408691
   0     -2.83     44.37     0.242    1467    0.7456408691
   0     -2.68     44.37     0.242    1383    0.6716074698
   0     -2.62     44.37     0.242    1349    0.6123128968
   0     -2.55     44.37     0.242    1348     0.543135895
   0     -2.48     44.37     0.242    1328    0.4739588931
   0     -2.42     44.37     0.242    1492    0.4146643202
   0     -2.42     44.62     0.242    1735    0.4123067093
   0     -2.48     44.62     0.242    2080     0.453205233
   0     -2.62     44.62     0.242    2634    0.5280940138
   0     -2.77     44.62     0.242    2280    0.5662761939
   0      -2.9     44.62     0.242    2262    0.5811964535
   0     -4.33     45.48     0.242    4173    0.7543788285
   0     -4.35     45.58     0.242    3943    0.6828151146
   0     -4.42     45.62     0.242    4102    0.6972639723
   0     -4.48     45.62     0.242    4231    0.7386401851
   0     -4.63     45.62     0.242    4360    0.8300959169
   0     -4.77     45.62     0.242    4390    0.8866168813
   0     -4.83     45.62     0.242    4414    0.9125366692
  10     -5.75     48.25      0.15     124    0.9227330417
   4     -6.25     48.25      0.15     140    0.6897133782
   6      -6.8     48.23      0.15     160    0.5694966125
   8     -7.27     48.25      0.15     173    0.3675161731
   1     -7.78     48.25      0.15     184     0.125272891
  20     -8.23     48.27      0.15     206    0.0116115201
  20     -8.75     48.23      0.15     192    0.0877499994
  26     -9.25     48.25      0.15     256    0.1541522464
  22     -9.75     48.23      0.15     875    0.4171101972
  76    -10.27     48.25      0.15    3000    0.5375307864
  27    -10.77     48.25      0.15    2700    0.7401461724
  22    -11.25     48.25      0.15    4000    0.8738787306
  27    -11.73     48.23      0.15    2800     1.153986225
  51    -12.73     48.73      0.15    1600    1.6615010232
  21    -12.25     48.73      0.15    1300    1.2186671569
  29    -11.87     48.73      0.15    1400    0.8894291719
  81    -11.62     48.75      0.15    2600    0.6857238367
  65    -11.37     48.75      0.15    1600    0.5293366535
  53    -11.12     48.73      0.15    1000    0.3914502223
  61    -10.87     48.75      0.15    1100    0.2694969977
  28    -10.62     48.73      0.15    1400    0.2622425884
  23    -10.37     48.73      0.15     162    0.1179571067
  30    -10.07     48.75      0.15     425    0.0157432556
  45     -9.85     48.75      0.15     193     0.025306547
  43      -9.6     48.75      0.15     176     0.119013888
  51     -9.37     48.75      0.15     175    0.2067097778
  17     -9.13     48.75      0.15     165     0.329107796
   2     -8.87     48.73      0.15     170    0.3988589266
   6     -8.63     48.73      0.15     146    0.4214112795
  10     -8.23     48.75      0.15     153    0.4462789579
   5     -7.75     48.75      0.15     158    0.6032527898
  23     -7.23     48.75      0.15     153    0.8105988408
  10     -6.77     48.75      0.15     145      1.07455991
   3     -6.27     48.75      0.15     137    1.1596222107
   3     -6.73     49.25      0.15     126    1.4904395437
   7     -7.22     49.25      0.15     135    1.2631817058
  10     -7.72     49.25      0.15     138    1.0828571439
  15     -8.25     49.25      0.15     127    0.9321593173
  18     -8.63     49.25      0.15     152    0.9390438736
   4     -8.88     49.25      0.15     154     0.887526452
   1     -9.12     49.25      0.15     151    0.7629012711
  13     -9.38     49.25      0.15     153    0.6526316781
   9     -9.65     49.23      0.15     154    0.5327590071
  10     -9.88     49.25      0.15     132     0.499412191
  26    -10.15     49.25      0.15     162    0.3409200497
  33    -10.37     49.25      0.15     136    0.2579881877
  40    -10.63     49.25      0.15     163    0.2573486008
  40    -10.88     49.23      0.15     177    0.1556872456
  57    -11.08     49.25      0.15     187    0.0439525458
  41    -11.37     49.23      0.15     350    0.2139985205
  37    -11.73     49.25      0.15     830    0.5692610808
  32    -12.25     49.25      0.15    1100       1.0889355
  65    -12.75     49.25      0.15    1200      1.58882337
  55    -12.73     49.73      0.15    2200    1.6324447345
  57    -12.25     49.73      0.15    1200    1.1787450108
  59    -11.75     49.75      0.15     780    0.7545504793
  59    -11.33     49.73      0.15     450    0.3765631369
  43     -11.1     49.73      0.15     320    0.1580631431
  30    -10.83     49.75      0.15     175    0.0496397168
  24     -10.6     49.75      0.15     147    0.1960982386
  12    -10.35     49.75      0.15     143    0.4192402102
  45    -10.08     49.75      0.15     141    0.6815092188
  40     -9.85     49.73      0.15     131    0.9051489134
  33     -9.62     49.73      0.15     158     1.030266117
  21     -9.33     49.73      0.15     118    1.1073121557
   4     -9.12     49.75      0.15     142    1.2157726486
   6     -8.85     49.75      0.15     126    1.3310745619
   3      -8.6     49.75      0.15     136    1.4377120457
   0      -8.2     49.75      0.15     138    1.4344814083
   3     -7.72     49.73      0.15     133    1.5342936614
   5     -7.73     50.25      0.15     103    2.0258870905
   0     -8.23     50.25      0.15     129    1.9293283944
  35     -8.73     50.25      0.15     126      1.83164208
  60     -9.23     50.25      0.15     140    1.4756778835
  46     -9.75     50.25      0.15     126    0.9658144214
  19     -10.1     50.23      0.15     137    0.6396928347
   9    -10.37     50.23      0.15     147    0.3980227523
   8    -10.62     50.23      0.15     164    0.1488695265
   5    -10.87     50.23      0.15     208     0.052118011
  56    -11.13     50.23      0.15     550    0.2380892631
  34    -11.38     50.23      0.15    1100    0.4783387553
  31    -11.73     50.25      0.15    1600    0.8269773693
  31    -12.25     50.25      0.15    2300    1.3439881777
  28    -12.75     50.23      0.15    2400    1.8415845371
  30    -12.75     50.75      0.15    2300    1.7069865241
  35    -12.25     50.73      0.15    2100    1.2328347626
  32    -11.87     50.75      0.15    2000    0.8753594353
  23     -11.6     50.73      0.15    1800    0.6519149416
  42    -11.37     50.75      0.15    1300    0.4537462859
  31    -11.12     50.73      0.15     330    0.2509354095
  28    -10.85     50.73      0.15     181    0.0218147471
  18    -10.55     50.73      0.15     166    0.2238786395
  10    -10.23     50.75      0.15     140    0.5199212285
   1     -9.72     50.75      0.15     123    1.0027754255
  20     -9.23     50.73      0.15     130    1.4775126146
  18     -8.73     50.73      0.15     105    1.9723065827
   0     -8.25     50.73      0.15     118    2.4071712185
   1     -8.25     51.25      0.15     101    2.5543951991
   2     -8.75     51.23      0.15     105    2.0742430025
   1     -9.25     51.25      0.15     110    1.6168704133
   2     -9.75     51.25      0.15     108    1.1773817858
   6    -10.25     51.23      0.15     136    0.7520089099
  12    -10.63     51.23      0.15     165    0.4279458648
  17    -10.88     51.25      0.15     178    0.2277079773
  58     -11.3     51.23      0.15     211    0.1922920227
  35    -11.65     51.23      0.15     610    0.5422920227
  78    -11.88     51.25      0.15    1100    0.7722920227
  35    -12.12     51.25      0.15    1200    1.0122920227
  29    -12.38     51.25      0.15    1300    1.2688766929
  69    -12.75     51.23      0.15    1500    1.5990568921
  14    -14.73     51.75      0.15     470    1.6548107034
  22    -14.35     51.73      0.15     375    1.5657689486
  29    -14.08     51.73      0.15     380     1.532859219
  54    -13.85     51.75      0.15     450    1.5175848241
  84    -13.62     51.73      0.15     590    1.5535554343
  46    -13.37     51.73      0.15     790    1.5816110843
  67    -13.12     51.73      0.15    1000    1.4404489283
  66    -12.88     51.73      0.15    1200    1.2811801214
  53    -12.62     51.75      0.15    1300    1.1187514574
  50    -12.38     51.75      0.15    1300    1.0045057132
  37    -12.12     51.73      0.15    1000    0.7998552594
  43    -11.88     51.75      0.15     710    0.5857323005
  36    -11.62     51.73      0.15     350     0.370994499
  46    -11.35     51.73      0.15     210    0.1361486674
  15     -11.1     51.75      0.15     183    0.0911692879
  24    -10.73     51.75      0.15     133    0.4272116779
   1     -8.23     51.23     0.222     101    2.5676033688
   0     -8.73     51.25     0.222     104    2.0998144654
   0     -9.23     51.25     0.222     105    1.6351086518
   2     -9.73     51.25     0.222     102    1.1945402346
  14    -10.23     51.25     0.222     133    0.7790685494
  23    -10.72     51.25     0.222     171    0.3647873377
  12    -11.22     51.23     0.222     203    0.1122920227
  78    -11.72     51.25     0.222    1161    0.6122920227
  90    -12.18     51.25     0.222    1645    1.0722920227
  77     -12.7     51.23     0.222    1885    1.5566992609
  33     -14.7     51.73     0.222     481    1.6622082011
  63    -14.28     51.75     0.222     377    1.5330746563
  46    -13.78     51.73     0.222     509    1.5429208526
  77    -13.28     51.73     0.222     912    1.5581619751
  41    -12.78     51.73     0.222    1085    1.2226241241
  73    -12.28     51.75     0.222     996    0.9253871512
  56    -11.77     51.73     0.222     523    0.5014644054
  16    -11.28     51.73     0.222     195     0.075262711
  11    -10.78     51.75     0.222     138    0.3802056563
   1     -10.3     51.75     0.222      68    0.8359314684
   3     -10.7     52.23     0.222     117    0.6759237064
   2     -11.2     52.25     0.222     139    0.2560062999
  26     -11.7     52.23     0.222     306     0.173370354
  66    -12.23     52.25     0.222     749    0.4886975668
  52    -12.72     52.25     0.222     689    0.7945254084
  76     -13.2     52.25     0.222     572    1.1112249253
  51     -13.7     52.23     0.222     387    1.0493084143
  37     -14.2     52.23     0.222     334    1.0464541627
  43    -14.68     52.25     0.222     436    1.1825388824
  73    -14.78     52.73     0.222     559    0.8758632468
  24    -14.28     52.75     0.222     332    0.5597685007
  47    -13.78     52.75     0.222     234    0.5251552118
  53    -13.28     52.73     0.222     347    0.6369096605
  57     -12.8     52.73     0.222     487    0.6172359158
  41    -12.28     52.75     0.222     337    0.1365312648
   1    -11.78     52.75     0.222     168    0.1552061145
   4    -11.27     52.73     0.222     130    0.4742166919
   1    -11.25     53.22     0.222     143    0.6238153661
   4    -11.68     53.25     0.222     182    0.3563155715
  10     -12.2     53.23     0.222     308    0.0094483276
  33     -12.7     53.25     0.222     338    0.4857796791
  32     -13.2     53.23     0.222     236    0.2341923116
  50    -13.68     53.23     0.222     171     0.055459343
  47     -14.2     53.25     0.222     208    0.1030179144
  53     -14.7     53.25     0.222     949    0.5948816159
 113    -14.75     53.72     0.222    1500    0.7661087538
  27    -14.28     53.75     0.222    1358     0.424455318
  44     -13.8     53.75     0.222     418     0.190150346
  37    -13.28     53.75     0.222     303    0.1171377308
  19    -12.77     53.73     0.222     283    0.4356861554
  21    -12.28     53.73     0.222     333    0.3994867254
  17    -11.77     53.73     0.222     295     0.082376324
  13     -11.3     53.75     0.222     205    0.1453250442
  21    -10.78     53.75     0.222     148    0.4245712891
   0    -10.27     53.75     0.222      74    0.6732115349
   0    -10.27     53.27     0.222     107    1.0927471177
   1    -10.72      53.2     0.222     119     0.928502272
   0    -10.75     52.75     0.222     123    0.9018657597
   2    -10.27     52.73     0.222     100    1.2989060637
   0    -10.22     54.25     0.222      92    0.3025224369
  58    -10.72     54.25     0.222     214    0.0175548779
  62     -11.2     54.25     0.222     340    0.2170562127
  34    -11.73     54.25     0.222    1172    0.5054580427
  49     -12.2     54.23     0.222    2030    0.7386455951
  24     -12.7     54.23     0.222    2200    0.8056505265
  16    -12.27     54.73     0.222    2000    1.1964552871
  10    -11.78     54.75     0.222    2873    0.9476267756
  13     -11.3     54.75     0.222    2900    0.7019613094
  43    -10.78     54.75     0.222    1750    0.4384999915
  42    -10.28     54.75     0.222     120    0.0331852024
  27     -9.73     54.73     0.222     100    0.3910312491
   0     -9.28     54.75     0.222      91    0.6847841032
   0     -8.25     55.23     0.222      32    1.0048444388
   1      -8.7     55.23     0.222      93    0.6555891558
   8     -9.22     55.25     0.222     109    0.2937613708
  27     -9.72     55.25     0.222     128    0.0318280984
  32    -10.18     55.23     0.222    1291    0.2745622713
  19    -10.65     55.23     0.222    2480    0.6214210959
  42    -10.27     55.73     0.222    2300    0.6827875476
  40     -9.78     55.75     0.222    1785    0.4312701461
  30     -9.27     55.73     0.222     140    0.1006893438
   7     -8.77     55.73     0.222     106    0.3264369189
   3     -8.28     55.73     0.222     104    0.8082043186
   1     -8.23     56.25     0.222     167    0.6484234451
   2     -8.72     56.23     0.222     131    0.2596481931
  27      -9.2     56.25     0.222     711    0.1909803452
 106     -9.68     56.23     0.222    1595    0.6582792428
  22     -9.77     56.75     0.222    1795    0.6673417415
  70     -9.28     56.75     0.222    1260    0.3117001765
   7     -8.78     56.73     0.222     117    0.0667314293
   4     -8.73     57.23     0.222     134       0.3715636
  37     -9.18     57.23     0.222     174     0.082095629
  37      -9.7     57.25     0.222    1940    0.3656364314
  49     -9.77     57.75     0.222    1088    0.3920204076
   3     -9.27     57.75     0.222     144    0.0873164841
   6     -8.78     57.75     0.222     146    0.3649476732
   3     -8.28     57.73     0.222     111    0.5841104348
   0     -8.23     57.27     0.222     133    0.8039890158
   1     -8.25     56.77     0.222     132    0.5666356525
   1     -5.75     48.25      0.19     120    0.9227330417
  10     -6.23     48.25      0.19     139    0.6979037904
  10     -6.75     48.25      0.19     155    0.5913811675
   3     -7.25     48.25      0.19     172    0.3821368185
  10     -7.75     48.25      0.19     177    0.1366381488
  16      -7.8     48.25      0.19     178    0.1176960525
  11     -8.25     48.25      0.19     205    0.0357602758
 140     -8.75     48.25      0.19     200    0.0678410534
 122      -8.8     48.25      0.19     190    0.0726067329
  24     -9.23     48.25      0.19     560    0.1428536135
  14     -9.25     48.25      0.19     520    0.1541522464
  41     -9.75     48.25      0.19     775    0.3980794402
  22     -9.77     48.25      0.19     906    0.4042300733
  75    -10.23     48.25      0.19    2000    0.5276628557
  44    -10.23     48.25      0.19    2000    0.5276628557
  87    -10.73     48.25      0.19    1800     0.716775297
  57    -10.75     48.25      0.19    1900    0.7282798854
  72    -11.25     48.25      0.19    3600    0.8738787306
  38    -11.73     48.25      0.19    3000    1.1389598949
   8    -10.77     47.75      0.19    2000    1.1697823809
   7    -10.25     47.75      0.19    2000    1.0242227267
   5     -9.75     47.75      0.19    4450    0.8437334395
  15     -9.23     47.75      0.19    3500    0.6031973168
  34     -8.73     47.75      0.19    2500    0.5636584298
  50     -8.27     47.75      0.19    1560    0.5147870417
  72     -8.23     47.75      0.19    1350    0.5011520938
  52     -7.77     47.75      0.19     900    0.3179692814
  50     -7.75     47.75      0.19     770    0.3066726343
  12     -7.25     47.75      0.19     176     0.017722909
   6     -6.75     47.75      0.19     167    0.0946135503
   2     -6.27     47.75      0.19     152    0.2511573132
   2     -5.75     47.75      0.19     134    0.4343295287
   0     -5.23     47.75      0.19     125    0.6090585857
   1     -4.75     47.75      0.19     100    0.9150075624
   2     -4.73     47.25      0.19     122    0.6274923845
   0     -5.25     47.25      0.19     140    0.2096735702
   0     -5.73     47.25      0.19     200    0.0660511438
   0     -5.75     47.25      0.19     224    0.0650769018
   9     -6.23     47.25      0.19    1314    0.1760703613
  28     -6.28     47.25      0.19    1115    0.1976657775
 133     -6.75     47.25      0.19    2300    0.4015999827
  53     -7.25     47.25      0.19    2500    0.4495331185
  40     -7.32     47.25      0.19    2500    0.4666956288
  15     -7.73     47.25      0.19    3000    0.6422156047
  18     -8.25     47.25      0.19    3000    0.9769665902
   0     -8.73     47.25      0.19    3000    1.0610467794
   4     -9.25     47.25      0.19    3000    1.0987471315
   3     -9.73     47.25      0.19    3000    1.2601360266
   2    -10.23     47.25      0.19    3000    1.5026291941
  10     -8.25     46.75      0.19    4000    1.3529101749
   1     -7.75     46.75      0.19    4000    1.0985733569
  23     -7.23     46.75      0.19    4000    0.9368645063
  24     -6.73     46.75      0.19    2000    0.8485726594
  19     -6.25     46.75      0.19    2000     0.640094992
   5     -6.18     46.75      0.19    2000    0.6146546419
  20     -5.78     46.75      0.19    2000    0.4411973968
  38     -5.73     46.75      0.19    2000    0.4229382941
  75     -5.22     46.75      0.19    1000    0.0233752426
  32     -5.12     46.75      0.19     800    0.0759745133
   4     -4.77     46.75      0.19     163    0.1286761194
   3     -4.72     46.75      0.19     158     0.143225385
   4      -4.3     46.75      0.19     138    0.3103123518
   5     -4.27     46.75      0.19     142    0.3297899965
   3     -3.78     46.75      0.19     128    0.5394946831
   4     -3.75     46.75      0.19     124    0.5581890288
   0     -3.25     46.75      0.19     107    0.9027161972
   0     -2.75     46.75      0.19      65    1.2198092629
   0     -2.25     46.75      0.19       8    1.5348188074
   0     -2.23     46.25      0.19      41    1.1403006583
   0     -2.75     46.25      0.19      99    0.8375665995
   3     -3.25     46.25      0.19     123    0.5206434916
   1     -3.72     46.25      0.19     133     0.221112692
   0     -3.75     46.27      0.19     136    0.2225501754
  22     -4.22     46.25      0.19     205    0.1181120828
  21     -4.28     46.25      0.19     237     0.129990854
   0     -4.73     46.25      0.19    2500    0.3346500536
   5     -4.25     47.23      0.19     114     0.718192582
   3     -3.72     47.25      0.19     100    0.9944669778
   0     -3.25     47.25      0.19      64    1.2639918431
;

data Temp;
   set Mackerel;
   Density = Egg_Count/Net_Area;
run;

%let off0 = offsetmin=0 offsetmax=0
            linearopts=(thresholdmin=0 thresholdmax=0);
%let off0 = xaxisopts=(&off0) yaxisopts=(&off0);
proc template;
   define statgraph surface;
      dynamic _title _z;
      begingraph / designwidth=defaultDesignHeight;
         entrytitle _title;
         layout overlay / &off0;
            contourplotparm z=_z y=Latitude x=Longitude / gridded=FALSE;
         endlayout;
      endgraph;
   end;
run;

proc sgrender data=Temp template=surface;
   dynamic _title='Mackerel Egg Density' _z='Density';
run;

data Mackerel;
   set Mackerel;
   Log_Net_Area = log(Net_Area);
run;

ods graphics on;
proc gampl data=Mackerel plots seed=2345;
   model Egg_Count = spline(Depth) spline(Distance)
                     spline(Longitude Latitude/maxdf=40)
                      / offset=Log_Net_Area dist=poisson;
   output out=MackerelOut1 pred=Pred1 pearson=Pea1;
run;

proc gampl data=Mackerel plots seed=2345;
   model Egg_Count = spline(Depth) spline(Distance)
                     spline(Longitude Latitude/maxdf=40)
                      / offset=Log_Net_Area dist=negbin;
   output out=MackerelOut2 pred=Pred2 pearson=Pea2;
run;

data SurfaceFit;
   set Mackerel(keep=Longitude Latitude Net_Area Egg_Count);
   set MackerelOut1(keep=Pred1);
   set MackerelOut2(keep=Pred2);
   Pred1 = Pred1 / Net_Area;
   Pred2 = Pred2 / Net_Area;
run;

%let eopt = location=outside valign=top textattrs=graphlabeltext;
proc template;
   define statgraph surfaces;
      begingraph / designheight=360px;
         layout lattice / columns=2;
            layout overlay / &off0;
               entry "Poisson Model" / &eopt;
               contourplotparm z=Pred1 y=Latitude x=Longitude /
                               gridded=FALSE;
            endlayout;
            layout overlay / &off0;
               entry "Negative Binomial Model" / &eopt;
               contourplotparm z=Pred2 y=Latitude x=Longitude /
                               gridded=FALSE;
            endlayout;
         endlayout;
      endgraph;
   end;
run;

proc sgrender data=SurfaceFit template=surfaces;
run;