R version 3.3.2 (2016-10-31) -- "Sincere Pumpkin Patch"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-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 <- c(1.386
+ ,1.388
+ ,1.391
+ ,1.394
+ ,1.402
+ ,1.414
+ ,1.419
+ ,1.419
+ ,1.42
+ ,1.424
+ ,1.431
+ ,1.438
+ ,1.44
+ ,1.437
+ ,1.437
+ ,1.439
+ ,1.439
+ ,1.438
+ ,1.437
+ ,1.437
+ ,1.43
+ ,1.43
+ ,1.432
+ ,1.427
+ ,1.422
+ ,1.417
+ ,1.416
+ ,1.414
+ ,1.414
+ ,1.415
+ ,1.414
+ ,1.414
+ ,1.415
+ ,1.415
+ ,1.414
+ ,1.414
+ ,1.415
+ ,1.42
+ ,1.425
+ ,1.425
+ ,1.425
+ ,1.424
+ ,1.426
+ ,1.428
+ ,1.429
+ ,1.433
+ ,1.435
+ ,1.44
+ ,1.442
+ ,1.442
+ ,1.444
+ ,1.444
+ ,1.445
+ ,1.448
+ ,1.447
+ ,1.445
+ ,1.445
+ ,1.447
+ ,1.447
+ ,1.447
+ ,1.446
+ ,1.447
+ ,1.446
+ ,1.445
+ ,1.447
+ ,1.447
+ ,1.447
+ ,1.445
+ ,1.448
+ ,1.45
+ ,1.449
+ ,1.452
+ ,1.454
+ ,1.455
+ ,1.459
+ ,1.465
+ ,1.463
+ ,1.463
+ ,1.465
+ ,1.467
+ ,1.466
+ ,1.466
+ ,1.466
+ ,1.462
+ ,1.463
+ ,1.464
+ ,1.464
+ ,1.464
+ ,1.461
+ ,1.462
+ ,1.458
+ ,1.458
+ ,1.459
+ ,1.459
+ ,1.458
+ ,1.457
+ ,1.458
+ ,1.456
+ ,1.456
+ ,1.456
+ ,1.457
+ ,1.457
+ ,1.455
+ ,1.456
+ ,1.455
+ ,1.455
+ ,1.457
+ ,1.457
+ ,1.456
+ ,1.451
+ ,1.449
+ ,1.445
+ ,1.445
+ ,1.447
+ ,1.446
+ ,1.446
+ ,1.444
+ ,1.446
+ ,1.445
+ ,1.445
+ ,1.444
+ ,1.446
+ ,1.446
+ ,1.445
+ ,1.446
+ ,1.442
+ ,1.442
+ ,1.443
+ ,1.439
+ ,1.434
+ ,1.429
+ ,1.427
+ ,1.425
+ ,1.425
+ ,1.426
+ ,1.426
+ ,1.426
+ ,1.424
+ ,1.425
+ ,1.423
+ ,1.423
+ ,1.425
+ ,1.423
+ ,1.422
+ ,1.418
+ ,1.417
+ ,1.415
+ ,1.415
+ ,1.415
+ ,1.417
+ ,1.416
+ ,1.416
+ ,1.416
+ ,1.407
+ ,1.408
+ ,1.409
+ ,1.403
+ ,1.4
+ ,1.397
+ ,1.397
+ ,1.391
+ ,1.391
+ ,1.392
+ ,1.388
+ ,1.386
+ ,1.384
+ ,1.381
+ ,1.376
+ ,1.375
+ ,1.377
+ ,1.376
+ ,1.374
+ ,1.371
+ ,1.369
+ ,1.366
+ ,1.367
+ ,1.367
+ ,1.367
+ ,1.367
+ ,1.367
+ ,1.367
+ ,1.365
+ ,1.366
+ ,1.367
+ ,1.367
+ ,1.371
+ ,1.375
+ ,1.38
+ ,1.385
+ ,1.385
+ ,1.387
+ ,1.391
+ ,1.394
+ ,1.396
+ ,1.397
+ ,1.402
+ ,1.402
+ ,1.404
+ ,1.408
+ ,1.415
+ ,1.419
+ ,1.429
+ ,1.434
+ ,1.435
+ ,1.437
+ ,1.437
+ ,1.437
+ ,1.436
+ ,1.438
+ ,1.436
+ ,1.436
+ ,1.437
+ ,1.437
+ ,1.437
+ ,1.436
+ ,1.437
+ ,1.435
+ ,1.435
+ ,1.437
+ ,1.437
+ ,1.44
+ ,1.445
+ ,1.451
+ ,1.45
+ ,1.45
+ ,1.452
+ ,1.454
+ ,1.457
+ ,1.459
+ ,1.466
+ ,1.469
+ ,1.469
+ ,1.47
+ ,1.474
+ ,1.475
+ ,1.475
+ ,1.475
+ ,1.473
+ ,1.473
+ ,1.475
+ ,1.476
+ ,1.476
+ ,1.476
+ ,1.476
+ ,1.475
+ ,1.474
+ ,1.476
+ ,1.476
+ ,1.477
+ ,1.476
+ ,1.478
+ ,1.478
+ ,1.478
+ ,1.479
+ ,1.477
+ ,1.476
+ ,1.475
+ ,1.476
+ ,1.473
+ ,1.473
+ ,1.475
+ ,1.474
+ ,1.474
+ ,1.472
+ ,1.472
+ ,1.471
+ ,1.471
+ ,1.472
+ ,1.472
+ ,1.473
+ ,1.469
+ ,1.468
+ ,1.467
+ ,1.467
+ ,1.483
+ ,1.487
+ ,1.486
+ ,1.484
+ ,1.485
+ ,1.484
+ ,1.484
+ ,1.486
+ ,1.486
+ ,1.486
+ ,1.493
+ ,1.496
+ ,1.498
+ ,1.499
+ ,1.5
+ ,1.501
+ ,1.502
+ ,1.5
+ ,1.497
+ ,1.494
+ ,1.494
+ ,1.496
+ ,1.493
+ ,1.492
+ ,1.489
+ ,1.489
+ ,1.484
+ ,1.485
+ ,1.486
+ ,1.486
+ ,1.485
+ ,1.483
+ ,1.485
+ ,1.481
+ ,1.482
+ ,1.483
+ ,1.483
+ ,1.483
+ ,1.48
+ ,1.479
+ ,1.475
+ ,1.475
+ ,1.477
+ ,1.475
+ ,1.474
+ ,1.472
+ ,1.473
+ ,1.471
+ ,1.471
+ ,1.473
+ ,1.475
+ ,1.474
+ ,1.472
+ ,1.474
+ ,1.472
+ ,1.472
+ ,1.474
+ ,1.474
+ ,1.473
+ ,1.471
+ ,1.472
+ ,1.468
+ ,1.468
+ ,1.47
+ ,1.468
+ ,1.466
+ ,1.463
+ ,1.46
+ ,1.452
+ ,1.452
+ ,1.454
+ ,1.449
+ ,1.446
+ ,1.445
+ ,1.447
+ ,1.443
+ ,1.443
+ ,1.445
+ ,1.445
+ ,1.442
+ ,1.437
+ ,1.437
+ ,1.434
+ ,1.434
+ ,1.436
+ ,1.434
+ ,1.434
+ ,1.433
+ ,1.434
+ ,1.433
+ ,1.433
+ ,1.435
+ ,1.441
+ ,1.443
+ ,1.442
+ ,1.444
+ ,1.443
+ ,1.443
+ ,1.444
+ ,1.445
+ ,1.449
+ ,1.449
+ ,1.453
+ ,1.452
+ ,1.452
+ ,1.455
+ ,1.454
+ ,1.454
+ ,1.453
+ ,1.453
+ ,1.451
+ ,1.452
+ ,1.455
+ ,1.454
+ ,1.455
+ ,1.453
+ ,1.457
+ ,1.456
+ ,1.456
+ ,1.458
+ ,1.459
+ ,1.458
+ ,1.457
+ ,1.461
+ ,1.46
+ ,1.46
+ ,1.462
+ ,1.463
+ ,1.462
+ ,1.462
+ ,1.464
+ ,1.464
+ ,1.464
+ ,1.466
+ ,1.469
+ ,1.474
+ ,1.473
+ ,1.478
+ ,1.477
+ ,1.478
+ ,1.48
+ ,1.482
+ ,1.482
+ ,1.48
+ ,1.482
+ ,1.48
+ ,1.48
+ ,1.482
+ ,1.481
+ ,1.479
+ ,1.474
+ ,1.475
+ ,1.468
+ ,1.467
+ ,1.468
+ ,1.466
+ ,1.464
+ ,1.462
+ ,1.464
+ ,1.46
+ ,1.46
+ ,1.461
+ ,1.458
+ ,1.455
+ ,1.453
+ ,1.45
+ ,1.444
+ ,1.444
+ ,1.446
+ ,1.444
+ ,1.443
+ ,1.437
+ ,1.435
+ ,1.432
+ ,1.43
+ ,1.434
+ ,1.434
+ ,1.433
+ ,1.432
+ ,1.434
+ ,1.435
+ ,1.436
+ ,1.436
+ ,1.438
+ ,1.442
+ ,1.44
+ ,1.442
+ ,1.438
+ ,1.437
+ ,1.44
+ ,1.439
+ ,1.439
+ ,1.437
+ ,1.433
+ ,1.422
+ ,1.422
+ ,1.424
+ ,1.418
+ ,1.409
+ ,1.402
+ ,1.395
+ ,1.391
+ ,1.391
+ ,1.392
+ ,1.392
+ ,1.391
+ ,1.389
+ ,1.391
+ ,1.39
+ ,1.391
+ ,1.393
+ ,1.391
+ ,1.392
+ ,1.387
+ ,1.384
+ ,1.381
+ ,1.382
+ ,1.383
+ ,1.384
+ ,1.383
+ ,1.381
+ ,1.384
+ ,1.382
+ ,1.382
+ ,1.384
+ ,1.384
+ ,1.384
+ ,1.382
+ ,1.389
+ ,1.394
+ ,1.397
+ ,1.396
+ ,1.398
+ ,1.399
+ ,1.398
+ ,1.4
+ ,1.396
+ ,1.397
+ ,1.399
+ ,1.398
+ ,1.399
+ ,1.399
+ ,1.401
+ ,1.398
+ ,1.399
+ ,1.401
+ ,1.399
+ ,1.395
+ ,1.392
+ ,1.393
+ ,1.39
+ ,1.39
+ ,1.391
+ ,1.392
+ ,1.392
+ ,1.39
+ ,1.391
+ ,1.39
+ ,1.391
+ ,1.393
+ ,1.394
+ ,1.393
+ ,1.392
+ ,1.396
+ ,1.394
+ ,1.394
+ ,1.396
+ ,1.397
+ ,1.396
+ ,1.395
+ ,1.396
+ ,1.395
+ ,1.396
+ ,1.397
+ ,1.399
+ ,1.396
+ ,1.398
+ ,1.404
+ ,1.407
+ ,1.408
+ ,1.411
+ ,1.411
+ ,1.412
+ ,1.416
+ ,1.419
+ ,1.418
+ ,1.419
+ ,1.421
+ ,1.421
+ ,1.421
+ ,1.42
+ ,1.427
+ ,1.426
+ ,1.429
+ ,1.429
+ ,1.429
+ ,1.429
+ ,1.427
+ ,1.429
+ ,1.426
+ ,1.428
+ ,1.428
+ ,1.429
+ ,1.428
+ ,1.427
+ ,1.429
+ ,1.427
+ ,1.428
+ ,1.43
+ ,1.43
+ ,1.429
+ ,1.426
+ ,1.421
+ ,1.418
+ ,1.418
+ ,1.419
+ ,1.419
+ ,1.419
+ ,1.417
+ ,1.418
+ ,1.417
+ ,1.419
+ ,1.421
+ ,1.424
+ ,1.429
+ ,1.427
+ ,1.43
+ ,1.427
+ ,1.428
+ ,1.43
+ ,1.429
+ ,1.433
+ ,1.437
+ ,1.442
+ ,1.441
+ ,1.443
+ ,1.444
+ ,1.445
+ ,1.444
+ ,1.442
+ ,1.444
+ ,1.443
+ ,1.444
+ ,1.446
+ ,1.445
+ ,1.444
+ ,1.442
+ ,1.444
+ ,1.442
+ ,1.442
+ ,1.444
+ ,1.443
+ ,1.442
+ ,1.439
+ ,1.441
+ ,1.435
+ ,1.436
+ ,1.438
+ ,1.434
+ ,1.426
+ ,1.421
+ ,1.418
+ ,1.412
+ ,1.413
+ ,1.414
+ ,1.413
+ ,1.412
+ ,1.411
+ ,1.414
+ ,1.41
+ ,1.411
+ ,1.413
+ ,1.413
+ ,1.412
+ ,1.411
+ ,1.413
+ ,1.411
+ ,1.413
+ ,1.414
+ ,1.41
+ ,1.409
+ ,1.407
+ ,1.408
+ ,1.405
+ ,1.405
+ ,1.407
+ ,1.407
+ ,1.407
+ ,1.404
+ ,1.401
+ ,1.392
+ ,1.393
+ ,1.394
+ ,1.392
+ ,1.391
+ ,1.388
+ ,1.39
+ ,1.386
+ ,1.388
+ ,1.39
+ ,1.389
+ ,1.386
+ ,1.382
+ ,1.38
+ ,1.376
+ ,1.378
+ ,1.379
+ ,1.379
+ ,1.378
+ ,1.378
+ ,1.381
+ ,1.379
+ ,1.38
+ ,1.382
+ ,1.384
+ ,1.384
+ ,1.383
+ ,1.391
+ ,1.392
+ ,1.394
+ ,1.396
+ ,1.398
+ ,1.402
+ ,1.403
+ ,1.408
+ ,1.405
+ ,1.406
+ ,1.408
+ ,1.408
+ ,1.407
+ ,1.405
+ ,1.409
+ ,1.406
+ ,1.407
+ ,1.408
+ ,1.408
+ ,1.407
+ ,1.405
+ ,1.407
+ ,1.403
+ ,1.404
+ ,1.405
+ ,1.404
+ ,1.399
+ ,1.396
+ ,1.398
+ ,1.399
+ ,1.4
+ ,1.401
+ ,1.402
+ ,1.402
+ ,1.4
+ ,1.404
+ ,1.403
+ ,1.404
+ ,1.404
+ ,1.404
+ ,1.44
+ ,1.438
+ ,1.441
+ ,1.438
+ ,1.439
+ ,1.44
+ ,1.44
+ ,1.437
+ ,1.434
+ ,1.437
+ ,1.435
+ ,1.435
+ ,1.437
+ ,1.434
+ ,1.431
+ ,1.426
+ ,1.426
+ ,1.423
+ ,1.425
+ ,1.426
+ ,1.427
+ ,1.424
+ ,1.422
+ ,1.427
+ ,1.424
+ ,1.426
+ ,1.427
+ ,1.429
+ ,1.427
+ ,1.425
+ ,1.427
+ ,1.425
+ ,1.426
+ ,1.428
+ ,1.427
+ ,1.427
+ ,1.425
+ ,1.428
+ ,1.426
+ ,1.428
+ ,1.429
+ ,1.434
+ ,1.434
+ ,1.435
+ ,1.439
+ ,1.437
+ ,1.438
+ ,1.439
+ ,1.439
+ ,1.439
+ ,1.438
+ ,1.442
+ ,1.441
+ ,1.441
+ ,1.442
+ ,1.442
+ ,1.44
+ ,1.438
+ ,1.439
+ ,1.435
+ ,1.435
+ ,1.437
+ ,1.438
+ ,1.434
+ ,1.429
+ ,1.43
+ ,1.424
+ ,1.424
+ ,1.425
+ ,1.423
+ ,1.42
+ ,1.415
+ ,1.413
+ ,1.409
+ ,1.409
+ ,1.411
+ ,1.41
+ ,1.406
+ ,1.4
+ ,1.4
+ ,1.4
+ ,1.401
+ ,1.402
+ ,1.407
+ ,1.408
+ ,1.407
+ ,1.41
+ ,1.41
+ ,1.411
+ ,1.412
+ ,1.418
+ ,1.419
+ ,1.416
+ ,1.419
+ ,1.417
+ ,1.418
+ ,1.42
+ ,1.419
+ ,1.419
+ ,1.417
+ ,1.419
+ ,1.417
+ ,1.418
+ ,1.42
+ ,1.421
+ ,1.421
+ ,1.425
+ ,1.429
+ ,1.428
+ ,1.428
+ ,1.429
+ ,1.43
+ ,1.429
+ ,1.428
+ ,1.431
+ ,1.429
+ ,1.43
+ ,1.431
+ ,1.431
+ ,1.429
+ ,1.428
+ ,1.428
+ ,1.424
+ ,1.424
+ ,1.425
+ ,1.422
+ ,1.42
+ ,1.418
+ ,1.42
+ ,1.418
+ ,1.419
+ ,1.421
+ ,1.421
+ ,1.421
+ ,1.423
+ ,1.428
+ ,1.427
+ ,1.428
+ ,1.43
+ ,1.431
+ ,1.429
+ ,1.428
+ ,1.43
+ ,1.428
+ ,1.429
+ ,1.43
+ ,1.43
+ ,1.429
+ ,1.427
+ ,1.43
+ ,1.428
+ ,1.429
+ ,1.43
+ ,1.431
+ ,1.427
+ ,1.424
+ ,1.422
+ ,1.419
+ ,1.42
+ ,1.421
+ ,1.421
+ ,1.421
+ ,1.42
+ ,1.423
+ ,1.427
+ ,1.428
+ ,1.43
+ ,1.43
+ ,1.429
+ ,1.434
+ ,1.442
+ ,1.445
+ ,1.447
+ ,1.448
+ ,1.448
+ ,1.448
+ ,1.447
+ ,1.449
+ ,1.447
+ ,1.448
+ ,1.449
+ ,1.449
+ ,1.444
+ ,1.439
+ ,1.439
+ ,1.432
+ ,1.433
+ ,1.435
+ ,1.433
+ ,1.427
+ ,1.424
+ ,1.424
+ ,1.42
+ ,1.421
+ ,1.422
+ ,1.42
+ ,1.42
+ ,1.419
+ ,1.42
+ ,1.414
+ ,1.415
+ ,1.416
+ ,1.416
+ ,1.418
+ ,1.418
+ ,1.419
+ ,1.417
+ ,1.418
+ ,1.42
+ ,1.419
+ ,1.422
+ ,1.425
+ ,1.427
+ ,1.425
+ ,1.426
+ ,1.427
+ ,1.425
+ ,1.421
+ ,1.418
+ ,1.42
+ ,1.419
+ ,1.421
+ ,1.421
+ ,1.421
+ ,1.419
+ ,1.417
+ ,1.419
+ ,1.418
+ ,1.418
+ ,1.419
+ ,1.415
+ ,1.412
+ ,1.408
+ ,1.41
+ ,1.409
+ ,1.409
+ ,1.411
+ ,1.411
+ ,1.414
+ ,1.413
+ ,1.421
+ ,1.421
+ ,1.421
+ ,1.424
+ ,1.425
+ ,1.426
+ ,1.422
+ ,1.424
+ ,1.427
+ ,1.426
+ ,1.428
+ ,1.427
+ ,1.427
+ ,1.422
+ ,1.419
+ ,1.415
+ ,1.415
+ ,1.416
+ ,1.413
+ ,1.409
+ ,1.408
+ ,1.409
+ ,1.408
+ ,1.407
+ ,1.41
+ ,1.409
+ ,1.407
+ ,1.404
+ ,1.405
+ ,1.406
+ ,1.406
+ ,1.408
+ ,1.408
+ ,1.412
+ ,1.414
+ ,1.411
+ ,1.406
+ ,1.406
+ ,1.409
+ ,1.405
+ ,1.4
+ ,1.393
+ ,1.387
+ ,1.378
+ ,1.378
+ ,1.378
+ ,1.374
+ ,1.37
+ ,1.361
+ ,1.355
+ ,1.351
+ ,1.35
+ ,1.352
+ ,1.35
+ ,1.349
+ ,1.351
+ ,1.356
+ ,1.355
+ ,1.355
+ ,1.357
+ ,1.357
+ ,1.357
+ ,1.355
+ ,1.358
+ ,1.357
+ ,1.358
+ ,1.361
+ ,1.365
+ ,1.366
+ ,1.364
+ ,1.366
+ ,1.364
+ ,1.364
+ ,1.367
+ ,1.372
+ ,1.373
+ ,1.371
+ ,1.374
+ ,1.365
+ ,1.365
+ ,1.368
+ ,1.362
+ ,1.36
+ ,1.356
+ ,1.357
+ ,1.352
+ ,1.352
+ ,1.354
+ ,1.347
+ ,1.344
+ ,1.342
+ ,1.342
+ ,1.33
+ ,1.33
+ ,1.332
+ ,1.323
+ ,1.319
+ ,1.313
+ ,1.313
+ ,1.304
+ ,1.304
+ ,1.305
+ ,1.299
+ ,1.291
+ ,1.283
+ ,1.28
+ ,1.27
+ ,1.27
+ ,1.271
+ ,1.261
+ ,1.252
+ ,1.24
+ ,1.236
+ ,1.228
+ ,1.228
+ ,1.229
+ ,1.227
+ ,1.225
+ ,1.221
+ ,1.223
+ ,1.222
+ ,1.222
+ ,1.223
+ ,1.217
+ ,1.212
+ ,1.207
+ ,1.208
+ ,1.206
+ ,1.206
+ ,1.208
+ ,1.207
+ ,1.206
+ ,1.203
+ ,1.205
+ ,1.198
+ ,1.2
+ ,1.201
+ ,1.197
+ ,1.191
+ ,1.186
+ ,1.187
+ ,1.184
+ ,1.185
+ ,1.187
+ ,1.187
+ ,1.186
+ ,1.185
+ ,1.186
+ ,1.184
+ ,1.184
+ ,1.187
+ ,1.188
+ ,1.189
+ ,1.188
+ ,1.19
+ ,1.19
+ ,1.192
+ ,1.194
+ ,1.201
+ ,1.21
+ ,1.219
+ ,1.222
+ ,1.225
+ ,1.225
+ ,1.229
+ ,1.237
+ ,1.249
+ ,1.255
+ ,1.26
+ ,1.267
+ ,1.268
+ ,1.271
+ ,1.276
+ ,1.279
+ ,1.28
+ ,1.282
+ ,1.28
+ ,1.28
+ ,1.282
+ ,1.283
+ ,1.283
+ ,1.282
+ ,1.289
+ ,1.291
+ ,1.291
+ ,1.293
+ ,1.293
+ ,1.293
+ ,1.291
+ ,1.293
+ ,1.292
+ ,1.292
+ ,1.294
+ ,1.294
+ ,1.294
+ ,1.292
+ ,1.293
+ ,1.292
+ ,1.291
+ ,1.293
+ ,1.292
+ ,1.287
+ ,1.283
+ ,1.284
+ ,1.278
+ ,1.278
+ ,1.28
+ ,1.279
+ ,1.28
+ ,1.277
+ ,1.278
+ ,1.276
+ ,1.277
+ ,1.278
+ ,1.277
+ ,1.275
+ ,1.273
+ ,1.27
+ ,1.268
+ ,1.268
+ ,1.269
+ ,1.27
+ ,1.269
+ ,1.268
+ ,1.268
+ ,1.261
+ ,1.262
+ ,1.263
+ ,1.264
+ ,1.269
+ ,1.272
+ ,1.282
+ ,1.288
+ ,1.289
+ ,1.291
+ ,1.294
+ ,1.299
+ ,1.3
+ ,1.304
+ ,1.31
+ ,1.31
+ ,1.309
+ ,1.312
+ ,1.312
+ ,1.312
+ ,1.31
+ ,1.306
+ ,1.306
+ ,1.307
+ ,1.306
+ ,1.307
+ ,1.312
+ ,1.313
+ ,1.312
+ ,1.312
+ ,1.313
+ ,1.31
+ ,1.309
+ ,1.309
+ ,1.311
+ ,1.31
+ ,1.31
+ ,1.312
+ ,1.313
+ ,1.313
+ ,1.308
+ ,1.309
+ ,1.308
+ ,1.309
+ ,1.309
+ ,1.31
+ ,1.313
+ ,1.312
+ ,1.31
+ ,1.306
+ ,1.306
+ ,1.304
+ ,1.309
+ ,1.31
+ ,1.31
+ ,1.309
+ ,1.303
+ ,1.302
+ ,1.304
+ ,1.302
+ ,1.301
+ ,1.302
+ ,1.302
+ ,1.301
+ ,1.302
+ ,1.304
+ ,1.305
+ ,1.305
+ ,1.304
+ ,1.305
+ ,1.299
+ ,1.299
+ ,1.301
+ ,1.297
+ ,1.296
+ ,1.295
+ ,1.295
+ ,1.293
+ ,1.293
+ ,1.295
+ ,1.295
+ ,1.295
+ ,1.293
+ ,1.294
+ ,1.294
+ ,1.293
+ ,1.295
+ ,1.294
+ ,1.289
+ ,1.285
+ ,1.284
+ ,1.28
+ ,1.28
+ ,1.281
+ ,1.28
+ ,1.275
+ ,1.267
+ ,1.262
+ ,1.254
+ ,1.254
+ ,1.255
+ ,1.251
+ ,1.245
+ ,1.242
+ ,1.242
+ ,1.239
+ ,1.24
+ ,1.239
+ ,1.233
+ ,1.225
+ ,1.22
+ ,1.22
+ ,1.218
+ ,1.219
+ ,1.22
+ ,1.218
+ ,1.215
+ ,1.21
+ ,1.209
+ ,1.203
+ ,1.203
+ ,1.204
+ ,1.201
+ ,1.197
+ ,1.194
+ ,1.194
+ ,1.193
+ ,1.193
+ ,1.194
+ ,1.194
+ ,1.194
+ ,1.192
+ ,1.19
+ ,1.181
+ ,1.181
+ ,1.181
+ ,1.172
+ ,1.163
+ ,1.157
+ ,1.156
+ ,1.155
+ ,1.155
+ ,1.157
+ ,1.165
+ ,1.174
+ ,1.175
+ ,1.181
+ ,1.186
+ ,1.188
+ ,1.19
+ ,1.196
+ ,1.197
+ ,1.199
+ ,1.2
+ ,1.199
+ ,1.199
+ ,1.201
+ ,1.2
+ ,1.198
+ ,1.198
+ ,1.199
+ ,1.193
+ ,1.194
+ ,1.196
+ ,1.194
+ ,1.193
+ ,1.193
+ ,1.189
+ ,1.185
+ ,1.184
+ ,1.186
+ ,1.185
+ ,1.182
+ ,1.181
+ ,1.184
+ ,1.181
+ ,1.181
+ ,1.185
+ ,1.184
+ ,1.187
+ ,1.19
+ ,1.193
+ ,1.196
+ ,1.197
+ ,1.2
+ ,1.201
+ ,1.199
+ ,1.195
+ ,1.193
+ ,1.187
+ ,1.188
+ ,1.19
+ ,1.186
+ ,1.179
+ ,1.172
+ ,1.172
+ ,1.17
+ ,1.17
+ ,1.172
+ ,1.172
+ ,1.17
+ ,1.169
+ ,1.175
+ ,1.173
+ ,1.173
+ ,1.175
+ ,1.18
+ ,1.187
+ ,1.196
+ ,1.199
+ ,1.198
+ ,1.198
+ ,1.202
+ ,1.206
+ ,1.207
+ ,1.208
+ ,1.208
+ ,1.203
+ ,1.204
+ ,1.206
+ ,1.199
+ ,1.195
+ ,1.191
+ ,1.191
+ ,1.186
+ ,1.186
+ ,1.188
+ ,1.184
+ ,1.181
+ ,1.181
+ ,1.182
+ ,1.18
+ ,1.18
+ ,1.182
+ ,1.177
+ ,1.172
+ ,1.168
+ ,1.162
+ ,1.157
+ ,1.157
+ ,1.159
+ ,1.158
+ ,1.151
+ ,1.143
+ ,1.139
+ ,1.129
+ ,1.129
+ ,1.131
+ ,1.122
+ ,1.113
+ ,1.107
+ ,1.101
+ ,1.094
+ ,1.094
+ ,1.096
+ ,1.087
+ ,1.081
+ ,1.078
+ ,1.071
+ ,1.07
+ ,1.07
+ ,1.071
+ ,1.068
+ ,1.064
+ ,1.067)
> par3 = 'additive'
> par2 = 'Single'
> par1 = '12'
> par3 <- 'additive'
> par2 <- 'Single'
> par1 <- '12'
> #'GNU S' R Code compiled by R2WASP v. 1.2.327 ()
> #Author: root
> #To cite this work: Wessa P., (2013), Exponential Smoothing (v1.0.5) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_exponentialsmoothing.wasp/
> #Source of accompanying publication: 
> #
> par1 <- as.numeric(par1)
> if (par2 == 'Single') K <- 1
> if (par2 == 'Double') K <- 2
> if (par2 == 'Triple') K <- par1
> nx <- length(x)
> nxmK <- nx - K
> x <- ts(x, frequency = par1)
> if (par2 == 'Single') fit <- HoltWinters(x, gamma=F, beta=F)
> if (par2 == 'Double') fit <- HoltWinters(x, gamma=F)
> if (par2 == 'Triple') fit <- HoltWinters(x, seasonal=par3)
> fit
Holt-Winters exponential smoothing without trend and without seasonal component.

