
SUMMARY OF SOME DATA-BASED DISCUSSION OF DESCRIPTIVE STATISTICS,
IN CONNECTION WITH THE POLYNESIAN DEFORESTATION DATA
==========================================================


(NOTES in connection with copy commands within SAS windows, using
select-mark-copy sequence of operations, all within EDIT menu).

[note could use word-processor search/replace commands to
 change specific missing indicators like X or NA to .]

libname mydir "."; 

data mydir.PolynesT;
  input ID Archip $  Islnd $ ForstRep  Deforstn  Area   NrArea  
  Isol_25 Isol_75  Elev   LatitS  Rain  Age PctMakat  Dust Tephra ;
datalines;
1  Bisma  NewBrL     2        2       35742  43784     30      95     2439    6.22   4156   2     0      550    3
2  Bisma  NewBrW     2        2       35742  43784     30      95     2439    6.22   6076   2     0      550    3
3  Bisma  NewIrld    2        2       7174   44482     30      30     2399    3      3302   3     0      550    3
4  Solom  Bougainv   2        2       8591   12471     50      200    2591    6.12   3009   2     0      550    3
5  Solom  Choiseul   2        2       2966   11698     45      45      970    7.05   3454   3     0      550    3
6  Solom  GuadalcL   2        4       5281    5856     50      50     2448    9.30   1500   3     0      550    3
7  Solom  GuadalcW   2        2       5281    5856     50      50     2448    9.30   4000   3     0      550    3
8  Solom  Malaita    2        2       4307    4745     45      45     1280    8.47   3096   3     0      550    3
9  Solom  SanCrist   2        2       3090    3168     45      45     1040    11     3332   3     0      550    3
10 Solom  SantaIsab  2        2       4095    4395     75      75     1250    8      3200   3     0      550    3
11 Solom  Tikopia    4        2          5       5    140     140      360    12     4000   2     0      370    3
12 NewCal NewCaled   2        3      16192   16192    400     1e5     1639    22     1105   3     0      370    3
13 Vanuat Aneityum   3      3.5        160     160     70      70      852    20.2   2290   2     0      370    3
14 Vanuat Efate      3        3        915     985     80     110      647    17.37  2293   2     0      370    3
15 Vanuat ErromanL   3        3        902    1482     45     100      886    18.3   1740   2     0      370    3
16 Vanuat ErromanW   3        3        902    1482     45     100      886    18.3   2822   2     0      370    3
17 Vanuat EspirSant  3        2       3900    6700     30     550     1879    15     3108   2     0      370    3
18 Vanuat Tanna      3        3        572    1482     45      45     1084    19     1643   2     0      370    3
19  Fiji  Lakeba     3        4         56      56     30      70      220    18.14  1930   3     0      260    3
20  Fiji  Rotuma     3.5      3         60      60    500     500      256    12     3533   3     0      260    3
21  Fiji  TaveuniL   3        2        431    5966     10      10     1241    16.50  2506   2     0      260    3
22  Fiji  TaveuniW   3        2        431    5966     10      10     1241    16.50  5784   2     0      260    3
23  Fiji  VanuLevL   3        4       5535    5966     65      65     1032    16.3   2445   3     0      260    3
24  Fiji  VanuLevW   3        3       5535    5966     65      65     1032    16.3   2404   3     0      260    3
25  Fiji  VitiLevL   3        4      10384   10384     65    1170     1323    17.5   1922   3     0      260    3
26  Fiji  VitiLevW   3        3      10384   10384     65    1170     1323    17.5   3158   3     0      260    3
27 Mcrones  Kosrae   4        2        109    109     550     550     629      5     5000   3     0      230    2
28 Mcrones  Yap      4        3        100    100     450     450     174     -9.3   3023   3     0      750    2
29 NewZea  NorthIsl  1        2     114453  266917     30      30     2797    38     1145   2     0      140    3
30 NewZea  SthIsl(E  1        4     150718  266917     30      30     3764    43.3    685   3     0      140    3
32  WPoly  Alofi     3        3       17.5    63.2    2.5     2.5      417    19      .     .     0      205    3
33  WPoly  Futuna    3        3       45.7    63.2    2.5     230      524    19     2700   3     0      205    3
34  WPoly  Niue      3        2       259      259    600     600       60    18     2009   .    50      150    2
35  Tonga  Tongatap  4        2       259      309     17     600       18    21.1   2000   .     0      205    3
36  Samoa  Savaii    3        2      1709     2823     20     400     1858    13     3000   2     0      150    2
37  Samoa  Upolu     3        2      1114     2823     20      20     1100    13     3000   2     0      150    2
38  Hawaii HawaiiL   3        4     10433    12317   3000    3175     4205    -19.3   500  1.5    0      420    1
39  Hawaii HawaiiW   3        3     10433    12317   3000    3175     4205    -19.3  3000  1.5    0      420    1
40  Hawaii KauaiL    3        4      1430     1610    116     116     1598    -22.1   600  2.5    0      420    1
41  Hawaii KauaiW    3        3      1430     1610    116     116     1598    -22.1  1500  2.5    0      420    1
42  Hawaii Lanai     3      3.5       364     3036     15      15     1026    -20.5   400   3     0      420    1
43  Hawaii MauiL     3        4      1884    13469     47      47     3055    -20.4   400   2     0      420    1
44  Hawaii MauiW     3        3      1884    13469     47      47     3055    -20.4  7000   2     0      420    1
45  Hawaii MolokaiL  3        4       673     4467     14      14     1514    -21.1   400   3     0      420    1
46  Hawaii MolokaiW  3        3       673     4467     14      14     1514    -21.1  1000  2.5    0      420    1
46.1 Hawaii  Necker  NA       5       0.2      0.2    290     290       85    -23.3   600   3     0      420    1
46.2 Hawaii  Nihoa   NA       5       0.7      0.7    209     209      278    -23     600   3     0      420    1
47  Hawaii  Niihau   5        5       180     1610     28      28      390    -21.5   600  2.5    0      420    1
48  Hawaii  OahuL    3        4      1546     2219     41      41     1220    -21.3   600  2.5    0      420    1
49  Hawaii  OahuW    3        3      1546     2219     41      41     1220    -21.3  1500  2.5    0      420    1
50  Cooks   Aitutaki 4        4        18       18    187     187      124     19    1894   3     0      210    1
51  Cooks   Atiu     4        4        28       48     22      22       71     20    1970   3    53      210    1
52  Cooks   Mangaia  4        4        70       70    165     165      169     21.5  1967   3    56      210    1
53  Cooks   Mauke    4        4        20       42     37      37       29     20.1  1626   3    75      210    1
54  Cooks   Mitiaro  4        4        22       70     22      22       15     19.5  1828   3    74      210    1
55  Cooks   Rarotong 4        3        67       67    165     165      653     21    2021   3     0      210    1
56  Austrls Raivavae 4        4        18       18    180     180      437     24    1841   3     0      168    1
57  Austrls Rapa     4        4        40       40    600     600      650     27.5  2738   3     0      156    1
58  Austrls Rimatara 4        4         9        9    150     150       83     23    1660   3    10      187    1
59  Austrls Rurutu   4        4        32       32    150     210      389     22    1899  2.5   10      181    1
60  Austrls Tubuai   4        4        45       45    180     210      422     23    2030   3     0      174    1
61  Societ  BoraBora 4        3        38      432     20      20      727     16.5  2248   3     0      182    1
62  Societ  Huahine  4        3        74      432     28      28      669     16.75 2987   3     0      182    1
63  Societ  Maupiti  4        3        12       50     47      47      372     16.5  1444   3     0      182    1
64  Societ  Moorea   4        3       132     1174     14      14     1212     17.5  2000   3     0      174    1
65  Societ  Raiatea  4        3       238      432      4     164     1033     16.5  2000   3     0      182    1
66  Societ  Tahaa    4        3        82      432      4       4      590     16.5  2055   3     0      182    1
67  Societ  Tahiti   4        3      1042     1174    164    2000     2241     17.5  1706  2.5    0      174    1
68  Tuamotu Makatea  3        2        31       31    240     240      113     15    1700   .    90      171    1
69  Marques Eiao     4        3        52      400      5      50      577      8    1000   3     0      138    1
70  Marques Fatuiva  4        3        80      476     35      35      960     10    1279   3     0      138    1
72  Marques HivaOa   4        3       320      476     40      65     1190     10    1156   3     0      138    1
73  Marques Nukuhiva 4        3       330      512     20      65     1208      9    1058   3     0      138    1
74  Marques Tahuata  4        3        61      476      4       4     1050     10    2403   3     0      138    1
75  Marques UaHuka   4      3.5        77      512     20      20      853      9     666   3     0      138    1
76  Marques UaPou    4      3.5       105      512     20      20     1252      9     673   3     0      138    1
77  SEPolyn Easter   5        5       166      166   1600    1600      510     27    1198  2.5    0        5    1
78  SEPolyn Pitcairn .        .         6        6    175     175      309     25    1663   2     0       99    1
79  SEPolyn MangarL  4        4        15       25      2     550      441     23    1349   3     0      118    1
80  SEPolyn MangarW  4        4        15       25      2     550      441     23    1932   3     0      118    1
         ;
  run;

NOTE: The data set MYDIR.POLYNEST has 80 observations and 16 variables.


proc plot data=PolynesT;
  plot rain * deforstn ;   * etc ;

proc gplot data=PolynesT;
  plot rain * deforstn ; 

[clearly shows decreasing tendency]

proc sort data=mydir.PolynesT out=PolynesS;
   by deforstn;
proc univariate data=polyness plot;
   var rain;
   by deforstn;
     run;



The UNIVARIATE Procedure
Variable:  Rain

Schematic Plots

         |
    8000 +
         |
         |                                    *
         |
    6000 +                        0
         |                        |
         |                        |
         |                        |
    4000 +                     +-----+
         |                     |  +  |        |
         |                     *-----*     +-----+
         |                        |        |     |        |           |
    2000 +                        |        *--+--*        |        *-----*
         |         *--+--*        |        +-----+     +-----+     |  +  |
         |                        |           |        |  +  |     |     |     +-----+
         |                                             *-----*     +-----+     *--+--*
       0 +
          ------------+-----------+-----------+-----------+-----------+-----------+-----------
Deforstn                 .           2           3         3.5           4           5



*  NEXT Histograms  ;

proc univariate data=polyness;
   var rain;
   histogram rain;
   run;

* This creates a high-quality graphical window with histogram.
  In future will try to control number of bars, locations, etc.
    using additional statement options.

* Next Tabulations ;

data TabPoly;
   set PolynesS;
   rainT = int(rain/1000);
proc freq data=TabPoly;
   table rainT;
      run;

The FREQ Procedure

                                  Cumulative    Cumulative
rainT    Frequency     Percent     Frequency      Percent
----------------------------------------------------------
    0          12       15.19            12        15.19
    1          29       36.71            41        51.90
    2          18       22.78            59        74.68
    3          13       16.46            72        91.14
    4           3        3.80            75        94.94
    5           2        2.53            77        97.47
    6           1        1.27            78        98.73
    7           1        1.27            79       100.00

Frequency Missing = 1
