R version 2.15.2 (2012-10-26) -- "Trick or Treat"
Copyright (C) 2012 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i686-pc-linux-gnu (32-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> x <- array(list(1418
+ ,210907
+ ,396
+ ,81
+ ,3
+ ,79
+ ,30
+ ,869
+ ,120982
+ ,297
+ ,55
+ ,4
+ ,58
+ ,28
+ ,1530
+ ,176508
+ ,559
+ ,50
+ ,12
+ ,60
+ ,38
+ ,2172
+ ,179321
+ ,967
+ ,125
+ ,2
+ ,108
+ ,30
+ ,901
+ ,123185
+ ,270
+ ,40
+ ,1
+ ,49
+ ,22
+ ,463
+ ,52746
+ ,143
+ ,37
+ ,3
+ ,0
+ ,26
+ ,3201
+ ,385534
+ ,1562
+ ,63
+ ,0
+ ,121
+ ,25
+ ,371
+ ,33170
+ ,109
+ ,44
+ ,0
+ ,1
+ ,18
+ ,1192
+ ,101645
+ ,371
+ ,88
+ ,0
+ ,20
+ ,11
+ ,1583
+ ,149061
+ ,656
+ ,66
+ ,5
+ ,43
+ ,26
+ ,1439
+ ,165446
+ ,511
+ ,57
+ ,0
+ ,69
+ ,25
+ ,1764
+ ,237213
+ ,655
+ ,74
+ ,0
+ ,78
+ ,38
+ ,1495
+ ,173326
+ ,465
+ ,49
+ ,7
+ ,86
+ ,44
+ ,1373
+ ,133131
+ ,525
+ ,52
+ ,7
+ ,44
+ ,30
+ ,2187
+ ,258873
+ ,885
+ ,88
+ ,3
+ ,104
+ ,40
+ ,1491
+ ,180083
+ ,497
+ ,36
+ ,9
+ ,63
+ ,34
+ ,4041
+ ,324799
+ ,1436
+ ,108
+ ,0
+ ,158
+ ,47
+ ,1706
+ ,230964
+ ,612
+ ,43
+ ,4
+ ,102
+ ,30
+ ,2152
+ ,236785
+ ,865
+ ,75
+ ,3
+ ,77
+ ,31
+ ,1036
+ ,135473
+ ,385
+ ,32
+ ,0
+ ,82
+ ,23
+ ,1882
+ ,202925
+ ,567
+ ,44
+ ,7
+ ,115
+ ,36
+ ,1929
+ ,215147
+ ,639
+ ,85
+ ,0
+ ,101
+ ,36
+ ,2242
+ ,344297
+ ,963
+ ,86
+ ,1
+ ,80
+ ,30
+ ,1220
+ ,153935
+ ,398
+ ,56
+ ,5
+ ,50
+ ,25
+ ,1289
+ ,132943
+ ,410
+ ,50
+ ,7
+ ,83
+ ,39
+ ,2515
+ ,174724
+ ,966
+ ,135
+ ,0
+ ,123
+ ,34
+ ,2147
+ ,174415
+ ,801
+ ,63
+ ,0
+ ,73
+ ,31
+ ,2352
+ ,225548
+ ,892
+ ,81
+ ,5
+ ,81
+ ,31
+ ,1638
+ ,223632
+ ,513
+ ,52
+ ,0
+ ,105
+ ,33
+ ,1222
+ ,124817
+ ,469
+ ,44
+ ,0
+ ,47
+ ,25
+ ,1812
+ ,221698
+ ,683
+ ,113
+ ,0
+ ,105
+ ,33
+ ,1677
+ ,210767
+ ,643
+ ,39
+ ,3
+ ,94
+ ,35
+ ,1579
+ ,170266
+ ,535
+ ,73
+ ,4
+ ,44
+ ,42
+ ,1731
+ ,260561
+ ,625
+ ,48
+ ,1
+ ,114
+ ,43
+ ,807
+ ,84853
+ ,264
+ ,33
+ ,4
+ ,38
+ ,30
+ ,2452
+ ,294424
+ ,992
+ ,59
+ ,2
+ ,107
+ ,33
+ ,829
+ ,101011
+ ,238
+ ,41
+ ,0
+ ,30
+ ,13
+ ,1940
+ ,215641
+ ,818
+ ,69
+ ,0
+ ,71
+ ,32
+ ,2662
+ ,325107
+ ,937
+ ,64
+ ,0
+ ,84
+ ,36
+ ,186
+ ,7176
+ ,70
+ ,1
+ ,0
+ ,0
+ ,0
+ ,1499
+ ,167542
+ ,507
+ ,59
+ ,2
+ ,59
+ ,28
+ ,865
+ ,106408
+ ,260
+ ,32
+ ,1
+ ,33
+ ,14
+ ,1793
+ ,96560
+ ,503
+ ,129
+ ,0
+ ,42
+ ,17
+ ,2527
+ ,265769
+ ,927
+ ,37
+ ,2
+ ,96
+ ,32
+ ,2747
+ ,269651
+ ,1269
+ ,31
+ ,10
+ ,106
+ ,30
+ ,1324
+ ,149112
+ ,537
+ ,65
+ ,6
+ ,56
+ ,35
+ ,2702
+ ,175824
+ ,910
+ ,107
+ ,0
+ ,57
+ ,20
+ ,1383
+ ,152871
+ ,532
+ ,74
+ ,5
+ ,59
+ ,28
+ ,1179
+ ,111665
+ ,345
+ ,54
+ ,4
+ ,39
+ ,28
+ ,2099
+ ,116408
+ ,918
+ ,76
+ ,1
+ ,34
+ ,39
+ ,4308
+ ,362301
+ ,1635
+ ,715
+ ,2
+ ,76
+ ,34
+ ,918
+ ,78800
+ ,330
+ ,57
+ ,2
+ ,20
+ ,26
+ ,1831
+ ,183167
+ ,557
+ ,66
+ ,0
+ ,91
+ ,39
+ ,3373
+ ,277965
+ ,1178
+ ,106
+ ,8
+ ,115
+ ,39
+ ,1713
+ ,150629
+ ,740
+ ,54
+ ,3
+ ,85
+ ,33
+ ,1438
+ ,168809
+ ,452
+ ,32
+ ,0
+ ,76
+ ,28
+ ,496
+ ,24188
+ ,218
+ ,20
+ ,0
+ ,8
+ ,4
+ ,2253
+ ,329267
+ ,764
+ ,71
+ ,8
+ ,79
+ ,39
+ ,744
+ ,65029
+ ,255
+ ,21
+ ,5
+ ,21
+ ,18
+ ,1161
+ ,101097
+ ,454
+ ,70
+ ,3
+ ,30
+ ,14
+ ,2352
+ ,218946
+ ,866
+ ,112
+ ,1
+ ,76
+ ,29
+ ,2144
+ ,244052
+ ,574
+ ,66
+ ,5
+ ,101
+ ,44
+ ,4691
+ ,341570
+ ,1276
+ ,190
+ ,1
+ ,94
+ ,21
+ ,1112
+ ,103597
+ ,379
+ ,66
+ ,1
+ ,27
+ ,16
+ ,2694
+ ,233328
+ ,825
+ ,165
+ ,5
+ ,92
+ ,28
+ ,1973
+ ,256462
+ ,798
+ ,56
+ ,0
+ ,123
+ ,35
+ ,1769
+ ,206161
+ ,663
+ ,61
+ ,12
+ ,75
+ ,28
+ ,3148
+ ,311473
+ ,1069
+ ,53
+ ,8
+ ,128
+ ,38
+ ,2474
+ ,235800
+ ,921
+ ,127
+ ,8
+ ,105
+ ,23
+ ,2084
+ ,177939
+ ,858
+ ,63
+ ,8
+ ,55
+ ,36
+ ,1954
+ ,207176
+ ,711
+ ,38
+ ,8
+ ,56
+ ,32
+ ,1226
+ ,196553
+ ,503
+ ,50
+ ,2
+ ,41
+ ,29
+ ,1389
+ ,174184
+ ,382
+ ,52
+ ,0
+ ,72
+ ,25
+ ,1496
+ ,143246
+ ,464
+ ,42
+ ,5
+ ,67
+ ,27
+ ,2269
+ ,187559
+ ,717
+ ,76
+ ,8
+ ,75
+ ,36
+ ,1833
+ ,187681
+ ,690
+ ,67
+ ,2
+ ,114
+ ,28
+ ,1268
+ ,119016
+ ,462
+ ,50
+ ,5
+ ,118
+ ,23
+ ,1943
+ ,182192
+ ,657
+ ,53
+ ,12
+ ,77
+ ,40
+ ,893
+ ,73566
+ ,385
+ ,39
+ ,6
+ ,22
+ ,23
+ ,1762
+ ,194979
+ ,577
+ ,50
+ ,7
+ ,66
+ ,40
+ ,1403
+ ,167488
+ ,619
+ ,77
+ ,2
+ ,69
+ ,28
+ ,1425
+ ,143756
+ ,479
+ ,57
+ ,0
+ ,105
+ ,34
+ ,1857
+ ,275541
+ ,817
+ ,73
+ ,4
+ ,116
+ ,33
+ ,1840
+ ,243199
+ ,752
+ ,34
+ ,3
+ ,88
+ ,28
+ ,1502
+ ,182999
+ ,430
+ ,39
+ ,6
+ ,73
+ ,34
+ ,1441
+ ,135649
+ ,451
+ ,46
+ ,2
+ ,99
+ ,30
+ ,1420
+ ,152299
+ ,537
+ ,63
+ ,0
+ ,62
+ ,33
+ ,1416
+ ,120221
+ ,519
+ ,35
+ ,1
+ ,53
+ ,22
+ ,2970
+ ,346485
+ ,1000
+ ,106
+ ,0
+ ,118
+ ,38
+ ,1317
+ ,145790
+ ,637
+ ,43
+ ,5
+ ,30
+ ,26
+ ,1644
+ ,193339
+ ,465
+ ,47
+ ,2
+ ,100
+ ,35
+ ,870
+ ,80953
+ ,437
+ ,31
+ ,0
+ ,49
+ ,8
+ ,1654
+ ,122774
+ ,711
+ ,162
+ ,0
+ ,24
+ ,24
+ ,1054
+ ,130585
+ ,299
+ ,57
+ ,5
+ ,67
+ ,29
+ ,937
+ ,112611
+ ,248
+ ,36
+ ,0
+ ,46
+ ,20
+ ,3004
+ ,286468
+ ,1162
+ ,263
+ ,1
+ ,57
+ ,29
+ ,2008
+ ,241066
+ ,714
+ ,78
+ ,0
+ ,75
+ ,45
+ ,2547
+ ,148446
+ ,905
+ ,63
+ ,1
+ ,135
+ ,37
+ ,1885
+ ,204713
+ ,649
+ ,54
+ ,1
+ ,68
+ ,33
+ ,1626
+ ,182079
+ ,512
+ ,63
+ ,2
+ ,124
+ ,33
+ ,1468
+ ,140344
+ ,472
+ ,77
+ ,6
+ ,33
+ ,25
+ ,2445
+ ,220516
+ ,905
+ ,79
+ ,1
+ ,98
+ ,32
+ ,1964
+ ,243060
+ ,786
+ ,110
+ ,4
+ ,58
+ ,29
+ ,1381
+ ,162765
+ ,489
+ ,56
+ ,2
+ ,68
+ ,28
+ ,1369
+ ,182613
+ ,479
+ ,56
+ ,3
+ ,81
+ ,28
+ ,1659
+ ,232138
+ ,617
+ ,43
+ ,0
+ ,131
+ ,31
+ ,2888
+ ,265318
+ ,925
+ ,111
+ ,10
+ ,110
+ ,52
+ ,1290
+ ,85574
+ ,351
+ ,71
+ ,0
+ ,37
+ ,21
+ ,2845
+ ,310839
+ ,1144
+ ,62
+ ,9
+ ,130
+ ,24
+ ,1982
+ ,225060
+ ,669
+ ,56
+ ,7
+ ,93
+ ,41
+ ,1904
+ ,232317
+ ,707
+ ,74
+ ,0
+ ,118
+ ,33
+ ,1391
+ ,144966
+ ,458
+ ,60
+ ,0
+ ,39
+ ,32
+ ,602
+ ,43287
+ ,214
+ ,43
+ ,4
+ ,13
+ ,19
+ ,1743
+ ,155754
+ ,599
+ ,68
+ ,4
+ ,74
+ ,20
+ ,1559
+ ,164709
+ ,572
+ ,53
+ ,0
+ ,81
+ ,31
+ ,2014
+ ,201940
+ ,897
+ ,87
+ ,0
+ ,109
+ ,31
+ ,2143
+ ,235454
+ ,819
+ ,46
+ ,0
+ ,151
+ ,32
+ ,2146
+ ,220801
+ ,720
+ ,105
+ ,1
+ ,51
+ ,18
+ ,874
+ ,99466
+ ,273
+ ,32
+ ,0
+ ,28
+ ,23
+ ,1590
+ ,92661
+ ,508
+ ,133
+ ,1
+ ,40
+ ,17
+ ,1590
+ ,133328
+ ,506
+ ,79
+ ,0
+ ,56
+ ,20
+ ,1210
+ ,61361
+ ,451
+ ,51
+ ,0
+ ,27
+ ,12
+ ,2072
+ ,125930
+ ,699
+ ,207
+ ,4
+ ,37
+ ,17
+ ,1281
+ ,100750
+ ,407
+ ,67
+ ,0
+ ,83
+ ,30
+ ,1401
+ ,224549
+ ,465
+ ,47
+ ,4
+ ,54
+ ,31
+ ,834
+ ,82316
+ ,245
+ ,34
+ ,4
+ ,27
+ ,10
+ ,1105
+ ,102010
+ ,370
+ ,66
+ ,3
+ ,28
+ ,13
+ ,1272
+ ,101523
+ ,316
+ ,76
+ ,0
+ ,59
+ ,22
+ ,1944
+ ,243511
+ ,603
+ ,65
+ ,0
+ ,133
+ ,42
+ ,391
+ ,22938
+ ,154
+ ,9
+ ,0
+ ,12
+ ,1
+ ,761
+ ,41566
+ ,229
+ ,42
+ ,5
+ ,0
+ ,9
+ ,1605
+ ,152474
+ ,577
+ ,45
+ ,0
+ ,106
+ ,32
+ ,530
+ ,61857
+ ,192
+ ,25
+ ,4
+ ,23
+ ,11
+ ,1988
+ ,99923
+ ,617
+ ,115
+ ,0
+ ,44
+ ,25
+ ,1386
+ ,132487
+ ,411
+ ,97
+ ,0
+ ,71
+ ,36
+ ,2395
+ ,317394
+ ,975
+ ,53
+ ,1
+ ,116
+ ,31
+ ,387
+ ,21054
+ ,146
+ ,2
+ ,0
+ ,4
+ ,0
+ ,1742
+ ,209641
+ ,705
+ ,52
+ ,5
+ ,62
+ ,24
+ ,620
+ ,22648
+ ,184
+ ,44
+ ,0
+ ,12
+ ,13
+ ,449
+ ,31414
+ ,200
+ ,22
+ ,0
+ ,18
+ ,8
+ ,800
+ ,46698
+ ,274
+ ,35
+ ,0
+ ,14
+ ,13
+ ,1684
+ ,131698
+ ,502
+ ,74
+ ,0
+ ,60
+ ,19
+ ,1050
+ ,91735
+ ,382
+ ,103
+ ,0
+ ,7
+ ,18
+ ,2699
+ ,244749
+ ,964
+ ,144
+ ,2
+ ,98
+ ,33
+ ,1606
+ ,184510
+ ,537
+ ,60
+ ,7
+ ,64
+ ,40
+ ,1502
+ ,79863
+ ,438
+ ,134
+ ,1
+ ,29
+ ,22
+ ,1204
+ ,128423
+ ,369
+ ,89
+ ,8
+ ,32
+ ,38
+ ,1138
+ ,97839
+ ,417
+ ,42
+ ,2
+ ,25
+ ,24
+ ,568
+ ,38214
+ ,276
+ ,52
+ ,0
+ ,16
+ ,8
+ ,1459
+ ,151101
+ ,514
+ ,98
+ ,2
+ ,48
+ ,35
+ ,2158
+ ,272458
+ ,822
+ ,99
+ ,0
+ ,100
+ ,43
+ ,1111
+ ,172494
+ ,389
+ ,52
+ ,0
+ ,46
+ ,43
+ ,1421
+ ,108043
+ ,466
+ ,29
+ ,1
+ ,45
+ ,14
+ ,2833
+ ,328107
+ ,1255
+ ,125
+ ,3
+ ,129
+ ,41
+ ,1955
+ ,250579
+ ,694
+ ,106
+ ,0
+ ,130
+ ,38
+ ,2922
+ ,351067
+ ,1024
+ ,95
+ ,3
+ ,136
+ ,45
+ ,1002
+ ,158015
+ ,400
+ ,40
+ ,0
+ ,59
+ ,31
+ ,1060
+ ,98866
+ ,397
+ ,140
+ ,0
+ ,25
+ ,13
+ ,956
+ ,85439
+ ,350
+ ,43
+ ,0
+ ,32
+ ,28
+ ,2186
+ ,229242
+ ,719
+ ,128
+ ,4
+ ,63
+ ,31
+ ,3604
+ ,351619
+ ,1277
+ ,142
+ ,4
+ ,95
+ ,40
+ ,1035
+ ,84207
+ ,356
+ ,73
+ ,11
+ ,14
+ ,30
+ ,1417
+ ,120445
+ ,457
+ ,72
+ ,0
+ ,36
+ ,16
+ ,3261
+ ,324598
+ ,1402
+ ,128
+ ,0
+ ,113
+ ,37
+ ,1587
+ ,131069
+ ,600
+ ,61
+ ,4
+ ,47
+ ,30
+ ,1424
+ ,204271
+ ,480
+ ,73
+ ,0
+ ,92
+ ,35
+ ,1701
+ ,165543
+ ,595
+ ,148
+ ,1
+ ,70
+ ,32
+ ,1249
+ ,141722
+ ,436
+ ,64
+ ,0
+ ,19
+ ,27
+ ,946
+ ,116048
+ ,230
+ ,45
+ ,0
+ ,50
+ ,20
+ ,1926
+ ,250047
+ ,651
+ ,58
+ ,0
+ ,41
+ ,18
+ ,3352
+ ,299775
+ ,1367
+ ,97
+ ,9
+ ,91
+ ,31
+ ,1641
+ ,195838
+ ,564
+ ,50
+ ,1
+ ,111
+ ,31
+ ,2035
+ ,173260
+ ,716
+ ,37
+ ,3
+ ,41
+ ,21
+ ,2312
+ ,254488
+ ,747
+ ,50
+ ,10
+ ,120
+ ,39
+ ,1369
+ ,104389
+ ,467
+ ,105
+ ,5
+ ,135
+ ,41
+ ,1577
+ ,136084
+ ,671
+ ,69
+ ,0
+ ,27
+ ,13
+ ,2201
+ ,199476
+ ,861
+ ,46
+ ,2
+ ,87
+ ,32
+ ,961
+ ,92499
+ ,319
+ ,57
+ ,0
+ ,25
+ ,18
+ ,1900
+ ,224330
+ ,612
+ ,52
+ ,1
+ ,131
+ ,39
+ ,1254
+ ,135781
+ ,433
+ ,98
+ ,2
+ ,45
+ ,14
+ ,1335
+ ,74408
+ ,434
+ ,61
+ ,4
+ ,29
+ ,7
+ ,1597
+ ,81240
+ ,503
+ ,89
+ ,0
+ ,58
+ ,17
+ ,207
+ ,14688
+ ,85
+ ,0
+ ,0
+ ,4
+ ,0
+ ,1645
+ ,181633
+ ,564
+ ,48
+ ,2
+ ,47
+ ,30
+ ,2429
+ ,271856
+ ,824
+ ,91
+ ,1
+ ,109
+ ,37
+ ,151
+ ,7199
+ ,74
+ ,0
+ ,0
+ ,7
+ ,0
+ ,474
+ ,46660
+ ,259
+ ,7
+ ,0
+ ,12
+ ,5
+ ,141
+ ,17547
+ ,69
+ ,3
+ ,0
+ ,0
+ ,1
+ ,1639
+ ,133368
+ ,535
+ ,54
+ ,1
+ ,37
+ ,16
+ ,872
+ ,95227
+ ,239
+ ,70
+ ,0
+ ,37
+ ,32
+ ,1318
+ ,152601
+ ,438
+ ,36
+ ,2
+ ,46
+ ,24
+ ,1018
+ ,98146
+ ,459
+ ,37
+ ,0
+ ,15
+ ,17
+ ,1383
+ ,79619
+ ,426
+ ,123
+ ,3
+ ,42
+ ,11
+ ,1314
+ ,59194
+ ,288
+ ,247
+ ,6
+ ,7
+ ,24
+ ,1335
+ ,139942
+ ,498
+ ,46
+ ,0
+ ,54
+ ,22
+ ,1403
+ ,118612
+ ,454
+ ,72
+ ,2
+ ,54
+ ,12
+ ,910
+ ,72880
+ ,376
+ ,41
+ ,0
+ ,14
+ ,19
+ ,616
+ ,65475
+ ,225
+ ,24
+ ,2
+ ,16
+ ,13
+ ,1407
+ ,99643
+ ,555
+ ,45
+ ,1
+ ,33
+ ,17
+ ,771
+ ,71965
+ ,252
+ ,33
+ ,1
+ ,32
+ ,15
+ ,766
+ ,77272
+ ,208
+ ,27
+ ,2
+ ,21
+ ,16
+ ,473
+ ,49289
+ ,130
+ ,36
+ ,1
+ ,15
+ ,24
+ ,1376
+ ,135131
+ ,481
+ ,87
+ ,0
+ ,38
+ ,15
+ ,1232
+ ,108446
+ ,389
+ ,90
+ ,1
+ ,22
+ ,17
+ ,1521
+ ,89746
+ ,565
+ ,114
+ ,3
+ ,28
+ ,18
+ ,572
+ ,44296
+ ,173
+ ,31
+ ,0
+ ,10
+ ,20
+ ,1059
+ ,77648
+ ,278
+ ,45
+ ,0
+ ,31
+ ,16
+ ,1544
+ ,181528
+ ,609
+ ,69
+ ,0
+ ,32
+ ,16
+ ,1230
+ ,134019
+ ,422
+ ,51
+ ,0
+ ,32
+ ,18
+ ,1206
+ ,124064
+ ,445
+ ,34
+ ,1
+ ,43
+ ,22
+ ,1205
+ ,92630
+ ,387
+ ,60
+ ,4
+ ,27
+ ,8
+ ,1255
+ ,121848
+ ,339
+ ,45
+ ,0
+ ,37
+ ,17
+ ,613
+ ,52915
+ ,181
+ ,54
+ ,0
+ ,20
+ ,18
+ ,721
+ ,81872
+ ,245
+ ,25
+ ,0
+ ,32
+ ,16
+ ,1109
+ ,58981
+ ,384
+ ,38
+ ,7
+ ,0
+ ,23
+ ,740
+ ,53515
+ ,212
+ ,52
+ ,2
+ ,5
+ ,22
+ ,1126
+ ,60812
+ ,399
+ ,67
+ ,0
+ ,26
+ ,13
+ ,728
+ ,56375
+ ,229
+ ,74
+ ,7
+ ,10
+ ,13
+ ,689
+ ,65490
+ ,224
+ ,38
+ ,3
+ ,27
+ ,16
+ ,592
+ ,80949
+ ,203
+ ,30
+ ,0
+ ,11
+ ,16
+ ,995
+ ,76302
+ ,333
+ ,26
+ ,0
+ ,29
+ ,20
+ ,1613
+ ,104011
+ ,384
+ ,67
+ ,6
+ ,25
+ ,22
+ ,2048
+ ,98104
+ ,636
+ ,132
+ ,2
+ ,55
+ ,17
+ ,705
+ ,67989
+ ,185
+ ,42
+ ,0
+ ,23
+ ,18
+ ,301
+ ,30989
+ ,93
+ ,35
+ ,0
+ ,5
+ ,17
+ ,1803
+ ,135458
+ ,581
+ ,118
+ ,3
+ ,43
+ ,12
+ ,799
+ ,73504
+ ,248
+ ,68
+ ,0
+ ,23
+ ,7
+ ,861
+ ,63123
+ ,304
+ ,43
+ ,1
+ ,34
+ ,17
+ ,1186
+ ,61254
+ ,344
+ ,76
+ ,1
+ ,36
+ ,14
+ ,1451
+ ,74914
+ ,407
+ ,64
+ ,0
+ ,35
+ ,23
+ ,628
+ ,31774
+ ,170
+ ,48
+ ,1
+ ,0
+ ,17
+ ,1161
+ ,81437
+ ,312
+ ,64
+ ,0
+ ,37
+ ,14
+ ,1463
+ ,87186
+ ,507
+ ,56
+ ,0
+ ,28
+ ,15
+ ,742
+ ,50090
+ ,224
+ ,71
+ ,0
+ ,16
+ ,17
+ ,979
+ ,65745
+ ,340
+ ,75
+ ,0
+ ,26
+ ,21
+ ,675
+ ,56653
+ ,168
+ ,39
+ ,0
+ ,38
+ ,18
+ ,1241
+ ,158399
+ ,443
+ ,42
+ ,0
+ ,23
+ ,18
+ ,676
+ ,46455
+ ,204
+ ,39
+ ,0
+ ,22
+ ,17
+ ,1049
+ ,73624
+ ,367
+ ,93
+ ,0
+ ,30
+ ,17
+ ,620
+ ,38395
+ ,210
+ ,38
+ ,0
+ ,16
+ ,16
+ ,1081
+ ,91899
+ ,335
+ ,60
+ ,0
+ ,18
+ ,15
+ ,1688
+ ,139526
+ ,364
+ ,71
+ ,0
+ ,28
+ ,21
+ ,736
+ ,52164
+ ,178
+ ,52
+ ,0
+ ,32
+ ,16
+ ,617
+ ,51567
+ ,206
+ ,27
+ ,2
+ ,21
+ ,14
+ ,812
+ ,70551
+ ,279
+ ,59
+ ,0
+ ,23
+ ,15
+ ,1051
+ ,84856
+ ,387
+ ,40
+ ,1
+ ,29
+ ,17
+ ,1656
+ ,102538
+ ,490
+ ,79
+ ,1
+ ,50
+ ,15
+ ,705
+ ,86678
+ ,238
+ ,44
+ ,0
+ ,12
+ ,15
+ ,945
+ ,85709
+ ,343
+ ,65
+ ,0
+ ,21
+ ,10
+ ,554
+ ,34662
+ ,232
+ ,10
+ ,0
+ ,18
+ ,6
+ ,1597
+ ,150580
+ ,530
+ ,124
+ ,0
+ ,27
+ ,22
+ ,982
+ ,99611
+ ,291
+ ,81
+ ,0
+ ,41
+ ,21
+ ,222
+ ,19349
+ ,67
+ ,15
+ ,0
+ ,13
+ ,1
+ ,1212
+ ,99373
+ ,397
+ ,92
+ ,1
+ ,12
+ ,18
+ ,1143
+ ,86230
+ ,467
+ ,42
+ ,0
+ ,21
+ ,17
+ ,435
+ ,30837
+ ,178
+ ,10
+ ,0
+ ,8
+ ,4
+ ,532
+ ,31706
+ ,175
+ ,24
+ ,0
+ ,26
+ ,10
+ ,882
+ ,89806
+ ,299
+ ,64
+ ,0
+ ,27
+ ,16
+ ,608
+ ,62088
+ ,154
+ ,45
+ ,1
+ ,13
+ ,16
+ ,459
+ ,40151
+ ,106
+ ,22
+ ,0
+ ,16
+ ,9
+ ,578
+ ,27634
+ ,189
+ ,56
+ ,0
+ ,2
+ ,16
+ ,826
+ ,76990
+ ,194
+ ,94
+ ,0
+ ,42
+ ,17
+ ,509
+ ,37460
+ ,135
+ ,19
+ ,0
+ ,5
+ ,7
+ ,717
+ ,54157
+ ,201
+ ,35
+ ,0
+ ,37
+ ,15
+ ,637
+ ,49862
+ ,207
+ ,32
+ ,0
+ ,17
+ ,14
+ ,857
+ ,84337
+ ,280
+ ,35
+ ,0
+ ,38
+ ,14
+ ,830
+ ,64175
+ ,260
+ ,48
+ ,0
+ ,37
+ ,18
+ ,652
+ ,59382
+ ,227
+ ,49
+ ,0
+ ,29
+ ,12
+ ,707
+ ,119308
+ ,239
+ ,48
+ ,0
+ ,32
+ ,16
+ ,954
+ ,76702
+ ,333
+ ,62
+ ,0
+ ,35
+ ,21
+ ,1461
+ ,103425
+ ,428
+ ,96
+ ,1
+ ,17
+ ,19
+ ,672
+ ,70344
+ ,230
+ ,45
+ ,0
+ ,20
+ ,16
+ ,778
+ ,43410
+ ,292
+ ,63
+ ,0
+ ,7
+ ,1
+ ,1141
+ ,104838
+ ,350
+ ,71
+ ,1
+ ,46
+ ,16
+ ,680
+ ,62215
+ ,186
+ ,26
+ ,0
+ ,24
+ ,10
+ ,1090
+ ,69304
+ ,326
+ ,48
+ ,6
+ ,40
+ ,19
+ ,616
+ ,53117
+ ,155
+ ,29
+ ,3
+ ,3
+ ,12
+ ,285
+ ,19764
+ ,75
+ ,19
+ ,1
+ ,10
+ ,2
+ ,1145
+ ,86680
+ ,361
+ ,45
+ ,2
+ ,37
+ ,14
+ ,733
+ ,84105
+ ,261
+ ,45
+ ,0
+ ,17
+ ,17
+ ,888
+ ,77945
+ ,299
+ ,67
+ ,0
+ ,28
+ ,19
+ ,849
+ ,89113
+ ,300
+ ,30
+ ,0
+ ,19
+ ,14
+ ,1182
+ ,91005
+ ,450
+ ,36
+ ,3
+ ,29
+ ,11
+ ,528
+ ,40248
+ ,183
+ ,34
+ ,1
+ ,8
+ ,4
+ ,642
+ ,64187
+ ,238
+ ,36
+ ,0
+ ,10
+ ,16
+ ,947
+ ,50857
+ ,165
+ ,34
+ ,0
+ ,15
+ ,20
+ ,819
+ ,56613
+ ,234
+ ,37
+ ,1
+ ,15
+ ,12
+ ,757
+ ,62792
+ ,176
+ ,46
+ ,0
+ ,28
+ ,15
+ ,894
+ ,72535
+ ,329
+ ,44
+ ,0
+ ,17
+ ,16)
+ ,dim=c(7
+ ,289)
+ ,dimnames=list(c('pageviews'
+ ,'time'
+ ,'compinfo'
+ ,'comppr'
+ ,'sharedcomp'
+ ,'bloggedcomp'
+ ,'compreviewed')
+ ,1:289))
>  y <- array(NA,dim=c(7,289),dimnames=list(c('pageviews','time','compinfo','comppr','sharedcomp','bloggedcomp','compreviewed'),1:289))
>  for (i in 1:dim(x)[1])
+  {
+  	for (j in 1:dim(x)[2])
+  	{
+  		y[i,j] <- as.numeric(x[i,j])
+  	}
+  }
> par4 = 'no'
> par3 = '3'
> par2 = 'none'
> par1 = '2'
> library(party)
Loading required package: survival
Loading required package: splines
Loading required package: grid
Loading required package: modeltools
Loading required package: stats4
Loading required package: coin
Loading required package: mvtnorm
Loading required package: zoo

Attaching package: 'zoo'

The following object(s) are masked from 'package:base':

    as.Date, as.Date.numeric

Loading required package: sandwich
Loading required package: strucchange
Loading required package: vcd
Loading required package: MASS
Loading required package: colorspace
> library(Hmisc)
Hmisc library by Frank E Harrell Jr

Type library(help='Hmisc'), ?Overview, or ?Hmisc.Overview')
to see overall documentation.

NOTE:Hmisc no longer redefines [.factor to drop unused levels when
subsetting.  To get the old behavior of Hmisc type dropUnusedLevels().


Attaching package: 'Hmisc'

The following object(s) are masked from 'package:survival':

    untangle.specials

The following object(s) are masked from 'package:base':

    format.pval, round.POSIXt, trunc.POSIXt, units

> par1 <- as.numeric(par1)
> par3 <- as.numeric(par3)
> x <- data.frame(t(y))
> is.data.frame(x)
[1] TRUE
> x <- x[!is.na(x[,par1]),]
> k <- length(x[1,])
> n <- length(x[,1])
> colnames(x)[par1]
[1] "time"
> x[,par1]
  [1] 210907 120982 176508 179321 123185  52746 385534  33170 101645 149061
 [11] 165446 237213 173326 133131 258873 180083 324799 230964 236785 135473
 [21] 202925 215147 344297 153935 132943 174724 174415 225548 223632 124817
 [31] 221698 210767 170266 260561  84853 294424 101011 215641 325107   7176
 [41] 167542 106408  96560 265769 269651 149112 175824 152871 111665 116408
 [51] 362301  78800 183167 277965 150629 168809  24188 329267  65029 101097
 [61] 218946 244052 341570 103597 233328 256462 206161 311473 235800 177939
 [71] 207176 196553 174184 143246 187559 187681 119016 182192  73566 194979
 [81] 167488 143756 275541 243199 182999 135649 152299 120221 346485 145790
 [91] 193339  80953 122774 130585 112611 286468 241066 148446 204713 182079
[101] 140344 220516 243060 162765 182613 232138 265318  85574 310839 225060
[111] 232317 144966  43287 155754 164709 201940 235454 220801  99466  92661
[121] 133328  61361 125930 100750 224549  82316 102010 101523 243511  22938
[131]  41566 152474  61857  99923 132487 317394  21054 209641  22648  31414
[141]  46698 131698  91735 244749 184510  79863 128423  97839  38214 151101
[151] 272458 172494 108043 328107 250579 351067 158015  98866  85439 229242
[161] 351619  84207 120445 324598 131069 204271 165543 141722 116048 250047
[171] 299775 195838 173260 254488 104389 136084 199476  92499 224330 135781
[181]  74408  81240  14688 181633 271856   7199  46660  17547 133368  95227
[191] 152601  98146  79619  59194 139942 118612  72880  65475  99643  71965
[201]  77272  49289 135131 108446  89746  44296  77648 181528 134019 124064
[211]  92630 121848  52915  81872  58981  53515  60812  56375  65490  80949
[221]  76302 104011  98104  67989  30989 135458  73504  63123  61254  74914
[231]  31774  81437  87186  50090  65745  56653 158399  46455  73624  38395
[241]  91899 139526  52164  51567  70551  84856 102538  86678  85709  34662
[251] 150580  99611  19349  99373  86230  30837  31706  89806  62088  40151
[261]  27634  76990  37460  54157  49862  84337  64175  59382 119308  76702
[271] 103425  70344  43410 104838  62215  69304  53117  19764  86680  84105
[281]  77945  89113  91005  40248  64187  50857  56613  62792  72535
> if (par2 == 'kmeans') {
+ cl <- kmeans(x[,par1], par3)
+ print(cl)
+ clm <- matrix(cbind(cl$centers,1:par3),ncol=2)
+ clm <- clm[sort.list(clm[,1]),]
+ for (i in 1:par3) {
+ cl$cluster[cl$cluster==clm[i,2]] <- paste('C',i,sep='')
+ }
+ cl$cluster <- as.factor(cl$cluster)
+ print(cl$cluster)
+ x[,par1] <- cl$cluster
+ }
> if (par2 == 'quantiles') {
+ x[,par1] <- cut2(x[,par1],g=par3)
+ }
> if (par2 == 'hclust') {
+ hc <- hclust(dist(x[,par1])^2, 'cen')
+ print(hc)
+ memb <- cutree(hc, k = par3)
+ dum <- c(mean(x[memb==1,par1]))
+ for (i in 2:par3) {
+ dum <- c(dum, mean(x[memb==i,par1]))
+ }
+ hcm <- matrix(cbind(dum,1:par3),ncol=2)
+ hcm <- hcm[sort.list(hcm[,1]),]
+ for (i in 1:par3) {
+ memb[memb==hcm[i,2]] <- paste('C',i,sep='')
+ }
+ memb <- as.factor(memb)
+ print(memb)
+ x[,par1] <- memb
+ }
> if (par2=='equal') {
+ ed <- cut(as.numeric(x[,par1]),par3,labels=paste('C',1:par3,sep=''))
+ x[,par1] <- as.factor(ed)
+ }
> table(x[,par1])

  7176   7199  14688  17547  19349  19764  21054  22648  22938  24188  27634 
     1      1      1      1      1      1      1      1      1      1      1 
 30837  30989  31414  31706  31774  33170  34662  37460  38214  38395  40151 
     1      1      1      1      1      1      1      1      1      1      1 
 40248  41566  43287  43410  44296  46455  46660  46698  49289  49862  50090 
     1      1      1      1      1      1      1      1      1      1      1 
 50857  51567  52164  52746  52915  53117  53515  54157  56375  56613  56653 
     1      1      1      1      1      1      1      1      1      1      1 
 58981  59194  59382  60812  61254  61361  61857  62088  62215  62792  63123 
     1      1      1      1      1      1      1      1      1      1      1 
 64175  64187  65029  65475  65490  65745  67989  69304  70344  70551  71965 
     1      1      1      1      1      1      1      1      1      1      1 
 72535  72880  73504  73566  73624  74408  74914  76302  76702  76990  77272 
     1      1      1      1      1      1      1      1      1      1      1 
 77648  77945  78800  79619  79863  80949  80953  81240  81437  81872  82316 
     1      1      1      1      1      1      1      1      1      1      1 
 84105  84207  84337  84853  84856  85439  85574  85709  86230  86678  86680 
     1      1      1      1      1      1      1      1      1      1      1 
 87186  89113  89746  89806  91005  91735  91899  92499  92630  92661  95227 
     1      1      1      1      1      1      1      1      1      1      1 
 96560  97839  98104  98146  98866  99373  99466  99611  99643  99923 100750 
     1      1      1      1      1      1      1      1      1      1      1 
101011 101097 101523 101645 102010 102538 103425 103597 104011 104389 104838 
     1      1      1      1      1      1      1      1      1      1      1 
106408 108043 108446 111665 112611 116048 116408 118612 119016 119308 120221 
     1      1      1      1      1      1      1      1      1      1      1 
120445 120982 121848 122774 123185 124064 124817 125930 128423 130585 131069 
     1      1      1      1      1      1      1      1      1      1      1 
131698 132487 132943 133131 133328 133368 134019 135131 135458 135473 135649 
     1      1      1      1      1      1      1      1      1      1      1 
135781 136084 139526 139942 140344 141722 143246 143756 144966 145790 148446 
     1      1      1      1      1      1      1      1      1      1      1 
149061 149112 150580 150629 151101 152299 152474 152601 152871 153935 155754 
     1      1      1      1      1      1      1      1      1      1      1 
158015 158399 162765 164709 165446 165543 167488 167542 168809 170266 172494 
     1      1      1      1      1      1      1      1      1      1      1 
173260 173326 174184 174415 174724 175824 176508 177939 179321 180083 181528 
     1      1      1      1      1      1      1      1      1      1      1 
181633 182079 182192 182613 182999 183167 184510 187559 187681 193339 194979 
     1      1      1      1      1      1      1      1      1      1      1 
195838 196553 199476 201940 202925 204271 204713 206161 207176 209641 210767 
     1      1      1      1      1      1      1      1      1      1      1 
210907 215147 215641 218946 220516 220801 221698 223632 224330 224549 225060 
     1      1      1      1      1      1      1      1      1      1      1 
225548 229242 230964 232138 232317 233328 235454 235800 236785 237213 241066 
     1      1      1      1      1      1      1      1      1      1      1 
243060 243199 243511 244052 244749 250047 250579 254488 256462 258873 260561 
     1      1      1      1      1      1      1      1      1      1      1 
265318 265769 269651 271856 272458 275541 277965 286468 294424 299775 310839 
     1      1      1      1      1      1      1      1      1      1      1 
311473 317394 324598 324799 325107 328107 329267 341570 344297 346485 351067 
     1      1      1      1      1      1      1      1      1      1      1 
351619 362301 385534 
     1      1      1 
> colnames(x)
[1] "pageviews"    "time"         "compinfo"     "comppr"       "sharedcomp"  
[6] "bloggedcomp"  "compreviewed"
> colnames(x)[par1]
[1] "time"
> x[,par1]
  [1] 210907 120982 176508 179321 123185  52746 385534  33170 101645 149061
 [11] 165446 237213 173326 133131 258873 180083 324799 230964 236785 135473
 [21] 202925 215147 344297 153935 132943 174724 174415 225548 223632 124817
 [31] 221698 210767 170266 260561  84853 294424 101011 215641 325107   7176
 [41] 167542 106408  96560 265769 269651 149112 175824 152871 111665 116408
 [51] 362301  78800 183167 277965 150629 168809  24188 329267  65029 101097
 [61] 218946 244052 341570 103597 233328 256462 206161 311473 235800 177939
 [71] 207176 196553 174184 143246 187559 187681 119016 182192  73566 194979
 [81] 167488 143756 275541 243199 182999 135649 152299 120221 346485 145790
 [91] 193339  80953 122774 130585 112611 286468 241066 148446 204713 182079
[101] 140344 220516 243060 162765 182613 232138 265318  85574 310839 225060
[111] 232317 144966  43287 155754 164709 201940 235454 220801  99466  92661
[121] 133328  61361 125930 100750 224549  82316 102010 101523 243511  22938
[131]  41566 152474  61857  99923 132487 317394  21054 209641  22648  31414
[141]  46698 131698  91735 244749 184510  79863 128423  97839  38214 151101
[151] 272458 172494 108043 328107 250579 351067 158015  98866  85439 229242
[161] 351619  84207 120445 324598 131069 204271 165543 141722 116048 250047
[171] 299775 195838 173260 254488 104389 136084 199476  92499 224330 135781
[181]  74408  81240  14688 181633 271856   7199  46660  17547 133368  95227
[191] 152601  98146  79619  59194 139942 118612  72880  65475  99643  71965
[201]  77272  49289 135131 108446  89746  44296  77648 181528 134019 124064
[211]  92630 121848  52915  81872  58981  53515  60812  56375  65490  80949
[221]  76302 104011  98104  67989  30989 135458  73504  63123  61254  74914
[231]  31774  81437  87186  50090  65745  56653 158399  46455  73624  38395
[241]  91899 139526  52164  51567  70551  84856 102538  86678  85709  34662
[251] 150580  99611  19349  99373  86230  30837  31706  89806  62088  40151
[261]  27634  76990  37460  54157  49862  84337  64175  59382 119308  76702
[271] 103425  70344  43410 104838  62215  69304  53117  19764  86680  84105
[281]  77945  89113  91005  40248  64187  50857  56613  62792  72535
> if (par2 == 'none') {
+ m <- ctree(as.formula(paste(colnames(x)[par1],' ~ .',sep='')),data = x)
+ }
> 
> #Note: the /var/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/wessaorg/rcomp/createtable")
> 
> if (par2 != 'none') {
+ m <- ctree(as.formula(paste('as.factor(',colnames(x)[par1],') ~ .',sep='')),data = x)
+ if (par4=='yes') {
+ a<-table.start()
+ a<-table.row.start(a)
+ a<-table.element(a,'10-Fold Cross Validation',3+2*par3,TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'',1,TRUE)
+ a<-table.element(a,'Prediction (training)',par3+1,TRUE)
+ a<-table.element(a,'Prediction (testing)',par3+1,TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'Actual',1,TRUE)
+ for (jjj in 1:par3) a<-table.element(a,paste('C',jjj,sep=''),1,TRUE)
+ a<-table.element(a,'CV',1,TRUE)
+ for (jjj in 1:par3) a<-table.element(a,paste('C',jjj,sep=''),1,TRUE)
+ a<-table.element(a,'CV',1,TRUE)
+ a<-table.row.end(a)
+ for (i in 1:10) {
+ ind <- sample(2, nrow(x), replace=T, prob=c(0.9,0.1))
+ m.ct <- ctree(as.formula(paste('as.factor(',colnames(x)[par1],') ~ .',sep='')),data =x[ind==1,])
+ if (i==1) {
+ m.ct.i.pred <- predict(m.ct, newdata=x[ind==1,])
+ m.ct.i.actu <- x[ind==1,par1]
+ m.ct.x.pred <- predict(m.ct, newdata=x[ind==2,])
+ m.ct.x.actu <- x[ind==2,par1]
+ } else {
+ m.ct.i.pred <- c(m.ct.i.pred,predict(m.ct, newdata=x[ind==1,]))
+ m.ct.i.actu <- c(m.ct.i.actu,x[ind==1,par1])
+ m.ct.x.pred <- c(m.ct.x.pred,predict(m.ct, newdata=x[ind==2,]))
+ m.ct.x.actu <- c(m.ct.x.actu,x[ind==2,par1])
+ }
+ }
+ print(m.ct.i.tab <- table(m.ct.i.actu,m.ct.i.pred))
+ numer <- 0
+ for (i in 1:par3) {
+ print(m.ct.i.tab[i,i] / sum(m.ct.i.tab[i,]))
+ numer <- numer + m.ct.i.tab[i,i]
+ }
+ print(m.ct.i.cp <- numer / sum(m.ct.i.tab))
+ print(m.ct.x.tab <- table(m.ct.x.actu,m.ct.x.pred))
+ numer <- 0
+ for (i in 1:par3) {
+ print(m.ct.x.tab[i,i] / sum(m.ct.x.tab[i,]))
+ numer <- numer + m.ct.x.tab[i,i]
+ }
+ print(m.ct.x.cp <- numer / sum(m.ct.x.tab))
+ for (i in 1:par3) {
+ a<-table.row.start(a)
+ a<-table.element(a,paste('C',i,sep=''),1,TRUE)
+ for (jjj in 1:par3) a<-table.element(a,m.ct.i.tab[i,jjj])
+ a<-table.element(a,round(m.ct.i.tab[i,i]/sum(m.ct.i.tab[i,]),4))
+ for (jjj in 1:par3) a<-table.element(a,m.ct.x.tab[i,jjj])
+ a<-table.element(a,round(m.ct.x.tab[i,i]/sum(m.ct.x.tab[i,]),4))
+ a<-table.row.end(a)
+ }
+ a<-table.row.start(a)
+ a<-table.element(a,'Overall',1,TRUE)
+ for (jjj in 1:par3) a<-table.element(a,'-')
+ a<-table.element(a,round(m.ct.i.cp,4))
+ for (jjj in 1:par3) a<-table.element(a,'-')
+ a<-table.element(a,round(m.ct.x.cp,4))
+ a<-table.row.end(a)
+ a<-table.end(a)
+ table.save(a,file="/var/wessaorg/rcomp/tmp/1sn531354894935.tab") 
+ }
+ }
> m

	 Conditional inference tree with 10 terminal nodes

Response:  time 
Inputs:  pageviews, compinfo, comppr, sharedcomp, bloggedcomp, compreviewed 
Number of observations:  289 

1) pageviews <= 1613; criterion = 1, statistic = 229.566
  2) pageviews <= 1212; criterion = 1, statistic = 126.863
    3) pageviews <= 800; criterion = 1, statistic = 72.955
      4) pageviews <= 578; criterion = 1, statistic = 31.889
        5) pageviews <= 391; criterion = 0.998, statistic = 13.203
          6)*  weights = 10 
        5) pageviews > 391
          7)*  weights = 15 
      4) pageviews > 578
        8)*  weights = 36 
    3) pageviews > 800
      9) bloggedcomp <= 42; criterion = 1, statistic = 21.332
        10)*  weights = 53 
      9) bloggedcomp > 42
        11)*  weights = 11 
  2) pageviews > 1212
    12) compreviewed <= 24; criterion = 1, statistic = 20.512
      13)*  weights = 31 
    12) compreviewed > 24
      14)*  weights = 39 
