FreeStatistics.org

Free Statistics

of Irreproducible Research!

Description of Statistical Computation

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_spectrum.wasp
Title produced by softwareSpectral Analysis
Date of computationFri, 24 Dec 2010 16:03:03 +0000
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2010/Dec/24/t1293206455bvn93xtsyt4iv5o.htm/, Retrieved Tue, 30 Apr 2024 03:02:55 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=115165, Retrieved Tue, 30 Apr 2024 03:02:55 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact159

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Multiple Regression] [Multiple Regressi...] [2010-12-17 13:46:43] [1251ac2db27b84d4a3ba43449388906b]
-   PD  [Multiple Regression] [Multiple Regressi...] [2010-12-17 15:41:32] [1251ac2db27b84d4a3ba43449388906b]
-   P     [Multiple Regression] [MR Paper (monthly...] [2010-12-17 16:35:58] [1251ac2db27b84d4a3ba43449388906b]
-   PD      [Multiple Regression] [MR Paper (month)] [2010-12-17 16:45:18] [1251ac2db27b84d4a3ba43449388906b]
-   P         [Multiple Regression] [MR Paper (trend)b] [2010-12-18 12:41:58] [1251ac2db27b84d4a3ba43449388906b]
- RMPD          [Classical Decomposition] [CD] [2010-12-18 14:47:46] [1251ac2db27b84d4a3ba43449388906b]
- RM              [Spectral Analysis] [Spectral Analysis...] [2010-12-18 18:17:41] [1251ac2db27b84d4a3ba43449388906b]
-   P                 [Spectral Analysis] [] [2010-12-24 16:03:03] [4f70e6cd0867f10d298e58e8e27859b5] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
14.458
13.594
17.814
20.235
21.811
21.439
21.393
19.831
20.468
21.080
21.600
17.390
17.848
19.592
21.092
20.899
25.890
24.965
22.225
20.977
22.897
22.785
22.769
19.637
20.203
20.450
23.083
21.738
26.766
25.280
22.574
22.729
21.378
22.902
24.989
21.116
15.169
15.846
20.927
18.273
22.538
15.596
14.034
11.366
14.861
15.149
13.577
13.026
13.190
13.196
15.826
14.733
16.307
15.703
14.589
12.043
15.057
14.053
12.698
10.888
10.045
11.549
13.767
12.434
13.116
14.211
12.266
12.602
15.714
13.742
12.745
10.491
10.057
10.900
11.771
11.992
11.933
14.504
11.727
11.477
13.578
11.555
11.846
11.397
10.066
10.269
14.279
13.870
13.695
14.420
11.424
9.704
12.464
14.301
13.464
9.893
11.572
12.380
16.692
16.052
16.459
14.761
13.654
13.480
18.068
16.560
14.530
10.650
11.651
13.735
13.360
17.818
20.613
16.231
13.862
12.004
17.734
15.034
12.609
12.320
10.833
11.350
13.648
14.890
16.325
18.045
15.616
11.926
16.855
15.083
12.520
12.355

Tables (Output of Computation)





Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time55 seconds
R Server'George Udny Yule' @ 72.249.76.132

\begin{tabular}{lllllllll}
\hline
Summary of computational transaction \tabularnewline
Raw Input & view raw input (R code)  \tabularnewline
Raw Output & view raw output of R engine  \tabularnewline
Computing time & 55 seconds \tabularnewline
R Server & 'George Udny Yule' @ 72.249.76.132 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=115165&T=0

[TABLE]
[ROW][C]Summary of computational transaction[/C][/ROW]
[ROW][C]Raw Input[/C][C]view raw input (R code) [/C][/ROW]
[ROW][C]Raw Output[/C][C]view raw output of R engine [/C][/ROW]
[ROW][C]Computing time[/C][C]55 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'George Udny Yule' @ 72.249.76.132[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=115165&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=115165&T=0

As an alternative you can also use a QR Code:  
QR Code


The GUIDs for individual cells are displayed in the table below:

Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time55 seconds
R Server'George Udny Yule' @ 72.249.76.132






Raw Periodogram
ParameterValue
Box-Cox transformation parameter (lambda)1
Degree of non-seasonal differencing (d)0
Degree of seasonal differencing (D)0
Seasonal Period (s)12
Frequency (Period)Spectrum
0.0074 (135)274.205195
0.0148 (67.5)123.904683
0.0222 (45)61.536764
0.0296 (33.75)10.29903
0.037 (27)5.067529
0.0444 (22.5)3.442247
0.0519 (19.2857)5.777284
0.0593 (16.875)11.269186
0.0667 (15)5.569656
0.0741 (13.5)5.489349
0.0815 (12.2727)44.886933
0.0889 (11.25)5.822954
0.0963 (10.3846)7.198393
0.1037 (9.6429)5.818261
0.1111 (9)2.508553
0.1185 (8.4375)0.010607
0.1259 (7.9412)2.853442
0.1333 (7.5)0.239436
0.1407 (7.1053)3.066153
0.1481 (6.75)0.836111
0.1556 (6.4286)18.394461
0.163 (6.1364)21.939773
0.1704 (5.8696)51.662073
0.1778 (5.625)7.282166
0.1852 (5.4)4.862474
0.1926 (5.1923)0.26397
0.2 (5)1.511483
0.2074 (4.8214)4.465788
0.2148 (4.6552)1.148821
0.2222 (4.5)2.951009
0.2296 (4.3548)3.127388
0.237 (4.2188)0.017582
0.2444 (4.0909)0.061233
0.2519 (3.9706)14.202995
0.2593 (3.8571)1.118618
0.2667 (3.75)0.751987
0.2741 (3.6486)4.462532
0.2815 (3.5526)2.095525
0.2889 (3.4615)0.803037
0.2963 (3.375)0.87938
0.3037 (3.2927)1.023119
0.3111 (3.2143)3.12649
0.3185 (3.1395)0.635779
0.3259 (3.0682)6.151054
0.3333 (3)6.750947
0.3407 (2.9348)0.276941
0.3481 (2.8723)1.221282
0.3556 (2.8125)0.157484
0.363 (2.7551)1.21057
0.3704 (2.7)0.047862
0.3778 (2.6471)1.146836
0.3852 (2.5962)0.726227
0.3926 (2.5472)0.967389
0.4 (2.5)0.082001
0.4074 (2.4545)2.805805
0.4148 (2.4107)0.130709
0.4222 (2.3684)1.760301
0.4296 (2.3276)3.66678
0.437 (2.2881)2.003551
0.4444 (2.25)0.842837
0.4519 (2.2131)0.228989
0.4593 (2.1774)0.475707
0.4667 (2.1429)0.278949
0.4741 (2.1094)0.141589
0.4815 (2.0769)0.538107
0.4889 (2.0455)2.723041
0.4963 (2.0149)7.609089

\begin{tabular}{lllllllll}
\hline
Raw Periodogram \tabularnewline
Parameter & Value \tabularnewline
Box-Cox transformation parameter (lambda) & 1 \tabularnewline
Degree of non-seasonal differencing (d) & 0 \tabularnewline
Degree of seasonal differencing (D) & 0 \tabularnewline
Seasonal Period (s) & 12 \tabularnewline
Frequency (Period) & Spectrum \tabularnewline
0.0074 (135) & 274.205195 \tabularnewline
0.0148 (67.5) & 123.904683 \tabularnewline
0.0222 (45) & 61.536764 \tabularnewline
0.0296 (33.75) & 10.29903 \tabularnewline
0.037 (27) & 5.067529 \tabularnewline
0.0444 (22.5) & 3.442247 \tabularnewline
0.0519 (19.2857) & 5.777284 \tabularnewline
0.0593 (16.875) & 11.269186 \tabularnewline
0.0667 (15) & 5.569656 \tabularnewline
0.0741 (13.5) & 5.489349 \tabularnewline
0.0815 (12.2727) & 44.886933 \tabularnewline
0.0889 (11.25) & 5.822954 \tabularnewline
0.0963 (10.3846) & 7.198393 \tabularnewline
0.1037 (9.6429) & 5.818261 \tabularnewline
0.1111 (9) & 2.508553 \tabularnewline
0.1185 (8.4375) & 0.010607 \tabularnewline
0.1259 (7.9412) & 2.853442 \tabularnewline
0.1333 (7.5) & 0.239436 \tabularnewline
0.1407 (7.1053) & 3.066153 \tabularnewline
0.1481 (6.75) & 0.836111 \tabularnewline
0.1556 (6.4286) & 18.394461 \tabularnewline
0.163 (6.1364) & 21.939773 \tabularnewline
0.1704 (5.8696) & 51.662073 \tabularnewline
0.1778 (5.625) & 7.282166 \tabularnewline
0.1852 (5.4) & 4.862474 \tabularnewline
0.1926 (5.1923) & 0.26397 \tabularnewline
0.2 (5) & 1.511483 \tabularnewline
0.2074 (4.8214) & 4.465788 \tabularnewline
0.2148 (4.6552) & 1.148821 \tabularnewline
0.2222 (4.5) & 2.951009 \tabularnewline
0.2296 (4.3548) & 3.127388 \tabularnewline
0.237 (4.2188) & 0.017582 \tabularnewline
0.2444 (4.0909) & 0.061233 \tabularnewline
0.2519 (3.9706) & 14.202995 \tabularnewline
0.2593 (3.8571) & 1.118618 \tabularnewline
0.2667 (3.75) & 0.751987 \tabularnewline
0.2741 (3.6486) & 4.462532 \tabularnewline
0.2815 (3.5526) & 2.095525 \tabularnewline
0.2889 (3.4615) & 0.803037 \tabularnewline
0.2963 (3.375) & 0.87938 \tabularnewline
0.3037 (3.2927) & 1.023119 \tabularnewline
0.3111 (3.2143) & 3.12649 \tabularnewline
0.3185 (3.1395) & 0.635779 \tabularnewline
0.3259 (3.0682) & 6.151054 \tabularnewline
0.3333 (3) & 6.750947 \tabularnewline
0.3407 (2.9348) & 0.276941 \tabularnewline
0.3481 (2.8723) & 1.221282 \tabularnewline
0.3556 (2.8125) & 0.157484 \tabularnewline
0.363 (2.7551) & 1.21057 \tabularnewline
0.3704 (2.7) & 0.047862 \tabularnewline
0.3778 (2.6471) & 1.146836 \tabularnewline
0.3852 (2.5962) & 0.726227 \tabularnewline
0.3926 (2.5472) & 0.967389 \tabularnewline
0.4 (2.5) & 0.082001 \tabularnewline
0.4074 (2.4545) & 2.805805 \tabularnewline
0.4148 (2.4107) & 0.130709 \tabularnewline
0.4222 (2.3684) & 1.760301 \tabularnewline
0.4296 (2.3276) & 3.66678 \tabularnewline
0.437 (2.2881) & 2.003551 \tabularnewline
0.4444 (2.25) & 0.842837 \tabularnewline
0.4519 (2.2131) & 0.228989 \tabularnewline
0.4593 (2.1774) & 0.475707 \tabularnewline
0.4667 (2.1429) & 0.278949 \tabularnewline
0.4741 (2.1094) & 0.141589 \tabularnewline
0.4815 (2.0769) & 0.538107 \tabularnewline
0.4889 (2.0455) & 2.723041 \tabularnewline
0.4963 (2.0149) & 7.609089 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=115165&T=1

[TABLE]
[ROW][C]Raw Periodogram[/C][/ROW]
[ROW][C]Parameter[/C][C]Value[/C][/ROW]
[ROW][C]Box-Cox transformation parameter (lambda)[/C][C]1[/C][/ROW]
[ROW][C]Degree of non-seasonal differencing (d)[/C][C]0[/C][/ROW]
[ROW][C]Degree of seasonal differencing (D)[/C][C]0[/C][/ROW]
[ROW][C]Seasonal Period (s)[/C][C]12[/C][/ROW]
[ROW][C]Frequency (Period)[/C][C]Spectrum[/C][/ROW]
[ROW][C]0.0074 (135)[/C][C]274.205195[/C][/ROW]
[ROW][C]0.0148 (67.5)[/C][C]123.904683[/C][/ROW]
[ROW][C]0.0222 (45)[/C][C]61.536764[/C][/ROW]
[ROW][C]0.0296 (33.75)[/C][C]10.29903[/C][/ROW]
[ROW][C]0.037 (27)[/C][C]5.067529[/C][/ROW]
[ROW][C]0.0444 (22.5)[/C][C]3.442247[/C][/ROW]
[ROW][C]0.0519 (19.2857)[/C][C]5.777284[/C][/ROW]
[ROW][C]0.0593 (16.875)[/C][C]11.269186[/C][/ROW]
[ROW][C]0.0667 (15)[/C][C]5.569656[/C][/ROW]
[ROW][C]0.0741 (13.5)[/C][C]5.489349[/C][/ROW]
[ROW][C]0.0815 (12.2727)[/C][C]44.886933[/C][/ROW]
[ROW][C]0.0889 (11.25)[/C][C]5.822954[/C][/ROW]
[ROW][C]0.0963 (10.3846)[/C][C]7.198393[/C][/ROW]
[ROW][C]0.1037 (9.6429)[/C][C]5.818261[/C][/ROW]
[ROW][C]0.1111 (9)[/C][C]2.508553[/C][/ROW]
[ROW][C]0.1185 (8.4375)[/C][C]0.010607[/C][/ROW]
[ROW][C]0.1259 (7.9412)[/C][C]2.853442[/C][/ROW]
[ROW][C]0.1333 (7.5)[/C][C]0.239436[/C][/ROW]
[ROW][C]0.1407 (7.1053)[/C][C]3.066153[/C][/ROW]
[ROW][C]0.1481 (6.75)[/C][C]0.836111[/C][/ROW]
[ROW][C]0.1556 (6.4286)[/C][C]18.394461[/C][/ROW]
[ROW][C]0.163 (6.1364)[/C][C]21.939773[/C][/ROW]
[ROW][C]0.1704 (5.8696)[/C][C]51.662073[/C][/ROW]
[ROW][C]0.1778 (5.625)[/C][C]7.282166[/C][/ROW]
[ROW][C]0.1852 (5.4)[/C][C]4.862474[/C][/ROW]
[ROW][C]0.1926 (5.1923)[/C][C]0.26397[/C][/ROW]
[ROW][C]0.2 (5)[/C][C]1.511483[/C][/ROW]
[ROW][C]0.2074 (4.8214)[/C][C]4.465788[/C][/ROW]
[ROW][C]0.2148 (4.6552)[/C][C]1.148821[/C][/ROW]
[ROW][C]0.2222 (4.5)[/C][C]2.951009[/C][/ROW]
[ROW][C]0.2296 (4.3548)[/C][C]3.127388[/C][/ROW]
[ROW][C]0.237 (4.2188)[/C][C]0.017582[/C][/ROW]
[ROW][C]0.2444 (4.0909)[/C][C]0.061233[/C][/ROW]
[ROW][C]0.2519 (3.9706)[/C][C]14.202995[/C][/ROW]
[ROW][C]0.2593 (3.8571)[/C][C]1.118618[/C][/ROW]
[ROW][C]0.2667 (3.75)[/C][C]0.751987[/C][/ROW]
[ROW][C]0.2741 (3.6486)[/C][C]4.462532[/C][/ROW]
[ROW][C]0.2815 (3.5526)[/C][C]2.095525[/C][/ROW]
[ROW][C]0.2889 (3.4615)[/C][C]0.803037[/C][/ROW]
[ROW][C]0.2963 (3.375)[/C][C]0.87938[/C][/ROW]
[ROW][C]0.3037 (3.2927)[/C][C]1.023119[/C][/ROW]
[ROW][C]0.3111 (3.2143)[/C][C]3.12649[/C][/ROW]
[ROW][C]0.3185 (3.1395)[/C][C]0.635779[/C][/ROW]
[ROW][C]0.3259 (3.0682)[/C][C]6.151054[/C][/ROW]
[ROW][C]0.3333 (3)[/C][C]6.750947[/C][/ROW]
[ROW][C]0.3407 (2.9348)[/C][C]0.276941[/C][/ROW]
[ROW][C]0.3481 (2.8723)[/C][C]1.221282[/C][/ROW]
[ROW][C]0.3556 (2.8125)[/C][C]0.157484[/C][/ROW]
[ROW][C]0.363 (2.7551)[/C][C]1.21057[/C][/ROW]
[ROW][C]0.3704 (2.7)[/C][C]0.047862[/C][/ROW]
[ROW][C]0.3778 (2.6471)[/C][C]1.146836[/C][/ROW]
[ROW][C]0.3852 (2.5962)[/C][C]0.726227[/C][/ROW]
[ROW][C]0.3926 (2.5472)[/C][C]0.967389[/C][/ROW]
[ROW][C]0.4 (2.5)[/C][C]0.082001[/C][/ROW]
[ROW][C]0.4074 (2.4545)[/C][C]2.805805[/C][/ROW]
[ROW][C]0.4148 (2.4107)[/C][C]0.130709[/C][/ROW]
[ROW][C]0.4222 (2.3684)[/C][C]1.760301[/C][/ROW]
[ROW][C]0.4296 (2.3276)[/C][C]3.66678[/C][/ROW]
[ROW][C]0.437 (2.2881)[/C][C]2.003551[/C][/ROW]
[ROW][C]0.4444 (2.25)[/C][C]0.842837[/C][/ROW]
[ROW][C]0.4519 (2.2131)[/C][C]0.228989[/C][/ROW]
[ROW][C]0.4593 (2.1774)[/C][C]0.475707[/C][/ROW]
[ROW][C]0.4667 (2.1429)[/C][C]0.278949[/C][/ROW]
[ROW][C]0.4741 (2.1094)[/C][C]0.141589[/C][/ROW]
[ROW][C]0.4815 (2.0769)[/C][C]0.538107[/C][/ROW]
[ROW][C]0.4889 (2.0455)[/C][C]2.723041[/C][/ROW]
[ROW][C]0.4963 (2.0149)[/C][C]7.609089[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=115165&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=115165&T=1

As an alternative you can also use a QR Code:  
QR Code


The GUIDs for individual cells are displayed in the table below:

Raw Periodogram
ParameterValue
Box-Cox transformation parameter (lambda)1
Degree of non-seasonal differencing (d)0
Degree of seasonal differencing (D)0
Seasonal Period (s)12
Frequency (Period)Spectrum
0.0074 (135)274.205195
0.0148 (67.5)123.904683
0.0222 (45)61.536764
0.0296 (33.75)10.29903
0.037 (27)5.067529
0.0444 (22.5)3.442247
0.0519 (19.2857)5.777284
0.0593 (16.875)11.269186
0.0667 (15)5.569656
0.0741 (13.5)5.489349
0.0815 (12.2727)44.886933
0.0889 (11.25)5.822954
0.0963 (10.3846)7.198393
0.1037 (9.6429)5.818261
0.1111 (9)2.508553
0.1185 (8.4375)0.010607
0.1259 (7.9412)2.853442
0.1333 (7.5)0.239436
0.1407 (7.1053)3.066153
0.1481 (6.75)0.836111
0.1556 (6.4286)18.394461
0.163 (6.1364)21.939773
0.1704 (5.8696)51.662073
0.1778 (5.625)7.282166
0.1852 (5.4)4.862474
0.1926 (5.1923)0.26397
0.2 (5)1.511483
0.2074 (4.8214)4.465788
0.2148 (4.6552)1.148821
0.2222 (4.5)2.951009
0.2296 (4.3548)3.127388
0.237 (4.2188)0.017582
0.2444 (4.0909)0.061233
0.2519 (3.9706)14.202995
0.2593 (3.8571)1.118618
0.2667 (3.75)0.751987
0.2741 (3.6486)4.462532
0.2815 (3.5526)2.095525
0.2889 (3.4615)0.803037
0.2963 (3.375)0.87938
0.3037 (3.2927)1.023119
0.3111 (3.2143)3.12649
0.3185 (3.1395)0.635779
0.3259 (3.0682)6.151054
0.3333 (3)6.750947
0.3407 (2.9348)0.276941
0.3481 (2.8723)1.221282
0.3556 (2.8125)0.157484
0.363 (2.7551)1.21057
0.3704 (2.7)0.047862
0.3778 (2.6471)1.146836
0.3852 (2.5962)0.726227
0.3926 (2.5472)0.967389
0.4 (2.5)0.082001
0.4074 (2.4545)2.805805
0.4148 (2.4107)0.130709
0.4222 (2.3684)1.760301
0.4296 (2.3276)3.66678
0.437 (2.2881)2.003551
0.4444 (2.25)0.842837
0.4519 (2.2131)0.228989
0.4593 (2.1774)0.475707
0.4667 (2.1429)0.278949
0.4741 (2.1094)0.141589
0.4815 (2.0769)0.538107
0.4889 (2.0455)2.723041
0.4963 (2.0149)7.609089


Figures (Output of Computation)

Input Parameters & R Code

Parameters (Session):
par1 = 1 ; par2 = 0 ; par3 = 0 ; par4 = 12 ;
Parameters (R input):
par1 = 1 ; par2 = 0 ; par3 = 0 ; par4 = 12 ; par5 = ; par6 = ; par7 = ; par8 = ; par9 = ; par10 = ; par11 = ; par12 = ; par13 = ; par14 = ; par15 = ; par16 = ; par17 = ; par18 = ; par19 = ; par20 = ;
R code (references can be found in the software module):
par1 <- as.numeric(par1)
par2 <- as.numeric(par2)
par3 <- as.numeric(par3)
par4 <- as.numeric(par4)
if (par1 == 0) {
x <- log(x)
} else {
x <- (x ^ par1 - 1) / par1
}
if (par2 > 0) x <- diff(x,lag=1,difference=par2)
if (par3 > 0) x <- diff(x,lag=par4,difference=par3)
bitmap(file='test1.png')
r <- spectrum(x,main='Raw Periodogram')
dev.off()
bitmap(file='test2.png')
cpgram(x,main='Cumulative Periodogram')
dev.off()
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Raw Periodogram',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,'Box-Cox transformation parameter (lambda)',header=TRUE)
a<-table.element(a,par1)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Degree of non-seasonal differencing (d)',header=TRUE)
a<-table.element(a,par2)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Degree of seasonal differencing (D)',header=TRUE)
a<-table.element(a,par3)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Seasonal Period (s)',header=TRUE)
a<-table.element(a,par4)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Frequency (Period)',header=TRUE)
a<-table.element(a,'Spectrum',header=TRUE)
a<-table.row.end(a)
for (i in 1:length(r$freq)) {
a<-table.row.start(a)
mylab <- round(r$freq[i],4)
mylab <- paste(mylab,' (',sep='')
mylab <- paste(mylab,round(1/r$freq[i],4),sep='')
mylab <- paste(mylab,')',sep='')
a<-table.element(a,mylab,header=TRUE)
a<-table.element(a,round(r$spec[i],6))
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable.tab')