Call:
HoltWinters(x = x, beta = F, gamma = F)

Smoothing parameters:
 alpha: 0.9999544
 beta : FALSE
 gamma: FALSE

Coefficients:
   [,1]
a 1.067
> myresid <- x - fit$fitted[,'xhat']
> postscript(file="/var/wessaorg/rcomp/tmp/1u8q21483909853.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> op <- par(mfrow=c(2,1))
> plot(fit,ylab='Observed (black) / Fitted (red)',main='Interpolation Fit of Exponential Smoothing')
> plot(myresid,ylab='Residuals',main='Interpolation Prediction Errors')
> par(op)
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/2rjre1483909853.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> p <- predict(fit, par1, prediction.interval=TRUE)
> np <- length(p[,1])
> plot(fit,p,ylab='Observed (black) / Fitted (red)',main='Extrapolation Fit of Exponential Smoothing')
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/3t30w1483909853.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> op <- par(mfrow = c(2,2))
> acf(as.numeric(myresid),lag.max = nx/2,main='Residual ACF')
> spectrum(myresid,main='Residals Periodogram')
> cpgram(myresid,main='Residal Cumulative Periodogram')
> qqnorm(myresid,main='Residual Normal QQ Plot')
> qqline(myresid)
> par(op)
> dev.off()
null device 
          1 
> 
> #Note: the /var/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/wessaorg/rcomp/createtable")
> 
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,'Estimated Parameters of Exponential Smoothing',2,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Parameter',header=TRUE)
> a<-table.element(a,'Value',header=TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'alpha',header=TRUE)
> a<-table.element(a,fit$alpha)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'beta',header=TRUE)
> a<-table.element(a,fit$beta)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'gamma',header=TRUE)
> a<-table.element(a,fit$gamma)
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/wessaorg/rcomp/tmp/4va1f1483909853.tab") 
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,'Interpolation Forecasts of Exponential Smoothing',4,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'t',header=TRUE)
> a<-table.element(a,'Observed',header=TRUE)
> a<-table.element(a,'Fitted',header=TRUE)
> a<-table.element(a,'Residuals',header=TRUE)
> a<-table.row.end(a)
> for (i in 1:nxmK) {
+ a<-table.row.start(a)
+ a<-table.element(a,i+K,header=TRUE)
+ a<-table.element(a,x[i+K])
+ a<-table.element(a,fit$fitted[i,'xhat'])
+ a<-table.element(a,myresid[i])
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/wessaorg/rcomp/tmp/56o8z1483909853.tab") 
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,'Extrapolation Forecasts of Exponential Smoothing',4,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'t',header=TRUE)
> a<-table.element(a,'Forecast',header=TRUE)
> a<-table.element(a,'95% Lower Bound',header=TRUE)
> a<-table.element(a,'95% Upper Bound',header=TRUE)
> a<-table.row.end(a)
> for (i in 1:np) {
+ a<-table.row.start(a)
+ a<-table.element(a,nx+i,header=TRUE)
+ a<-table.element(a,p[i,'fit'])
+ a<-table.element(a,p[i,'lwr'])
+ a<-table.element(a,p[i,'upr'])
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/wessaorg/rcomp/tmp/69q6s1483909853.tab") 
> 
> try(system("convert tmp/1u8q21483909853.ps tmp/1u8q21483909853.png",intern=TRUE))
character(0)
> try(system("convert tmp/2rjre1483909853.ps tmp/2rjre1483909853.png",intern=TRUE))
character(0)
> try(system("convert tmp/3t30w1483909853.ps tmp/3t30w1483909853.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
 44.012   0.169  44.300