1) pageviews > 1613
  15) compinfo <= 967; criterion = 1, statistic = 46.779
    16) compreviewed <= 25; criterion = 1, statistic = 22.209
      17)*  weights = 17 
    16) compreviewed > 25
      18)*  weights = 60 
  15) compinfo > 967
    19)*  weights = 17 
> postscript(file="/var/wessaorg/rcomp/tmp/2mwre1354894935.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> plot(m)
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/3z5gm1354894935.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> plot(x[,par1] ~ as.factor(where(m)),main='Response by Terminal Node',xlab='Terminal Node',ylab='Response')
> dev.off()
null device 
          1 
> if (par2 == 'none') {
+ forec <- predict(m)
+ result <- as.data.frame(cbind(x[,par1],forec,x[,par1]-forec))
+ colnames(result) <- c('Actuals','Forecasts','Residuals')
+ print(result)
+ }
    Actuals Forecasts     Residuals
1    210907 158293.33   52613.66667
2    120982 125386.18   -4404.18182
3    176508 158293.33   18214.66667
4    179321 222620.37  -43299.36667
5    123185 125386.18   -2201.18182
6     52746  39424.13   13321.86667
7    385534 322592.29   62941.70588
8     33170  19387.40   13782.60000
9    101645  84318.92   17326.07547
10   149061 158293.33   -9232.33333
11   165446 158293.33    7152.66667
12   237213 222620.37   14592.63333
13   173326 158293.33   15032.66667
14   133131 158293.33  -25162.33333
15   258873 222620.37   36252.63333
16   180083 158293.33   21789.66667
17   324799 322592.29    2206.70588
18   230964 222620.37    8343.63333
19   236785 222620.37   14164.63333
20   135473 125386.18   10086.81818
21   202925 222620.37  -19695.36667
22   215147 222620.37   -7473.36667
23   344297 222620.37  121676.63333
24   153935 158293.33   -4358.33333
25   132943 158293.33  -25350.33333
26   174724 222620.37  -47896.36667
27   174415 222620.37  -48205.36667
28   225548 222620.37    2927.63333
29   223632 222620.37    1011.63333
30   124817 158293.33  -33476.33333
31   221698 222620.37    -922.36667
32   210767 222620.37  -11853.36667
33   170266 158293.33   11972.66667
34   260561 222620.37   37940.63333
35    84853  84318.92     534.07547
36   294424 322592.29  -28168.29412
37   101011  84318.92   16692.07547
38   215641 222620.37   -6979.36667
39   325107 222620.37  102486.63333
40     7176  19387.40  -12211.40000
41   167542 158293.33    9248.66667
42   106408  84318.92   22089.07547
43    96560 153353.29  -56793.29412
44   265769 222620.37   43148.63333
45   269651 322592.29  -52941.29412
46   149112 158293.33   -9181.33333
47   175824 153353.29   22470.70588
48   152871 158293.33   -5422.33333
49   111665  84318.92   27346.07547
50   116408 222620.37 -106212.36667
51   362301 322592.29   39708.70588
52    78800  84318.92   -5518.92453
53   183167 222620.37  -39453.36667
54   277965 322592.29  -44627.29412
55   150629 222620.37  -71991.36667
56   168809 158293.33   10515.66667
57    24188  39424.13  -15236.13333
58   329267 222620.37  106646.63333
59    65029  60147.00    4882.00000
60   101097  84318.92   16778.07547
61   218946 222620.37   -3674.36667
62   244052 222620.37   21431.63333
63   341570 322592.29   18977.70588
64   103597  84318.92   19278.07547
65   233328 222620.37   10707.63333
66   256462 222620.37   33841.63333
67   206161 222620.37  -16459.36667
68   311473 322592.29  -11119.29412
69   235800 153353.29   82446.70588
70   177939 222620.37  -44681.36667
71   207176 222620.37  -15444.36667
72   196553 158293.33   38259.66667
73   174184 158293.33   15890.66667
74   143246 158293.33  -15047.33333
75   187559 222620.37  -35061.36667
76   187681 222620.37  -34939.36667
77   119016 112484.87    6531.12903
78   182192 222620.37  -40428.36667
79    73566  84318.92  -10752.92453
80   194979 222620.37  -27641.36667
81   167488 158293.33    9194.66667
82   143756 158293.33  -14537.33333
83   275541 222620.37   52920.63333
84   243199 222620.37   20578.63333
85   182999 158293.33   24705.66667
86   135649 158293.33  -22644.33333
87   152299 158293.33   -5994.33333
88   120221 112484.87    7736.12903
89   346485 322592.29   23892.70588
90   145790 158293.33  -12503.33333
91   193339 222620.37  -29281.36667
92    80953 125386.18  -44433.18182
93   122774 153353.29  -30579.29412
94   130585 125386.18    5198.81818
95   112611 125386.18  -12775.18182
96   286468 322592.29  -36124.29412
97   241066 222620.37   18445.63333
98   148446 222620.37  -74174.36667
99   204713 222620.37  -17907.36667
100  182079 222620.37  -40541.36667
101  140344 158293.33  -17949.33333
102  220516 222620.37   -2104.36667
103  243060 222620.37   20439.63333
104  162765 158293.33    4471.66667
105  182613 158293.33   24319.66667
106  232138 222620.37    9517.63333
107  265318 222620.37   42697.63333
108   85574 112484.87  -26910.87097
109  310839 322592.29  -11753.29412
110  225060 222620.37    2439.63333
111  232317 222620.37    9696.63333
112  144966 158293.33  -13327.33333
113   43287  60147.00  -16860.00000
114  155754 153353.29    2400.70588
115  164709 158293.33    6415.66667
116  201940 222620.37  -20680.36667
117  235454 222620.37   12833.63333
118  220801 153353.29   67447.70588
119   99466  84318.92   15147.07547
120   92661 112484.87  -19823.87097
121  133328 112484.87   20843.12903
122   61361  84318.92  -22957.92453
123  125930 153353.29  -27423.29412
124  100750 158293.33  -57543.33333
125  224549 158293.33   66255.66667
126   82316  84318.92   -2002.92453
127  102010  84318.92   17691.07547
128  101523 112484.87  -10961.87097
129  243511 222620.37   20890.63333
130   22938  19387.40    3550.60000
131   41566  60147.00  -18581.00000
132  152474 158293.33   -5819.33333
133   61857  39424.13   22432.86667
134   99923 153353.29  -53430.29412
135  132487 158293.33  -25806.33333
136  317394 322592.29   -5198.29412
137   21054  19387.40    1666.60000
138  209641 153353.29   56287.70588
139   22648  60147.00  -37499.00000
140   31414  39424.13   -8010.13333
141   46698  60147.00  -13449.00000
142  131698 153353.29  -21655.29412
143   91735  84318.92    7416.07547
144  244749 222620.37   22128.63333
145  184510 158293.33   26216.66667
146   79863 112484.87  -32621.87097
147  128423  84318.92   44104.07547
148   97839  84318.92   13520.07547
149   38214  39424.13   -1210.13333
150  151101 158293.33   -7192.33333
151  272458 222620.37   49837.63333
152  172494 125386.18   47107.81818
153  108043 112484.87   -4441.87097
154  328107 322592.29    5514.70588
155  250579 222620.37   27958.63333
156  351067 322592.29   28474.70588
157  158015 125386.18   32628.81818
158   98866  84318.92   14547.07547
159   85439  84318.92    1120.07547
160  229242 222620.37    6621.63333
161  351619 322592.29   29026.70588
162   84207  84318.92    -111.92453
163  120445 112484.87    7960.12903
164  324598 322592.29    2005.70588
165  131069 158293.33  -27224.33333
166  204271 158293.33   45977.66667
167  165543 222620.37  -57077.36667
168  141722 158293.33  -16571.33333
169  116048 125386.18   -9338.18182
170  250047 153353.29   96693.70588
171  299775 322592.29  -22817.29412
172  195838 222620.37  -26782.36667
173  173260 153353.29   19906.70588
174  254488 222620.37   31867.63333
175  104389 158293.33  -53904.33333
176  136084 112484.87   23599.12903
177  199476 222620.37  -23144.36667
178   92499  84318.92    8180.07547
179  224330 222620.37    1709.63333
180  135781 112484.87   23296.12903
181   74408 112484.87  -38076.87097
182   81240 112484.87  -31244.87097
183   14688  19387.40   -4699.40000
184  181633 222620.37  -40987.36667
185  271856 222620.37   49235.63333
186    7199  19387.40  -12188.40000
187   46660  39424.13    7235.86667
188   17547  19387.40   -1840.40000
189  133368 153353.29  -19985.29412
190   95227  84318.92   10908.07547
191  152601 112484.87   40116.12903
192   98146  84318.92   13827.07547
193   79619 112484.87  -32865.87097
194   59194 112484.87  -53290.87097
195  139942 112484.87   27457.12903
196  118612 112484.87    6127.12903
197   72880  84318.92  -11438.92453
198   65475  60147.00    5328.00000
199   99643 112484.87  -12841.87097
200   71965  60147.00   11818.00000
201   77272  60147.00   17125.00000
202   49289  39424.13    9864.86667
203  135131 112484.87   22646.12903
204  108446 112484.87   -4038.87097
205   89746 112484.87  -22738.87097
206   44296  39424.13    4871.86667
207   77648  84318.92   -6670.92453
208  181528 112484.87   69043.12903
209  134019 112484.87   21534.12903
210  124064 125386.18   -1322.18182
211   92630  84318.92    8311.07547
212  121848 112484.87    9363.12903
213   52915  60147.00   -7232.00000
214   81872  60147.00   21725.00000
215   58981  84318.92  -25337.92453
216   53515  60147.00   -6632.00000
217   60812  84318.92  -23506.92453
218   56375  60147.00   -3772.00000
219   65490  60147.00    5343.00000
220   80949  60147.00   20802.00000
221   76302  84318.92   -8016.92453
222  104011 112484.87   -8473.87097
223   98104 153353.29  -55249.29412
224   67989  60147.00    7842.00000
225   30989  19387.40   11601.60000
226  135458 153353.29  -17895.29412
227   73504  60147.00   13357.00000
228   63123  84318.92  -21195.92453
229   61254  84318.92  -23064.92453
230   74914 112484.87  -37570.87097
231   31774  60147.00  -28373.00000
232   81437  84318.92   -2881.92453
233   87186 112484.87  -25298.87097
234   50090  60147.00  -10057.00000
235   65745  84318.92  -18573.92453
236   56653  60147.00   -3494.00000
237  158399 112484.87   45914.12903
238   46455  60147.00  -13692.00000
239   73624  84318.92  -10694.92453
240   38395  60147.00  -21752.00000
241   91899  84318.92    7580.07547
242  139526 153353.29  -13827.29412
243   52164  60147.00   -7983.00000
244   51567  60147.00   -8580.00000
245   70551  84318.92  -13767.92453
246   84856  84318.92     537.07547
247  102538 153353.29  -50815.29412
248   86678  60147.00   26531.00000
249   85709  84318.92    1390.07547
250   34662  39424.13   -4762.13333
251  150580 112484.87   38095.12903
252   99611  84318.92   15292.07547
253   19349  19387.40     -38.40000
254   99373  84318.92   15054.07547
255   86230  84318.92    1911.07547
256   30837  39424.13   -8587.13333
257   31706  39424.13   -7718.13333
258   89806  84318.92    5487.07547
259   62088  60147.00    1941.00000
260   40151  39424.13     726.86667
261   27634  39424.13  -11790.13333
262   76990  84318.92   -7328.92453
263   37460  39424.13   -1964.13333
264   54157  60147.00   -5990.00000
265   49862  60147.00  -10285.00000
266   84337  84318.92      18.07547
267   64175  84318.92  -20143.92453
268   59382  60147.00    -765.00000
269  119308  60147.00   59161.00000
270   76702  84318.92   -7616.92453
271  103425 112484.87   -9059.87097
272   70344  60147.00   10197.00000
273   43410  60147.00  -16737.00000
274  104838 125386.18  -20548.18182
275   62215  60147.00    2068.00000
276   69304  84318.92  -15014.92453
277   53117  60147.00   -7030.00000
278   19764  19387.40     376.60000
279   86680  84318.92    2361.07547
280   84105  60147.00   23958.00000
281   77945  84318.92   -6373.92453
282   89113  84318.92    4794.07547
283   91005  84318.92    6686.07547
284   40248  39424.13     823.86667
285   64187  60147.00    4040.00000
286   50857  84318.92  -33461.92453
287   56613  84318.92  -27705.92453
288   62792  60147.00    2645.00000
289   72535  84318.92  -11783.92453
> if (par2 != 'none') {
+ print(cbind(as.factor(x[,par1]),predict(m)))
+ myt <- table(as.factor(x[,par1]),predict(m))
+ print(myt)
+ }
> postscript(file="/var/wessaorg/rcomp/tmp/4fel21354894935.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> if(par2=='none') {
+ op <- par(mfrow=c(2,2))
+ plot(density(result$Actuals),main='Kernel Density Plot of Actuals')
+ plot(density(result$Residuals),main='Kernel Density Plot of Residuals')
+ plot(result$Forecasts,result$Actuals,main='Actuals versus Predictions',xlab='Predictions',ylab='Actuals')
+ plot(density(result$Forecasts),main='Kernel Density Plot of Predictions')
+ par(op)
+ }
> if(par2!='none') {
+ plot(myt,main='Confusion Matrix',xlab='Actual',ylab='Predicted')
+ }
> dev.off()
null device 
          1 
> if (par2 == 'none') {
+ detcoef <- cor(result$Forecasts,result$Actuals)
+ a<-table.start()
+ a<-table.row.start(a)
+ a<-table.element(a,'Goodness of Fit',2,TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'Correlation',1,TRUE)
+ a<-table.element(a,round(detcoef,4))
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'R-squared',1,TRUE)
+ a<-table.element(a,round(detcoef*detcoef,4))
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'RMSE',1,TRUE)
+ a<-table.element(a,round(sqrt(mean((result$Residuals)^2)),4))
+ a<-table.row.end(a)
+ a<-table.end(a)
+ table.save(a,file="/var/wessaorg/rcomp/tmp/5i32h1354894935.tab") 
+ a<-table.start()
+ a<-table.row.start(a)
+ a<-table.element(a,'Actuals, Predictions, and Residuals',4,TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'#',header=TRUE)
+ a<-table.element(a,'Actuals',header=TRUE)
+ a<-table.element(a,'Forecasts',header=TRUE)
+ a<-table.element(a,'Residuals',header=TRUE)
+ a<-table.row.end(a)
+ for (i in 1:length(result$Actuals)) {
+ a<-table.row.start(a)
+ a<-table.element(a,i,header=TRUE)
+ a<-table.element(a,result$Actuals[i])
+ a<-table.element(a,result$Forecasts[i])
+ a<-table.element(a,result$Residuals[i])
+ a<-table.row.end(a)
+ }
+ a<-table.end(a)
+ table.save(a,file="/var/wessaorg/rcomp/tmp/6mf3f1354894935.tab") 
+ }
> if (par2 != 'none') {
+ a<-table.start()
+ a<-table.row.start(a)
+ a<-table.element(a,'Confusion Matrix (predicted in columns / actuals in rows)',par3+1,TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'',1,TRUE)
+ for (i in 1:par3) {
+ a<-table.element(a,paste('C',i,sep=''),1,TRUE)
+ }
+ a<-table.row.end(a)
+ for (i in 1:par3) {
+ a<-table.row.start(a)
+ a<-table.element(a,paste('C',i,sep=''),1,TRUE)
+ for (j in 1:par3) {
+ a<-table.element(a,myt[i,j])
+ }
+ a<-table.row.end(a)
+ }
+ a<-table.end(a)
+ table.save(a,file="/var/wessaorg/rcomp/tmp/79hd41354894935.tab") 
+ }
> 
> try(system("convert tmp/2mwre1354894935.ps tmp/2mwre1354894935.png",intern=TRUE))
character(0)
> try(system("convert tmp/3z5gm1354894935.ps tmp/3z5gm1354894935.png",intern=TRUE))
character(0)
> try(system("convert tmp/4fel21354894935.ps tmp/4fel21354894935.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  8.350   0.456   8.804